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