ceph/src/boost/boost/numeric/odeint/stepper/detail/adams_moulton_coefficients.hpp

   1 /*
   2  [auto_generated]
   3  boost/numeric/odeint/stepper/detail/adams_moulton_coefficients.hpp
   4
   5  [begin_description]
   6  Coefficients for the Adams Moulton method.
   7  [end_description]
   8
   9  Copyright 2011-2012 Karsten Ahnert
  10  Copyright 2011-2012 Mario Mulansky
  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_ADAMS_MOULTON_COEFFICIENTS_HPP_INCLUDED
  19 #define BOOST_NUMERIC_ODEINT_STEPPER_DETAIL_ADAMS_MOULTON_COEFFICIENTS_HPP_INCLUDED
  20
  21
  22 #include <boost/array.hpp>
  23
  24
  25 namespace boost {
  26 namespace numeric {
  27 namespace odeint {
  28 namespace detail {
  29
  30 template< class Value , size_t Steps >
  31 class adams_moulton_coefficients ;
  32
  33 template< class Value >
  34 class adams_moulton_coefficients< Value , 1 > : public boost::array< Value , 1 >
  35 {
  36 public:
  37     adams_moulton_coefficients( void )
  38     : boost::array< Value , 1 >()
  39       {
  40         (*this)[0] = static_cast< Value >( 1 );
  41       }
  42 };
  43
  44
  45 template< class Value >
  46 class adams_moulton_coefficients< Value , 2 > : public boost::array< Value , 2 >
  47 {
  48 public:
  49     adams_moulton_coefficients( void )
  50     : boost::array< Value , 2 >()
  51       {
  52         (*this)[0] = static_cast< Value >( 1 ) / static_cast< Value >( 2 );
  53         (*this)[1] = static_cast< Value >( 1 ) / static_cast< Value >( 2 );
  54       }
  55 };
  56
  57
  58 template< class Value >
  59 class adams_moulton_coefficients< Value , 3 > : public boost::array< Value , 3 >
  60 {
  61 public:
  62     adams_moulton_coefficients( void )
  63     : boost::array< Value , 3 >()
  64       {
  65         (*this)[0] = static_cast< Value >( 5 ) / static_cast< Value >( 12 );
  66         (*this)[1] = static_cast< Value >( 2 ) / static_cast< Value >( 3 );
  67         (*this)[2] = -static_cast< Value >( 1 ) / static_cast< Value >( 12 );
  68       }
  69 };
  70
  71
  72 template< class Value >
  73 class adams_moulton_coefficients< Value , 4 > : public boost::array< Value , 4 >
  74 {
  75 public:
  76     adams_moulton_coefficients( void )
  77     : boost::array< Value , 4 >()
  78       {
  79         (*this)[0] = static_cast< Value >( 3 ) / static_cast< Value >( 8 );
  80         (*this)[1] = static_cast< Value >( 19 ) / static_cast< Value >( 24 );
  81         (*this)[2] = -static_cast< Value >( 5 ) / static_cast< Value >( 24 );
  82         (*this)[3] = static_cast< Value >( 1 ) / static_cast< Value >( 24 );
  83       }
  84 };
  85
  86
  87 template< class Value >
  88 class adams_moulton_coefficients< Value , 5 > : public boost::array< Value , 5 >
  89 {
  90 public:
  91     adams_moulton_coefficients( void )
  92     : boost::array< Value , 5 >()
  93       {
  94         (*this)[0] = static_cast< Value >( 251 ) / static_cast< Value >( 720 );
  95         (*this)[1] = static_cast< Value >( 323 ) / static_cast< Value >( 360 );
  96         (*this)[2] = -static_cast< Value >( 11 ) / static_cast< Value >( 30 );
  97         (*this)[3] = static_cast< Value >( 53 ) / static_cast< Value >( 360 );
  98         (*this)[4] = -static_cast< Value >( 19 ) / static_cast< Value >( 720 );
  99       }
 100 };
 101
 102
 103 template< class Value >
 104 class adams_moulton_coefficients< Value , 6 > : public boost::array< Value , 6 >
 105 {
 106 public:
 107     adams_moulton_coefficients( void )
 108     : boost::array< Value , 6 >()
 109       {
 110         (*this)[0] = static_cast< Value >( 95 ) / static_cast< Value >( 288 );
 111         (*this)[1] = static_cast< Value >( 1427 ) / static_cast< Value >( 1440 );
 112         (*this)[2] = -static_cast< Value >( 133 ) / static_cast< Value >( 240 );
 113         (*this)[3] = static_cast< Value >( 241 ) / static_cast< Value >( 720 );
 114         (*this)[4] = -static_cast< Value >( 173 ) / static_cast< Value >( 1440 );
 115         (*this)[5] = static_cast< Value >( 3 ) / static_cast< Value >( 160 );
 116       }
 117 };
 118
 119 template< class Value >
 120 class adams_moulton_coefficients< Value , 7 > : public boost::array< Value , 7 >
 121 {
 122 public:
 123     adams_moulton_coefficients( void )
 124     : boost::array< Value , 7 >()
 125       {
 126         (*this)[0] = static_cast< Value >( 19087 ) / static_cast< Value >( 60480 );
 127         (*this)[1] = static_cast< Value >( 2713 ) / static_cast< Value >( 2520 );
 128         (*this)[2] = -static_cast< Value >( 15487 ) / static_cast< Value >( 20160 );
 129         (*this)[3] = static_cast< Value >( 586 ) / static_cast< Value >( 945 );
 130         (*this)[4] = -static_cast< Value >( 6737 ) / static_cast< Value >( 20160 );
 131         (*this)[5] = static_cast< Value >( 263 ) / static_cast< Value >( 2520 );
 132         (*this)[6] = -static_cast< Value >( 863 ) / static_cast< Value >( 60480 );
 133       }
 134 };
 135
 136
 137 template< class Value >
 138 class adams_moulton_coefficients< Value , 8 > : public boost::array< Value , 8 >
 139 {
 140 public:
 141     adams_moulton_coefficients( void )
 142     : boost::array< Value , 8 >()
 143       {
 144         (*this)[0] = static_cast< Value >( 5257 ) / static_cast< Value >( 17280 );
 145         (*this)[1] = static_cast< Value >( 139849 ) / static_cast< Value >( 120960 );
 146         (*this)[2] = -static_cast< Value >( 4511 ) / static_cast< Value >( 4480 );
 147         (*this)[3] = static_cast< Value >( 123133 ) / static_cast< Value >( 120960 );
 148         (*this)[4] = -static_cast< Value >( 88547 ) / static_cast< Value >( 120960 );
 149         (*this)[5] = static_cast< Value >( 1537 ) / static_cast< Value >( 4480 );
 150         (*this)[6] = -static_cast< Value >( 11351 ) / static_cast< Value >( 120960 );
 151         (*this)[7] = static_cast< Value >( 275 ) / static_cast< Value >( 24192 );
 152       }
 153 };
 154
 155
 156
 157
 158
 159
 160
 161 } // detail
 162 } // odeint
 163 } // numeric
 164 } // boost
 165
 166
 167
 168 #endif // BOOST_NUMERIC_ODEINT_STEPPER_DETAIL_ADAMS_MOULTON_COEFFICIENTS_HPP_INCLUDED