--- trunk/OOPSE/libmdtools/SimInfo.cpp 2004/04/15 16:18:26 1113 +++ trunk/OOPSE/libmdtools/SimInfo.cpp 2004/04/19 22:13:01 1125 @@ -326,8 +326,12 @@ int SimInfo::getNDF(){ for(int i = 0; i < integrableObjects.size(); i++){ ndf_local += 3; - if (integrableObjects[i]->isDirectional()) - ndf_local += 3; + if (integrableObjects[i]->isDirectional()) { + if (integrableObjects[i]->isLinear()) + ndf_local += 2; + else + ndf_local += 3; + } } // n_constraints is local, so subtract them on each processor: @@ -345,6 +349,8 @@ int SimInfo::getNDF(){ ndf = ndf - 3 - nZconstraints; + std::cerr << "ndf = " << ndf; + return ndf; } @@ -356,8 +362,12 @@ int SimInfo::getNDFraw() { for(int i = 0; i < integrableObjects.size(); i++){ ndfRaw_local += 3; - if (integrableObjects[i]->isDirectional()) - ndfRaw_local += 3; + if (integrableObjects[i]->isDirectional()) { + if (integrableObjects[i]->isLinear()) + ndfRaw_local += 2; + else + ndfRaw_local += 3; + } } #ifdef IS_MPI @@ -581,13 +591,3 @@ GenericData* SimInfo::getProperty(const string& propNa return NULL; } -vector SimInfo::getProperties(){ - - vector result; - map::iterator i; - - for(i = properties.begin(); i != properties.end(); i++) - result.push_back((*i).second); - - return result; -}