1 /*
2
3 [begin_description]
4 Test case for issue 147
5 [end_description]
6
7 Copyright 2011-2015 Karsten Ahnert
8 Copyright 2011-2015 Mario Mulansky
9
10 Distributed under the Boost Software License, Version 1.0.
11 (See accompanying file LICENSE_1_0.txt or
12 copy at http://www.boost.org/LICENSE_1_0.txt)
13 */
14
15
16 // disable checked iterator warning for msvc
17
18 #include <boost/config.hpp>
19 #ifdef BOOST_MSVC
20 #pragma warning(disable:4996)
21 #endif
22
23 #define BOOST_TEST_MODULE odeint_regression_147
24
25 #include <utility>
26
27 #include <boost/array.hpp>
28
29 #include <boost/test/unit_test.hpp>
30
31 #include <boost/mpl/vector.hpp>
32
33 #include <boost/numeric/odeint.hpp>
34
35 using namespace boost::unit_test;
36 using namespace boost::numeric::odeint;
37 namespace mpl = boost::mpl;
38
39 typedef double state_type;
40
rhs(const state_type & x,state_type & dxdt,const double t)41 void rhs( const state_type &x , state_type &dxdt , const double t )
42 {
43 dxdt = 1;
44 }
45
46
47 template<class Stepper, class InitStepper>
48 struct perform_init_test
49 {
operator ()perform_init_test50 void operator()( void )
51 {
52 double t = 0;
53 const double dt = 0.1;
54
55 state_type x = 0;
56
57 Stepper stepper;
58 InitStepper init_stepper;
59 stepper.initialize( init_stepper, rhs, x, t, dt );
60
61 // ab-stepper needs order-1 init steps: t and x should be (order-1)*dt
62 BOOST_CHECK_CLOSE( t , (stepper.order()-1)*dt , 1E-16 );
63 BOOST_CHECK_CLOSE( x, ( stepper.order() - 1 ) * dt, 2E-14 );
64 }
65 };
66
67 typedef mpl::vector<
68 euler< state_type > ,
69 modified_midpoint< state_type > ,
70 runge_kutta4< state_type > ,
71 runge_kutta4_classic< state_type > ,
72 runge_kutta_cash_karp54_classic< state_type > ,
73 runge_kutta_cash_karp54< state_type > ,
74 runge_kutta_dopri5< state_type > ,
75 runge_kutta_fehlberg78< state_type >
76 > runge_kutta_steppers;
77
78
79 BOOST_AUTO_TEST_SUITE( regression_147_test )
80
BOOST_AUTO_TEST_CASE_TEMPLATE(init_test,InitStepper,runge_kutta_steppers)81 BOOST_AUTO_TEST_CASE_TEMPLATE( init_test , InitStepper,
82 runge_kutta_steppers )
83 {
84 perform_init_test< adams_bashforth<4, state_type>, InitStepper > tester;
85 tester();
86 }
87
88 BOOST_AUTO_TEST_SUITE_END()
89