[Boost-bugs] [Boost C++ Libraries] #10498: Boost.Odeint: Infinite loop in integrate_adaptive with rosenbrock4 dense output stepper

Subject: [Boost-bugs] [Boost C++ Libraries] #10498: Boost.Odeint: Infinite loop in integrate_adaptive with rosenbrock4 dense output stepper
From: Boost C++ Libraries (noreply_at_[hidden])
Date: 2014-09-12 12:10:37

#10498: Boost.Odeint: Infinite loop in integrate_adaptive with rosenbrock4 dense
output stepper
-------------------------------------------------+---------------------
 Reporter:  vygis.d@…                            |      Owner:  karsten
     Type:  Bugs                                 |     Status:  new
Milestone:  To Be Determined                     |  Component:  odeint
  Version:  Boost 1.56.0                         |   Severity:  Problem
 Keywords:  infinite loop, integrator algorithm  |
-------------------------------------------------+---------------------
 '''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.

-- 
Ticket URL: <https://svn.boost.org/trac/boost/ticket/10498>
Boost C++ Libraries <http://www.boost.org/>
Boost provides free peer-reviewed portable C++ source libraries.