3 boost/numeric/odeint/stepper/detail/generic_rk_operations.hpp
6 Operations caller for the generic Runge Kutta method.
9 Copyright 2011 Mario Mulansky
10 Copyright 2011-2012 Karsten Ahnert
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)
18 #ifndef BOOST_NUMERIC_ODEINT_STEPPER_DETAIL_GENERIC_RK_OPERATIONS_HPP_INCLUDED
19 #define BOOST_NUMERIC_ODEINT_STEPPER_DETAIL_GENERIC_RK_OPERATIONS_HPP_INCLUDED
27 template< size_t StageNumber , class Operations , class Fac , class Time >
28 struct generic_rk_scale_sum;
30 template< class Operations , class Fac , class Time >
31 struct generic_rk_scale_sum< 1 , Operations , Fac , Time > : public Operations::template scale_sum2< Fac , Time >
33 generic_rk_scale_sum( const boost::array<Fac,1> &a , Time dt ) : Operations::template scale_sum2< Fac , Time >( 1.0 , a[0]*dt )
36 typedef void result_type;
40 template< class Operations , class Fac , class Time >
41 struct generic_rk_scale_sum< 2 , Operations , Fac , Time > : public Operations::template scale_sum3< Fac , Time >
43 generic_rk_scale_sum( const boost::array<Fac,2> &a , Time dt )
44 : Operations::template scale_sum3< Fac , Time >( 1.0 , a[0]*dt , a[1]*dt )
47 typedef void result_type;
50 template< class Operations , class Fac , class Time >
51 struct generic_rk_scale_sum< 3 , Operations , Fac , Time > : public Operations::template scale_sum4< Fac , Time >
53 generic_rk_scale_sum( const boost::array<Fac,3> &a , Time dt )
54 : Operations::template scale_sum4< Fac , Time >( 1.0 , a[0]*dt , a[1]*dt , a[2]*dt )
57 typedef void result_type;
60 template< class Operations , class Fac , class Time >
61 struct generic_rk_scale_sum< 4 , Operations , Fac , Time > : public Operations::template scale_sum5< Fac , Time >
63 generic_rk_scale_sum( const boost::array<Fac,4> &a , Time dt )
64 : Operations::template scale_sum5< Fac , Time >( 1.0 , a[0]*dt , a[1]*dt , a[2]*dt , a[3]*dt )
67 typedef void result_type;
70 template< class Operations , class Fac , class Time >
71 struct generic_rk_scale_sum< 5 , Operations , Fac , Time > : public Operations::template scale_sum6< Fac , Time >
73 generic_rk_scale_sum( const boost::array<Fac,5> &a , Time dt )
74 : Operations::template scale_sum6< Fac , Time >( 1.0 , a[0]*dt , a[1]*dt , a[2]*dt , a[3]*dt , a[4]*dt )
77 typedef void result_type;
80 template< class Operations , class Fac , class Time >
81 struct generic_rk_scale_sum< 6 , Operations , Fac , Time > : public Operations::template scale_sum7< Fac , Time >
83 generic_rk_scale_sum( const boost::array<Fac,6> &a , Time dt )
84 : Operations::template scale_sum7< Fac , Time >( 1.0 , a[0]*dt , a[1]*dt , a[2]*dt , a[3]*dt , a[4]*dt , a[5]*dt )
87 typedef void result_type;
90 template< class Operations , class Fac , class Time >
91 struct generic_rk_scale_sum< 7 , Operations , Fac , Time > : public Operations::template scale_sum8< Fac , Time >
93 generic_rk_scale_sum( const boost::array<Fac,7> &a , Time dt )
94 : Operations::template scale_sum8< Fac , Time >( 1.0 , a[0]*dt , a[1]*dt , a[2]*dt , a[3]*dt , a[4]*dt , a[5]*dt , a[6]*dt )
97 typedef void result_type;
100 template< class Operations , class Fac , class Time >
101 struct generic_rk_scale_sum< 8 , Operations , Fac , Time > : public Operations::template scale_sum9< Fac , Time >
103 generic_rk_scale_sum( const boost::array<Fac,8> &a , Time dt )
104 : Operations::template scale_sum9< Fac , Time >( 1.0 , a[0]*dt , a[1]*dt , a[2]*dt , a[3]*dt , a[4]*dt ,
105 a[5]*dt , a[6]*dt , a[7]*dt )
108 typedef void result_type;
111 template< class Operations , class Fac , class Time >
112 struct generic_rk_scale_sum< 9 , Operations , Fac , Time > : public Operations::template scale_sum10< Fac , Time >
114 generic_rk_scale_sum( const boost::array<Fac,9> &a , Time dt )
115 : Operations::template scale_sum10< Fac , Time >( 1.0 , a[0]*dt , a[1]*dt , a[2]*dt , a[3]*dt , a[4]*dt ,
116 a[5]*dt , a[6]*dt , a[7]*dt , a[8]*dt )
119 typedef void result_type;
122 template< class Operations , class Fac , class Time >
123 struct generic_rk_scale_sum< 10 , Operations , Fac , Time > : public Operations::template scale_sum11< Fac , Time >
125 generic_rk_scale_sum( const boost::array<Fac,10> &a , Time dt )
126 : Operations::template scale_sum11< Fac , Time >( 1.0 , a[0]*dt , a[1]*dt , a[2]*dt , a[3]*dt , a[4]*dt ,
127 a[5]*dt , a[6]*dt , a[7]*dt , a[8]*dt , a[9]*dt )
130 typedef void result_type;
133 template< class Operations , class Fac , class Time >
134 struct generic_rk_scale_sum< 11 , Operations , Fac , Time > : public Operations::template scale_sum12< Fac , Time >
136 generic_rk_scale_sum( const boost::array<Fac,11> &a , Time dt )
137 : Operations::template scale_sum12< Fac , Time >( 1.0 , a[0]*dt , a[1]*dt , a[2]*dt , a[3]*dt , a[4]*dt ,
138 a[5]*dt , a[6]*dt , a[7]*dt , a[8]*dt , a[9]*dt , a[10]*dt )
141 typedef void result_type;
144 template< class Operations , class Fac , class Time >
145 struct generic_rk_scale_sum< 12 , Operations , Fac , Time > : public Operations::template scale_sum13< Fac , Time >
147 generic_rk_scale_sum( const boost::array<Fac,12> &a , Time dt )
148 : Operations::template scale_sum13< Fac , Time >( 1.0 , a[0]*dt , a[1]*dt , a[2]*dt , a[3]*dt , a[4]*dt ,
149 a[5]*dt , a[6]*dt , a[7]*dt , a[8]*dt , a[9]*dt , a[10]*dt , a[11]*dt )
152 typedef void result_type;
155 template< class Operations , class Fac , class Time >
156 struct generic_rk_scale_sum< 13 , Operations , Fac , Time > : public Operations::template scale_sum14< Fac , Time >
158 generic_rk_scale_sum( const boost::array<Fac,13> &a , Time dt )
159 : Operations::template scale_sum14< Fac , Time >( 1.0 , a[0]*dt , a[1]*dt , a[2]*dt , a[3]*dt , a[4]*dt ,
160 a[5]*dt , a[6]*dt , a[7]*dt , a[8]*dt , a[9]*dt , a[10]*dt , a[11]*dt , a[12]*dt )
163 typedef void result_type;
167 // for error estimates
168 template< size_t StageNumber , class Operations , class Fac , class Time >
169 struct generic_rk_scale_sum_err;
171 template< class Operations , class Fac , class Time >
172 struct generic_rk_scale_sum_err< 1 , Operations , Fac , Time > : public Operations::template scale_sum1< Time >
174 generic_rk_scale_sum_err( const boost::array<Fac,1> &a , Time dt ) : Operations::template scale_sum1< Time >( a[0]*dt )
177 typedef void result_type;
181 template< class Operations , class Fac , class Time >
182 struct generic_rk_scale_sum_err< 2 , Operations , Fac , Time > : public Operations::template scale_sum2< Time >
184 generic_rk_scale_sum_err( const boost::array<Fac,2> &a , Time dt )
185 : Operations::template scale_sum2< Time >( a[0]*dt , a[1]*dt )
188 typedef void result_type;
191 template< class Operations , class Fac , class Time >
192 struct generic_rk_scale_sum_err< 3 , Operations , Fac , Time > : public Operations::template scale_sum3< Time >
194 generic_rk_scale_sum_err( const boost::array<Fac,3> &a , Time dt )
195 : Operations::template scale_sum3< Time >( a[0]*dt , a[1]*dt , a[2]*dt )
198 typedef void result_type;
201 template< class Operations , class Fac , class Time >
202 struct generic_rk_scale_sum_err< 4 , Operations , Fac , Time > : public Operations::template scale_sum4< Time >
204 generic_rk_scale_sum_err( const boost::array<Fac,4> &a , Time dt )
205 : Operations::template scale_sum4< Time >( a[0]*dt , a[1]*dt , a[2]*dt , a[3]*dt )
208 typedef void result_type;
211 template< class Operations , class Fac , class Time >
212 struct generic_rk_scale_sum_err< 5 , Operations , Fac , Time > : public Operations::template scale_sum5< Fac >
214 generic_rk_scale_sum_err( const boost::array<Fac,5> &a , Time dt )
215 : Operations::template scale_sum5< Time >( a[0]*dt , a[1]*dt , a[2]*dt , a[3]*dt , a[4]*dt )
218 typedef void result_type;
222 template< class Operations , class Fac , class Time >
223 struct generic_rk_scale_sum_err< 6 , Operations , Fac , Time > : public Operations::template scale_sum6< Time >
225 generic_rk_scale_sum_err( const boost::array<Fac,6> &a , Time dt )
226 : Operations::template scale_sum6< Time >( a[0]*dt , a[1]*dt , a[2]*dt , a[3]*dt , a[4]*dt , a[5]*dt )
229 typedef void result_type;
234 template< class Operations , class Fac , class Time >
235 struct generic_rk_scale_sum_err< 13 , Operations , Fac , Time > : public Operations::template scale_sum13< Time >
237 generic_rk_scale_sum_err( const boost::array<Fac,13> &a , Time dt )
238 : Operations::template scale_sum13< Time >( a[0]*dt , a[1]*dt , a[2]*dt , a[3]*dt , a[4]*dt , a[5]*dt ,
239 a[6]*dt , a[7]*dt , a[8]*dt , a[9]*dt , a[10]*dt , a[11]*dt , a[12]*dt )
242 typedef void result_type;
252 #endif // BOOST_NUMERIC_ODEINT_STEPPER_DETAIL_GENERIC_RK_OPERATIONS_HPP_INCLUDED