1 /*
2 [auto_generated]
3 boost/numeric/odeint/external/openmp/openmp_range_algebra.hpp
4
5 [begin_description]
6 Range algebra for OpenMP.
7 [end_description]
8
9 Copyright 2013 Karsten Ahnert
10 Copyright 2013 Mario Mulansky
11 Copyright 2013 Pascal Germroth
12
13 Distributed under the Boost Software License, Version 1.0.
14 (See accompanying file LICENSE_1_0.txt or
15 copy at http://www.boost.org/LICENSE_1_0.txt)
16 */
17
18
19 #ifndef BOOST_NUMERIC_ODEINT_EXTERNAL_OPENMP_OPENMP_RANGE_ALGEBRA_HPP_INCLUDED
20 #define BOOST_NUMERIC_ODEINT_EXTERNAL_OPENMP_OPENMP_RANGE_ALGEBRA_HPP_INCLUDED
21
22 #include <boost/assert.hpp>
23 #include <boost/range.hpp>
24 #include <boost/numeric/odeint/algebra/norm_result_type.hpp>
25 #include <boost/numeric/odeint/util/n_ary_helper.hpp>
26
27 namespace boost {
28 namespace numeric {
29 namespace odeint {
30
31 /** \brief OpenMP-parallelized range algebra.
32 *
33 * State must be a model of Random Access Range.
34 */
35 struct openmp_range_algebra
36 {
37
38 #if __cplusplus >= 201103L // C++11 supports _Pragma
39
40 #define BOOST_ODEINT_GEN_LOCAL(z, n, unused) \
41 BOOST_ASSERT_MSG( len == boost::size(s ## n), "All state ranges must have the same size." ); \
42 typename boost::range_iterator<S ## n>::type beg ## n = boost::begin(s ## n);
43 #define BOOST_ODEINT_GEN_BODY(n) \
44 const size_t len = boost::size(s0); \
45 BOOST_PP_REPEAT(n, BOOST_ODEINT_GEN_LOCAL, ~) \
46 _Pragma("omp parallel for schedule(runtime)") \
47 for( size_t i = 0 ; i < len ; i++ ) \
48 op( BOOST_PP_ENUM_BINARY_PARAMS(n, beg, [i] BOOST_PP_INTERCEPT) );
BOOST_ODEINT_GEN_FOR_EACHboost::numeric::odeint::openmp_range_algebra49 BOOST_ODEINT_GEN_FOR_EACH(BOOST_ODEINT_GEN_BODY)
50 #undef BOOST_ODEINT_GEN_BODY
51 #undef BOOST_ODEINT_GEN_LOCAL
52
53 #else
54
55 template< class S0 , class Op > static void for_each1 ( S0 &s0 , Op op ) {
56 const size_t len = boost::size(s0);
57 typename boost::range_iterator<S0>::type beg0 = boost::begin(s0);
58 #pragma omp parallel for schedule(runtime)
59 for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] );
60 }
61 template< class S0 , class S1 , class Op > static void for_each2 ( S0 &s0 , S1 &s1 , Op op ) {
62 const size_t len = boost::size(s0);
63 typename boost::range_iterator<S0>::type beg0 = boost::begin(s0);
64 typename boost::range_iterator<S1>::type beg1 = boost::begin(s1);
65 #pragma omp parallel for schedule(runtime)
66 for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] , beg1 [i] );
67 }
68 template< class S0 , class S1 , class S2 , class Op > static void for_each3 ( S0 &s0 , S1 &s1 , S2 &s2 , Op op ) {
69 const size_t len = boost::size(s0);
70 typename boost::range_iterator<S0>::type beg0 = boost::begin(s0);
71 typename boost::range_iterator<S1>::type beg1 = boost::begin(s1);
72 typename boost::range_iterator<S2>::type beg2 = boost::begin(s2);
73 #pragma omp parallel for schedule(runtime)
74 for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] , beg1 [i] , beg2 [i] );
75 }
76 template< class S0 , class S1 , class S2 , class S3 , class Op > static void for_each4 ( S0 &s0 , S1 &s1 , S2 &s2 , S3 &s3 , Op op ) {
77 const size_t len = boost::size(s0);
78 typename boost::range_iterator<S0>::type beg0 = boost::begin(s0);
79 typename boost::range_iterator<S1>::type beg1 = boost::begin(s1);
80 typename boost::range_iterator<S2>::type beg2 = boost::begin(s2);
81 typename boost::range_iterator<S3>::type beg3 = boost::begin(s3);
82 #pragma omp parallel for schedule(runtime)
83 for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] , beg1 [i] , beg2 [i] , beg3 [i] );
84 }
85 template< class S0 , class S1 , class S2 , class S3 , class S4 , class Op > static void for_each5 ( S0 &s0 , S1 &s1 , S2 &s2 , S3 &s3 , S4 &s4 , Op op ) {
86 const size_t len = boost::size(s0);
87 typename boost::range_iterator<S0>::type beg0 = boost::begin(s0);
88 typename boost::range_iterator<S1>::type beg1 = boost::begin(s1);
89 typename boost::range_iterator<S2>::type beg2 = boost::begin(s2);
90 typename boost::range_iterator<S3>::type beg3 = boost::begin(s3);
91 typename boost::range_iterator<S4>::type beg4 = boost::begin(s4);
92 #pragma omp parallel for schedule(runtime)
93 for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] , beg1 [i] , beg2 [i] , beg3 [i] , beg4 [i] );
94 }
95 template< class S0 , class S1 , class S2 , class S3 , class S4 , class S5 , class Op > static void for_each6 ( S0 &s0 , S1 &s1 , S2 &s2 , S3 &s3 , S4 &s4 , S5 &s5 , Op op ) {
96 const size_t len = boost::size(s0);
97 typename boost::range_iterator<S0>::type beg0 = boost::begin(s0);
98 typename boost::range_iterator<S1>::type beg1 = boost::begin(s1);
99 typename boost::range_iterator<S2>::type beg2 = boost::begin(s2);
100 typename boost::range_iterator<S3>::type beg3 = boost::begin(s3);
101 typename boost::range_iterator<S4>::type beg4 = boost::begin(s4);
102 typename boost::range_iterator<S5>::type beg5 = boost::begin(s5);
103 #pragma omp parallel for schedule(runtime)
104 for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] , beg1 [i] , beg2 [i] , beg3 [i] , beg4 [i] , beg5 [i] );
105 }
106 template< class S0 , class S1 , class S2 , class S3 , class S4 , class S5 , class S6 , class Op > static void for_each7 ( S0 &s0 , S1 &s1 , S2 &s2 , S3 &s3 , S4 &s4 , S5 &s5 , S6 &s6 , Op op ) {
107 const size_t len = boost::size(s0);
108 typename boost::range_iterator<S0>::type beg0 = boost::begin(s0);
109 typename boost::range_iterator<S1>::type beg1 = boost::begin(s1);
110 typename boost::range_iterator<S2>::type beg2 = boost::begin(s2);
111 typename boost::range_iterator<S3>::type beg3 = boost::begin(s3);
112 typename boost::range_iterator<S4>::type beg4 = boost::begin(s4);
113 typename boost::range_iterator<S5>::type beg5 = boost::begin(s5);
114 typename boost::range_iterator<S6>::type beg6 = boost::begin(s6);
115 #pragma omp parallel for schedule(runtime)
116 for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] , beg1 [i] , beg2 [i] , beg3 [i] , beg4 [i] , beg5 [i] , beg6 [i] );
117 }
118 template< class S0 , class S1 , class S2 , class S3 , class S4 , class S5 , class S6 , class S7 , class Op > static void for_each8 ( S0 &s0 , S1 &s1 , S2 &s2 , S3 &s3 , S4 &s4 , S5 &s5 , S6 &s6 , S7 &s7 , Op op ) {
119 const size_t len = boost::size(s0);
120 typename boost::range_iterator<S0>::type beg0 = boost::begin(s0);
121 typename boost::range_iterator<S1>::type beg1 = boost::begin(s1);
122 typename boost::range_iterator<S2>::type beg2 = boost::begin(s2);
123 typename boost::range_iterator<S3>::type beg3 = boost::begin(s3);
124 typename boost::range_iterator<S4>::type beg4 = boost::begin(s4);
125 typename boost::range_iterator<S5>::type beg5 = boost::begin(s5);
126 typename boost::range_iterator<S6>::type beg6 = boost::begin(s6);
127 typename boost::range_iterator<S7>::type beg7 = boost::begin(s7);
128 #pragma omp parallel for schedule(runtime)
129 for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] , beg1 [i] , beg2 [i] , beg3 [i] , beg4 [i] , beg5 [i] , beg6 [i] , beg7 [i] );
130 }
131 template< class S0 , class S1 , class S2 , class S3 , class S4 , class S5 , class S6 , class S7 , class S8 , class Op > static void for_each9 ( S0 &s0 , S1 &s1 , S2 &s2 , S3 &s3 , S4 &s4 , S5 &s5 , S6 &s6 , S7 &s7 , S8 &s8 , Op op ) {
132 const size_t len = boost::size(s0);
133 typename boost::range_iterator<S0>::type beg0 = boost::begin(s0);
134 typename boost::range_iterator<S1>::type beg1 = boost::begin(s1);
135 typename boost::range_iterator<S2>::type beg2 = boost::begin(s2);
136 typename boost::range_iterator<S3>::type beg3 = boost::begin(s3);
137 typename boost::range_iterator<S4>::type beg4 = boost::begin(s4);
138 typename boost::range_iterator<S5>::type beg5 = boost::begin(s5);
139 typename boost::range_iterator<S6>::type beg6 = boost::begin(s6);
140 typename boost::range_iterator<S7>::type beg7 = boost::begin(s7);
141 typename boost::range_iterator<S8>::type beg8 = boost::begin(s8);
142 #pragma omp parallel for schedule(runtime)
143 for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] , beg1 [i] , beg2 [i] , beg3 [i] , beg4 [i] , beg5 [i] , beg6 [i] , beg7 [i] , beg8 [i] );
144 }
145 template< class S0 , class S1 , class S2 , class S3 , class S4 , class S5 , class S6 , class S7 , class S8 , class S9 , class Op > static void for_each10 ( S0 &s0 , S1 &s1 , S2 &s2 , S3 &s3 , S4 &s4 , S5 &s5 , S6 &s6 , S7 &s7 , S8 &s8 , S9 &s9 , Op op ) {
146 const size_t len = boost::size(s0);
147 typename boost::range_iterator<S0>::type beg0 = boost::begin(s0);
148 typename boost::range_iterator<S1>::type beg1 = boost::begin(s1);
149 typename boost::range_iterator<S2>::type beg2 = boost::begin(s2);
150 typename boost::range_iterator<S3>::type beg3 = boost::begin(s3);
151 typename boost::range_iterator<S4>::type beg4 = boost::begin(s4);
152 typename boost::range_iterator<S5>::type beg5 = boost::begin(s5);
153 typename boost::range_iterator<S6>::type beg6 = boost::begin(s6);
154 typename boost::range_iterator<S7>::type beg7 = boost::begin(s7);
155 typename boost::range_iterator<S8>::type beg8 = boost::begin(s8);
156 typename boost::range_iterator<S9>::type beg9 = boost::begin(s9);
157 #pragma omp parallel for schedule(runtime)
158 for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] , beg1 [i] , beg2 [i] , beg3 [i] , beg4 [i] , beg5 [i] , beg6 [i] , beg7 [i] , beg8 [i] , beg9 [i] );
159 }
160 template< class S0 , class S1 , class S2 , class S3 , class S4 , class S5 , class S6 , class S7 , class S8 , class S9 , class S10 , class Op > static void for_each11 ( S0 &s0 , S1 &s1 , S2 &s2 , S3 &s3 , S4 &s4 , S5 &s5 , S6 &s6 , S7 &s7 , S8 &s8 , S9 &s9 , S10 &s10 , Op op ) {
161 const size_t len = boost::size(s0);
162 typename boost::range_iterator<S0>::type beg0 = boost::begin(s0);
163 typename boost::range_iterator<S1>::type beg1 = boost::begin(s1);
164 typename boost::range_iterator<S2>::type beg2 = boost::begin(s2);
165 typename boost::range_iterator<S3>::type beg3 = boost::begin(s3);
166 typename boost::range_iterator<S4>::type beg4 = boost::begin(s4);
167 typename boost::range_iterator<S5>::type beg5 = boost::begin(s5);
168 typename boost::range_iterator<S6>::type beg6 = boost::begin(s6);
169 typename boost::range_iterator<S7>::type beg7 = boost::begin(s7);
170 typename boost::range_iterator<S8>::type beg8 = boost::begin(s8);
171 typename boost::range_iterator<S9>::type beg9 = boost::begin(s9);
172 typename boost::range_iterator<S10>::type beg10 = boost::begin(s10);
173 #pragma omp parallel for schedule(runtime)
174 for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] , beg1 [i] , beg2 [i] , beg3 [i] , beg4 [i] , beg5 [i] , beg6 [i] , beg7 [i] , beg8 [i] , beg9 [i] , beg10 [i] );
175 }
176 template< class S0 , class S1 , class S2 , class S3 , class S4 , class S5 , class S6 , class S7 , class S8 , class S9 , class S10 , class S11 , class Op > static void for_each12 ( S0 &s0 , S1 &s1 , S2 &s2 , S3 &s3 , S4 &s4 , S5 &s5 , S6 &s6 , S7 &s7 , S8 &s8 , S9 &s9 , S10 &s10 , S11 &s11 , Op op ) {
177 const size_t len = boost::size(s0);
178 typename boost::range_iterator<S0>::type beg0 = boost::begin(s0);
179 typename boost::range_iterator<S1>::type beg1 = boost::begin(s1);
180 typename boost::range_iterator<S2>::type beg2 = boost::begin(s2);
181 typename boost::range_iterator<S3>::type beg3 = boost::begin(s3);
182 typename boost::range_iterator<S4>::type beg4 = boost::begin(s4);
183 typename boost::range_iterator<S5>::type beg5 = boost::begin(s5);
184 typename boost::range_iterator<S6>::type beg6 = boost::begin(s6);
185 typename boost::range_iterator<S7>::type beg7 = boost::begin(s7);
186 typename boost::range_iterator<S8>::type beg8 = boost::begin(s8);
187 typename boost::range_iterator<S9>::type beg9 = boost::begin(s9);
188 typename boost::range_iterator<S10>::type beg10 = boost::begin(s10);
189 typename boost::range_iterator<S11>::type beg11 = boost::begin(s11);
190 #pragma omp parallel for schedule(runtime)
191 for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] , beg1 [i] , beg2 [i] , beg3 [i] , beg4 [i] , beg5 [i] , beg6 [i] , beg7 [i] , beg8 [i] , beg9 [i] , beg10 [i] , beg11 [i] );
192 }
193 template< class S0 , class S1 , class S2 , class S3 , class S4 , class S5 , class S6 , class S7 , class S8 , class S9 , class S10 , class S11 , class S12 , class Op > static void for_each13 ( S0 &s0 , S1 &s1 , S2 &s2 , S3 &s3 , S4 &s4 , S5 &s5 , S6 &s6 , S7 &s7 , S8 &s8 , S9 &s9 , S10 &s10 , S11 &s11 , S12 &s12 , Op op ) {
194 const size_t len = boost::size(s0);
195 typename boost::range_iterator<S0>::type beg0 = boost::begin(s0);
196 typename boost::range_iterator<S1>::type beg1 = boost::begin(s1);
197 typename boost::range_iterator<S2>::type beg2 = boost::begin(s2);
198 typename boost::range_iterator<S3>::type beg3 = boost::begin(s3);
199 typename boost::range_iterator<S4>::type beg4 = boost::begin(s4);
200 typename boost::range_iterator<S5>::type beg5 = boost::begin(s5);
201 typename boost::range_iterator<S6>::type beg6 = boost::begin(s6);
202 typename boost::range_iterator<S7>::type beg7 = boost::begin(s7);
203 typename boost::range_iterator<S8>::type beg8 = boost::begin(s8);
204 typename boost::range_iterator<S9>::type beg9 = boost::begin(s9);
205 typename boost::range_iterator<S10>::type beg10 = boost::begin(s10);
206 typename boost::range_iterator<S11>::type beg11 = boost::begin(s11);
207 typename boost::range_iterator<S12>::type beg12 = boost::begin(s12);
208 #pragma omp parallel for schedule(runtime)
209 for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] , beg1 [i] , beg2 [i] , beg3 [i] , beg4 [i] , beg5 [i] , beg6 [i] , beg7 [i] , beg8 [i] , beg9 [i] , beg10 [i] , beg11 [i] , beg12 [i] );
210 }
211 template< class S0 , class S1 , class S2 , class S3 , class S4 , class S5 , class S6 , class S7 , class S8 , class S9 , class S10 , class S11 , class S12 , class S13 , class Op > static void for_each14 ( S0 &s0 , S1 &s1 , S2 &s2 , S3 &s3 , S4 &s4 , S5 &s5 , S6 &s6 , S7 &s7 , S8 &s8 , S9 &s9 , S10 &s10 , S11 &s11 , S12 &s12 , S13 &s13 , Op op ) {
212 const size_t len = boost::size(s0);
213 typename boost::range_iterator<S0>::type beg0 = boost::begin(s0);
214 typename boost::range_iterator<S1>::type beg1 = boost::begin(s1);
215 typename boost::range_iterator<S2>::type beg2 = boost::begin(s2);
216 typename boost::range_iterator<S3>::type beg3 = boost::begin(s3);
217 typename boost::range_iterator<S4>::type beg4 = boost::begin(s4);
218 typename boost::range_iterator<S5>::type beg5 = boost::begin(s5);
219 typename boost::range_iterator<S6>::type beg6 = boost::begin(s6);
220 typename boost::range_iterator<S7>::type beg7 = boost::begin(s7);
221 typename boost::range_iterator<S8>::type beg8 = boost::begin(s8);
222 typename boost::range_iterator<S9>::type beg9 = boost::begin(s9);
223 typename boost::range_iterator<S10>::type beg10 = boost::begin(s10);
224 typename boost::range_iterator<S11>::type beg11 = boost::begin(s11);
225 typename boost::range_iterator<S12>::type beg12 = boost::begin(s12);
226 typename boost::range_iterator<S13>::type beg13 = boost::begin(s13);
227 #pragma omp parallel for schedule(runtime)
228 for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] , beg1 [i] , beg2 [i] , beg3 [i] , beg4 [i] , beg5 [i] , beg6 [i] , beg7 [i] , beg8 [i] , beg9 [i] , beg10 [i] , beg11 [i] , beg12 [i] , beg13 [i] );
229 }
230 template< class S0 , class S1 , class S2 , class S3 , class S4 , class S5 , class S6 , class S7 , class S8 , class S9 , class S10 , class S11 , class S12 , class S13 , class S14 , class Op > static void for_each15 ( S0 &s0 , S1 &s1 , S2 &s2 , S3 &s3 , S4 &s4 , S5 &s5 , S6 &s6 , S7 &s7 , S8 &s8 , S9 &s9 , S10 &s10 , S11 &s11 , S12 &s12 , S13 &s13 , S14 &s14 , Op op ) {
231 const size_t len = boost::size(s0);
232 typename boost::range_iterator<S0>::type beg0 = boost::begin(s0);
233 typename boost::range_iterator<S1>::type beg1 = boost::begin(s1);
234 typename boost::range_iterator<S2>::type beg2 = boost::begin(s2);
235 typename boost::range_iterator<S3>::type beg3 = boost::begin(s3);
236 typename boost::range_iterator<S4>::type beg4 = boost::begin(s4);
237 typename boost::range_iterator<S5>::type beg5 = boost::begin(s5);
238 typename boost::range_iterator<S6>::type beg6 = boost::begin(s6);
239 typename boost::range_iterator<S7>::type beg7 = boost::begin(s7);
240 typename boost::range_iterator<S8>::type beg8 = boost::begin(s8);
241 typename boost::range_iterator<S9>::type beg9 = boost::begin(s9);
242 typename boost::range_iterator<S10>::type beg10 = boost::begin(s10);
243 typename boost::range_iterator<S11>::type beg11 = boost::begin(s11);
244 typename boost::range_iterator<S12>::type beg12 = boost::begin(s12);
245 typename boost::range_iterator<S13>::type beg13 = boost::begin(s13);
246 typename boost::range_iterator<S14>::type beg14 = boost::begin(s14);
247 #pragma omp parallel for schedule(runtime)
248 for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] , beg1 [i] , beg2 [i] , beg3 [i] , beg4 [i] , beg5 [i] , beg6 [i] , beg7 [i] , beg8 [i] , beg9 [i] , beg10 [i] , beg11 [i] , beg12 [i] , beg13 [i] , beg14 [i] );
249 }
250
251 #endif
252
253
254 template< class S >
255 static typename norm_result_type< S >::type norm_inf( const S &s )
256 {
257 using std::max;
258 using std::abs;
259 typedef typename norm_result_type< S >::type result_type;
260 result_type init = static_cast< result_type >( 0 );
261 const size_t len = boost::size(s);
262 typename boost::range_iterator<const S>::type beg = boost::begin(s);
263 # pragma omp parallel for reduction(max: init) schedule(dynamic)
264 for( size_t i = 0 ; i < len ; ++i )
265 init = max( init , abs( beg[i] ) );
266 return init;
267 }
268
269 };
270
271
272 }
273 }
274 }
275
276 #endif
277