From 63fade416ca69b2b5dcaa08d5783c416fe6b2d2f Mon Sep 17 00:00:00 2001 From: Mario Mulansky Date: Wed, 8 Feb 2012 19:05:40 +0100 Subject: [PATCH] removed is_resizable from state_wrapper --- .../odeint/external/gsl/gsl_wrapper.hpp | 1 - .../odeint/stepper/adams_bashforth.hpp | 3 +- .../numeric/odeint/stepper/adams_moulton.hpp | 3 +- .../base/explicit_error_stepper_base.hpp | 2 +- ...xplicit_stepper_and_error_stepper_base.hpp | 3 +- ...it_stepper_and_error_stepper_fsal_base.hpp | 3 +- .../stepper/base/explicit_stepper_base.hpp | 2 +- .../base/symplectic_rkn_stepper_base.hpp | 4 +- .../numeric/odeint/stepper/bulirsch_stoer.hpp | 9 +- .../stepper/bulirsch_stoer_dense_out.hpp | 19 +-- .../odeint/stepper/controlled_runge_kutta.hpp | 15 +- .../stepper/dense_output_runge_kutta.hpp | 13 +- .../stepper/explicit_error_generic_rk.hpp | 5 +- .../odeint/stepper/explicit_generic_rk.hpp | 5 +- .../numeric/odeint/stepper/implicit_euler.hpp | 11 +- .../odeint/stepper/modified_midpoint.hpp | 11 +- boost/numeric/odeint/stepper/rosenbrock4.hpp | 30 ++-- .../odeint/stepper/rosenbrock4_controller.hpp | 5 +- .../stepper/rosenbrock4_dense_output.hpp | 5 +- .../odeint/stepper/runge_kutta4_classic.hpp | 9 +- .../stepper/runge_kutta_cash_karp54.hpp | 2 + .../runge_kutta_cash_karp54_classic.hpp | 13 +- .../odeint/stepper/runge_kutta_dopri5.hpp | 17 +- .../odeint/stepper/runge_kutta_fehlberg78.hpp | 2 + boost/numeric/odeint/util/state_wrapper.hpp | 23 +-- boost/numeric/odeint/util/ublas_wrapper.hpp | 21 +-- .../concepts.html | 3 +- .../concepts/controlled_stepper.html | 27 ++-- .../concepts/dense_output_stepper.html | 33 ++-- .../concepts/error_stepper.html | 27 ++-- .../concepts/implicit_system.html | 23 +-- .../concepts/observer.html | 21 +-- .../concepts/simple_symplectic_system.html | 21 +-- .../concepts/state_algebra_operations.html | 66 ++++---- .../concepts/state_wrapper.html | 21 +-- .../concepts/stepper.html | 33 ++-- .../concepts/symplectic_system.html | 21 +-- .../concepts/system.html | 21 +-- .../extend_odeint.html | 3 +- .../adapt_your_own_operations.html | 3 +- .../adapt_your_own_state_types.html | 91 ++++++----- .../extend_odeint/write_own_steppers.html | 3 +- .../getting_started.html | 3 +- .../getting_started/overview.html | 3 +- .../getting_started/short_example.html | 14 +- .../usage__compilation__headers.html | 6 +- .../odeint_in_detail.html | 3 +- .../algebras_and_operations.html | 3 +- .../generation_functions.html | 3 +- .../odeint_in_detail/integrate_functions.html | 38 ++--- .../odeint_in_detail/steppers.html | 57 +++---- .../odeint_in_detail/using_boost__range.html | 3 +- .../odeint_in_detail/using_boost__ref.html | 3 +- .../old_concepts.html | 3 +- .../old_concepts/basic_stepper.html | 3 +- .../old_concepts/compositestepper.html | 3 +- .../old_concepts/controlled_stepper.html | 3 +- .../old_concepts/dense_ouput_stepper.html | 3 +- .../old_concepts/error_stepper.html | 3 +- .../old_concepts/size_adjusting_stepper.html | 3 +- .../old_reference.html | 3 +- .../old_reference/algebras.html | 3 +- .../old_reference/integration_functions.html | 3 +- .../old_reference/operations.html | 3 +- .../old_reference/resizing.html | 3 +- .../old_reference/stepper_classes.html | 3 +- .../tutorial.html | 3 +- .../tutorial/all_examples.html | 3 +- ...haotic_systems_and_lyapunov_exponents.html | 25 ++- .../tutorial/harmonic_oscillator.html | 12 +- .../tutorial/references.html | 3 +- .../tutorial/solar_system.html | 31 ++-- .../tutorial/special_topics.html | 93 +++++------ .../tutorial/stiff_systems.html | 7 +- .../tutorial/using_cuda_and_thrust.html | 149 ++++++++---------- doc/index.html | 7 +- doc/odeint/reference.html | 3 +- libs/numeric/odeint/test/stepper_copying.cpp | 10 +- 78 files changed, 556 insertions(+), 618 deletions(-) diff --git a/boost/numeric/odeint/external/gsl/gsl_wrapper.hpp b/boost/numeric/odeint/external/gsl/gsl_wrapper.hpp index 8c91cdec..999b06ab 100644 --- a/boost/numeric/odeint/external/gsl/gsl_wrapper.hpp +++ b/boost/numeric/odeint/external/gsl/gsl_wrapper.hpp @@ -199,7 +199,6 @@ struct state_wrapper< gsl_vector* > typedef double value_type; typedef gsl_vector* state_type; typedef state_wrapper< gsl_vector* > state_wrapper_type; - typedef boost::true_type is_resizeable; state_type m_v; diff --git a/boost/numeric/odeint/stepper/adams_bashforth.hpp b/boost/numeric/odeint/stepper/adams_bashforth.hpp index 7ef6b106..f8bf953d 100644 --- a/boost/numeric/odeint/stepper/adams_bashforth.hpp +++ b/boost/numeric/odeint/stepper/adams_bashforth.hpp @@ -27,6 +27,7 @@ #include #include +#include #include #include @@ -290,7 +291,7 @@ private: bool resized( false ); for( size_t i=0 ; i::type() ); } return resized; } diff --git a/boost/numeric/odeint/stepper/adams_moulton.hpp b/boost/numeric/odeint/stepper/adams_moulton.hpp index ccf1df64..793f3486 100644 --- a/boost/numeric/odeint/stepper/adams_moulton.hpp +++ b/boost/numeric/odeint/stepper/adams_moulton.hpp @@ -27,6 +27,7 @@ #include #include +#include #include #include @@ -169,7 +170,7 @@ private: template< class StateIn > bool resize_impl( const StateIn &x ) { - return adjust_size_by_resizeability( m_dxdt , x , typename wrapped_deriv_type::is_resizeable() ); + return adjust_size_by_resizeability( m_dxdt , x , typename is_resizeable::type() ); } diff --git a/boost/numeric/odeint/stepper/base/explicit_error_stepper_base.hpp b/boost/numeric/odeint/stepper/base/explicit_error_stepper_base.hpp index 45cd3b2d..89f886c2 100644 --- a/boost/numeric/odeint/stepper/base/explicit_error_stepper_base.hpp +++ b/boost/numeric/odeint/stepper/base/explicit_error_stepper_base.hpp @@ -155,7 +155,7 @@ private: template< class StateIn > bool resize_impl( const StateIn &x ) { - return adjust_size_by_resizeability( m_dxdt , x , typename wrapped_deriv_type::is_resizeable() ); + return adjust_size_by_resizeability( m_dxdt , x , typename is_resizeable::type() ); } diff --git a/boost/numeric/odeint/stepper/base/explicit_stepper_and_error_stepper_base.hpp b/boost/numeric/odeint/stepper/base/explicit_stepper_and_error_stepper_base.hpp index 2b0337f7..3bd87a9b 100644 --- a/boost/numeric/odeint/stepper/base/explicit_stepper_and_error_stepper_base.hpp +++ b/boost/numeric/odeint/stepper/base/explicit_stepper_and_error_stepper_base.hpp @@ -23,6 +23,7 @@ #include #include +#include #include #include @@ -244,7 +245,7 @@ private: template< class StateIn > bool resize_impl( const StateIn &x ) { - return adjust_size_by_resizeability( m_dxdt , x , typename wrapped_deriv_type::is_resizeable() ); + return adjust_size_by_resizeability( m_dxdt , x , typename is_resizeable::type() ); } stepper_type& stepper( void ) diff --git a/boost/numeric/odeint/stepper/base/explicit_stepper_and_error_stepper_fsal_base.hpp b/boost/numeric/odeint/stepper/base/explicit_stepper_and_error_stepper_fsal_base.hpp index b8bf41af..7e3c1ff2 100644 --- a/boost/numeric/odeint/stepper/base/explicit_stepper_and_error_stepper_fsal_base.hpp +++ b/boost/numeric/odeint/stepper/base/explicit_stepper_and_error_stepper_fsal_base.hpp @@ -23,6 +23,7 @@ #include #include +#include #include #include @@ -283,7 +284,7 @@ private: template< class StateIn > bool resize_impl( const StateIn &x ) { - return adjust_size_by_resizeability( m_dxdt , x , typename wrapped_deriv_type::is_resizeable() ); + return adjust_size_by_resizeability( m_dxdt , x , typename is_resizeable::type() ); } diff --git a/boost/numeric/odeint/stepper/base/explicit_stepper_base.hpp b/boost/numeric/odeint/stepper/base/explicit_stepper_base.hpp index 19257078..1532d8c2 100644 --- a/boost/numeric/odeint/stepper/base/explicit_stepper_base.hpp +++ b/boost/numeric/odeint/stepper/base/explicit_stepper_base.hpp @@ -159,7 +159,7 @@ private: template< class StateIn > bool resize_impl( const StateIn &x ) { - return adjust_size_by_resizeability( m_dxdt , x , typename wrapped_deriv_type::is_resizeable() ); + return adjust_size_by_resizeability( m_dxdt , x , typename is_resizeable::type() ); } diff --git a/boost/numeric/odeint/stepper/base/symplectic_rkn_stepper_base.hpp b/boost/numeric/odeint/stepper/base/symplectic_rkn_stepper_base.hpp index 6cbbbc12..185b924e 100644 --- a/boost/numeric/odeint/stepper/base/symplectic_rkn_stepper_base.hpp +++ b/boost/numeric/odeint/stepper/base/symplectic_rkn_stepper_base.hpp @@ -277,13 +277,13 @@ private: template< class StateIn > bool resize_dqdt( const StateIn &x ) { - return adjust_size_by_resizeability( m_dqdt , x , typename wrapped_coor_deriv_type::is_resizeable() ); + return adjust_size_by_resizeability( m_dqdt , x , typename is_resizeable::type() ); } template< class StateIn > bool resize_dpdt( const StateIn &x ) { - return adjust_size_by_resizeability( m_dpdt , x , typename wrapped_momentum_deriv_type::is_resizeable() ); + return adjust_size_by_resizeability( m_dpdt , x , typename is_resizeable::type() ); } diff --git a/boost/numeric/odeint/stepper/bulirsch_stoer.hpp b/boost/numeric/odeint/stepper/bulirsch_stoer.hpp index abe68d78..2c7eefbe 100644 --- a/boost/numeric/odeint/stepper/bulirsch_stoer.hpp +++ b/boost/numeric/odeint/stepper/bulirsch_stoer.hpp @@ -37,6 +37,7 @@ #include #include +#include #include namespace boost { @@ -341,13 +342,13 @@ private: template< class StateIn > bool resize_m_dxdt( const StateIn &x ) { - return adjust_size_by_resizeability( m_dxdt , x , typename wrapped_deriv_type::is_resizeable() ); + return adjust_size_by_resizeability( m_dxdt , x , typename is_resizeable::type() ); } template< class StateIn > bool resize_m_xnew( const StateIn &x ) { - return adjust_size_by_resizeability( m_xnew , x , typename wrapped_state_type::is_resizeable() ); + return adjust_size_by_resizeability( m_xnew , x , typename is_resizeable::type() ); } template< class StateIn > @@ -355,8 +356,8 @@ private: { bool resized( false ); for( size_t i = 0 ; i < m_k_max ; ++i ) - resized |= adjust_size_by_resizeability( m_table[i] , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_err , x , typename wrapped_state_type::is_resizeable() ); + resized |= adjust_size_by_resizeability( m_table[i] , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_err , x , typename is_resizeable::type() ); return resized; } diff --git a/boost/numeric/odeint/stepper/bulirsch_stoer_dense_out.hpp b/boost/numeric/odeint/stepper/bulirsch_stoer_dense_out.hpp index 78fc2f82..1a347e2f 100644 --- a/boost/numeric/odeint/stepper/bulirsch_stoer_dense_out.hpp +++ b/boost/numeric/odeint/stepper/bulirsch_stoer_dense_out.hpp @@ -36,6 +36,7 @@ #include #include +#include #include #include @@ -762,22 +763,22 @@ private: { bool resized( false ); - resized |= adjust_size_by_resizeability( m_x1 , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_x2 , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_dxdt1 , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_dxdt2 , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_err , x , typename wrapped_state_type::is_resizeable() ); + resized |= adjust_size_by_resizeability( m_x1 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_x2 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_dxdt1 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_dxdt2 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_err , x , typename is_resizeable::type() ); for( size_t i = 0 ; i < m_k_max ; ++i ) - resized |= adjust_size_by_resizeability( m_table[i] , x , typename wrapped_state_type::is_resizeable() ); + resized |= adjust_size_by_resizeability( m_table[i] , x , typename is_resizeable::type() ); for( size_t i = 0 ; i < m_k_max+1 ; ++i ) - resized |= adjust_size_by_resizeability( m_mp_states[i] , x , typename wrapped_state_type::is_resizeable() ); + resized |= adjust_size_by_resizeability( m_mp_states[i] , x , typename is_resizeable::type() ); for( size_t i = 0 ; i < m_k_max+1 ; ++i ) for( size_t j = 0 ; j < m_derivs[i].size() ; ++j ) - resized |= adjust_size_by_resizeability( m_derivs[i][j] , x , typename wrapped_deriv_type::is_resizeable() ); + resized |= adjust_size_by_resizeability( m_derivs[i][j] , x , typename is_resizeable::type() ); for( size_t i = 0 ; i < 2*m_k_max+1 ; ++i ) for( size_t j = 0 ; j < m_diffs[i].size() ; ++j ) - resized |= adjust_size_by_resizeability( m_diffs[i][j] , x , typename wrapped_deriv_type::is_resizeable() ); + resized |= adjust_size_by_resizeability( m_diffs[i][j] , x , typename is_resizeable::type() ); return resized; } diff --git a/boost/numeric/odeint/stepper/controlled_runge_kutta.hpp b/boost/numeric/odeint/stepper/controlled_runge_kutta.hpp index 9ab98681..3e340685 100644 --- a/boost/numeric/odeint/stepper/controlled_runge_kutta.hpp +++ b/boost/numeric/odeint/stepper/controlled_runge_kutta.hpp @@ -27,6 +27,7 @@ #include #include +#include #include #include @@ -295,19 +296,19 @@ private: template< class StateIn > bool resize_m_xerr_impl( const StateIn &x ) { - return adjust_size_by_resizeability( m_xerr , x , typename wrapped_state_type::is_resizeable() ); + return adjust_size_by_resizeability( m_xerr , x , typename is_resizeable::type() ); } template< class StateIn > bool resize_m_dxdt_impl( const StateIn &x ) { - return adjust_size_by_resizeability( m_dxdt , x , typename wrapped_deriv_type::is_resizeable() ); + return adjust_size_by_resizeability( m_dxdt , x , typename is_resizeable::type() ); } template< class StateIn > bool resize_m_xnew_impl( const StateIn &x ) { - return adjust_size_by_resizeability( m_xnew , x , typename wrapped_state_type::is_resizeable() ); + return adjust_size_by_resizeability( m_xnew , x , typename is_resizeable::type() ); } @@ -506,25 +507,25 @@ private: template< class StateIn > bool resize_m_xerr_impl( const StateIn &x ) { - return adjust_size_by_resizeability( m_xerr , x , typename wrapped_state_type::is_resizeable() ); + return adjust_size_by_resizeability( m_xerr , x , typename is_resizeable::type() ); } template< class StateIn > bool resize_m_dxdt_impl( const StateIn &x ) { - return adjust_size_by_resizeability( m_dxdt , x , typename wrapped_deriv_type::is_resizeable() ); + return adjust_size_by_resizeability( m_dxdt , x , typename is_resizeable::type() ); } template< class StateIn > bool resize_m_dxdt_new_impl( const StateIn &x ) { - return adjust_size_by_resizeability( m_dxdtnew , x , typename wrapped_deriv_type::is_resizeable() ); + return adjust_size_by_resizeability( m_dxdtnew , x , typename is_resizeable::type() ); } template< class StateIn > bool resize_m_xnew_impl( const StateIn &x ) { - return adjust_size_by_resizeability( m_xnew , x , typename wrapped_state_type::is_resizeable() ); + return adjust_size_by_resizeability( m_xnew , x , typename is_resizeable::type() ); } diff --git a/boost/numeric/odeint/stepper/dense_output_runge_kutta.hpp b/boost/numeric/odeint/stepper/dense_output_runge_kutta.hpp index f44bef0a..3eca6ec9 100644 --- a/boost/numeric/odeint/stepper/dense_output_runge_kutta.hpp +++ b/boost/numeric/odeint/stepper/dense_output_runge_kutta.hpp @@ -29,6 +29,7 @@ #include #include +#include #include #include @@ -183,8 +184,8 @@ private: bool resize_impl( const StateIn &x ) { bool resized = false; - resized |= adjust_size_by_resizeability( m_x1 , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_x2 , x , typename wrapped_state_type::is_resizeable() ); + resized |= adjust_size_by_resizeability( m_x1 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_x2 , x , typename is_resizeable::type() ); return resized; } @@ -346,10 +347,10 @@ public: bool resize( const StateIn &x ) { bool resized = false; - resized |= adjust_size_by_resizeability( m_x1 , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_x2 , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_dxdt1 , x , typename wrapped_deriv_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_dxdt2 , x , typename wrapped_deriv_type::is_resizeable() ); + resized |= adjust_size_by_resizeability( m_x1 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_x2 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_dxdt1 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_dxdt2 , x , typename is_resizeable::type() ); return resized; } diff --git a/boost/numeric/odeint/stepper/explicit_error_generic_rk.hpp b/boost/numeric/odeint/stepper/explicit_error_generic_rk.hpp index d46e9e02..cd75bb51 100644 --- a/boost/numeric/odeint/stepper/explicit_error_generic_rk.hpp +++ b/boost/numeric/odeint/stepper/explicit_error_generic_rk.hpp @@ -27,6 +27,7 @@ #include #include +#include #include @@ -143,10 +144,10 @@ private: bool resize_impl( const StateIn &x ) { bool resized( false ); - resized |= adjust_size_by_resizeability( m_x_tmp , x , typename wrapped_state_type::is_resizeable() ); + resized |= adjust_size_by_resizeability( m_x_tmp , x , typename is_resizeable::type() ); for( size_t i = 0 ; i < StageCount-1 ; ++i ) { - resized |= adjust_size_by_resizeability( m_F[i] , x , typename wrapped_deriv_type::is_resizeable() ); + resized |= adjust_size_by_resizeability( m_F[i] , x , typename is_resizeable::type() ); } return resized; } diff --git a/boost/numeric/odeint/stepper/explicit_generic_rk.hpp b/boost/numeric/odeint/stepper/explicit_generic_rk.hpp index 1ca982cb..cc8b48a6 100644 --- a/boost/numeric/odeint/stepper/explicit_generic_rk.hpp +++ b/boost/numeric/odeint/stepper/explicit_generic_rk.hpp @@ -29,6 +29,7 @@ #include #include +#include #include namespace mpl = boost::mpl; @@ -180,10 +181,10 @@ private: bool resize_impl( const StateIn &x ) { bool resized( false ); - resized |= adjust_size_by_resizeability( m_x_tmp , x , typename wrapped_state_type::is_resizeable() ); + resized |= adjust_size_by_resizeability( m_x_tmp , x , typename is_resizeable::type() ); for( size_t i = 0 ; i < StageCount-1 ; ++i ) { - resized |= adjust_size_by_resizeability( m_F[i] , x , typename wrapped_deriv_type::is_resizeable() ); + resized |= adjust_size_by_resizeability( m_F[i] , x , typename is_resizeable::type() ); } return resized; } diff --git a/boost/numeric/odeint/stepper/implicit_euler.hpp b/boost/numeric/odeint/stepper/implicit_euler.hpp index 3158ae4e..d0e7074a 100644 --- a/boost/numeric/odeint/stepper/implicit_euler.hpp +++ b/boost/numeric/odeint/stepper/implicit_euler.hpp @@ -27,6 +27,7 @@ #include #include +#include #include #include @@ -131,11 +132,11 @@ private: bool resize_impl( const StateIn &x ) { bool resized = false; - resized |= adjust_size_by_resizeability( m_dxdt , x , typename wrapped_deriv_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_x , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_b , x , typename wrapped_deriv_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_jacobi , x , typename wrapped_matrix_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_pm , x , typename wrapped_pmatrix_type::is_resizeable() ); + resized |= adjust_size_by_resizeability( m_dxdt , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_x , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_b , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_jacobi , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_pm , x , typename is_resizeable::type() ); return resized; } diff --git a/boost/numeric/odeint/stepper/modified_midpoint.hpp b/boost/numeric/odeint/stepper/modified_midpoint.hpp index 41dcc745..dabe0746 100644 --- a/boost/numeric/odeint/stepper/modified_midpoint.hpp +++ b/boost/numeric/odeint/stepper/modified_midpoint.hpp @@ -22,6 +22,7 @@ #include #include +#include #include #include #include @@ -119,9 +120,9 @@ private: bool resize_impl( const StateIn &x ) { bool resized( false ); - resized |= adjust_size_by_resizeability( m_x0 , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_x1 , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_dxdt , x , typename wrapped_state_type::is_resizeable() ); + resized |= adjust_size_by_resizeability( m_x0 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_x1 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_dxdt , x , typename is_resizeable::type() ); return resized; } @@ -242,8 +243,8 @@ public : bool resize( const StateIn &x ) { bool resized( false ); - resized |= adjust_size_by_resizeability( m_x0 , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_x1 , x , typename wrapped_state_type::is_resizeable() ); + resized |= adjust_size_by_resizeability( m_x0 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_x1 , x , typename is_resizeable::type() ); return resized; } diff --git a/boost/numeric/odeint/stepper/rosenbrock4.hpp b/boost/numeric/odeint/stepper/rosenbrock4.hpp index 0fa0fd6b..37b58b7f 100644 --- a/boost/numeric/odeint/stepper/rosenbrock4.hpp +++ b/boost/numeric/odeint/stepper/rosenbrock4.hpp @@ -30,7 +30,7 @@ #include #include -//#include +#include #include #include @@ -269,26 +269,26 @@ protected: bool resize_impl( const StateIn &x ) { bool resized = false; - resized |= adjust_size_by_resizeability( m_dxdt , x , typename wrapped_deriv_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_dfdt , x , typename wrapped_deriv_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_dxdtnew , x , typename wrapped_deriv_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_xtmp , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_g1 , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_g2 , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_g3 , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_g4 , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_g5 , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_cont3 , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_cont4 , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_jac , x , typename wrapped_matrix_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_pm , x , typename wrapped_pmatrix_type::is_resizeable() ); + resized |= adjust_size_by_resizeability( m_dxdt , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_dfdt , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_dxdtnew , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_xtmp , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_g1 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_g2 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_g3 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_g4 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_g5 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_cont3 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_cont4 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_jac , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_pm , x , typename is_resizeable::type() ); return resized; } template< class StateIn > bool resize_x_err( const StateIn &x ) { - return adjust_size_by_resizeability( m_x_err , x , typename wrapped_state_type::is_resizeable() ); + return adjust_size_by_resizeability( m_x_err , x , typename is_resizeable::type() ); } private: diff --git a/boost/numeric/odeint/stepper/rosenbrock4_controller.hpp b/boost/numeric/odeint/stepper/rosenbrock4_controller.hpp index 40cca00e..150e40b5 100644 --- a/boost/numeric/odeint/stepper/rosenbrock4_controller.hpp +++ b/boost/numeric/odeint/stepper/rosenbrock4_controller.hpp @@ -26,6 +26,7 @@ #include #include +#include #include @@ -174,13 +175,13 @@ private: template< class StateIn > bool resize_m_xerr( const StateIn &x ) { - return adjust_size_by_resizeability( m_xerr , x , typename wrapped_state_type::is_resizeable() ); + return adjust_size_by_resizeability( m_xerr , x , typename is_resizeable::type() ); } template< class StateIn > bool resize_m_xnew( const StateIn &x ) { - return adjust_size_by_resizeability( m_xnew , x , typename wrapped_state_type::is_resizeable() ); + return adjust_size_by_resizeability( m_xnew , x , typename is_resizeable::type() ); } diff --git a/boost/numeric/odeint/stepper/rosenbrock4_dense_output.hpp b/boost/numeric/odeint/stepper/rosenbrock4_dense_output.hpp index f6302471..0a6e94fa 100644 --- a/boost/numeric/odeint/stepper/rosenbrock4_dense_output.hpp +++ b/boost/numeric/odeint/stepper/rosenbrock4_dense_output.hpp @@ -25,6 +25,7 @@ #include #include +#include namespace boost { namespace numeric { @@ -177,8 +178,8 @@ private: bool resize_impl( const StateIn &x ) { bool resized = false; - resized |= adjust_size_by_resizeability( m_x1 , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_x2 , x , typename wrapped_state_type::is_resizeable() ); + resized |= adjust_size_by_resizeability( m_x1 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_x2 , x , typename is_resizeable::type() ); return resized; } diff --git a/boost/numeric/odeint/stepper/runge_kutta4_classic.hpp b/boost/numeric/odeint/stepper/runge_kutta4_classic.hpp index 1d103790..ba1cbd35 100644 --- a/boost/numeric/odeint/stepper/runge_kutta4_classic.hpp +++ b/boost/numeric/odeint/stepper/runge_kutta4_classic.hpp @@ -27,6 +27,7 @@ #include #include +#include #include namespace boost { @@ -117,10 +118,10 @@ private: bool resize_impl( const StateIn &x ) { bool resized = false; - resized |= adjust_size_by_resizeability( m_x_tmp , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_dxm , x , typename wrapped_deriv_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_dxt , x , typename wrapped_deriv_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_dxh , x , typename wrapped_deriv_type::is_resizeable() ); + resized |= adjust_size_by_resizeability( m_x_tmp , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_dxm , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_dxt , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_dxh , x , typename is_resizeable::type() ); return resized; } diff --git a/boost/numeric/odeint/stepper/runge_kutta_cash_karp54.hpp b/boost/numeric/odeint/stepper/runge_kutta_cash_karp54.hpp index c9861516..dc6833a6 100644 --- a/boost/numeric/odeint/stepper/runge_kutta_cash_karp54.hpp +++ b/boost/numeric/odeint/stepper/runge_kutta_cash_karp54.hpp @@ -24,6 +24,8 @@ #include #include +#include +#include #include #include diff --git a/boost/numeric/odeint/stepper/runge_kutta_cash_karp54_classic.hpp b/boost/numeric/odeint/stepper/runge_kutta_cash_karp54_classic.hpp index 2e068d46..30d4159e 100644 --- a/boost/numeric/odeint/stepper/runge_kutta_cash_karp54_classic.hpp +++ b/boost/numeric/odeint/stepper/runge_kutta_cash_karp54_classic.hpp @@ -28,6 +28,7 @@ #include #include #include +#include #include namespace boost { @@ -166,12 +167,12 @@ private: bool resize_impl( const StateIn &x ) { bool resized = false; - resized |= adjust_size_by_resizeability( m_x_tmp , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_k2 , x , typename wrapped_deriv_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_k3 , x , typename wrapped_deriv_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_k4 , x , typename wrapped_deriv_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_k5 , x , typename wrapped_deriv_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_k6 , x , typename wrapped_deriv_type::is_resizeable() ); + resized |= adjust_size_by_resizeability( m_x_tmp , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_k2 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_k3 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_k4 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_k5 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_k6 , x , typename is_resizeable::type() ); return resized; } diff --git a/boost/numeric/odeint/stepper/runge_kutta_dopri5.hpp b/boost/numeric/odeint/stepper/runge_kutta_dopri5.hpp index b7724c26..cdf27c6d 100644 --- a/boost/numeric/odeint/stepper/runge_kutta_dopri5.hpp +++ b/boost/numeric/odeint/stepper/runge_kutta_dopri5.hpp @@ -28,6 +28,11 @@ #include #include +#include +#include +#include + + namespace boost { namespace numeric { namespace odeint { @@ -238,12 +243,12 @@ private: bool resize_impl( const StateIn &x ) { bool resized = false; - resized |= adjust_size_by_resizeability( m_x_tmp , x , typename wrapped_state_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_k2 , x , typename wrapped_deriv_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_k3 , x , typename wrapped_deriv_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_k4 , x , typename wrapped_deriv_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_k5 , x , typename wrapped_deriv_type::is_resizeable() ); - resized |= adjust_size_by_resizeability( m_k6 , x , typename wrapped_deriv_type::is_resizeable() ); + resized |= adjust_size_by_resizeability( m_x_tmp , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_k2 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_k3 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_k4 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_k5 , x , typename is_resizeable::type() ); + resized |= adjust_size_by_resizeability( m_k6 , x , typename is_resizeable::type() ); return resized; } diff --git a/boost/numeric/odeint/stepper/runge_kutta_fehlberg78.hpp b/boost/numeric/odeint/stepper/runge_kutta_fehlberg78.hpp index 778df40d..b3edf97e 100644 --- a/boost/numeric/odeint/stepper/runge_kutta_fehlberg78.hpp +++ b/boost/numeric/odeint/stepper/runge_kutta_fehlberg78.hpp @@ -27,6 +27,8 @@ #include +#include +#include #include diff --git a/boost/numeric/odeint/util/state_wrapper.hpp b/boost/numeric/odeint/util/state_wrapper.hpp index 4d3e28f5..3c2bfe0b 100644 --- a/boost/numeric/odeint/util/state_wrapper.hpp +++ b/boost/numeric/odeint/util/state_wrapper.hpp @@ -32,34 +32,15 @@ namespace numeric { namespace odeint { - - -template< class V , bool resizeable = is_resizeable< V >::value > -struct state_wrapper; - - - -//two standard implementations, with and without resizing depending on is_resizeable< StateType > - template< class V > -struct state_wrapper< V , true > // with resizing +struct state_wrapper { - typedef state_wrapper< V , true > state_wrapper_type; - typedef boost::true_type is_resizeable; + typedef state_wrapper< V > state_wrapper_type; V m_v; }; -template< class V > -struct state_wrapper< V , false > // without resizing -{ - typedef state_wrapper< V , false > state_wrapper_type; - typedef boost::false_type is_resizeable; - - V m_v; -}; - } } } diff --git a/boost/numeric/odeint/util/ublas_wrapper.hpp b/boost/numeric/odeint/util/ublas_wrapper.hpp index 7c7e044d..64833274 100644 --- a/boost/numeric/odeint/util/ublas_wrapper.hpp +++ b/boost/numeric/odeint/util/ublas_wrapper.hpp @@ -140,35 +140,16 @@ struct resize_impl< boost::numeric::ublas::vector< T_V , A_V > , boost::numeric: template< class T , class A > -struct state_wrapper< boost::numeric::ublas::permutation_matrix< T , A > , true > // with resizing +struct state_wrapper< boost::numeric::ublas::permutation_matrix< T , A > > // with resizing { typedef boost::numeric::ublas::permutation_matrix< T , A > state_type; typedef state_wrapper< state_type > state_wrapper_type; - //typedef typename V::value_type value_type; - typedef boost::true_type is_resizeable; state_type m_v; state_wrapper() : m_v( 1 ) // permutation matrix constructor requires a size, choose 1 as default { } - template< class T_V , class A_V > - bool same_size( const boost::numeric::ublas::vector< T_V , A_V > &x ) - { - return boost::numeric::odeint::same_size_impl< state_type , boost::numeric::ublas::vector< T_V , A_V > >::same_size( m_v , x ); - } - - template< class T_V , class A_V > - bool resize( const boost::numeric::ublas::vector< T_V , A_V > &x ) - { - //standard resizing done like for std::vector - if( !same_size( x ) ) - { - boost::numeric::odeint::resize_impl< state_type , boost::numeric::ublas::vector< T_V , A_V > >::resize( m_v , x ); - return true; - } else - return false; - } }; diff --git a/doc/boost_sandbox_numeric_odeint/concepts.html b/doc/boost_sandbox_numeric_odeint/concepts.html index 4b916cc1..0fd682f2 100644 --- a/doc/boost_sandbox_numeric_odeint/concepts.html +++ b/doc/boost_sandbox_numeric_odeint/concepts.html @@ -44,7 +44,8 @@ -
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/concepts/controlled_stepper.html b/doc/boost_sandbox_numeric_odeint/concepts/controlled_stepper.html index 6b56a819..153c5230 100644 --- a/doc/boost_sandbox_numeric_odeint/concepts/controlled_stepper.html +++ b/doc/boost_sandbox_numeric_odeint/concepts/controlled_stepper.html @@ -25,9 +25,9 @@ to be used within integrate functions.

-
- - Description +
+ + Description

A controlled stepper following this Controlled Stepper concept provides the @@ -36,9 +36,9 @@ with a given step-size dt. Depending on an error estimate of the solution the step might be rejected and a smaller step-size is suggested.

-
- - Notation +
+ + Notation

@@ -77,9 +77,9 @@

-
- - Valid +
+ + Valid Expressions
@@ -178,9 +178,9 @@
-
- - Models +
+ + Models
  • @@ -202,7 +202,8 @@
-
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/concepts/dense_output_stepper.html b/doc/boost_sandbox_numeric_odeint/concepts/dense_output_stepper.html index 86e56619..112ee1f9 100644 --- a/doc/boost_sandbox_numeric_odeint/concepts/dense_output_stepper.html +++ b/doc/boost_sandbox_numeric_odeint/concepts/dense_output_stepper.html @@ -25,9 +25,9 @@ to be used within integrate functions.

-
- - Description +
+ + Description

A dense ouput stepper following this Dense Output Stepper concept provides @@ -37,9 +37,9 @@ to calculate the state x(t') at any point t <= t' <= t+dt.

-
- - Associated +
+ + Associated types
@@ -55,9 +55,9 @@

-
- - Notation +
+ + Notation

@@ -92,9 +92,9 @@

-
- - Valid +
+ + Valid Expressions
@@ -284,9 +284,9 @@
-
- - Models +
+ + Models
  • @@ -300,7 +300,8 @@
-
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/concepts/error_stepper.html b/doc/boost_sandbox_numeric_odeint/concepts/error_stepper.html index 204902ea..6ccef802 100644 --- a/doc/boost_sandbox_numeric_odeint/concepts/error_stepper.html +++ b/doc/boost_sandbox_numeric_odeint/concepts/error_stepper.html @@ -24,9 +24,9 @@ This concepts specifies the interface an error stepper has to fulfill to be used within a ControlledErrorStepper.

-
- - Description +
+ + Description

An error stepper following this Error Stepper concept is capable of doing @@ -40,9 +40,9 @@ Symplectic System or Implicit System.

-
- - Notation +
+ + Notation

@@ -90,9 +90,9 @@

-
- - Valid +
+ + Valid Expressions
@@ -187,9 +187,9 @@
-
- - Models +
+ + Models
  • @@ -205,7 +205,8 @@
-
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/concepts/implicit_system.html b/doc/boost_sandbox_numeric_odeint/concepts/implicit_system.html index 4c70061f..497e6ac2 100644 --- a/doc/boost_sandbox_numeric_odeint/concepts/implicit_system.html +++ b/doc/boost_sandbox_numeric_odeint/concepts/implicit_system.html @@ -20,22 +20,22 @@ Implicit System
-
- - Description +
+ + Description

This concept describes how to define a ODE that can be solved by an implicit routine. Implicit routines need not only the function f(x,t) - but also the Jacobian /dfdx = A(x,t). A + but also the Jacobian df/dx = A(x,t). A is a matrix and implicit routines need to solve the linear problem Ax = b. In odeint this is implemented with use of Boost.UBlas, therefore, the state_type implicit routines is ublas::vector and the matrix is defined as ublas::matrix.

-
- - Notation +
+ + Notation

@@ -70,9 +70,9 @@

-
- - Valid +
+ + Valid Expressions
@@ -162,7 +162,8 @@
-
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/concepts/observer.html b/doc/boost_sandbox_numeric_odeint/concepts/observer.html index ca335835..8dc8c4d2 100644 --- a/doc/boost_sandbox_numeric_odeint/concepts/observer.html +++ b/doc/boost_sandbox_numeric_odeint/concepts/observer.html @@ -19,9 +19,9 @@
-
- - Description +
+ + Description

The Observer concept is used in integrate @@ -32,9 +32,9 @@ or at user defined time points depending on which integrate function is used.

-
- - Notation +
+ + Notation

@@ -65,9 +65,9 @@

-
- - Valid +
+ + Valid Expressions
@@ -128,7 +128,8 @@
-
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/concepts/simple_symplectic_system.html b/doc/boost_sandbox_numeric_odeint/concepts/simple_symplectic_system.html index 46f9c21b..d9556262 100644 --- a/doc/boost_sandbox_numeric_odeint/concepts/simple_symplectic_system.html +++ b/doc/boost_sandbox_numeric_odeint/concepts/simple_symplectic_system.html @@ -20,9 +20,9 @@ Simple Symplectic System
-
- - Description +
+ + Description

In most Hamiltonian systems the kinetic term is a quadratic term in the momentum @@ -44,9 +44,9 @@ equation for p' has to be provided to the symplectic stepper. We call this concept Simple_Symplectic_System

-
- - Notation +
+ + Notation

@@ -77,9 +77,9 @@

-
- - Valid +
+ + Valid Expressions
@@ -166,7 +166,8 @@
-
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/concepts/state_algebra_operations.html b/doc/boost_sandbox_numeric_odeint/concepts/state_algebra_operations.html index 3f581c47..9bdbb1c6 100644 --- a/doc/boost_sandbox_numeric_odeint/concepts/state_algebra_operations.html +++ b/doc/boost_sandbox_numeric_odeint/concepts/state_algebra_operations.html @@ -40,9 +40,9 @@ and no algebra/operations are used.

-
- - Description +
+ + Description

The State, Algebra and Operations @@ -82,9 +82,9 @@

-
- - Notation +
+ + Notation

@@ -111,7 +111,7 @@ Type that represents a general scale_sum operation, N should be replaced by a number from 1 to 14.

-
scale_sumN
+
scale_sum/N/

Object of type ScaleSumN, N should be replaced by a @@ -142,9 +142,9 @@

-
- - Valid +
+ + Valid Expressions
@@ -258,21 +258,18 @@ @@ -284,13 +281,13 @@

- Operations::scale_sumN< Value1 - , ... , ValueN > + Operations::scale_sum/N/< Value1 , ... , ValueN >

- ScaleSumN + ScaleSum/N/

- Get the ScaleSumN type - from Operations, - N should be replaced - by a number from 1 to 14. + Get the ScaleSum/N/ type from Operations, N + should be replaced by a number from 1 to 14.

- ScaleSumN< Value1 , ... , ValueN - >( a1 , ... , aN ) + ScaleSum/N/< Value1 , ... , ValueN >( a1 , ... + , aN )

- ScaleSumN + ScaleSum/N/

@@ -310,7 +307,7 @@

- scale_sumN( y , x1 , ... , xN ) + scale_sum/N/( y , x1 , ... , xN )

@@ -417,9 +414,9 @@
-
- - Notation +
+ + Notation

@@ -432,7 +429,7 @@

The algebra type

-
OperationN
+
Operation/N/

An N-ary operation type, N should be a number from @@ -442,7 +439,7 @@

Object of type Algebra

-
operationN
+
operation/N/

Object of type Operation/N

@@ -455,9 +452,9 @@

-
- - Valid +
+ + Valid Expressions
@@ -557,8 +554,8 @@

- algebra.for_eachN( y , x1 , ... - , xN , operationN ) + algebra.for_each/N/( y , x1 , ... , xN , operation/N/ + )

@@ -568,8 +565,8 @@

- Calls operationN( y_i , x1_i , - ... , xN_i ) for each element y_i + Calls operation/N/( y_i , x1_i , ... , xN_i ) + for each element y_i of y and x1_i of x1 and so on. N should be replaced by a number between 1 and 14. @@ -817,7 +814,8 @@ -
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/concepts/state_wrapper.html b/doc/boost_sandbox_numeric_odeint/concepts/state_wrapper.html index 8846787a..91821c5e 100644 --- a/doc/boost_sandbox_numeric_odeint/concepts/state_wrapper.html +++ b/doc/boost_sandbox_numeric_odeint/concepts/state_wrapper.html @@ -20,9 +20,9 @@ State Wrapper
-
- - Description +
+ + Description

The State Wrapper @@ -30,9 +30,9 @@ intermediate results within the stepper's do_step methods.

-
- - Notation +
+ + Notation

@@ -56,9 +56,9 @@

-
- - Valid +
+ + Valid Expressions
@@ -267,7 +267,8 @@
-
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/concepts/stepper.html b/doc/boost_sandbox_numeric_odeint/concepts/stepper.html index e274001f..7e1a48ea 100644 --- a/doc/boost_sandbox_numeric_odeint/concepts/stepper.html +++ b/doc/boost_sandbox_numeric_odeint/concepts/stepper.html @@ -24,9 +24,9 @@ be used within the integrate functions.

-
- - Description +
+ + Description

The basic stepper concept. A basic stepper following this Stepper concept @@ -40,9 +40,9 @@ Symplectic System or Implicit System. Note that all error steppers are also basic steppers.

-
- - Refinement +
+ + Refinement of
    @@ -53,9 +53,9 @@ CopyConstructable
-
- - Notation +
+ + Notation

@@ -95,9 +95,9 @@

-
- - Valid +
+ + Valid Expressions
@@ -187,9 +187,9 @@
-
- - Models +
+ + Models
  • @@ -211,7 +211,8 @@
-
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/concepts/symplectic_system.html b/doc/boost_sandbox_numeric_odeint/concepts/symplectic_system.html index e999880b..a0940292 100644 --- a/doc/boost_sandbox_numeric_odeint/concepts/symplectic_system.html +++ b/doc/boost_sandbox_numeric_odeint/concepts/symplectic_system.html @@ -20,9 +20,9 @@ Symplectic System
-
- - Description +
+ + Description

This concept describes how to define a symplectic system written with generalized @@ -58,9 +58,9 @@ std::pair of functions or a functors. Symplectic systems are used in symplectic steppers like symplectic_rkn_sb3a_mclachlan.

-
- - Notation +
+ + Notation

@@ -107,9 +107,9 @@

-
- - Valid +
+ + Valid expressions
@@ -220,7 +220,8 @@
-
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/concepts/system.html b/doc/boost_sandbox_numeric_odeint/concepts/system.html index b14bf5bb..a6598de0 100644 --- a/doc/boost_sandbox_numeric_odeint/concepts/system.html +++ b/doc/boost_sandbox_numeric_odeint/concepts/system.html @@ -19,9 +19,9 @@
-
- - Description +
+ + Description

The System concept models the algorithmic implementation of the rhs. of the @@ -34,9 +34,9 @@ System and Implicit System.

-
- - Notation +
+ + Notation

@@ -75,9 +75,9 @@

-
- - Valid +
+ + Valid expressions
@@ -138,7 +138,8 @@
-
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/extend_odeint.html b/doc/boost_sandbox_numeric_odeint/extend_odeint.html index 561f0252..b8db00cd 100644 --- a/doc/boost_sandbox_numeric_odeint/extend_odeint.html +++ b/doc/boost_sandbox_numeric_odeint/extend_odeint.html @@ -33,7 +33,8 @@
-
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/extend_odeint/adapt_your_own_operations.html b/doc/boost_sandbox_numeric_odeint/extend_odeint/adapt_your_own_operations.html index 36e2bd8a..d53b04b9 100644 --- a/doc/boost_sandbox_numeric_odeint/extend_odeint/adapt_your_own_operations.html +++ b/doc/boost_sandbox_numeric_odeint/extend_odeint/adapt_your_own_operations.html @@ -29,7 +29,8 @@
- @@ -1040,7 +1044,6 @@ so all that has to be done in this case is to declare resizability:

- 

typedef boost::numeric::ublas::vector< double > state_type;
 
@@ -1064,7 +1067,6 @@
             code shows the corresponding definitions:
           

 

-            
 

 
int main()
 {
@@ -1097,9 +1099,8 @@
             looks as follows:
           

 

-            
 

-
class point3D :     
+
class point3D :
     boost::additive1< point3D ,
     boost::additive2< point3D , double ,
     boost::multiplicative2< point3D , double > > >
@@ -1111,7 +1112,7 @@
     point3D()
         : x( 0.0 ) , y( 0.0 ) , z( 0.0 )
     { }
-    
+
     point3D( const double val )
         : x( val ) , y( val ) , z( val )
     { }
@@ -1150,10 +1151,9 @@
             function calculating the element-wise absolute value for the state type:
           

 

-            
 

-
// only required for steppers with error control
-point3D operator/( const point3D &p1 , const point3D &p2 )
+
// only required for steppers with error control
+point3D operator/( const point3D &p1 , const point3D &p2 )
 {
     return point3D( p1.x/p2.x , p1.y/p2.y , p1.z/p1.z );
 }
@@ -1170,11 +1170,10 @@
             implementing a reduction over the state type:
           

 

-            
 

 
namespace boost { namespace numeric { namespace odeint {
-// specialization of vector_space_reduce, only required for steppers with error control
-template<>
+// specialization of vector_space_reduce, only required for steppers with error control
+template<>
 struct vector_space_reduce< point3D >
 {
     template< class Value , class Op >
@@ -1186,7 +1185,7 @@
         return init;
     }
 };
-} } } 
+} } }
 

 

           

@@ -1198,7 +1197,6 @@
             perform the integration on a lorenz system by using the vector_space_algebra again:
           

 

-            
 

 
const double sigma = 10.0;
 const double R = 28.0;
@@ -1217,11 +1215,11 @@
 {
 
     point3D x( 10.0 , 5.0 , 5.0 );
-    //typedef runge_kutta4< point3D , double , point3D , 
-    //                      double , vector_space_algebra > stepper;
-    typedef runge_kutta_dopri5< point3D , double , point3D , 
+    //typedef runge_kutta4< point3D , double , point3D , 
+    //                      double , vector_space_algebra > stepper;
+    typedef runge_kutta_dopri5< point3D , double , point3D ,
                                 double , vector_space_algebra > stepper;
-    integrate_adaptive( make_controlled<stepper>( 1E-8 , 1E-8 ) , lorenz , x , 
+    integrate_adaptive( make_controlled<stepper>( 1E-8 , 1E-8 ) , lorenz , x ,
                         0.0 , 10.0 , 0.1 );
     std::cout << x << std::endl;
 }
@@ -1236,7 +1234,8 @@
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/extend_odeint/adapt_your_own_state_types.html b/doc/boost_sandbox_numeric_odeint/extend_odeint/adapt_your_own_state_types.html index 8be7929f..e2c64286 100644 --- a/doc/boost_sandbox_numeric_odeint/extend_odeint/adapt_your_own_state_types.html +++ b/doc/boost_sandbox_numeric_odeint/extend_odeint/adapt_your_own_state_types.html @@ -381,7 +381,6 @@ as resizeable as it is a dynamically sized array:

- 

class my_vector : public std::vector< double >
 {
@@ -394,11 +393,11 @@
         : std::vector< double >()
     { }
 
-// ...
-};
+// ...
+};
+
+// define my_vector as reizeable
 
-// define my_vector as reizeable
-
 namespace boost { namespace numeric { namespace odeint {
 
 template<>
@@ -515,7 +514,6 @@
             to work with odeint. The following code shows the required template specializations:
           

 

-            
 

 
typedef std::list< double > state_type;
 
@@ -523,15 +521,15 @@
 
 template< >
 struct is_resizeable< state_type >
-{ // declare resizeablility
-    typedef boost::true_type type;
+{ // declare resizeablility
+    typedef boost::true_type type;
     const static bool value = type::value;
 };
 
-template< > 
+template< >
 struct same_size_impl< state_type , state_type >
-{ // define how to check size
-    static bool same_size( const state_type &v1 ,
+{ // define how to check size
+    static bool same_size( const state_type &v1 ,
                            const state_type &v2 )
     {
         return v1.size() == v2.size();
@@ -540,15 +538,15 @@
 
 template< >
 struct resize_impl< state_type , state_type >
-{ // define how to resize
-    static void resize( state_type &v1 ,
+{ // define how to resize
+    static void resize( state_type &v1 ,
                         const state_type &v2 )
     {
         v1.resize( v2.size() );
     }
 };
 
-} } } 
+} } }
 

 

           

@@ -617,7 +615,6 @@
             instances of state types:
           

 

-            
 

 
template<>
 struct state_wrapper< gsl_vector* >
@@ -654,7 +651,7 @@
             This state_wrapper specialization
             tells odeint how gsl_vectors are created, copied and destroyed. Next
             we need resizing, this is required because gsl_vectors are dynamically
-            sized objects: 
+            sized objects:
 

 
template<>
 struct is_resizeable< gsl_vector* >
@@ -689,13 +686,12 @@
             support iterators, so we first implement a gsl iterator:
           

 

-            
 

 
/*
  * defines an iterator for gsl_vector
  */
-class gsl_vector_iterator 
-      : public boost::iterator_facade< gsl_vector_iterator , double , 
+class gsl_vector_iterator
+      : public boost::iterator_facade< gsl_vector_iterator , double ,
                                        boost::random_access_traversal_tag >
 {
 public :
@@ -724,7 +720,6 @@
             A similar class exists for the const
             version of the iterator. Then we have a function returning the end iterator
             (similarily for const again):
-            
 

 
gsl_vector_iterator end_iterator( gsl_vector *x )
 {
@@ -737,16 +732,16 @@
           

 

             Finally, the bindings for Boost.Range
-            are added: 
+            are added:
 

-
// template<>
-inline gsl_vector_iterator range_begin( gsl_vector *x )
+
// template<>
+inline gsl_vector_iterator range_begin( gsl_vector *x )
 {
     return gsl_vector_iterator( x );
 }
 
-// template<>
-inline gsl_vector_iterator range_end( gsl_vector *x )
+// template<>
+inline gsl_vector_iterator range_end( gsl_vector *x )
 {
     return end_iterator( x );
 }
@@ -1012,15 +1007,24 @@
                     

                       Performs the operation operator
                       for subsequently each element of state
-                      and returns the aggregate value. E.g. 
 init =
+                      and returns the aggregate value. E.g.
+                    

+                    

+                      init =
                       operator(
                       init ,
                       state[0]
-                      ); 
 init =
+                      );
+                    

+                    

+                      init =
                       operator(
                       init ,
                       state[1]
-                      ) 
 ...
+                      )
+                    

+                    

+                      ...
-
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/extend_odeint/write_own_steppers.html b/doc/boost_sandbox_numeric_odeint/extend_odeint/write_own_steppers.html index b2fec6f7..e3803f1d 100644 --- a/doc/boost_sandbox_numeric_odeint/extend_odeint/write_own_steppers.html +++ b/doc/boost_sandbox_numeric_odeint/extend_odeint/write_own_steppers.html @@ -21,7 +21,8 @@
-
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/getting_started.html b/doc/boost_sandbox_numeric_odeint/getting_started.html index a92d20fa..7c940689 100644 --- a/doc/boost_sandbox_numeric_odeint/getting_started.html +++ b/doc/boost_sandbox_numeric_odeint/getting_started.html @@ -29,7 +29,8 @@
-
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/getting_started/overview.html b/doc/boost_sandbox_numeric_odeint/getting_started/overview.html index 1f1578e5..4dbabbd9 100644 --- a/doc/boost_sandbox_numeric_odeint/getting_started/overview.html +++ b/doc/boost_sandbox_numeric_odeint/getting_started/overview.html @@ -1095,7 +1095,8 @@
-
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/getting_started/short_example.html b/doc/boost_sandbox_numeric_odeint/getting_started/short_example.html index b79da97e..922c3eb8 100644 --- a/doc/boost_sandbox_numeric_odeint/getting_started/short_example.html +++ b/doc/boost_sandbox_numeric_odeint/getting_started/short_example.html @@ -30,7 +30,6 @@ = (p,-q+γ p):

- 

/* The type of container used to hold the state vector */
 typedef std::vector< double > state_type;
@@ -60,11 +59,10 @@
         start:
       

 

-        
 

 
state_type x(2);
-x[0] = 1.0; // start at x=1.0, p=0.0
-x[1] = 0.0;
+x[0] = 1.0; // start at x=1.0, p=0.0
+x[1] = 0.0;
 

 

       

@@ -75,7 +73,6 @@
         stepper (5th order) and uses adaptive stepsize.
       

 

-        
 

 
size_t steps = integrate( harmonic_oscillator ,
         x , 0.0 , 10.0 , 0.1 );
@@ -96,7 +93,6 @@
         rhs must then be implemented as a functor having defined the ()-operator:
       

 

-        
 

 
/* The rhs of x' = f(x) defined as a class */
 class harm_osc {
@@ -119,7 +115,6 @@
         which can be used via
       

 

-        
 

 
harm_osc ho(0.15);
 steps = integrate( ho ,
@@ -134,7 +129,6 @@
         to provide a reasonable observer. An example is
       

 

-        
 

 
struct push_back_state_and_time
 {
@@ -158,7 +152,6 @@
         pass this container to the integration function:
       

 

-        
 

 
vector<state_type> x_vec;
 vector<double> times;
@@ -186,7 +179,8 @@
-
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/getting_started/usage__compilation__headers.html b/doc/boost_sandbox_numeric_odeint/getting_started/usage__compilation__headers.html index 3cf5f598..b685263e 100644 --- a/doc/boost_sandbox_numeric_odeint/getting_started/usage__compilation__headers.html +++ b/doc/boost_sandbox_numeric_odeint/getting_started/usage__compilation__headers.html @@ -25,13 +25,12 @@ library. It can be include by

- 

#include <boost/numeric/odeint.hpp>

which includes all headers in the library. All functions and classes from - odeint live in the namespace + odeint live in the namespace 

using namespace boost::numeric::odeint;
@@ -40,7 +39,8 @@
-
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/odeint_in_detail.html b/doc/boost_sandbox_numeric_odeint/odeint_in_detail.html index 4d673c88..ba34acd0 100644 --- a/doc/boost_sandbox_numeric_odeint/odeint_in_detail.html +++ b/doc/boost_sandbox_numeric_odeint/odeint_in_detail.html @@ -36,7 +36,8 @@
-
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/odeint_in_detail/algebras_and_operations.html b/doc/boost_sandbox_numeric_odeint/odeint_in_detail/algebras_and_operations.html index 4014b28e..6654a535 100644 --- a/doc/boost_sandbox_numeric_odeint/odeint_in_detail/algebras_and_operations.html +++ b/doc/boost_sandbox_numeric_odeint/odeint_in_detail/algebras_and_operations.html @@ -26,7 +26,8 @@
-
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/odeint_in_detail/generation_functions.html b/doc/boost_sandbox_numeric_odeint/odeint_in_detail/generation_functions.html index 06a0263d..66a743a1 100644 --- a/doc/boost_sandbox_numeric_odeint/odeint_in_detail/generation_functions.html +++ b/doc/boost_sandbox_numeric_odeint/odeint_in_detail/generation_functions.html @@ -202,7 +202,8 @@
-
Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

+

Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

diff --git a/doc/boost_sandbox_numeric_odeint/odeint_in_detail/integrate_functions.html b/doc/boost_sandbox_numeric_odeint/odeint_in_detail/integrate_functions.html index e626a582..4ee57d8c 100644 --- a/doc/boost_sandbox_numeric_odeint/odeint_in_detail/integrate_functions.html +++ b/doc/boost_sandbox_numeric_odeint/odeint_in_detail/integrate_functions.html @@ -35,9 +35,9 @@ Output Stepper. Depending on the abilities of the stepper, the integrate functions make use of step-size control or dense output.

-
- - Equidistant +
+ + Equidistant observer calls

@@ -101,16 +101,16 @@ is considerably larger than typical step sizes used by the stepper.

  • - If stepper is a __dense_out_stepper + If stepper is a __dense_outout_stepper then dt is the initial step size. The actual step size will be adjusted during integration due to error control. If an observer is provided dense output is used to calculate x(t) at t = t0 + n dt.
    • -
    - - Observer +
    + + Observer calls at each step

    @@ -170,15 +170,16 @@ the first step at t0).

  • - If stepper is a __dense_out_stepper then dt + If stepper is a Dense + Output Stepper then dt is the initial step size and integrate_adaptive behaves just like for Controlled Stepper above. No dense output is used.
    • -
    - - Observer +
    + + Observer calls at given time points

    @@ -251,9 +252,9 @@ at the time points from the sequence. -

    - - Integrate +
    + + Integrate a given number of steps

    @@ -308,9 +309,9 @@ This function returns the end time t0 + n*dt.

    -
    - - Convenience +
    + + Convenience integrate function

    @@ -344,7 +345,8 @@ -
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

    +

    Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

    Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

    diff --git a/doc/boost_sandbox_numeric_odeint/odeint_in_detail/steppers.html b/doc/boost_sandbox_numeric_odeint/odeint_in_detail/steppers.html index c2520504..5b6f8d7c 100644 --- a/doc/boost_sandbox_numeric_odeint/odeint_in_detail/steppers.html +++ b/doc/boost_sandbox_numeric_odeint/odeint_in_detail/steppers.html @@ -130,15 +130,14 @@ is the Runge Kutta stepper of fourth order:

    - 

    runge_kutta4< state_type > rk;
-rk.do_step( sys1 , inout , t , dt );               // In-place transformation of inout
-rk.do_step( sys2 , inout , t , dt );               // Ok
-rk.do_step( sys1 , in , t , out , dt );            // Out-of-place transformation
-rk.do_step( sys1 , inout , dxdtin , t , dt );      // In-place tranformation of inout
-rk.do_step( sys1 , in , dxdtin , t , out , dt );   // Out-of-place transformation
-
    +rk.do_step( sys1 , inout , t , dt ); // In-place transformation of inout +rk.do_step( sys2 , inout , t , dt ); // Ok +rk.do_step( sys1 , in , t , out , dt ); // Out-of-place transformation +rk.do_step( sys1 , inout , dxdtin , t , dt ); // In-place tranformation of inout +rk.do_step( sys1 , in , dxdtin , t , out , dt ); // Out-of-place transformation +

    @@ -175,21 +174,20 @@ how the FSAL steppers can be used is

    - 

    runge_kutta_dopri5< state_type > rk;
 rk.do_step( sys1 , in , t , out , dt );
-rk.do_step( sys2 , in , t , out , dt );         // DONT do this, sys1 is assumed
-
+rk.do_step( sys2 , in , t , out , dt );         // DONT do this, sys1 is assumed
+
 rk.do_step( sys2 , in2 , t , out , dt );
-rk.do_step( sys2 , in3 , t , out , dt );        // DONT do this, in2 is assumed
-
+rk.do_step( sys2 , in3 , t , out , dt );        // DONT do this, in2 is assumed
+
 rk.do_step( sys1 , inout , dxdtinout , t , dt );
-rk.do_step( sys2 , inout , dxdtinout , t , dt );           // Ok, internal derivative is not used, dxdtinout is updated
-
+rk.do_step( sys2 , inout , dxdtinout , t , dt );           // Ok, internal derivative is not used, dxdtinout is updated
+
 rk.do_step( sys1 , in , dxdtin , t , out , dxdtout , dt );
-rk.do_step( sys2 , in , dxdtin , t , out , dxdtout , dt ); // Ok, internal derivative is not used
-
    +rk.do_step( sys2 , in , dxdtin , t , out , dxdtout , dt ); // Ok, internal derivative is not used +

    @@ -241,7 +239,6 @@ again simple C-functions of functors. An example is the harmonic oscillator:

    - 

    typedef boost::array< double , 1 > vector_type;
 
@@ -263,7 +260,6 @@
           example for the harmonic oscillator is now:
         
    
 
    
-          
 
    
 
    pair< vector_type , vector_type > x;
 x.first[0] = 1.0; x.second[0] = 0.0;
@@ -277,7 +273,6 @@
           two public method:
         
    
 
    
-          
 
    
 
    struct harm_osc
 {
@@ -295,7 +290,6 @@
 
    
         
    
 
    
-          
 
    
 
    harm_osc h;
 rkn.do_step( make_pair( boost::bind( &harm_osc::f1 , h , _1 , _2 ) , boost::bind( &harm_osc::f2 , h , _1 , _2 ) ) ,
@@ -315,7 +309,6 @@
           can be written as
         
    
 
    
-          
 
    
 
    pair< vector_type , vector_type > x;
 x.first[0] = 1.0; x.second[0] = 0.0;
@@ -339,7 +332,6 @@
           takes q and p without combining them into a pair:
         
    
 
    
-          
 
    
 
    rkn.do_step( harm_osc_f1 , make_pair( boost::ref( q ) , boost::ref( p ) ) , t , dt );
 rkn.do_step( harm_osc_f1 , q , p , t , dt );
@@ -405,7 +397,6 @@
           steps is instantiated and initialized;
         
    
 
    
-          
 
    
 
    adams_bashforth_moulton< 5 , state_type > abm;
 abm.initialize( sys , inout , t , dt );
@@ -423,7 +414,6 @@
           state of the Adams-Bashforth-Stepper:
         
    
 
    
-          
 
    
 
    abm.initialize( runge_kutta_fehlberg78< state_type >() , sys , inout , t , dt );
 
    
@@ -686,7 +676,6 @@
           An example is
         
    
 
    
-          
 
    
 
    dense_output_runge_kutta< controlled_runge_kutta< runge_kutta_dopri5< state_type > > > dense;
 dense.initialize( in , t , dt );
@@ -712,7 +701,6 @@
           exist:
         
    
 
    
-          
 
    
 
    result_of::make_dense_output< runge_kutta_dopri5< state_type > >::type dense2 = make_dense_output( 1.0e-6 , 1.0e-6 , runge_kutta_dopri5< state_type >() );
 
    
@@ -730,7 +718,6 @@
           for easy use with the integrate functions:
         
    
 
    
-          
 
    
 
    integrate_const( make_dense_output( 1.0e-6 , 1.0e-6 , runge_kutta_dopri5< state_type >() ) , sys , inout , t_start , t_end , dt );
 
    
@@ -799,15 +786,14 @@
           4 internal derivatives of the solution at times (t-dt,t-2dt,t-3dt,t-4dt).
         
    
 
    
-          
 
    
 
    adams_bashforth_moulton< 4 , state_type > stepper;
-stepper.do_step( sys , x , t , dt );   // make one step with the classical Runge-Kutta stepper and initialize the first internal state at time t-4dt
-stepper.do_step( sys , x , t , dt );   // make one step with the classical Runge-Kutta stepper and initialize the first internal state at time t-3dt
-stepper.do_step( sys , x , t , dt );   // make one step with the classical Runge-Kutta stepper and initialize the first internal state at time t-2dt
-stepper.do_step( sys , x , t , dt );   // make one step with the classical Runge-Kutta stepper and initialize the first internal state at time t-dt
-stepper.do_step( sys , x , t , dt );   // make one step with Adam-Bashforth-Moulton, the internal array of states is now rotated
-
    
+stepper.do_step( sys , x , t , dt );   // make one step with the classical Runge-Kutta stepper and initialize the first internal state at time t-4dt
+stepper.do_step( sys , x , t , dt );   // make one step with the classical Runge-Kutta stepper and initialize the first internal state at time t-3dt
+stepper.do_step( sys , x , t , dt );   // make one step with the classical Runge-Kutta stepper and initialize the first internal state at time t-2dt
+stepper.do_step( sys , x , t , dt );   // make one step with the classical Runge-Kutta stepper and initialize the first internal state at time t-dt
+stepper.do_step( sys , x , t , dt );   // make one step with Adam-Bashforth-Moulton, the internal array of states is now rotated
+
    
 
    
         
    
 
    
@@ -1900,7 +1886,8 @@
    -
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

    +

    Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

    Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

    diff --git a/doc/boost_sandbox_numeric_odeint/odeint_in_detail/using_boost__range.html b/doc/boost_sandbox_numeric_odeint/odeint_in_detail/using_boost__range.html index 657ea5af..25f9c671 100644 --- a/doc/boost_sandbox_numeric_odeint/odeint_in_detail/using_boost__range.html +++ b/doc/boost_sandbox_numeric_odeint/odeint_in_detail/using_boost__range.html @@ -21,7 +21,8 @@
    -
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

    +

    Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

    Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

    diff --git a/doc/boost_sandbox_numeric_odeint/odeint_in_detail/using_boost__ref.html b/doc/boost_sandbox_numeric_odeint/odeint_in_detail/using_boost__ref.html index e93491f2..b4b2521c 100644 --- a/doc/boost_sandbox_numeric_odeint/odeint_in_detail/using_boost__ref.html +++ b/doc/boost_sandbox_numeric_odeint/odeint_in_detail/using_boost__ref.html @@ -21,7 +21,8 @@
    -
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

    +

    Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

    Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

    diff --git a/doc/boost_sandbox_numeric_odeint/old_concepts.html b/doc/boost_sandbox_numeric_odeint/old_concepts.html index 0e3192e9..2f0e8e3e 100644 --- a/doc/boost_sandbox_numeric_odeint/old_concepts.html +++ b/doc/boost_sandbox_numeric_odeint/old_concepts.html @@ -41,7 +41,8 @@
    -
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

    +

    Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

    Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

    diff --git a/doc/boost_sandbox_numeric_odeint/old_concepts/basic_stepper.html b/doc/boost_sandbox_numeric_odeint/old_concepts/basic_stepper.html index 8fb795ec..1b03d138 100644 --- a/doc/boost_sandbox_numeric_odeint/old_concepts/basic_stepper.html +++ b/doc/boost_sandbox_numeric_odeint/old_concepts/basic_stepper.html @@ -294,7 +294,8 @@
    -
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

    +

    Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

    Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

    diff --git a/doc/boost_sandbox_numeric_odeint/old_concepts/compositestepper.html b/doc/boost_sandbox_numeric_odeint/old_concepts/compositestepper.html index 638b9bb4..38755a13 100644 --- a/doc/boost_sandbox_numeric_odeint/old_concepts/compositestepper.html +++ b/doc/boost_sandbox_numeric_odeint/old_concepts/compositestepper.html @@ -20,7 +20,8 @@
    -
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

    +

    Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

    Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

    diff --git a/doc/boost_sandbox_numeric_odeint/old_concepts/controlled_stepper.html b/doc/boost_sandbox_numeric_odeint/old_concepts/controlled_stepper.html index 3bb7951e..3d8d5c0a 100644 --- a/doc/boost_sandbox_numeric_odeint/old_concepts/controlled_stepper.html +++ b/doc/boost_sandbox_numeric_odeint/old_concepts/controlled_stepper.html @@ -223,7 +223,8 @@
    -
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

    +

    Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

    Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

    diff --git a/doc/boost_sandbox_numeric_odeint/old_concepts/dense_ouput_stepper.html b/doc/boost_sandbox_numeric_odeint/old_concepts/dense_ouput_stepper.html index 125b75c4..6afec565 100644 --- a/doc/boost_sandbox_numeric_odeint/old_concepts/dense_ouput_stepper.html +++ b/doc/boost_sandbox_numeric_odeint/old_concepts/dense_ouput_stepper.html @@ -21,7 +21,8 @@
    -
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

    +

    Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

    Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

    diff --git a/doc/boost_sandbox_numeric_odeint/old_concepts/error_stepper.html b/doc/boost_sandbox_numeric_odeint/old_concepts/error_stepper.html index bce3b512..95fd84c4 100644 --- a/doc/boost_sandbox_numeric_odeint/old_concepts/error_stepper.html +++ b/doc/boost_sandbox_numeric_odeint/old_concepts/error_stepper.html @@ -232,7 +232,8 @@
    -
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

    +

    Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

    Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

    diff --git a/doc/boost_sandbox_numeric_odeint/old_concepts/size_adjusting_stepper.html b/doc/boost_sandbox_numeric_odeint/old_concepts/size_adjusting_stepper.html index c54bc218..b227fd1b 100644 --- a/doc/boost_sandbox_numeric_odeint/old_concepts/size_adjusting_stepper.html +++ b/doc/boost_sandbox_numeric_odeint/old_concepts/size_adjusting_stepper.html @@ -21,7 +21,8 @@
    -
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

    +

    Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

    Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

    diff --git a/doc/boost_sandbox_numeric_odeint/old_reference.html b/doc/boost_sandbox_numeric_odeint/old_reference.html index bc0874da..b95c5cc8 100644 --- a/doc/boost_sandbox_numeric_odeint/old_reference.html +++ b/doc/boost_sandbox_numeric_odeint/old_reference.html @@ -31,7 +31,8 @@
    -
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

    +

    Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

    Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

    diff --git a/doc/boost_sandbox_numeric_odeint/old_reference/algebras.html b/doc/boost_sandbox_numeric_odeint/old_reference/algebras.html index 33abc8dc..89cc2426 100644 --- a/doc/boost_sandbox_numeric_odeint/old_reference/algebras.html +++ b/doc/boost_sandbox_numeric_odeint/old_reference/algebras.html @@ -20,7 +20,8 @@
    -
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

    +

    Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

    Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

    diff --git a/doc/boost_sandbox_numeric_odeint/old_reference/integration_functions.html b/doc/boost_sandbox_numeric_odeint/old_reference/integration_functions.html index 6c4b3705..249ad967 100644 --- a/doc/boost_sandbox_numeric_odeint/old_reference/integration_functions.html +++ b/doc/boost_sandbox_numeric_odeint/old_reference/integration_functions.html @@ -37,7 +37,8 @@
    -
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

    +

    Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

    Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

    diff --git a/doc/boost_sandbox_numeric_odeint/old_reference/operations.html b/doc/boost_sandbox_numeric_odeint/old_reference/operations.html index d9d8d5ca..7b71d6ff 100644 --- a/doc/boost_sandbox_numeric_odeint/old_reference/operations.html +++ b/doc/boost_sandbox_numeric_odeint/old_reference/operations.html @@ -20,7 +20,8 @@
    -
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

    +

    Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

    Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

    diff --git a/doc/boost_sandbox_numeric_odeint/old_reference/resizing.html b/doc/boost_sandbox_numeric_odeint/old_reference/resizing.html index 9a88f763..61681d77 100644 --- a/doc/boost_sandbox_numeric_odeint/old_reference/resizing.html +++ b/doc/boost_sandbox_numeric_odeint/old_reference/resizing.html @@ -19,7 +19,8 @@
    -
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

    +

    Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

    Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

    diff --git a/doc/boost_sandbox_numeric_odeint/old_reference/stepper_classes.html b/doc/boost_sandbox_numeric_odeint/old_reference/stepper_classes.html index 0775a70a..04ebacf7 100644 --- a/doc/boost_sandbox_numeric_odeint/old_reference/stepper_classes.html +++ b/doc/boost_sandbox_numeric_odeint/old_reference/stepper_classes.html @@ -191,7 +191,8 @@
    -
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

    +

    Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

    Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

    diff --git a/doc/boost_sandbox_numeric_odeint/tutorial.html b/doc/boost_sandbox_numeric_odeint/tutorial.html index 0ea093c8..7709f593 100644 --- a/doc/boost_sandbox_numeric_odeint/tutorial.html +++ b/doc/boost_sandbox_numeric_odeint/tutorial.html @@ -39,7 +39,8 @@
    -
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

    +

    Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

    Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

    diff --git a/doc/boost_sandbox_numeric_odeint/tutorial/all_examples.html b/doc/boost_sandbox_numeric_odeint/tutorial/all_examples.html index 4ab2eadd..9e6b8072 100644 --- a/doc/boost_sandbox_numeric_odeint/tutorial/all_examples.html +++ b/doc/boost_sandbox_numeric_odeint/tutorial/all_examples.html @@ -236,7 +236,8 @@
    -
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky

    +

    Copyright © 2009-2011 Karsten Ahnert + and Mario Mulansky

    Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

    diff --git a/doc/boost_sandbox_numeric_odeint/tutorial/chaotic_systems_and_lyapunov_exponents.html b/doc/boost_sandbox_numeric_odeint/tutorial/chaotic_systems_and_lyapunov_exponents.html index 8494961e..ec030360 100644 --- a/doc/boost_sandbox_numeric_odeint/tutorial/chaotic_systems_and_lyapunov_exponents.html +++ b/doc/boost_sandbox_numeric_odeint/tutorial/chaotic_systems_and_lyapunov_exponents.html @@ -61,7 +61,6 @@ The implementation of the Lorenz system is

    - 

    const double sigma = 10.0;
 const double R = 28.0;
@@ -84,7 +83,6 @@
         the perturbation could look like:
       
    
 
    
-        
 
    
 
    const size_t n = 3;
 const size_t num_of_lyap = 3;
@@ -116,7 +114,6 @@
         inside a functor with templatized arguments:
       
    
 
    
-        
 
    
 
    struct lorenz
 {
@@ -138,11 +135,11 @@
 
    
 
    
         This works fine and lorenz_with_lyap
-        can be used for example via 
+        can be used for example via
 
    
 
    state_type x;
-// initialize x
-explicit_rk4< state_type > rk4;
+// initialize x
+explicit_rk4< state_type > rk4;
 integrate_n_steps( rk4 , lorenz_with_lyap , x , 0.0 , 0.01 , 1000 );
 
    
 
    
@@ -154,11 +151,10 @@
         current state lies on the attractor, hence it would look like
       
    
 
    
-        
 
    
 
    state_type x;
-// initialize x
-explicit_rk4< state_type > rk4;
+// initialize x
+explicit_rk4< state_type > rk4;
 integrate_n_steps( rk4 , lorenz , x , 0.0 , 0.01 , 1000 );
 
    
 
    
@@ -170,7 +166,6 @@
         by pointer to the state and an appropriate length, something similar to
       
    
 
    
-        
 
    
 
    void lorenz( double* x , double *dxdt , double t, void* params )
 {
@@ -188,7 +183,6 @@
         the system has to by changed:
       
    
 
    
-        
 
    
 
    struct lorenz
 {
@@ -211,10 +205,9 @@
         with an suited range:
       
    
 
    
-        
 
    
-
    // perform 10000 transient steps
-integrate_n_steps( rk4 , lorenz() , std::make_pair( x.begin() , x.begin() + n ) , 0.0 , dt , 10000 );
+
    // perform 10000 transient steps
+integrate_n_steps( rk4 , lorenz() , std::make_pair( x.begin() , x.begin() + n ) , 0.0 , dt , 10000 );
 
    
 
    
       
    
@@ -242,7 +235,6 @@

    - 

    fill( x.begin()+n , x.end() , 0.0 );
 for( size_t i=0 ; i<num_of_lyap ; ++i ) x[n+n*i+i] = 1.0;
@@ -273,7 +265,8 @@
 
 
-
    
 
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky
    
+
    Copyright © 2009-2011 Karsten Ahnert
+      and Mario Mulansky
    
         Distributed under the Boost Software License, Version 1.0. (See accompanying
         file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
       
    
diff --git a/doc/boost_sandbox_numeric_odeint/tutorial/harmonic_oscillator.html b/doc/boost_sandbox_numeric_odeint/tutorial/harmonic_oscillator.html
index 26036355..33a6ebf6 100644
--- a/doc/boost_sandbox_numeric_odeint/tutorial/harmonic_oscillator.html
+++ b/doc/boost_sandbox_numeric_odeint/tutorial/harmonic_oscillator.html
@@ -52,7 +52,6 @@
           to just define a function, e.g:
         
    
 
    
-          
 
    
 
    /* The type of container used to hold the state vector */
 typedef std::vector< double > state_type;
@@ -79,7 +78,6 @@
           parameter structure as above:
         
    
 
    
-          
 
    
 
    /* The rhs of x' = f(x) defined as a class */
 class harm_osc {
@@ -1138,7 +1136,6 @@
           ) function from odeint:
         
    
 
    
-          
 
    
 
    runge_kutta4< state_type > stepper;
 integrate_const( stepper , harmonic_oscillator , x , 0.0 , 10.0 , 0.01 );
@@ -1155,7 +1152,6 @@
           method which can used directly. So, you write down the above example as
         
    
 
    
-          
 
    
 
    const double dt = 0.01;
 for( double t=0.0 ; t<10.0 ; t+= dt )
@@ -1179,7 +1175,6 @@
           with 4th order error estimation and coefficients introduced by Cash-Karp.
         
    
 
    
-          
 
    
 
    typedef runge_kutta_cash_karp54< state_type > error_stepper_type;
 
    
@@ -1192,7 +1187,6 @@
           stepper exists which can be used via
         
    
 
    
-          
 
    
 
    typedef controlled_runge_kutta< error_stepper_type > controlled_stepper_type;
 controlled_stepper_type controlled_stepper;
@@ -1219,7 +1213,6 @@
           by using make_controlled:
         
    
 
    
-          
 
    
 
    integrate_adaptive( make_controlled< error_stepper_type >( 1.0e-10 , 1.0e-6 ) , harmonic_oscillator , x , 0.0 , 10.0 , 0.01 );
 
    
@@ -1235,7 +1228,6 @@
           is
         
    
 
    
-          
 
    
 
    integrate_adaptive( make_controlled( 1.0e-10 , 1.0e-6 , error_stepper_type() ) , harmonic_oscillator , x , 0.0 , 10.0 , 0.01 );
 
    
@@ -1251,7 +1243,6 @@
           with all parameters is therefore
         
    
 
    
-          
 
    
 
    double abs_err = 1.0e-10 , rel_err = 1.0e-6 , a_x = 1.0 , a_dxdt = 1.0;
 controlled_stepper_type controlled_stepper( default_error_checker< double >( abs_err , rel_err , a_x , a_dxdt ) );
@@ -1453,7 +1444,8 @@
 
    
 
-
    
 
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky
    
+
    Copyright © 2009-2011 Karsten Ahnert
+      and Mario Mulansky
    
         Distributed under the Boost Software License, Version 1.0. (See accompanying
         file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
       
    
diff --git a/doc/boost_sandbox_numeric_odeint/tutorial/references.html b/doc/boost_sandbox_numeric_odeint/tutorial/references.html
index 03ce8a67..d12d0ede 100644
--- a/doc/boost_sandbox_numeric_odeint/tutorial/references.html
+++ b/doc/boost_sandbox_numeric_odeint/tutorial/references.html
@@ -74,7 +74,8 @@
 
    
 
-
    
 
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky
    
+
    Copyright © 2009-2011 Karsten Ahnert
+      and Mario Mulansky
    
         Distributed under the Boost Software License, Version 1.0. (See accompanying
         file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
       
    
diff --git a/doc/boost_sandbox_numeric_odeint/tutorial/solar_system.html b/doc/boost_sandbox_numeric_odeint/tutorial/solar_system.html
index d7e89781..b51df30a 100644
--- a/doc/boost_sandbox_numeric_odeint/tutorial/solar_system.html
+++ b/doc/boost_sandbox_numeric_odeint/tutorial/solar_system.html
@@ -98,7 +98,6 @@
           space as well as the velocity. Therefore, we use the operators from Boost.Operators:
         
    
 
    
-          
 
    
 
    /*the point type */
 template< class T , size_t Dim >
@@ -114,20 +113,20 @@
         typedef T value_type;
         typedef point< value_type , dim > point_type;
 
-        // ...
-        // constructors
-
-        // ...
-        // operators
-
+        // ...
+        // constructors
+
+        // ...
+        // operators
+
     private:
 
         T m_val[dim];
     };
 
-    //...
-    // more operators
-    
    
+    //...
+    // more operators
+    
 
    
         
    
 
    
@@ -136,10 +135,9 @@
           type we use std::tr1::array
         
    
 
    
-          
 
    
-
    // we simulate 5 planets and the sun
-const size_t n = 6;
+
    // we simulate 5 planets and the sun
+const size_t n = 6;
 
 typedef point< double , 3 > point_type;
 typedef boost::array< point_type , n > container_type;
@@ -159,7 +157,6 @@
           As system function we have to provide f(p) and f(q):
         
    
 
    
-          
 
    
 
    const double gravitational_constant = 2.95912208286e-4;
 
@@ -179,7 +176,6 @@
 
    
         
    
 
    
-          
 
    
 
    struct solar_system_momentum
 {
@@ -221,7 +217,6 @@
           apply here:
         
    
 
    
-          
 
    
 
    typedef symplectic_rkn_sb3a_mclachlan< container_type > stepper_type;
 const double dt = 100.0;
@@ -246,7 +241,6 @@
           The reference wrapper is also passed, but this is not a problem at all:
         
    
 
    
-          
 
    
 
    struct streaming_observer
 {
@@ -273,7 +267,8 @@
 
    
 
-
    
 
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky
    
+
    Copyright © 2009-2011 Karsten Ahnert
+      and Mario Mulansky
    
         Distributed under the Boost Software License, Version 1.0. (See accompanying
         file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
       
    
diff --git a/doc/boost_sandbox_numeric_odeint/tutorial/special_topics.html b/doc/boost_sandbox_numeric_odeint/tutorial/special_topics.html
index 2ec48199..a4934780 100644
--- a/doc/boost_sandbox_numeric_odeint/tutorial/special_topics.html
+++ b/doc/boost_sandbox_numeric_odeint/tutorial/special_topics.html
@@ -63,7 +63,6 @@
           code is very simple
         
    
 
    
-          
 
    
 
    typedef boost::array< complex< double > , 1 > state_type;
 
@@ -89,7 +88,6 @@
           function. In this cast the Stuart-Landau oscillator looks like
         
    
 
    
-          
 
    
 
    double eta = 1.0;
 double alpha = 1.0;
@@ -111,7 +109,6 @@
           Integration is also very easy:
         
    
 
    
-          
 
    
 
    state_type x = {{ complex< double >( 1.0 , 0.0 ) }};
 
@@ -156,7 +153,6 @@
           components. Here, is how this function looks like
         
    
 
    
-          
 
    
 
    typedef vector< double > container_type;
 
@@ -166,8 +162,8 @@
 
     fpu( const double beta = 1.0 ) : m_beta( beta ) { }
 
-    // system function defining the ODE
-    void operator()( const container_type &q , container_type &dpdt ) const
+    // system function defining the ODE
+    void operator()( const container_type &q , container_type &dpdt ) const
     {
         size_t n = q.size();
         double tmp = q[0] - 0.0;
@@ -185,17 +181,17 @@
         dpdt[n-1] += tmp2;
     }
 
-    // calculates the energy of the system
-    double energy( const container_type &q , const container_type &p ) const
+    // calculates the energy of the system
+    double energy( const container_type &q , const container_type &p ) const
     {
-        // ...
-    }
+        // ...
+    }
 
-    // calculates the local energy of the system
-    void local_energy( const container_type &q , const container_type &p , container_type &e ) const
+    // calculates the local energy of the system
+    void local_energy( const container_type &q , const container_type &p , container_type &e ) const
     {
-        // ...
-    }
+        // ...
+    }
 };
 
    
 
    
@@ -208,7 +204,6 @@
           All this can be easily done with the following piece of code:
         
    
 
    
-          
 
    
 
    const size_t n = 64;
 container_type q( n , 0.0 ) , p( n , 0.0 );
@@ -235,7 +230,6 @@
           The observer is more or less trivial
         
    
 
    
-          
 
    
 
    struct streaming_observer
 {
@@ -304,7 +298,6 @@
           also need to store the individual frequencies of each oscillator.
         
    
 
    
-          
 
    
 
    typedef vector< double > container_type;
 
@@ -367,7 +360,6 @@
           and we record Z for different values of ε.
         
    
 
    
-          
 
    
 
    struct statistics_observer
 {
@@ -400,7 +392,6 @@
           2.
         
    
 
    
-          
 
    
 
    const size_t n = 16384;
 const double dt = 0.1;
@@ -411,8 +402,8 @@
 boost::uniform_real<> unif( 0.0 , 2.0 * M_PI );
 boost::variate_generator< boost::mt19937&, boost::uniform_real<> > gen( rng , unif );
 
-// gamma = 1, the phase transition occurs at epsilon = 2
-phase_ensemble ensemble( n , 1.0 );
+// gamma = 1, the phase transition occurs at epsilon = 2
+phase_ensemble ensemble( n , 1.0 );
 statistics_observer obs;
 
 for( double epsilon = 0.0 ; epsilon < 5.0 ; epsilon += 0.1 )
@@ -420,14 +411,14 @@
     ensemble.set_epsilon( epsilon );
     obs.reset();
 
-    // start with random initial conditions
-    generate( x.begin() , x.end() , gen );
+    // start with random initial conditions
+    generate( x.begin() , x.end() , gen );
 
-    // calculate some transients steps
-    integrate_const( runge_kutta4< container_type >() , boost::ref( ensemble ) , x , 0.0 , 10.0 , dt );
+    // calculate some transients steps
+    integrate_const( runge_kutta4< container_type >() , boost::ref( ensemble ) , x , 0.0 , 10.0 , dt );
 
-    // integrate and compute the statistics
-    integrate_const( runge_kutta4< container_type >() , boost::ref( ensemble ) , x , 0.0 , 100.0 , dt , boost::ref( obs ) );
+    // integrate and compute the statistics
+    integrate_const( runge_kutta4< container_type >() , boost::ref( ensemble ) , x , 0.0 , 100.0 , dt , boost::ref( obs ) );
     cout << epsilon << "\t" << obs.get_K_mean() << endl;
 }
 
@@ -458,7 +449,6 @@
           all quantities
         
    
 
    
-          
 
    
 
    #include <boost/numeric/odeint.hpp>
 #include <boost/numeric/odeint/algebra/fusion_algebra.hpp>
@@ -496,7 +486,6 @@
           equivalent to the example in the first section of this tutorial
         
    
 
    
-          
 
    
 
    struct oscillator
 {
@@ -520,7 +509,6 @@
           Furthermore, the basic calculations are now performed by the fusion_algebra which must also be given.
         
    
 
    
-          
 
    
 
    typedef runge_kutta_dopri5< state_type , double , deriv_type , time_type , fusion_algebra > stepper_type;
 
@@ -535,7 +523,6 @@
           the observer is defined a bit different
         
    
 
    
-          
 
    
 
    struct streaming_observer
 {
@@ -607,7 +594,6 @@
           The definition of the system is
         
    
 
    
-          
 
    
 
    typedef boost::numeric::ublas::matrix< double > state_type;
 
@@ -701,7 +687,6 @@
           Here is a simple example:
         
    
 
    
-          
 
    
 
    typedef mpf_class value_type;
 typedef boost::array< value_type , 3 > state_type;
@@ -726,7 +711,6 @@
           which can be used
         
    
 
    
-          
 
    
 
    const int precision = 1024;
 mpf_set_default_prec( precision );
@@ -788,7 +772,6 @@
           is kept which should be changed whenever the states' size change.
         
    
 
    
-          
 
    
 
    typedef vector< double > coord_type;
 typedef pair< coord_type , coord_type > state_type;
@@ -804,8 +787,8 @@
         : m_max_N( max_N ) , m_beta( beta ) , m_pot_start_index( pot_start_index ) , m_pot( max_N )
     {
         srand( time( NULL ) );
-        // fill random potential with iid values from [0,1]
-        boost::mt19937 rng;
+        // fill random potential with iid values from [0,1]
+        boost::mt19937 rng;
         boost::uniform_real<> unif( 0.0 , 1.0 );
         boost::variate_generator< boost::mt19937&, boost::uniform_real<> > gen( rng , unif );
         generate( m_pot.begin() , m_pot.end() , gen );
@@ -813,9 +796,9 @@
 
     void operator()( const coord_type &q , coord_type &dpdt )
     {
-        // calculate dpdt = -dH/dq of this hamiltonian system
-        // dp_i/dt = - V_i * q_i^3 - beta*(q_i - q_{i-1})^3 + beta*(q_{i+1} - q_i)^3
-        const int N = q.size();
+        // calculate dpdt = -dH/dq of this hamiltonian system
+        // dp_i/dt = - V_i * q_i^3 - beta*(q_i - q_{i-1})^3 + beta*(q_{i+1} - q_i)^3
+        const int N = q.size();
         double diff = q[0] - q[N-1];
         for( int i=0 ; i<N ; ++i )
         {
@@ -828,8 +811,8 @@
 
     void energy_distribution( const coord_type &q , const coord_type &p , coord_type &energies )
     {
-        // computes the energy per lattice site normalized by total energy
-        const size_t N = q.size();
+        // computes the energy per lattice site normalized by total energy
+        const size_t N = q.size();
         double en = 0.0;
         for( size_t i=0 ; i<N ; i++ )
         {
@@ -848,8 +831,8 @@
 
     double energy( const coord_type &q , const coord_type &p )
     {
-        // calculates the total energy of the excitation
-        const size_t N = q.size();
+        // calculates the total energy of the excitation
+        const size_t N = q.size();
         double en = 0.0;
         for( size_t i=0 ; i<N ; i++ )
         {
@@ -874,25 +857,24 @@
           of 60.
         
    
 
    
-          
 
    
-
    //start with 60 sites
-const int N_start = 60;
+
    //start with 60 sites
+const int N_start = 60;
 coord_type q( N_start , 0.0 );
 q.reserve( max_N );
 coord_type p( N_start , 0.0 );
 p.reserve( max_N );
-// start with uniform momentum distribution over 20 sites
-fill( p.begin()+20 , p.end()-20 , 1.0/sqrt(20) );
+// start with uniform momentum distribution over 20 sites
+fill( p.begin()+20 , p.end()-20 , 1.0/sqrt(20) );
 
 coord_type distr( N_start , 0.0 );
 distr.reserve( max_N );
 
-// create the system
-compacton_lattice lattice( max_N , beta , (max_N-N_start)/2 );
+// create the system
+compacton_lattice lattice( max_N , beta , (max_N-N_start)/2 );
 
-//create the stepper, note that we use an always_resizer because state size might change during steps
-typedef symplectic_rkn_sb3a_mclachlan< coord_type , coord_type , double , coord_type , coord_type , double ,
+//create the stepper, note that we use an always_resizer because state size might change during steps
+typedef symplectic_rkn_sb3a_mclachlan< coord_type , coord_type , double , coord_type , coord_type , double ,
         range_algebra , default_operations , always_resizer > hamiltonian_stepper;
 hamiltonian_stepper stepper;
 hamiltonian_stepper::state_type state = make_pair( q , p );
@@ -910,7 +892,6 @@
           the start index of the potential changes in this case.
         
    
 
    
-          
 
    
 
    double t = 0.0;
 const double dt = 0.1;
@@ -943,7 +924,6 @@
           , p and distr.
         
    
 
    
-          
 
    
 
    void do_resize( coord_type &q , coord_type &p , coord_type &distr , const int N )
 {
@@ -961,7 +941,8 @@
 
    
 
-
    
 
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky
    
+
    Copyright © 2009-2011 Karsten Ahnert
+      and Mario Mulansky
    
         Distributed under the Boost Software License, Version 1.0. (See accompanying
         file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
       
    
diff --git a/doc/boost_sandbox_numeric_odeint/tutorial/stiff_systems.html b/doc/boost_sandbox_numeric_odeint/tutorial/stiff_systems.html
index 6732f950..4fcca818 100644
--- a/doc/boost_sandbox_numeric_odeint/tutorial/stiff_systems.html
+++ b/doc/boost_sandbox_numeric_odeint/tutorial/stiff_systems.html
@@ -45,7 +45,6 @@
         is needed. Here is the definition of the above example
       
    
 
    
-        
 
    
 
    typedef boost::numeric::ublas::vector< double > vector_type;
 typedef boost::numeric::ublas::matrix< double > matrix_type;
@@ -82,7 +81,6 @@
         just templatize the operator():
       
    
 
    
-        
 
    
 
    typedef boost::numeric::ublas::vector< double > vector_type;
 typedef boost::numeric::ublas::matrix< double > matrix_type;
@@ -118,7 +116,6 @@
         all the other stepper:
       
    
 
    
-        
 
    
 
    vector_type x( 3 , 1.0 );
 
@@ -135,7 +132,6 @@
         5 method with step size control and dense output yields 1531 steps.
       
    
 
    
-        
 
    
 
    vector_type x2( 3 , 1.0 );
 
@@ -155,7 +151,8 @@
 
    
 
-
    
 
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky
    
+
    Copyright © 2009-2011 Karsten Ahnert
+      and Mario Mulansky
    
         Distributed under the Boost Software License, Version 1.0. (See accompanying
         file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
       
    
diff --git a/doc/boost_sandbox_numeric_odeint/tutorial/using_cuda_and_thrust.html b/doc/boost_sandbox_numeric_odeint/tutorial/using_cuda_and_thrust.html
index 57a56968..fc0ebba5 100644
--- a/doc/boost_sandbox_numeric_odeint/tutorial/using_cuda_and_thrust.html
+++ b/doc/boost_sandbox_numeric_odeint/tutorial/using_cuda_and_thrust.html
@@ -87,15 +87,14 @@
           website.
         
    
 
    
-          
 
    
-
    //change this to float if your device does not support double computation
-typedef double value_type;
+
    //change this to float if your device does not support double computation
+typedef double value_type;
 
-//change this to host_vector< ... > of you want to run on CPU
-typedef thrust::device_vector< value_type > state_type;
-// typedef thrust::host_vector< value_type > state_type;
-
    
+//change this to host_vector< ... > of you want to run on CPU
+typedef thrust::device_vector< value_type > state_type;
+// typedef thrust::host_vector< value_type > state_type;
+
    
 
    
         
    
 
    
@@ -103,7 +102,7 @@
           a calculation on the GPU you usually have to call a global function like
           thrust::for_each, thrust::reduce,
           ... with an appropriate local functor which performs the basic operation.
-          An example is 
+          An example is
 
    
 
    struct add_two
 {
@@ -115,8 +114,8 @@
     }
 };
 
-// ...
-
+// ...
+
 thrust::for_each( x.begin() , x.end() , add_two() );
 
    
 
    
@@ -138,7 +137,6 @@
           The mean field is calculated in a class mean_field_calculator
         
    
 
    
-          
 
    
 
    struct mean_field_calculator
 {
@@ -192,7 +190,6 @@
           GPU as well as on the CPU. The line
         
    
 
    
-          
 
    
 
    value_type sin_sum = thrust::reduce(
         thrust::make_transform_iterator( x.begin() , sin_functor() ) ,
@@ -208,7 +205,6 @@
           The system function is defined via
         
    
 
    
-          
 
    
 
    class phase_oscillator_ensemble
 {
@@ -230,8 +226,8 @@
         }
     };
 
-    // ...
-
+    // ...
+
     void operator() ( const state_type &x , state_type &dxdt , const value_type dt ) const
     {
         std::pair< value_type , value_type > mean_field = mean_field_calculator::get_mean( x );
@@ -243,8 +239,8 @@
                 );
     }
 
-    // ...
-};
+    // ...
+};
 
    
 
    
         
    
@@ -264,7 +260,6 @@
           thrust algebra:
         
    
 
    
-          
 
    
 
    typedef runge_kutta4< state_type , value_type , state_type , value_type , thrust_algebra , thrust_operations > stepper_type;
 
    
@@ -274,7 +269,6 @@
           Of course, you can also use a controlled or dense output stepper, e.g.
         
    
 
    
-          
 
    
 
    typedef runge_kutta_dopri5< state_type , value_type , state_type , value_type , thrust_algebra , thrust_operations > stepper_type;
 
    
@@ -285,14 +279,12 @@
           an instance of the rhs class and using an integration function:
         
    
 
    
-          
 
    
 
    phase_oscillator_ensemble ensemble( N , 1.0 );
 
    
 
    
         
    
 
    
-          
 
    
 
    size_t steps1 = integrate_const( make_controlled( 1.0e-6 , 1.0e-6 , stepper_type() ) , boost::ref( ensemble ) , x , 0.0 , t_transients , dt );
 
    
@@ -334,14 +326,13 @@
           in front. The full system class is:
         
    
 
    
-          
 
    
-
    //change this to host_vector< ... > of you want to run on CPU
-typedef thrust::device_vector< value_type > state_type;
+
    //change this to host_vector< ... > of you want to run on CPU
+typedef thrust::device_vector< value_type > state_type;
 typedef thrust::device_vector< size_t > index_vector_type;
-//typedef thrust::host_vector< value_type > state_type;
-//typedef thrust::host_vector< size_t > index_vector_type;
-
+//typedef thrust::host_vector< value_type > state_type;
+//typedef thrust::host_vector< size_t > index_vector_type;
+
 class phase_oscillators
 {
 
@@ -351,29 +342,29 @@
     {
         template< class Tuple >
         __host__ __device__
-        void operator()( Tuple t )  // this functor works on tuples of values
-        {
-            // first, unpack the tuple into value, neighbors and omega
-            const value_type phi = thrust::get<0>(t);
-            const value_type phi_left = thrust::get<1>(t);  // left neighbor
-            const value_type phi_right = thrust::get<2>(t); // right neighbor
-            const value_type omega = thrust::get<3>(t);
-            // the dynamical equation
-            thrust::get<4>(t) = omega + sin( phi_right - phi ) + sin( phi - phi_left );
+        void operator()( Tuple t )  // this functor works on tuples of values
+        {
+            // first, unpack the tuple into value, neighbors and omega
+            const value_type phi = thrust::get<0>(t);
+            const value_type phi_left = thrust::get<1>(t);  // left neighbor
+            const value_type phi_right = thrust::get<2>(t); // right neighbor
+            const value_type omega = thrust::get<3>(t);
+            // the dynamical equation
+            thrust::get<4>(t) = omega + sin( phi_right - phi ) + sin( phi - phi_left );
         }
     };
 
     phase_oscillators( const state_type &omega )
         : m_omega( omega ) , m_N( omega.size() ) , m_prev( omega.size() ) , m_next( omega.size() )
     {
-        // build indices pointing to left and right neighbours
-        thrust::counting_iterator<size_t> c( 0 );
+        // build indices pointing to left and right neighbours
+        thrust::counting_iterator<size_t> c( 0 );
         thrust::copy( c , c+m_N-1 , m_prev.begin()+1 );
-        m_prev[0] = 0; // m_prev = { 0 , 0 , 1 , 2 , 3 , ... , N-1 }
-
+        m_prev[0] = 0; // m_prev = { 0 , 0 , 1 , 2 , 3 , ... , N-1 }
+
         thrust::copy( c+1 , c+m_N , m_next.begin() );
-        m_next[m_N-1] = m_N-1; // m_next = { 1 , 2 , 3 , ... , N-1 , N-1 }
-    }
+        m_next[m_N-1] = m_N-1; // m_next = { 1 , 2 , 3 , ... , N-1 , N-1 }
+    }
 
     void operator() ( const state_type &x , state_type &dxdt , const value_type dt )
     {
@@ -423,34 +414,28 @@
           current state, hence the phase of each oscillator.
         
    
 
    
-          
 
    
-
    // create random number generators
-boost::mt19937 rng;
-boost::uniform_real< value_type > unif( 0.0 , 2.0 * pi );
-boost::variate_generator< boost::mt19937&, boost::uniform_real< value_type > > gen( rng , unif );
-
-// create initial conditions and omegas on host:
-vector< value_type > x_host( N );
+
    // create initial conditions and omegas on host:
+vector< value_type > x_host( N );
 vector< value_type > omega_host( N );
 for( size_t i=0 ; i<N ; ++i )
 {
-    x_host[i] = gen();
-    omega_host[i] = ( N - i ) * epsilon; // decreasing frequencies
-}
+    x_host[i] = 2.0 * pi * drand48();
+    omega_host[i] = ( N - i ) * epsilon; // decreasing frequencies
+}
 
-// copy to device
-state_type x = x_host;
+// copy to device
+state_type x = x_host;
 state_type omega = omega_host;
 
-// create stepper
-runge_kutta4< state_type , value_type , state_type , value_type , thrust_algebra , thrust_operations > stepper;
+// create stepper
+runge_kutta4< state_type , value_type , state_type , value_type , thrust_algebra , thrust_operations > stepper;
 
-// create phase oscillator system function
-phase_oscillators sys( omega );
+// create phase oscillator system function
+phase_oscillators sys( omega );
 
-// integrate
-integrate_const( stepper , sys , x , 0.0 , 10.0 , dt );
+// integrate
+integrate_const( stepper , sys , x , 0.0 , 10.0 , dt );
 
 thrust::copy( x.begin() , x.end() , std::ostream_iterator< value_type >( std::cout , "\n" ) );
 std::cout << std::endl;
@@ -491,13 +476,8 @@
           >.
         
    
 
    
-          
 
    
-
    boost::mt19937 rng;
-boost::uniform_real< value_type > unif( 0.0 , 1.0 );
-boost::variate_generator< boost::mt19937&, boost::uniform_real< value_type > > gen( rng , unif );
-
-vector< value_type > beta_host( N );
+
    vector< value_type > beta_host( N );
 const value_type beta_min = 0.0 , beta_max = 56.0;
 for( size_t i=0 ; i<N ; ++i )
     beta_host[i] = beta_min + value_type( i ) * ( beta_max - beta_min ) / value_type( N - 1 );
@@ -513,7 +493,6 @@
           calculates one particular realization of the Lorenz ensemble
         
    
 
    
-          
 
    
 
    struct lorenz_system
 {
@@ -523,8 +502,8 @@
         __host__ __device__
         void operator()( T t ) const
         {
-            // unpack the parameter we want to vary and the Lorenz variables
-            value_type R = thrust::get< 3 >( t );
+            // unpack the parameter we want to vary and the Lorenz variables
+            value_type R = thrust::get< 3 >( t );
             value_type x = thrust::get< 0 >( t );
             value_type y = thrust::get< 1 >( t );
             value_type z = thrust::get< 2 >( t );
@@ -560,7 +539,6 @@
                         boost::begin( dxdt ) + 3 * m_N  ) ) ,
                 lorenz_functor() );
     }
-
     size_t m_N;
     const state_type &m_beta;
 };
@@ -596,7 +574,6 @@
           which is used within this observer is defined
         
    
 
    
-          
 
    
 
    struct lyap_functor
 {
@@ -636,33 +613,32 @@
           and written to stdout.
         
    
 
    
-          
 
    
 
    state_type x( 6 * N );
 
-// initialize x,y,z
-thrust::fill( x.begin() , x.begin() + 3 * N , 10.0 );
+// initialize x,y,z
+thrust::fill( x.begin() , x.begin() + 3 * N , 10.0 );
 
-// initial dx
-thrust::fill( x.begin() + 3 * N , x.begin() + 4 * N , 1.0 );
+// initial dx
+thrust::fill( x.begin() + 3 * N , x.begin() + 4 * N , 1.0 );
 
-// initialize dy,dz
-thrust::fill( x.begin() + 4 * N , x.end() , 0.0 );
+// initialize dy,dz
+thrust::fill( x.begin() + 4 * N , x.end() , 0.0 );
 
 
-// create error stepper, can be used with make_controlled or make_dense_output
-typedef runge_kutta_dopri5< state_type , value_type , state_type , value_type , thrust_algebra , thrust_operations > stepper_type;
+// create error stepper, can be used with make_controlled or make_dense_output
+typedef runge_kutta_dopri5< state_type , value_type , state_type , value_type , thrust_algebra , thrust_operations > stepper_type;
 
 
 lorenz_system lorenz( N , beta );
 lorenz_perturbation_system lorenz_perturbation( N , beta );
 lyap_observer obs( N , 1 );
 
-// calculate transients
-integrate_adaptive( make_controlled( 1.0e-6 , 1.0e-6 , stepper_type() ) , lorenz , std::make_pair( x.begin() , x.begin() + 3 * N ) , 0.0 , 10.0 , dt );
+// calculate transients
+integrate_adaptive( make_controlled( 1.0e-6 , 1.0e-6 , stepper_type() ) , lorenz , std::make_pair( x.begin() , x.begin() + 3 * N ) , 0.0 , 10.0 , dt );
 
-// calculate the Lyapunov exponents -- the main loop
-double t = 0.0;
+// calculate the Lyapunov exponents -- the main loop
+double t = 0.0;
 while( t < 10000.0 )
 {
     integrate_adaptive( make_controlled( 1.0e-6 , 1.0e-6 , stepper_type() ) , lorenz_perturbation , x , t , t + 1.0 , 0.1 );
@@ -685,7 +661,8 @@
 
    
 
-
    
 
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky
    
+
    Copyright © 2009-2011 Karsten Ahnert
+      and Mario Mulansky
    
         Distributed under the Boost Software License, Version 1.0. (See accompanying
         file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
       
    
diff --git a/doc/index.html b/doc/index.html
index d5266d49..16726404 100644
--- a/doc/index.html
+++ b/doc/index.html
@@ -21,9 +21,10 @@
 
    
 Mario Mulansky
 
    
-
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky
    
+
    Copyright © 2009-2011 Karsten Ahnert
+      and Mario Mulansky
    
 
    
-
    
+
    
         Distributed under the Boost Software License, Version 1.0. (See accompanying
         file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
       
    
@@ -133,7 +134,7 @@
 
    
 
    
 
-
+
    Last revised: November 30, 2011 at 13:14:22 GMT
    Last revised: February 08, 2012 at 15:38:05 GMT
    		
 
    		
 
    
 
    
diff --git a/doc/odeint/reference.html b/doc/odeint/reference.html
index 06c2dd53..b742c23f 100644
--- a/doc/odeint/reference.html
+++ b/doc/odeint/reference.html
@@ -19,7 +19,8 @@
 Reference
    
 
-
    
 
    Copyright © 2009 -2011 Karsten Ahnert and Mario Mulansky
    
+
    Copyright © 2009-2011 Karsten Ahnert
+      and Mario Mulansky
    
         Distributed under the Boost Software License, Version 1.0. (See accompanying
         file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
       
    
diff --git a/libs/numeric/odeint/test/stepper_copying.cpp b/libs/numeric/odeint/test/stepper_copying.cpp
index 8928cc85..f4fa34c1 100644
--- a/libs/numeric/odeint/test/stepper_copying.cpp
+++ b/libs/numeric/odeint/test/stepper_copying.cpp
@@ -172,12 +172,11 @@ struct copy_impl< test_array2< T , Dim > , test_array2< T , Dim > >
 
 //provide the state_wrapper
 template< class T , size_t Dim >
-struct state_wrapper< test_array< T , Dim > , false >
+struct state_wrapper< test_array< T , Dim > >
 {
-    typedef state_wrapper< test_array< T , Dim > , false > state_wrapper_type;
+    typedef state_wrapper< test_array< T , Dim > > state_wrapper_type;
     typedef test_array< T , Dim > state_type;
     typedef T value_type;
-    typedef boost::false_type is_resizeable;
 
     state_type m_v;
 
@@ -212,12 +211,11 @@ struct state_wrapper< test_array< T , Dim > , false >
 
 //provide the state_wrapper
 template< class T , size_t Dim >
-struct state_wrapper< test_array2< T , Dim > , false >
+struct state_wrapper< test_array2< T , Dim > >
 {
-    typedef state_wrapper< test_array2< T , Dim > , false > state_wrapper_type;
+    typedef state_wrapper< test_array2< T , Dim > > state_wrapper_type;
     typedef test_array2< T , Dim > state_type;
     typedef T value_type;
-    typedef boost::false_type is_resizeable;
 
     state_type m_v;