| 37 |  | vector<ConstraintElement*> tempConsElem; | 
| 38 |  |  | 
| 39 |  | foundResult = iterInfoMap.find(mol); | 
| 40 | < | if (foundResult != iterInfoMap.end()){ | 
| 40 | > | if (foundResult == iterInfoMap.end()){ | 
| 41 |  |  | 
| 42 |  | //to improve efficiency, we will not record the molecule without constraint | 
| 43 | < | if (mol->getNConstrains() == 0){ | 
| 43 | > | if (mol->getNConstrains() > 0){ | 
| 44 |  |  | 
| 45 |  | //if a molecule has constraint, all of its myIntegrableOjects must be constraint elements | 
| 46 |  |  | 
| 47 |  | integrableObjects = mol->getIntegrableObjects(); | 
| 48 |  |  | 
| 49 |  | for(sdIter = integrableObjects.begin(); sdIter != integrableObjects.end(); sdIter++){ | 
| 50 | < |  | 
| 50 | > | sd = *sdIter; | 
| 51 |  | if(sd->isRigidBody()) | 
| 52 |  | ce = new ConstraintRigidBody((RigidBody*)sd, 0); | 
| 53 |  | else |