RKF.H

Go to the documentation of this file.

 
#ifndef _RKF_H_
#define _RKF_H_
 
#include "Types.H"
 
namespace criptic {
 
  namespace rkf {
    
    static constexpr Real ah[] =
      { 1.0 / 4.0, 3.0 / 8.0, 12.0 / 13.0,
    1.0, 1.0 / 2.0 }; 
    static constexpr Real b3[] =
      { 3.0 / 32.0, 9.0 / 32.0 }; 
    static constexpr Real b4[] =
      { 1932.0 / 2197.0, -7200.0 / 2197.0,
    7296.0 / 2197.0 }; 
    static constexpr Real b5[] =
      { 8341.0 / 4104.0, -32832.0 / 4104.0, 29440.0 / 4104.0,
    -845.0 / 4104.0 }; 
    static constexpr Real b6[] =
      { -6080.0 / 20520.0, 41040.0 / 20520.0, -28352.0 / 20520.0,
    9295.0 / 20520.0, -5643.0 / 20520.0
      }; 
    static constexpr Real c1
      = 902880.0 / 7618050.0; 
    static constexpr Real c3
      = 3953664.0 / 7618050.0; 
    static constexpr Real c4
      = 3855735.0 / 7618050.0; 
    static constexpr Real c5
      = -1371249.0 / 7618050.0; 
    static constexpr Real c6
      = 277020.0 / 7618050.0; 
    static constexpr Real ec[] = { 0.0,
      1.0 / 360.0,
      0.0,
      -128.0 / 4275.0,
      -2197.0 / 75240.0,
      1.0 / 50.0,
      2.0 / 55.0
    }; 
  }
 
}
 
#endif
// _RKF_H_