123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263 |
- #ifndef BOOST_NUMERIC_ODEINT_STEPPER_DETAIL_GENERIC_RK_CALL_ALGEBRA_HPP_INCLUDED
- #define BOOST_NUMERIC_ODEINT_STEPPER_DETAIL_GENERIC_RK_CALL_ALGEBRA_HPP_INCLUDED
- namespace boost {
- namespace numeric {
- namespace odeint {
- namespace detail {
- template< size_t StageNumber , class Algebra >
- struct generic_rk_call_algebra
- template< class Algebra >
- struct generic_rk_call_algebra< 1 , Algebra >
- {
- typedef Algebra algebra_type
- template< class S1 , class S2 , class S3 , class S4 , class Op>
- void operator()( algebra_type &algebra , S1 &s1 , S2 &s2 , S3 &s3 , S4 * , Op op ) const
- {
- algebra.for_each3( s1 , s2 , s3 , op )
- }
- template< class S1 , class S2 , class S4 , class Op>
- void operator()( algebra_type &algebra , S1 &s1 , S2 &s2 , S4 * , Op op ) const
- {
- algebra.for_each2( s1 , s2 , op )
- }
- }
- template< class Algebra >
- struct generic_rk_call_algebra< 2 , Algebra >
- {
- template< class S1 , class S2 , class S3 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S3 &s3 , S4 s4_array[1] , Op op ) const
- {
- algebra.for_each4( s1 , s2 , s3 , s4_array[0].m_v , op )
- }
- template< class S1 , class S2 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S4 s4_array[1] , Op op ) const
- {
- algebra.for_each3( s1 , s2 , s4_array[0].m_v , op )
- }
- }
- template< class Algebra >
- struct generic_rk_call_algebra< 3 , Algebra >
- {
- template< class S1 , class S2 , class S3 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S3 &s3 , S4 s4_array[2] , Op op ) const
- {
- algebra.for_each5( s1 , s2 , s3 , s4_array[0].m_v , s4_array[1].m_v , op )
- }
- template< class S1 , class S2 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S4 s4_array[2] , Op op ) const
- {
- algebra.for_each4( s1 , s2 , s4_array[0].m_v , s4_array[1].m_v , op )
- }
- }
- template< class Algebra >
- struct generic_rk_call_algebra< 4 , Algebra >
- {
- template< class S1 , class S2 , class S3 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S3 &s3 , S4 s4_array[3] , Op op ) const
- {
- algebra.for_each6( s1 , s2 , s3 , s4_array[0].m_v , s4_array[1].m_v , s4_array[2].m_v , op )
- }
- template< class S1 , class S2 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S4 s4_array[3] , Op op ) const
- {
- algebra.for_each5( s1 , s2 , s4_array[0].m_v , s4_array[1].m_v , s4_array[2].m_v , op )
- }
- }
- template< class Algebra >
- struct generic_rk_call_algebra< 5 , Algebra >
- {
- template< class S1 , class S2 , class S3 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S3 &s3 , S4 s4_array[4] , Op op ) const
- {
- algebra.for_each7( s1 , s2 , s3 , s4_array[0].m_v , s4_array[1].m_v , s4_array[2].m_v , s4_array[3].m_v , op )
- }
- template< class S1 , class S2 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S4 s4_array[4] , Op op ) const
- {
- algebra.for_each6( s1 , s2 , s4_array[0].m_v , s4_array[1].m_v , s4_array[2].m_v , s4_array[3].m_v , op )
- }
- }
- template< class Algebra >
- struct generic_rk_call_algebra< 6 , Algebra >
- {
- template< class S1 , class S2 , class S3 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S3 &s3 , S4 s4_array[5] , Op op ) const
- {
- algebra.for_each8( s1 , s2 , s3 , s4_array[0].m_v , s4_array[1].m_v , s4_array[2].m_v , s4_array[3].m_v , s4_array[4].m_v , op )
- }
- template< class S1 , class S2 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S4 s4_array[5] , Op op ) const
- {
- algebra.for_each7( s1 , s2 , s4_array[0].m_v , s4_array[1].m_v , s4_array[2].m_v , s4_array[3].m_v , s4_array[4].m_v , op )
- }
- }
- template< class Algebra >
- struct generic_rk_call_algebra< 7 , Algebra >
- {
- template< class S1 , class S2 , class S3 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S3 &s3 , S4 s4_array[6] , Op op ) const
- {
- algebra.for_each9( s1 , s2 , s3 , s4_array[0].m_v , s4_array[1].m_v , s4_array[2].m_v , s4_array[3].m_v , s4_array[4].m_v ,
- s4_array[5].m_v , op )
- }
- template< class S1 , class S2 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S4 s4_array[6] , Op op ) const
- {
- algebra.for_each8( s1 , s2 , s4_array[0].m_v , s4_array[1].m_v , s4_array[2].m_v , s4_array[3].m_v , s4_array[4].m_v ,
- s4_array[5].m_v , op )
- }
- }
- template< class Algebra >
- struct generic_rk_call_algebra< 8 , Algebra >
- {
- template< class S1 , class S2 , class S3 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S3 &s3 , S4 s4_array[7] , Op op ) const
- {
- algebra.for_each10( s1 , s2 , s3 , s4_array[0].m_v , s4_array[1].m_v , s4_array[2].m_v , s4_array[3].m_v , s4_array[4].m_v ,
- s4_array[5].m_v , s4_array[6].m_v , op )
- }
- template< class S1 , class S2 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S4 s4_array[7] , Op op ) const
- {
- algebra.for_each9( s1 , s2 , s4_array[0].m_v , s4_array[1].m_v , s4_array[2].m_v , s4_array[3].m_v , s4_array[4].m_v ,
- s4_array[5].m_v , s4_array[6].m_v , op )
- }
- }
- template< class Algebra >
- struct generic_rk_call_algebra< 9 , Algebra >
- {
- template< class S1 , class S2 , class S3 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S3 &s3 , S4 s4_array[8] , Op op ) const
- {
- algebra.for_each11( s1 , s2 , s3 , s4_array[0].m_v , s4_array[1].m_v , s4_array[2].m_v , s4_array[3].m_v , s4_array[4].m_v ,
- s4_array[5].m_v , s4_array[6].m_v , s4_array[7].m_v , op )
- }
- template< class S1 , class S2 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S4 s4_array[8] , Op op ) const
- {
- algebra.for_each10( s1 , s2 , s4_array[0].m_v , s4_array[1].m_v , s4_array[2].m_v , s4_array[3].m_v , s4_array[4].m_v ,
- s4_array[5].m_v , s4_array[6].m_v , s4_array[7].m_v , op )
- }
- }
- template< class Algebra >
- struct generic_rk_call_algebra< 10 , Algebra >
- {
- template< class S1 , class S2 , class S3 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S3 &s3 , S4 s4_array[9] , Op op ) const
- {
- algebra.for_each12( s1 , s2 , s3 , s4_array[0].m_v , s4_array[1].m_v , s4_array[2].m_v , s4_array[3].m_v , s4_array[4].m_v ,
- s4_array[5].m_v , s4_array[6].m_v , s4_array[7].m_v , s4_array[8].m_v , op )
- }
- template< class S1 , class S2 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S4 s4_array[9] , Op op ) const
- {
- algebra.for_each11( s1 , s2 , s4_array[0].m_v , s4_array[1].m_v , s4_array[2].m_v , s4_array[3].m_v , s4_array[4].m_v ,
- s4_array[5].m_v , s4_array[6].m_v , s4_array[7].m_v , s4_array[8].m_v , op )
- }
- }
- template< class Algebra >
- struct generic_rk_call_algebra< 11 , Algebra >
- {
- template< class S1 , class S2 , class S3 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S3 &s3 , S4 s4_array[10] , Op op ) const
- {
- algebra.for_each13( s1 , s2 , s3 , s4_array[0].m_v , s4_array[1].m_v , s4_array[2].m_v , s4_array[3].m_v , s4_array[4].m_v ,
- s4_array[5].m_v , s4_array[6].m_v , s4_array[7].m_v , s4_array[8].m_v , s4_array[9].m_v , op )
- }
- template< class S1 , class S2 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S4 s4_array[10] , Op op ) const
- {
- algebra.for_each12( s1 , s2 , s4_array[0].m_v , s4_array[1].m_v , s4_array[2].m_v , s4_array[3].m_v , s4_array[4].m_v ,
- s4_array[5].m_v , s4_array[6].m_v , s4_array[7].m_v , s4_array[8].m_v , s4_array[9].m_v , op )
- }
- }
- template< class Algebra >
- struct generic_rk_call_algebra< 12 , Algebra >
- {
- template< class S1 , class S2 , class S3 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S3 &s3 , S4 s4_array[11] , Op op ) const
- {
- algebra.for_each14( s1 , s2 , s3 , s4_array[0].m_v , s4_array[1].m_v , s4_array[2].m_v , s4_array[3].m_v , s4_array[4].m_v ,
- s4_array[5].m_v , s4_array[6].m_v , s4_array[7].m_v , s4_array[8].m_v , s4_array[9].m_v , s4_array[10].m_v , op )
- }
- template< class S1 , class S2 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S4 s4_array[11] , Op op ) const
- {
- algebra.for_each13( s1 , s2 , s4_array[0].m_v , s4_array[1].m_v , s4_array[2].m_v , s4_array[3].m_v , s4_array[4].m_v ,
- s4_array[5].m_v , s4_array[6].m_v , s4_array[7].m_v , s4_array[8].m_v , s4_array[9].m_v , s4_array[10].m_v , op )
- }
- }
- template< class Algebra >
- struct generic_rk_call_algebra< 13 , Algebra >
- {
- template< class S1 , class S2 , class S3 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S3 &s3 , S4 s4_array[12] , Op op ) const
- {
- algebra.for_each15( s1 , s2 , s3 , s4_array[0].m_v , s4_array[1].m_v , s4_array[2].m_v , s4_array[3].m_v , s4_array[4].m_v ,
- s4_array[5].m_v , s4_array[6].m_v , s4_array[7].m_v , s4_array[8].m_v , s4_array[9].m_v , s4_array[10].m_v , s4_array[11].m_v , op )
- }
- template< class S1 , class S2 , class S4 , class Op>
- void operator()( Algebra &algebra , S1 &s1 , S2 &s2 , S4 s4_array[12] , Op op ) const
- {
- algebra.for_each14( s1 , s2 , s4_array[0].m_v , s4_array[1].m_v , s4_array[2].m_v , s4_array[3].m_v , s4_array[4].m_v ,
- s4_array[5].m_v , s4_array[6].m_v , s4_array[7].m_v , s4_array[8].m_v , s4_array[9].m_v , s4_array[10].m_v , s4_array[11].m_v , op )
- }
- }
- }
- }
- }
- }
- #endif // BOOST_NUMERIC_ODEINT_STEPPER_DETAIL_GENERIC_RK_CALL_ALGEBRA_HPP_INCLUDED
|