karsten — Tue, 11 Nov 2014 17:08:44 GMT

Does this problem still exists for you with boost version > 1.56 ? There are checks for numeric limits within less_eq_with_sign .

Replying to vygis.d@…:

Problem description:

integrate_adaptive goes into infinite loop with rosenbrock4 dense stepper.

After the call to

st.initialize( st.current_state() , st.current_time() , static_cast<Time>(end_time - st.current_time()) );

Loop condition

while( less_eq_with_sign( static_cast<Time>(st.current_time() + st.current_time_step()) ,

end_time , st.current_time_step() ) )

evaluates to false, when t1-t2 < std::numeric_limits<T>::epsilon(); is evaluated.

Expected behaviour:

No infinite loops for any well behaving (albeit stiff) system.

status changed; resolution set

mariomulansky — Wed, 10 May 2017 04:19:05 GMT

status new → closed
resolution → fixed

Boost C++ Libraries: Ticket #10498: Boost.Odeint: Infinite loop in integrate_adaptive with rosenbrock4 dense output stepper

status changed; resolution set