ode_sdirk_43.h

Go to the documentation of this file.

#ifndef ODE_SDIRK_43_H_
#define ODE_SDIRK_43_H_
 
 
#include "ode_embedded.h"
#include "ode_rk.h"
#include "ode_newton_bridge.h"
 
namespace ode {
 
//forward declaration to set up Newton class
class OdeSDIRK43;
 
class NewtonSDIRK43 : public OdeNewtonSDIRK<OdeSDIRK43> {
    public:
 
        NewtonSDIRK43 (unsigned long neq, OdeSDIRK43 *integrator) : OdeNewtonSDIRK (neq, integrator) {};
    private:
        void f_Newton (double *x, double *y);
        void J_Newton (double *x, double **J);
};
 
 
class OdeSDIRK43 : public OdeEmbedded, private OdeRK {
    //friends!
    friend class OdeNewtonBridge<OdeSDIRK43>;
    friend class OdeNewtonSDIRK<OdeSDIRK43>;
 
    public:
 
        OdeSDIRK43 (unsigned long neq);
 
        ~OdeSDIRK43 ();
 
        NewtonSDIRK43 *get_newton () { return(newton_); }
 
    private:
        double gam;
        double **a;
        double *b;
        double *d;
        NewtonSDIRK43 *newton_;
        void step_ (double dt);
};
 
} // namespace ode
 
#endif