213 lines
7.2 KiB
C++
213 lines
7.2 KiB
C++
/*
|
|
[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
|
|
{
|
|
using std::abs;
|
|
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
|
|
{
|
|
using std::abs;
|
|
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()
|