1 /*
2  [auto_generated]
3  boost/numeric/odeint/stepper/detail/generic_rk_operations.hpp
4 
5  [begin_description]
6  Operations caller for the generic Runge Kutta method.
7  [end_description]
8 
9  Copyright 2011 Mario Mulansky
10  Copyright 2011-2012 Karsten Ahnert
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 #ifndef BOOST_NUMERIC_ODEINT_STEPPER_DETAIL_GENERIC_RK_OPERATIONS_HPP_INCLUDED
19 #define BOOST_NUMERIC_ODEINT_STEPPER_DETAIL_GENERIC_RK_OPERATIONS_HPP_INCLUDED
20 
21 
22 namespace boost {
23 namespace numeric {
24 namespace odeint {
25 namespace detail {
26 
27 template< size_t StageNumber , class Operations , class Fac , class Time >
28 struct generic_rk_scale_sum;
29 
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 >
32 {
33     generic_rk_scale_sum( const boost::array<Fac,1> &a , Time dt ) : Operations::template scale_sum2< Fac , Time >( 1.0 , a[0]*dt )
34                 { }
35 
36     typedef void result_type;
37 };
38 
39 
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 >
42 {
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 )
45                       { }
46 
47     typedef void result_type;
48 };
49 
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 >
52 {
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 )
55                       { }
56 
57     typedef void result_type;
58 };
59 
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 >
62 {
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 )
65                       { }
66 
67     typedef void result_type;
68 };
69 
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 >
72 {
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 )
75                       { }
76 
77     typedef void result_type;
78 };
79 
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 >
82 {
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 )
85                       { }
86 
87     typedef void result_type;
88 };
89 
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 >
92 {
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 )
95                       { }
96 
97     typedef void result_type;
98 };
99 
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 >
102 {
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 )
106                       { }
107 
108     typedef void result_type;
109 };
110 
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 >
113 {
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 )
117                       { }
118 
119     typedef void result_type;
120 };
121 
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 >
124 {
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 )
128                       { }
129 
130     typedef void result_type;
131 };
132 
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 >
135 {
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 )
139                       { }
140 
141     typedef void result_type;
142 };
143 
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 >
146 {
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 )
150                       { }
151 
152     typedef void result_type;
153 };
154 
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 >
157 {
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  )
161                       { }
162 
163     typedef void result_type;
164 };
165 
166 
167 // for error estimates
168 template< size_t StageNumber , class Operations , class Fac , class Time >
169 struct generic_rk_scale_sum_err;
170 
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 >
173 {
174     generic_rk_scale_sum_err( const boost::array<Fac,1> &a , Time dt ) : Operations::template scale_sum1< Time >( a[0]*dt )
175                 { }
176 
177     typedef void result_type;
178 };
179 
180 
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 >
183 {
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 )
186                       { }
187 
188     typedef void result_type;
189 };
190 
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 >
193 {
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 )
196                       { }
197 
198     typedef void result_type;
199 };
200 
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 >
203 {
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 )
206                       { }
207 
208     typedef void result_type;
209 };
210 
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 >
213 {
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 )
216                       { }
217 
218     typedef void result_type;
219 };
220 
221 
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 >
224 {
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 )
227                       { }
228 
229     typedef void result_type;
230 };
231 
232 
233 // for rk87
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 >
236 {
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 )
240                       { }
241 
242     typedef void result_type;
243 };
244 
245 
246 }
247 }
248 }
249 }
250 
251 
252 #endif // BOOST_NUMERIC_ODEINT_STEPPER_DETAIL_GENERIC_RK_OPERATIONS_HPP_INCLUDED
253