size_t integrate_adaptive(
state_type &start_state,
const time_type start_time,
const time_type end_time,
const time_type dt_init)
{
ode::Solver<ode::Steppers::RKDP5,CSystem> solver(std::ref(*this));
last_observed_time=start_time;
last_observed_intermediates.zeros();
next_sudden_change_time=end_time;
size_t return_val=
solver.integrate_adaptive(
start_state,// is manipulated
start_time ,
end_time ,
dt_init);
results_finalize();
post_solve();
// std::cout<<"Results("<<files::now_print<<") SHA1: ";
// std::cout<<sha1sum(results).signature()<<std::endl;
return return_val;
}
void PostSolve(b2Contact* contact, const b2ContactImpulse* impulse) {
post_solve(object, contact, impulse);
}
