1 /*
2 [auto_generated]
3 libs/numeric/odeint/test/fusion_algebra.cpp
4
5 [begin_description]
6 This file tests the Fusion algebra.
7 [end_description]
8
9 Copyright 2011-2012 Karsten Ahnert
10 Copyright 2011-2013 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_fusion_algebra
19
20 // using fusion vectors as state types requires increased macro variables
21 #define BOOST_FUSION_INVOKE_MAX_ARITY 15
22 #define BOOST_RESULT_OF_NUM_ARGS 15
23
24 #include <cmath>
25 #include <complex>
26 #include <utility>
27 #include <functional>
28
29 #include <boost/numeric/odeint/config.hpp>
30
31 #include <boost/test/unit_test.hpp>
32 #include <boost/test/floating_point_comparison.hpp>
33
34 #include <boost/units/systems/si/length.hpp>
35 #include <boost/units/systems/si/time.hpp>
36 #include <boost/units/systems/si/velocity.hpp>
37 #include <boost/units/systems/si/acceleration.hpp>
38 #include <boost/units/systems/si/io.hpp>
39
40 #include <boost/fusion/container.hpp>
41
42 #include <boost/numeric/odeint/algebra/default_operations.hpp>
43 #include <boost/numeric/odeint/algebra/fusion_algebra.hpp>
44
45 namespace units = boost::units;
46 namespace si = boost::units::si;
47 namespace fusion = boost::fusion;
48
49 using boost::numeric::odeint::default_operations;
50 using boost::numeric::odeint::fusion_algebra;
51
52 typedef double value_type;
53 typedef units::quantity< si::time , value_type > time_type;
54 typedef units::quantity< si::length , value_type > length_type;
55 typedef units::quantity< si::velocity , value_type > velocity_type;
56 typedef units::quantity< si::acceleration , value_type > acceleration_type;
57 typedef fusion::vector< length_type , velocity_type > state_type;
58 typedef fusion::vector< velocity_type , acceleration_type > deriv_type;
59
60 const time_type dt = 0.1 * si::second;
61
62
63 struct fusion_fixture
64 {
fusion_fixturefusion_fixture65 fusion_fixture( void )
66 : res( 0.0 * si::meter , 0.0 * si::meter_per_second ) ,
67 x( 1.0 * si::meter , 1.0 * si::meter_per_second ) ,
68 k1( 1.0 * si::meter_per_second , 1.0 * si::meter_per_second_squared ) ,
69 k2( 2.0 * si::meter_per_second , 2.0 * si::meter_per_second_squared ) ,
70 k3( 3.0 * si::meter_per_second , 3.0 * si::meter_per_second_squared ) ,
71 k4( 4.0 * si::meter_per_second , 4.0 * si::meter_per_second_squared ) ,
72 k5( 5.0 * si::meter_per_second , 5.0 * si::meter_per_second_squared ) ,
73 k6( 6.0 * si::meter_per_second , 6.0 * si::meter_per_second_squared ) { }
74
~fusion_fixturefusion_fixture75 ~fusion_fixture( void )
76 {
77 BOOST_CHECK_CLOSE( fusion::at_c< 0 >( x ).value() , 1.0 , 1.0e-10 );
78 BOOST_CHECK_CLOSE( fusion::at_c< 1 >( x ).value() , 1.0 , 1.0e-10 );
79 BOOST_CHECK_CLOSE( fusion::at_c< 0 >( k1 ).value() , 1.0 , 1.0e-10 );
80 BOOST_CHECK_CLOSE( fusion::at_c< 1 >( k1 ).value() , 1.0 , 1.0e-10 );
81 BOOST_CHECK_CLOSE( fusion::at_c< 0 >( k2 ).value() , 2.0 , 1.0e-10 );
82 BOOST_CHECK_CLOSE( fusion::at_c< 1 >( k2 ).value() , 2.0 , 1.0e-10 );
83 BOOST_CHECK_CLOSE( fusion::at_c< 0 >( k3 ).value() , 3.0 , 1.0e-10 );
84 BOOST_CHECK_CLOSE( fusion::at_c< 1 >( k3 ).value() , 3.0 , 1.0e-10 );
85 BOOST_CHECK_CLOSE( fusion::at_c< 0 >( k4 ).value() , 4.0 , 1.0e-10 );
86 BOOST_CHECK_CLOSE( fusion::at_c< 1 >( k4 ).value() , 4.0 , 1.0e-10 );
87 BOOST_CHECK_CLOSE( fusion::at_c< 0 >( k5 ).value() , 5.0 , 1.0e-10 );
88 BOOST_CHECK_CLOSE( fusion::at_c< 1 >( k5 ).value() , 5.0 , 1.0e-10 );
89 BOOST_CHECK_CLOSE( fusion::at_c< 0 >( k6 ).value() , 6.0 , 1.0e-10 );
90 BOOST_CHECK_CLOSE( fusion::at_c< 1 >( k6 ).value() , 6.0 , 1.0e-10 );
91 }
92
93 state_type res , x;
94 deriv_type k1 , k2 , k3 , k4 , k5 , k6 ;
95 };
96
97
98
99 BOOST_AUTO_TEST_SUITE( fusion_algebra_test )
100
101 fusion_algebra algebra;
102
BOOST_AUTO_TEST_CASE(for_each2)103 BOOST_AUTO_TEST_CASE( for_each2 )
104 {
105 fusion_fixture f;
106 algebra.for_each2( f.res , f.k1 ,
107 default_operations::scale_sum1< time_type >( dt ) );
108 BOOST_CHECK_CLOSE( fusion::at_c< 0 >( f.res ).value() , 0.1 , 1.0e-10 );
109 BOOST_CHECK_CLOSE( fusion::at_c< 1 >( f.res ).value() , 0.1 , 1.0e-10 );
110 }
111
BOOST_AUTO_TEST_CASE(for_each3)112 BOOST_AUTO_TEST_CASE( for_each3 )
113 {
114 fusion_fixture f;
115 algebra.for_each3( f.res , f.x , f.k1 ,
116 default_operations::scale_sum2< value_type , time_type >( 1.0 , dt ) );
117 BOOST_CHECK_CLOSE( fusion::at_c< 0 >( f.res ).value() , 1.1 , 1.0e-10 );
118 BOOST_CHECK_CLOSE( fusion::at_c< 1 >( f.res ).value() , 1.1 , 1.0e-10 );
119 }
120
BOOST_AUTO_TEST_CASE(for_each4)121 BOOST_AUTO_TEST_CASE( for_each4 )
122 {
123 fusion_fixture f;
124 algebra.for_each4( f.res , f.x , f.k1 , f.k2 ,
125 default_operations::scale_sum3< value_type , time_type , time_type >( 1.0 , dt , dt ) );
126 BOOST_CHECK_CLOSE( fusion::at_c< 0 >( f.res ).value() , 1.3 , 1.0e-10 );
127 BOOST_CHECK_CLOSE( fusion::at_c< 1 >( f.res ).value() , 1.3 , 1.0e-10 );
128 }
129
BOOST_AUTO_TEST_CASE(for_each5)130 BOOST_AUTO_TEST_CASE( for_each5 )
131 {
132 fusion_fixture f;
133 algebra.for_each5( f.res , f.x , f.k1 , f.k2 , f.k3 ,
134 default_operations::scale_sum4< value_type , time_type , time_type , time_type >( 1.0 , dt , dt , dt ) );
135 BOOST_CHECK_CLOSE( fusion::at_c< 0 >( f.res ).value() , 1.6 , 1.0e-10 );
136 BOOST_CHECK_CLOSE( fusion::at_c< 1 >( f.res ).value() , 1.6 , 1.0e-10 );
137 }
138
BOOST_AUTO_TEST_CASE(for_each6)139 BOOST_AUTO_TEST_CASE( for_each6 )
140 {
141 fusion_fixture f;
142 algebra.for_each6( f.res , f.x , f.k1 , f.k2 , f.k3 , f.k4 ,
143 default_operations::scale_sum5< value_type , time_type , time_type , time_type , time_type >( 1.0 , dt , dt , dt , dt ) );
144 BOOST_CHECK_CLOSE( fusion::at_c< 0 >( f.res ).value() , 2.0 , 1.0e-10 );
145 BOOST_CHECK_CLOSE( fusion::at_c< 1 >( f.res ).value() , 2.0 , 1.0e-10 );
146 }
147
BOOST_AUTO_TEST_CASE(for_each7)148 BOOST_AUTO_TEST_CASE( for_each7 )
149 {
150 fusion_fixture f;
151 algebra.for_each7( f.res , f.x , f.k1 , f.k2 , f.k3 , f.k4 , f.k5 ,
152 default_operations::scale_sum6< value_type , time_type , time_type , time_type , time_type , time_type >( 1.0 , dt , dt , dt , dt , dt ) );
153 BOOST_CHECK_CLOSE( fusion::at_c< 0 >( f.res ).value() , 2.5 , 1.0e-10 );
154 BOOST_CHECK_CLOSE( fusion::at_c< 1 >( f.res ).value() , 2.5 , 1.0e-10 );
155 }
156
BOOST_AUTO_TEST_CASE(for_each8)157 BOOST_AUTO_TEST_CASE( for_each8 )
158 {
159 fusion_fixture f;
160 algebra.for_each8( f.res , f.x , f.k1 , f.k2 , f.k3 , f.k4 , f.k5 , f.k6 ,
161 default_operations::scale_sum7< value_type , time_type , time_type , time_type , time_type , time_type , time_type >( 1.0 , dt , dt , dt , dt , dt , dt ) );
162 BOOST_CHECK_CLOSE( fusion::at_c< 0 >( f.res ).value() , 3.1 , 1.0e-10 );
163 BOOST_CHECK_CLOSE( fusion::at_c< 1 >( f.res ).value() , 3.1 , 1.0e-10 );
164 }
165
BOOST_AUTO_TEST_CASE(for_each15)166 BOOST_AUTO_TEST_CASE( for_each15 )
167 {
168 fusion_fixture f;
169 algebra.for_each15( f.res , f.x , f.k1 , f.k2 , f.k3 , f.k4 , f.k5 , f.k6 , f.k1 , f.k2 , f.k3 , f.k4 , f.k5 , f.k6 , f.k1 ,
170 default_operations::scale_sum14< value_type , time_type , time_type , time_type , time_type , time_type , time_type ,
171 time_type , time_type , time_type , time_type , time_type , time_type , time_type >( 1.0 , dt , dt , dt , dt , dt , dt , dt , dt , dt , dt , dt , dt , dt ) );
172 //BOOST_CHECK_CLOSE( fusion::at_c< 0 >( f.res ).value() , 3.1 , 1.0e-10 );
173 //BOOST_CHECK_CLOSE( fusion::at_c< 1 >( f.res ).value() , 3.1 , 1.0e-10 );
174 }
175
176
BOOST_AUTO_TEST_CASE(norm_inf)177 BOOST_AUTO_TEST_CASE( norm_inf )
178 {
179 double nrm = algebra.norm_inf( fusion::make_vector( 1.0 , 2.0 , 3.0 ) );
180 BOOST_CHECK_CLOSE( nrm , 3.0 , 1.0e-10 );
181 }
182
183
184 BOOST_AUTO_TEST_SUITE_END()
185