| 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 |