void System::update( float dt ) {
computeGravitation();
integrateSystem( dt );
Vector3f r, v;
for( size_t i = 0; i < _nBodies; ++i ){
// printf("start of for\n");
// printf("X: pos = %g, vel = %g \n", _body[i].position().x(), _body[i].velocity().x());
r = _body[i].position();
v = _body[i].velocity();
_body[i].oldPosition() = r;
_body[i].oldVelocity() = v;
// printf("X old: pos = %g, vel = %g \n", _body[i].oldPosition().x(), _body[i].oldVelocity().x());
// printf("Y old: pos = %g \n", _body[i].oldPosition().y());
// printf("Z old: pos = %g \n", _body[i].oldPosition().z());
// printf("end of for\n");
}
integrateSystem( dt );
average();
}
void System::update( float dt ) {
computeGravitation();
integrateSystem( dt );
}
