[Boost-commit] svn:boost r65646 - in sandbox/odeint/branches/karsten: boost/numeric/odeint/stepper libs/numeric/odeint/test

Subject: [Boost-commit] svn:boost r65646 - in sandbox/odeint/branches/karsten: boost/numeric/odeint/stepper libs/numeric/odeint/test
From: mario.mulansky_at_[hidden]
Date: 2010-09-28 05:49:02

Author: mariomulansky
Date: 2010-09-28 05:49:00 EDT (Tue, 28 Sep 2010)
New Revision: 65646
URL: http://svn.boost.org/trac/boost/changeset/65646

Log:
fixed small bug in error_checker
controlled_stepper_tests working now
Text files modified: 
   sandbox/odeint/branches/karsten/boost/numeric/odeint/stepper/error_checker.hpp      |     2 +-                                      
   sandbox/odeint/branches/karsten/libs/numeric/odeint/test/check_stepper_concepts.cpp |     4 ++--                                    
   2 files changed, 3 insertions(+), 3 deletions(-)

Modified: sandbox/odeint/branches/karsten/boost/numeric/odeint/stepper/error_checker.hpp
==============================================================================
--- sandbox/odeint/branches/karsten/boost/numeric/odeint/stepper/error_checker.hpp	(original)
+++ sandbox/odeint/branches/karsten/boost/numeric/odeint/stepper/error_checker.hpp	2010-09-28 05:49:00 EDT (Tue, 28 Sep 2010)
@@ -43,7 +43,7 @@
                 algebra_type::for_each3( x_old , dxdt_old , x_err ,
                                                      typename operations_type::rel_error( m_eps_abs , m_eps_rel , m_a_x , m_a_dxdt*dt ) );
 
-		return algebra_type::template reduce< time_type >( x_err , typename operations_type::maximum() , 100.0 );
+		return algebra_type::template reduce< time_type >( x_err , typename operations_type::maximum() , 0.0 );
         }
 
 private:

Modified: sandbox/odeint/branches/karsten/libs/numeric/odeint/test/check_stepper_concepts.cpp
==============================================================================
--- sandbox/odeint/branches/karsten/libs/numeric/odeint/test/check_stepper_concepts.cpp	(original)
+++ sandbox/odeint/branches/karsten/libs/numeric/odeint/test/check_stepper_concepts.cpp	2010-09-28 05:49:00 EDT (Tue, 28 Sep 2010)
@@ -5,7 +5,6 @@
  
  This file tests the use of the all different steppers with several state types:
  std::vector< double >
- gsl_vector
  vector_space_1d< double >  (see vector_space_1d.hpp)
  std::tr1::array< double , 1 >
   
@@ -98,7 +97,8 @@
     typedef typename stepper_type::time_type time_type;
 
     time_type t = 0.0 , dt = 0.1;
-    stepper.try_step( system , x , t , dt );
+    controlled_step_result step_result = stepper.try_step( system , x , t , dt );
+    BOOST_CHECK( step_result == success_step_size_increased ); // error = 0 for constant system -> step size is always too small
 }