--- trunk/OOPSE/libmdtools/ConstraintManager.cpp 2004/06/03 21:51:55 1232 +++ trunk/OOPSE/libmdtools/ConstraintManager.cpp 2004/08/23 15:11:36 1452 @@ -37,17 +37,17 @@ void ConstraintManager::addConstraints(Molecule* mol){ vector tempConsElem; foundResult = iterInfoMap.find(mol); - if (foundResult != iterInfoMap.end()){ + if (foundResult == iterInfoMap.end()){ //to improve efficiency, we will not record the molecule without constraint - if (mol->getNConstrains() == 0){ + if (mol->getNConstrains() > 0){ //if a molecule has constraint, all of its myIntegrableOjects must be constraint elements integrableObjects = mol->getIntegrableObjects(); for(sdIter = integrableObjects.begin(); sdIter != integrableObjects.end(); sdIter++){ - + sd = *sdIter; if(sd->isRigidBody()) ce = new ConstraintRigidBody((RigidBody*)sd, 0); else @@ -85,11 +85,11 @@ void ConstraintManager::removeConstraints(Molecule* mo } -ConstraintElementIterator* ConstraintManager::creatElementIterator(){ +ConstraintElementIterator* ConstraintManager::createElementIterator(){ return new ConstraintElementIterator(consElements); } -ConstraintPairIterator* ConstraintManager::creatPairIterator(){ +ConstraintPairIterator* ConstraintManager::createPairIterator(){ return new ConstraintPairIterator(consPairs); }