ceph/src/boost/boost/numeric/odeint/external/openmp/openmp_range_algebra.hpp

   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) );
  49 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