| 178 | 
  | 
  dt2 = 0.5 * dt; | 
| 179 | 
  | 
 | 
| 180 | 
  | 
  readyCheck(); | 
| 181 | 
+ | 
 | 
| 182 | 
+ | 
  // remove center of mass drift velocity (in case we passed in a configuration | 
| 183 | 
+ | 
  // that was drifting | 
| 184 | 
+ | 
  tStats->removeCOMdrift(); | 
| 185 | 
  | 
 | 
| 186 | 
  | 
  // initialize the forces before the first step | 
| 187 | 
  | 
 | 
| 346 | 
  | 
    integrableObjects[i]->getVel(vel); | 
| 347 | 
  | 
    integrableObjects[i]->getPos(pos); | 
| 348 | 
  | 
    integrableObjects[i]->getFrc(frc); | 
| 349 | 
+ | 
 | 
| 350 | 
+ | 
    std::cerr << "i =\t" << i << "\t" << frc[0] << "\t" << frc[1]<< "\t" << frc[2] << "\n"; | 
| 351 | 
  | 
     | 
| 352 | 
  | 
    mass = integrableObjects[i]->getMass(); | 
| 353 | 
  | 
 |