ode_grk4a.h

Go to the documentation of this file.

#ifndef ODE_GRK4A_H_
#define ODE_GRK4A_H_
 
 
#include "ode_linalg.h"
#include "ode_embedded.h"
#include "ode_rosenbrock.h"
 
namespace ode {
 
 
class OdeGRK4A : public OdeEmbedded, private OdeRosenbrock {
 
    public:
 
        OdeGRK4A (unsigned long neq);
 
    private:
        //function for taking a single step
        void step_ (double dt);
        //coefficients
        double gam21,
               gam31, gam32,
               gam41, gam42, gam43;
        double g21,
               g31, g32,
               g41, g42, g43;
        double alp21,
               alp31, alp32;
        double b1, b2, b3, b4,
               d1, d2, d3;
 
};
 
} // namespace ode
 
#endif