example stiff_system.cpp performance

headmyshoulder / odeint-v2

odeint - solving ordinary differential equations in c++ v2

http://headmyshoulder.github.com/odeint-v2/

Other

341 stars 101 forks source link

example stiff_system.cpp performance #189

Closed xlai closed 8 years ago

xlai commented 8 years ago

Hi, I followed the example without any change and compiled it with g++ 5.2.1. Using the rosenbrock4 stepper, it finishes in 40030 steps rather than 71 steps, as described in the document. The dopri5 stepper finishes in 1531 steps which agrees with the document. I'm not sure where it went wrong.

Thank you!

mariomulansky commented 8 years ago

I can confirm this observation. This is definitely some bug, I will look into this.

bjodah commented 8 years ago

@mariomulansky do you know when this fix will be in boost? (it's not in 1.62 released end of august..) EDIT: and not in 1.63