[Boost-commit] svn:boost r68478 - sandbox/odeint/boost/numeric/odeint

Subject: [Boost-commit] svn:boost r68478 - sandbox/odeint/boost/numeric/odeint
From: karsten.ahnert_at_[hidden]
Date: 2011-01-27 12:53:48

Author: karsten
Date: 2011-01-27 12:53:46 EST (Thu, 27 Jan 2011)
New Revision: 68478
URL: http://svn.boost.org/trac/boost/changeset/68478

Log:
* halp stepper modified
Text files modified: 
   sandbox/odeint/boost/numeric/odeint/stepper_half_step.hpp |     6 ++++--                                  
   1 files changed, 4 insertions(+), 2 deletions(-)

Modified: sandbox/odeint/boost/numeric/odeint/stepper_half_step.hpp
==============================================================================
--- sandbox/odeint/boost/numeric/odeint/stepper_half_step.hpp	(original)
+++ sandbox/odeint/boost/numeric/odeint/stepper_half_step.hpp	2011-01-27 12:53:46 EST (Thu, 27 Jan 2011)
@@ -51,6 +51,7 @@
         container_type m_dxdt;
         container_type m_xtemp;
         stepper_type m_stepper;
+        const time_time m_fac;
         
 
         //
@@ -82,6 +83,7 @@
 
         // standard constructor
         stepper_half_step( void )
+        : m_fac( 1.0 / std::pow( 2.0 , m_stepper.order_step() ) - 1.0 )
         {
         }
 
@@ -146,9 +148,9 @@
 
             detail::it_algebra::scale_sum( traits_type::begin(xerr) ,
                                            traits_type::end(xerr) ,
-                                           static_cast< value_type >(1.0),
+                                           static_cast< value_type >( m_fac ),
                                            traits_type::begin(m_xtemp) ,
-                                           static_cast< value_type >(-1.0),
+                                           static_cast< value_type >(- m_fac ),
                                            traits_type::begin(x) );
         }