| /* |
| [auto_generated] |
| libs/numeric/odeint/test/runge_kutta_controlled_concepts.cpp |
| |
| [begin_description] |
| This file tests the Stepper concepts of odeint with all Controlled Runge-Kutta steppers. |
| It's one of the main tests of odeint. |
| [end_description] |
| |
| Copyright 2012 Karsten Ahnert |
| Copyright 2012-2013 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) |
| */ |
| |
| // disable checked iterator warning for msvc |
| #include <boost/config.hpp> |
| #ifdef BOOST_MSVC |
| #pragma warning(disable:4996) |
| #endif |
| |
| #define BOOST_TEST_MODULE odeint_runge_kutta_controlled_concepts |
| |
| #include <vector> |
| #include <cmath> |
| #include <iostream> |
| |
| #include <boost/numeric/odeint/config.hpp> |
| |
| #include <boost/array.hpp> |
| |
| #include <boost/test/unit_test.hpp> |
| |
| #include <boost/ref.hpp> |
| #include <boost/bind.hpp> |
| #include <boost/utility.hpp> |
| #include <boost/type_traits/add_reference.hpp> |
| |
| #include <boost/mpl/vector.hpp> |
| #include <boost/mpl/for_each.hpp> |
| #include <boost/mpl/insert_range.hpp> |
| #include <boost/mpl/end.hpp> |
| #include <boost/mpl/copy.hpp> |
| #include <boost/mpl/placeholders.hpp> |
| #include <boost/mpl/inserter.hpp> |
| |
| #include <boost/numeric/odeint/stepper/runge_kutta_cash_karp54_classic.hpp> |
| #include <boost/numeric/odeint/stepper/runge_kutta_cash_karp54.hpp> |
| #include <boost/numeric/odeint/stepper/runge_kutta_dopri5.hpp> |
| #include <boost/numeric/odeint/stepper/runge_kutta_fehlberg78.hpp> |
| #include <boost/numeric/odeint/stepper/controlled_runge_kutta.hpp> |
| #include <boost/numeric/odeint/stepper/bulirsch_stoer.hpp> |
| |
| #include "prepare_stepper_testing.hpp" |
| #include "dummy_odes.hpp" |
| |
| using std::vector; |
| |
| using namespace boost::unit_test; |
| using namespace boost::numeric::odeint; |
| namespace mpl = boost::mpl; |
| |
| const double result = 2.2; // two steps |
| |
| const double eps = 1.0e-14; |
| |
| template< class Stepper , class System > |
| void check_controlled_stepper_concept( Stepper &stepper , System system , typename Stepper::state_type &x ) |
| { |
| typedef Stepper stepper_type; |
| typedef typename stepper_type::deriv_type container_type; |
| //typedef typename stepper_type::order_type order_type; controlled_error_stepper don't necessarily have a order (burlish-stoer) |
| typedef typename stepper_type::time_type time_type; |
| |
| time_type t = 0.0 , dt = 0.1; |
| controlled_step_result step_result = stepper.try_step( system , x , t , dt ); |
| |
| BOOST_CHECK_MESSAGE( step_result == success , "step result: " << step_result ); // error = 0 for constant system -> step size is always too small |
| } |
| |
| |
| template< class ControlledStepper , typename T > |
| struct perform_controlled_stepper_test |
| { |
| typedef T vector_space_type; |
| void operator()( void ) const |
| { |
| vector_space_type x; |
| x = 2.0; |
| ControlledStepper controlled_stepper; |
| constant_system_functor_vector_space sys; |
| #ifndef _MSC_VER |
| // dont run this for MSVC due to compiler bug 697006 |
| check_controlled_stepper_concept( controlled_stepper , |
| constant_system_vector_space< vector_space_type , vector_space_type , typename ControlledStepper::time_type > |
| , x ); |
| #else |
| check_controlled_stepper_concept( controlled_stepper , boost::cref( sys ) , x ); |
| #endif |
| check_controlled_stepper_concept( controlled_stepper , boost::cref( sys ) , x ); |
| BOOST_CHECK( (abs( x - result )) < eps ); |
| } |
| }; |
| |
| template< class ControlledStepper , typename T > |
| struct perform_controlled_stepper_test< ControlledStepper , std::vector<T> > |
| { |
| typedef std::vector<T> vector_type; |
| void operator()( void ) |
| { |
| using std::abs; |
| vector_type x( 1 , 2.0 ); |
| ControlledStepper controlled_stepper; |
| constant_system_functor_standard sys; |
| #ifndef _MSC_VER |
| // dont run this for MSVC due to compiler bug 697006 |
| |
| check_controlled_stepper_concept( controlled_stepper , |
| constant_system_standard< vector_type , vector_type , typename ControlledStepper::time_type > , |
| x ); |
| #else |
| check_controlled_stepper_concept( controlled_stepper , boost::cref( sys ) , x ); |
| #endif |
| check_controlled_stepper_concept( controlled_stepper , boost::cref( sys ) , x ); |
| BOOST_CHECK( (abs( x[0] - result )) < eps ); |
| } |
| }; |
| |
| |
| template< class ControlledStepper > |
| struct perform_controlled_stepper_test< ControlledStepper , vector_space_type > |
| { |
| void operator()( void ) const |
| { |
| vector_space_type x; |
| x = 2.0; |
| ControlledStepper controlled_stepper; |
| constant_system_functor_vector_space sys; |
| #ifndef _MSC_VER |
| // dont run this for MSVC due to compiler bug 697006 |
| check_controlled_stepper_concept( controlled_stepper , |
| constant_system_vector_space< vector_space_type , vector_space_type , typename ControlledStepper::time_type > |
| , x ); |
| #else |
| check_controlled_stepper_concept( controlled_stepper , boost::cref( sys ) , x ); |
| #endif |
| check_controlled_stepper_concept( controlled_stepper , boost::cref( sys ) , x ); |
| BOOST_CHECK( (abs( x - result )) < eps ); |
| } |
| }; |
| |
| template< class ControlledStepper , typename T > |
| struct perform_controlled_stepper_test< ControlledStepper , boost::array<T,1> > |
| { |
| typedef boost::array<T,1> array_type; |
| void operator()( void ) |
| { |
| using std::abs; |
| array_type x; |
| x[0] = 2.0; |
| ControlledStepper controlled_stepper; |
| constant_system_functor_standard sys; |
| #ifndef _MSC_VER |
| // dont run this for MSVC due to compiler bug 697006 |
| check_controlled_stepper_concept( controlled_stepper , constant_system_standard< array_type , array_type , typename ControlledStepper::time_type > , x ); |
| #else |
| check_controlled_stepper_concept( controlled_stepper , boost::cref( sys ) , x ); |
| #endif |
| check_controlled_stepper_concept( controlled_stepper , boost::cref( sys ) , x ); |
| BOOST_CHECK( (abs( x[0] - result )) < eps ); |
| } |
| }; |
| |
| |
| template< class State > class controlled_stepper_methods : public mpl::vector< |
| controlled_runge_kutta< runge_kutta_cash_karp54_classic< State , typename detail::extract_value_type<State>::type > > , |
| controlled_runge_kutta< runge_kutta_dopri5< State , typename detail::extract_value_type<State>::type > > , |
| controlled_runge_kutta< runge_kutta_fehlberg78< State , typename detail::extract_value_type<State>::type > > , |
| bulirsch_stoer< State , typename detail::extract_value_type<State>::type > |
| > { }; |
| |
| typedef mpl::copy |
| < |
| container_types , |
| mpl::inserter |
| < |
| mpl::vector0<> , |
| mpl::insert_range |
| < |
| mpl::_1 , |
| mpl::end< mpl::_1 > , |
| controlled_stepper_methods< mpl::_2 > |
| > |
| > |
| >::type all_controlled_stepper_methods; |
| |
| |
| |
| |
| BOOST_AUTO_TEST_SUITE( controlled_runge_kutta_concept_test ) |
| |
| BOOST_AUTO_TEST_CASE_TEMPLATE( controlled_stepper_test , ControlledStepper , all_controlled_stepper_methods ) |
| { |
| perform_controlled_stepper_test< ControlledStepper , typename ControlledStepper::state_type > tester; |
| tester(); |
| } |
| |
| BOOST_AUTO_TEST_SUITE_END() |