blob: 6432692d7e4ad7061f5f9c89da4bc582b232e9e2 [file] [log] [blame]
Brian Silverman7c33ab22018-08-04 17:14:51 -07001/*
2 [auto_generated]
3 libs/numeric/odeint/test/generation.cpp
4
5 [begin_description]
6 This file tests the generation functions.
7 [end_description]
8
9 Copyright 2011-2012 Karsten Ahnert
10 Copyright 2011-2012 Mario Mulansky
11
12 Distributed under the Boost Software License, Version 1.0.
13 (See accompanying file LICENSE_1_0.txt or
14 copy at http://www.boost.org/LICENSE_1_0.txt)
15 */
16
17
18#define BOOST_TEST_MODULE odeint_generation
19
20#include <vector>
21
22#include <boost/test/unit_test.hpp>
23#include <boost/static_assert.hpp>
24#include <boost/type_traits/is_same.hpp>
25
26#include <boost/numeric/odeint/config.hpp>
27#include <boost/numeric/odeint/stepper/generation.hpp>
28
29using namespace boost::numeric::odeint;
30
31template< class Stepper1 , class Stepper2 >
32void check_stepper_type( const Stepper1 &s1 , const Stepper2 &s2 )
33{
34 BOOST_STATIC_ASSERT(( boost::is_same< Stepper1 , Stepper2 >::value ));
35}
36
37BOOST_AUTO_TEST_SUITE( generation_test )
38
39typedef std::vector< double > state_type;
40
41typedef runge_kutta_cash_karp54< state_type > rk54_type;
42typedef runge_kutta_cash_karp54_classic< state_type > rk54_classic_type;
43typedef runge_kutta_dopri5< state_type > dopri5_type;
44typedef runge_kutta_fehlberg78< state_type > fehlberg78_type;
45typedef rosenbrock4< double > rosenbrock4_type;
46
47
48BOOST_AUTO_TEST_CASE( test_generation_dopri5 )
49{
50 check_stepper_type( make_controlled( 1.0e-6 , 1.0e-10 , dopri5_type() ) , controlled_runge_kutta< dopri5_type >() );
51 check_stepper_type( make_controlled< dopri5_type >( 1.0e-6 , 1.0e-10 ) , controlled_runge_kutta< dopri5_type >() );
52
53 check_stepper_type( make_dense_output( 1.0e-6 , 1.0e-10 , dopri5_type() ) , dense_output_runge_kutta< controlled_runge_kutta< dopri5_type > >() );
54 check_stepper_type( make_dense_output< dopri5_type >( 1.0e-6 , 1.0e-10 ) , dense_output_runge_kutta< controlled_runge_kutta< dopri5_type > >() );
55}
56
57BOOST_AUTO_TEST_CASE( test_generation_cash_karp54 )
58{
59 check_stepper_type( make_controlled( 1.0e-6 , 1.0e-10 , rk54_type() ) , controlled_runge_kutta< rk54_type >() );
60 check_stepper_type( make_controlled< rk54_type >( 1.0e-6 , 1.0e-10 ) , controlled_runge_kutta< rk54_type >() );
61}
62
63BOOST_AUTO_TEST_CASE( test_generation_cash_karp54_classic )
64{
65 check_stepper_type( make_controlled( 1.0e-6 , 1.0e-10 , rk54_classic_type() ) , controlled_runge_kutta< rk54_classic_type >() );
66 check_stepper_type( make_controlled< rk54_classic_type >( 1.0e-6 , 1.0e-10 ) , controlled_runge_kutta< rk54_classic_type >() );
67}
68
69
70BOOST_AUTO_TEST_CASE( test_generation_fehlberg78 )
71{
72 check_stepper_type( make_controlled( 1.0e-6 , 1.0e-10 , fehlberg78_type() ) , controlled_runge_kutta< fehlberg78_type >() );
73 check_stepper_type( make_controlled< fehlberg78_type >( 1.0e-6 , 1.0e-10 ) , controlled_runge_kutta< fehlberg78_type >() );
74}
75
76
77BOOST_AUTO_TEST_CASE( test_generation_rosenbrock4 )
78{
79 check_stepper_type( make_controlled( 1.0e-6 , 1.0e-10 , rosenbrock4_type() ) , rosenbrock4_controller< rosenbrock4_type >() );
80 check_stepper_type( make_controlled< rosenbrock4_type >( 1.0e-6 , 1.0e-10 ) , rosenbrock4_controller< rosenbrock4_type >() );
81
82 check_stepper_type( make_dense_output( 1.0e-6 , 1.0e-10 , rosenbrock4_type() ) , rosenbrock4_dense_output< rosenbrock4_controller< rosenbrock4_type > >() );
83 check_stepper_type( make_dense_output< rosenbrock4_type >( 1.0e-6 , 1.0e-10 ) , rosenbrock4_dense_output< rosenbrock4_controller< rosenbrock4_type > >() );
84}
85
86
87BOOST_AUTO_TEST_SUITE_END()