1 /*
2 [auto_generated]
3 libs/numeric/odeint/test/adams_bashforth_moulton.cpp
4
5 [begin_description]
6 This file tests the use of the Adams-Bashforth-Moulton.
7 [end_description]
8
9 Copyright 2011-2012 Karsten Ahnert
10 Copyright 2011 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 #define BOOST_TEST_MODULE odeint_adams_bashforth_moulton
18
19 #include <utility>
20
21 #include <boost/test/unit_test.hpp>
22
23 #include <boost/mpl/list.hpp>
24 #include <boost/mpl/size_t.hpp>
25 #include <boost/mpl/range_c.hpp>
26
27 #include <boost/numeric/odeint/stepper/adams_bashforth_moulton.hpp>
28
29 using namespace boost::unit_test;
30 using namespace boost::numeric::odeint;
31
32 typedef double value_type;
33
34 struct lorenz
35 {
36 template< class State , class Deriv , class Value >
operator ()lorenz37 void operator()( const State &_x , Deriv &_dxdt , const Value &dt ) const
38 {
39 const value_type sigma = 10.0;
40 const value_type R = 28.0;
41 const value_type b = 8.0 / 3.0;
42
43 typename boost::range_iterator< const State >::type x = boost::begin( _x );
44 typename boost::range_iterator< Deriv >::type dxdt = boost::begin( _dxdt );
45
46 dxdt[0] = sigma * ( x[1] - x[0] );
47 dxdt[1] = R * x[0] - x[1] - x[0] * x[2];
48 dxdt[2] = x[0]*x[1] - b * x[2];
49 }
50 };
51
52 BOOST_AUTO_TEST_SUITE( adams_bashforth_moulton_test )
53
54 typedef boost::mpl::range_c< size_t , 1 , 6 > vector_of_steps;
55
BOOST_AUTO_TEST_CASE_TEMPLATE(test_init_and_steps,step_type,vector_of_steps)56 BOOST_AUTO_TEST_CASE_TEMPLATE( test_init_and_steps , step_type , vector_of_steps )
57 {
58 const static size_t steps = step_type::value;
59 typedef boost::array< value_type , 3 > state_type;
60
61 adams_bashforth_moulton< steps , state_type > stepper;
62 state_type x = {{ 10.0 , 10.0 , 10.0 }};
63 const value_type dt = 0.01;
64 value_type t = 0.0;
65
66 stepper.initialize( lorenz() , x , t , dt );
67 BOOST_CHECK_CLOSE( t , value_type( steps - 1 ) * dt , 1.0e-14 );
68
69 stepper.do_step( lorenz() , x , t , dt );
70 }
71
72
BOOST_AUTO_TEST_CASE(test_copying)73 BOOST_AUTO_TEST_CASE( test_copying )
74 {
75 typedef boost::array< double , 1 > state_type;
76 typedef adams_bashforth_moulton< 2 , state_type > stepper_type;
77
78 stepper_type s1;
79
80 stepper_type s2( s1 );
81
82 stepper_type s3;
83 s3 = s1;
84 }
85
86
BOOST_AUTO_TEST_CASE(test_instantiation)87 BOOST_AUTO_TEST_CASE( test_instantiation )
88 {
89 typedef boost::array< double , 3 > state_type;
90 adams_bashforth_moulton< 1 , state_type > s1;
91 adams_bashforth_moulton< 2 , state_type > s2;
92 adams_bashforth_moulton< 3 , state_type > s3;
93 adams_bashforth_moulton< 4 , state_type > s4;
94 adams_bashforth_moulton< 5 , state_type > s5;
95 adams_bashforth_moulton< 6 , state_type > s6;
96 adams_bashforth_moulton< 7 , state_type > s7;
97 adams_bashforth_moulton< 8 , state_type > s8;
98
99 state_type x = {{ 10.0 , 10.0 , 10.0 }};
100 value_type t = 0.0 , dt = 0.01;
101 s1.do_step( lorenz() , x , t , dt );
102 s2.do_step( lorenz() , x , t , dt );
103 s3.do_step( lorenz() , x , t , dt );
104 s4.do_step( lorenz() , x , t , dt );
105 s5.do_step( lorenz() , x , t , dt );
106 s6.do_step( lorenz() , x , t , dt );
107 // s7.do_step( lorenz() , x , t , dt );
108 // s8.do_step( lorenz() , x , t , dt );
109 }
110
111
112 BOOST_AUTO_TEST_SUITE_END()
113