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