363 |
|
integrableObjects[i]->getVel(vel); |
364 |
|
integrableObjects[i]->getPos(pos); |
365 |
|
integrableObjects[i]->getFrc(frc); |
366 |
< |
std::cerr << "f = " << frc[0] << "\t" << frc[1] << "\t" << frc[2] << "\n"; |
366 |
> |
// std::cerr << "f = " << frc[0] << "\t" << frc[1] << "\t" << frc[2] << "\n"; |
367 |
|
|
368 |
|
mass = integrableObjects[i]->getMass(); |
369 |
|
|
384 |
|
|
385 |
|
integrableObjects[i]->getTrq(Tb); |
386 |
|
|
387 |
< |
std::cerr << "t = " << Tb[0] << "\t" << Tb[1] << "\t" << Tb[2] << "\n"; |
387 |
> |
// std::cerr << "t = " << Tb[0] << "\t" << Tb[1] << "\t" << Tb[2] << "\n"; |
388 |
|
integrableObjects[i]->lab2Body(Tb); |
389 |
|
|
390 |
|
// get the angular momentum, and propagate a half step |