headmyshoulder / odeint-v2

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

problem about usging integrate_const when end_time euqals to dt #249

Open CooperBruce opened 4 years ago

CooperBruce commented 4 years ago

i am testing OdeInt under Ubuntu 18.04. using the offical example code harmonic_oscillator,when i using integrate_const with end_time euqals to dt,coredump occurs,here is the code:

  `state_type x(2);
x[0] = 1.0; // start at x=1.0, p=0.0
x[1] = 0.0;
//[ define_adapt_stepper

double abs_err = 1.0e-10, rel_err = 1.0e-6;
integrate_const(
        make_dense_output(abs_err , rel_err, runge_kutta_dopri5< state_type >()),
        harmonicoscillator::harmonic_oscillator, x, 0.0, 0.001,
        0.001);`

when i use integrate_adaptive instead,it is OK.

is the way i use integrate_const not right?