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

Subject: Re: [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-11-11 17:08:44

#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       |  Component:  odeint
  Determined             |   Severity:  Problem
   Version:  Boost       |   Keywords:  infinite loop, integrator algorithm
  1.56.0                 |
Resolution:              |
-------------------------+-------------------------------------------------

Comment (by karsten):

 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 [ticket:10498 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.
 >

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