search

headmyshoulder / odeint-v2

odeint - solving ordinary differential equations in c++ v2
http://headmyshoulder.github.com/odeint-v2/
Other
344 stars 101 forks source link

integrate_adaptive wrong behaviour #40

Closed Tohiko closed 11 years ago

Tohiko commented 12 years ago

It seems that calling integrate_adaptive with the final time equal to the time-step makes it nonadaptive. That is, the method is applied for a single time step regardless of the error.

I tried this with both methods: runge_kutta_cash_karp54 runge_kutta_dopri5

mariomulansky commented 12 years ago

Hrm I can not reproduce that. Playing around with the lorenz_point.cpp in examples I can see a clear dependence of the steps taken when changing the error bounds. Also, choosing dt = t_end doesn't cause any problems at all, the method behaves as it should, decreasing the step size before getting to work. Can you add a minimal example showing the observed behavior?

mariomulansky commented 11 years ago

no answer in >6 month, close