ode_radau_iia_5.h

Go to the documentation of this file.

#ifndef ODE_RADAU_IIA_5_H_
#define ODE_RADAU_IIA_5_H_
 
 
#include "ode_adaptive.h"
#include "ode_irk.h"
#include "ode_newton_bridge.h"
 
namespace ode {
 
//forward declaration to set up Newton class
class OdeRadauIIA5;
 
class NewtonRadauIIA5 : public OdeNewtonIRK<OdeRadauIIA5> {
    public:
 
        NewtonRadauIIA5 (unsigned long neq, unsigned long nnew, OdeRadauIIA5 *integrator) : OdeNewtonIRK (neq, nnew, integrator) {};
    private:
        void f_Newton (double *x, double *y);
        void J_Newton (double *x, double **J);
};
 
 
class OdeRadauIIA5 : public OdeAdaptive, private OdeIRK {
    //friends!
    friend class OdeNewtonBridge<OdeRadauIIA5>;
    friend class OdeNewtonIRK<OdeRadauIIA5>;
 
    public:
 
        OdeRadauIIA5 (unsigned long neq);
 
        ~OdeRadauIIA5 ();
 
        NewtonRadauIIA5 *get_newton () { return(newton_); }
 
    private:
        double **a;
        double *b;
        NewtonRadauIIA5 *newton_;
        void step_ (double dt);
};
 
} // namespace ode
 
#endif