diff options
Diffstat (limited to 'boost/numeric/odeint/external/openmp/openmp_range_algebra.hpp')
-rw-r--r-- | boost/numeric/odeint/external/openmp/openmp_range_algebra.hpp | 276 |
1 files changed, 276 insertions, 0 deletions
diff --git a/boost/numeric/odeint/external/openmp/openmp_range_algebra.hpp b/boost/numeric/odeint/external/openmp/openmp_range_algebra.hpp new file mode 100644 index 0000000000..f5668e4d0b --- /dev/null +++ b/boost/numeric/odeint/external/openmp/openmp_range_algebra.hpp @@ -0,0 +1,276 @@ +/* + [auto_generated] + boost/numeric/odeint/external/openmp/openmp_range_algebra.hpp + + [begin_description] + Range algebra for OpenMP. + [end_description] + + Copyright 2013 Karsten Ahnert + Copyright 2013 Mario Mulansky + Copyright 2013 Pascal Germroth + + Distributed under the Boost Software License, Version 1.0. + (See accompanying file LICENSE_1_0.txt or + copy at http://www.boost.org/LICENSE_1_0.txt) + */ + + +#ifndef BOOST_NUMERIC_ODEINT_EXTERNAL_OPENMP_OPENMP_RANGE_ALGEBRA_HPP_INCLUDED +#define BOOST_NUMERIC_ODEINT_EXTERNAL_OPENMP_OPENMP_RANGE_ALGEBRA_HPP_INCLUDED + +#include <boost/assert.hpp> +#include <boost/range.hpp> +#include <boost/numeric/odeint/algebra/norm_result_type.hpp> +#include <boost/numeric/odeint/util/n_ary_helper.hpp> + +namespace boost { +namespace numeric { +namespace odeint { + +/** \brief OpenMP-parallelized range algebra. + * + * State must be a model of Random Access Range. + */ +struct openmp_range_algebra +{ + +#if __cplusplus >= 201103L // C++11 supports _Pragma + +#define BOOST_ODEINT_GEN_LOCAL(z, n, unused) \ + BOOST_ASSERT_MSG( len == boost::size(s ## n), "All state ranges must have the same size." ); \ + typename boost::range_iterator<S ## n>::type beg ## n = boost::begin(s ## n); +#define BOOST_ODEINT_GEN_BODY(n) \ + const size_t len = boost::size(s0); \ + BOOST_PP_REPEAT(n, BOOST_ODEINT_GEN_LOCAL, ~) \ + _Pragma("omp parallel for schedule(runtime)") \ + for( size_t i = 0 ; i < len ; i++ ) \ + op( BOOST_PP_ENUM_BINARY_PARAMS(n, beg, [i] BOOST_PP_INTERCEPT) ); +BOOST_ODEINT_GEN_FOR_EACH(BOOST_ODEINT_GEN_BODY) +#undef BOOST_ODEINT_GEN_BODY +#undef BOOST_ODEINT_GEN_LOCAL + +#else + + template< class S0 , class Op > static void for_each1 ( S0 &s0 , Op op ) { + const size_t len = boost::size(s0); + typename boost::range_iterator<S0>::type beg0 = boost::begin(s0); + #pragma omp parallel for schedule(runtime) + for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] ); + } + template< class S0 , class S1 , class Op > static void for_each2 ( S0 &s0 , S1 &s1 , Op op ) { + const size_t len = boost::size(s0); + typename boost::range_iterator<S0>::type beg0 = boost::begin(s0); + typename boost::range_iterator<S1>::type beg1 = boost::begin(s1); + #pragma omp parallel for schedule(runtime) + for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] , beg1 [i] ); + } + template< class S0 , class S1 , class S2 , class Op > static void for_each3 ( S0 &s0 , S1 &s1 , S2 &s2 , Op op ) { + const size_t len = boost::size(s0); + typename boost::range_iterator<S0>::type beg0 = boost::begin(s0); + typename boost::range_iterator<S1>::type beg1 = boost::begin(s1); + typename boost::range_iterator<S2>::type beg2 = boost::begin(s2); + #pragma omp parallel for schedule(runtime) + for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] , beg1 [i] , beg2 [i] ); + } + 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 ) { + const size_t len = boost::size(s0); + typename boost::range_iterator<S0>::type beg0 = boost::begin(s0); + typename boost::range_iterator<S1>::type beg1 = boost::begin(s1); + typename boost::range_iterator<S2>::type beg2 = boost::begin(s2); + typename boost::range_iterator<S3>::type beg3 = boost::begin(s3); + #pragma omp parallel for schedule(runtime) + for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] , beg1 [i] , beg2 [i] , beg3 [i] ); + } + 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 ) { + const size_t len = boost::size(s0); + typename boost::range_iterator<S0>::type beg0 = boost::begin(s0); + typename boost::range_iterator<S1>::type beg1 = boost::begin(s1); + typename boost::range_iterator<S2>::type beg2 = boost::begin(s2); + typename boost::range_iterator<S3>::type beg3 = boost::begin(s3); + typename boost::range_iterator<S4>::type beg4 = boost::begin(s4); + #pragma omp parallel for schedule(runtime) + for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] , beg1 [i] , beg2 [i] , beg3 [i] , beg4 [i] ); + } + 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 ) { + const size_t len = boost::size(s0); + typename boost::range_iterator<S0>::type beg0 = boost::begin(s0); + typename boost::range_iterator<S1>::type beg1 = boost::begin(s1); + typename boost::range_iterator<S2>::type beg2 = boost::begin(s2); + typename boost::range_iterator<S3>::type beg3 = boost::begin(s3); + typename boost::range_iterator<S4>::type beg4 = boost::begin(s4); + typename boost::range_iterator<S5>::type beg5 = boost::begin(s5); + #pragma omp parallel for schedule(runtime) + for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] , beg1 [i] , beg2 [i] , beg3 [i] , beg4 [i] , beg5 [i] ); + } + 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 ) { + const size_t len = boost::size(s0); + typename boost::range_iterator<S0>::type beg0 = boost::begin(s0); + typename boost::range_iterator<S1>::type beg1 = boost::begin(s1); + typename boost::range_iterator<S2>::type beg2 = boost::begin(s2); + typename boost::range_iterator<S3>::type beg3 = boost::begin(s3); + typename boost::range_iterator<S4>::type beg4 = boost::begin(s4); + typename boost::range_iterator<S5>::type beg5 = boost::begin(s5); + typename boost::range_iterator<S6>::type beg6 = boost::begin(s6); + #pragma omp parallel for schedule(runtime) + for( size_t i = 0 ; i < len ; i++ ) op( beg0 [i] , beg1 [i] , beg2 [i] , beg3 [i] , beg4 [i] , beg5 [i] , beg6 [i] ); + } + 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 ) { + const size_t len = boost::size(s0); + typename boost::range_iterator<S0>::type beg0 = boost::begin(s0); + typename boost::range_iterator<S1>::type beg1 = boost::begin(s1); + typename boost::range_iterator<S2>::type beg2 = boost::begin(s2); + typename boost::range_iterator<S3>::type beg3 = boost::begin(s3); + typename boost::range_iterator<S4>::type beg4 = boost::begin(s4); + typename boost::range_iterator<S5>::type beg5 = boost::begin(s5); + typename boost::range_iterator<S6>::type beg6 = boost::begin(s6); + typename boost::range_iterator<S7>::type beg7 = boost::begin(s7); + #pragma omp parallel for schedule(runtime) + 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] ); + } + 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 ) { + const size_t len = boost::size(s0); + typename boost::range_iterator<S0>::type beg0 = boost::begin(s0); + typename boost::range_iterator<S1>::type beg1 = boost::begin(s1); + typename boost::range_iterator<S2>::type beg2 = boost::begin(s2); + typename boost::range_iterator<S3>::type beg3 = boost::begin(s3); + typename boost::range_iterator<S4>::type beg4 = boost::begin(s4); + typename boost::range_iterator<S5>::type beg5 = boost::begin(s5); + typename boost::range_iterator<S6>::type beg6 = boost::begin(s6); + typename boost::range_iterator<S7>::type beg7 = boost::begin(s7); + typename boost::range_iterator<S8>::type beg8 = boost::begin(s8); + #pragma omp parallel for schedule(runtime) + 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] ); + } + 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 ) { + const size_t len = boost::size(s0); + typename boost::range_iterator<S0>::type beg0 = boost::begin(s0); + typename boost::range_iterator<S1>::type beg1 = boost::begin(s1); + typename boost::range_iterator<S2>::type beg2 = boost::begin(s2); + typename boost::range_iterator<S3>::type beg3 = boost::begin(s3); + typename boost::range_iterator<S4>::type beg4 = boost::begin(s4); + typename boost::range_iterator<S5>::type beg5 = boost::begin(s5); + typename boost::range_iterator<S6>::type beg6 = boost::begin(s6); + typename boost::range_iterator<S7>::type beg7 = boost::begin(s7); + typename boost::range_iterator<S8>::type beg8 = boost::begin(s8); + typename boost::range_iterator<S9>::type beg9 = boost::begin(s9); + #pragma omp parallel for schedule(runtime) + 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] ); + } + 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 ) { + const size_t len = boost::size(s0); + typename boost::range_iterator<S0>::type beg0 = boost::begin(s0); + typename boost::range_iterator<S1>::type beg1 = boost::begin(s1); + typename boost::range_iterator<S2>::type beg2 = boost::begin(s2); + typename boost::range_iterator<S3>::type beg3 = boost::begin(s3); + typename boost::range_iterator<S4>::type beg4 = boost::begin(s4); + typename boost::range_iterator<S5>::type beg5 = boost::begin(s5); + typename boost::range_iterator<S6>::type beg6 = boost::begin(s6); + typename boost::range_iterator<S7>::type beg7 = boost::begin(s7); + typename boost::range_iterator<S8>::type beg8 = boost::begin(s8); + typename boost::range_iterator<S9>::type beg9 = boost::begin(s9); + typename boost::range_iterator<S10>::type beg10 = boost::begin(s10); + #pragma omp parallel for schedule(runtime) + 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] ); + } + 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 ) { + const size_t len = boost::size(s0); + typename boost::range_iterator<S0>::type beg0 = boost::begin(s0); + typename boost::range_iterator<S1>::type beg1 = boost::begin(s1); + typename boost::range_iterator<S2>::type beg2 = boost::begin(s2); + typename boost::range_iterator<S3>::type beg3 = boost::begin(s3); + typename boost::range_iterator<S4>::type beg4 = boost::begin(s4); + typename boost::range_iterator<S5>::type beg5 = boost::begin(s5); + typename boost::range_iterator<S6>::type beg6 = boost::begin(s6); + typename boost::range_iterator<S7>::type beg7 = boost::begin(s7); + typename boost::range_iterator<S8>::type beg8 = boost::begin(s8); + typename boost::range_iterator<S9>::type beg9 = boost::begin(s9); + typename boost::range_iterator<S10>::type beg10 = boost::begin(s10); + typename boost::range_iterator<S11>::type beg11 = boost::begin(s11); + #pragma omp parallel for schedule(runtime) + 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] ); + } + 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 ) { + const size_t len = boost::size(s0); + typename boost::range_iterator<S0>::type beg0 = boost::begin(s0); + typename boost::range_iterator<S1>::type beg1 = boost::begin(s1); + typename boost::range_iterator<S2>::type beg2 = boost::begin(s2); + typename boost::range_iterator<S3>::type beg3 = boost::begin(s3); + typename boost::range_iterator<S4>::type beg4 = boost::begin(s4); + typename boost::range_iterator<S5>::type beg5 = boost::begin(s5); + typename boost::range_iterator<S6>::type beg6 = boost::begin(s6); + typename boost::range_iterator<S7>::type beg7 = boost::begin(s7); + typename boost::range_iterator<S8>::type beg8 = boost::begin(s8); + typename boost::range_iterator<S9>::type beg9 = boost::begin(s9); + typename boost::range_iterator<S10>::type beg10 = boost::begin(s10); + typename boost::range_iterator<S11>::type beg11 = boost::begin(s11); + typename boost::range_iterator<S12>::type beg12 = boost::begin(s12); + #pragma omp parallel for schedule(runtime) + 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] ); + } + 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 ) { + const size_t len = boost::size(s0); + typename boost::range_iterator<S0>::type beg0 = boost::begin(s0); + typename boost::range_iterator<S1>::type beg1 = boost::begin(s1); + typename boost::range_iterator<S2>::type beg2 = boost::begin(s2); + typename boost::range_iterator<S3>::type beg3 = boost::begin(s3); + typename boost::range_iterator<S4>::type beg4 = boost::begin(s4); + typename boost::range_iterator<S5>::type beg5 = boost::begin(s5); + typename boost::range_iterator<S6>::type beg6 = boost::begin(s6); + typename boost::range_iterator<S7>::type beg7 = boost::begin(s7); + typename boost::range_iterator<S8>::type beg8 = boost::begin(s8); + typename boost::range_iterator<S9>::type beg9 = boost::begin(s9); + typename boost::range_iterator<S10>::type beg10 = boost::begin(s10); + typename boost::range_iterator<S11>::type beg11 = boost::begin(s11); + typename boost::range_iterator<S12>::type beg12 = boost::begin(s12); + typename boost::range_iterator<S13>::type beg13 = boost::begin(s13); + #pragma omp parallel for schedule(runtime) + 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] ); + } + 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 ) { + const size_t len = boost::size(s0); + typename boost::range_iterator<S0>::type beg0 = boost::begin(s0); + typename boost::range_iterator<S1>::type beg1 = boost::begin(s1); + typename boost::range_iterator<S2>::type beg2 = boost::begin(s2); + typename boost::range_iterator<S3>::type beg3 = boost::begin(s3); + typename boost::range_iterator<S4>::type beg4 = boost::begin(s4); + typename boost::range_iterator<S5>::type beg5 = boost::begin(s5); + typename boost::range_iterator<S6>::type beg6 = boost::begin(s6); + typename boost::range_iterator<S7>::type beg7 = boost::begin(s7); + typename boost::range_iterator<S8>::type beg8 = boost::begin(s8); + typename boost::range_iterator<S9>::type beg9 = boost::begin(s9); + typename boost::range_iterator<S10>::type beg10 = boost::begin(s10); + typename boost::range_iterator<S11>::type beg11 = boost::begin(s11); + typename boost::range_iterator<S12>::type beg12 = boost::begin(s12); + typename boost::range_iterator<S13>::type beg13 = boost::begin(s13); + typename boost::range_iterator<S14>::type beg14 = boost::begin(s14); + #pragma omp parallel for schedule(runtime) + 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] ); + } + +#endif + + + template< class S > + static typename norm_result_type< S >::type norm_inf( const S &s ) + { + using std::max; + using std::abs; + typedef typename norm_result_type< S >::type result_type; + result_type init = static_cast< result_type >( 0 ); + const size_t len = boost::size(s); + typename boost::range_iterator<S>::type beg = boost::begin(s); +# pragma omp parallel for reduction(max: init) schedule(dynamic) + for( size_t i = 0 ; i < len ; ++i ) + init = max( init , abs( beg[i] ) ); + return init; + } + +}; + + +} +} +} + +#endif |