ode_rosenbrock.h

Go to the documentation of this file.

#ifndef ODE_ROSENBROCK_H_
#define ODE_ROSENBROCK_H_
 
 
#include "ode_linalg.h"
 
namespace ode {
 
class OdeRosenbrock {
 
    public:
 
        OdeRosenbrock (unsigned long neq, int nk);
        ~OdeRosenbrock ();
 
    protected:
        double gam;
        int *p_;
        double *rhs_;
        double *soltemp_;
        double **k_;
 
        void prep_jac (double **Jac, unsigned long n, double dt, int *p);
 
    private:
        //number of stages
        int nk_;
};
 
} // namespace ode
 
#endif