ode_adaptive.h

Go to the documentation of this file.

#ifndef ODE_ADAPTIVE_H_
#define ODE_ADAPTIVE_H_
 
 
#include <cstdio>
#include <cmath>
#include <vector>
#include <cstring>
#include <string>
 
#include "ode_io.h"
#include "ode_util.h"
#include "ode_base.h"
 
namespace ode {
 
 
class OdeAdaptive : public OdeBase {
 
    public:
 
        OdeAdaptive (unsigned long neq, bool need_jac);
        virtual ~OdeAdaptive ();
 
        //-------------------
        //getters and setters
 
        long unsigned get_nrej () { return(nrej_); }
        double get_abstol () { return(abstol_); }
        double get_reltol () { return(reltol_); }
        double get_dtmax () { return(dtmax_); }
 
        void set_abstol (double tol) { abstol_ = tol; };
        void set_reltol (double tol) { reltol_ = tol; };
        void set_tol (double tol) { abstol_ = tol; reltol_ = tol; };
        void set_dtmax (double dtmax) { dtmax_ = dtmax; }
 
        //---------------------------------------
        //integration with adaptive time stepping
 
 
        void solve_adaptive (double tint, double dt0, bool extras=true);
 
 
        void solve_adaptive (double tint, double dt0, const char *dirout, int inter);
 
 
        void solve_adaptive (double tint, double dt0, unsigned long nsnap, const char *dirout);
 
 
        void solve_adaptive (double dt0, double *tsnap, unsigned long nsnap, const char *dirout);
 
    protected:
 
 
        void solve_adaptive_ (double tint, double dt0, bool extra=true);
 
        //------------------------
        //basic adaptive variables
 
        long unsigned nrej_;
        double abstol_;
        double reltol_;
        double dtmax_;
 
        virtual void adapt (double abstol, double reltol);
        virtual bool is_rejected ();
        virtual double dt_adapt ();
 
        bool solve_done_adaptive (double tend);
 
        //wrappers
        bool step_adaptive_ (double dt, bool extra=true);
        double dt_adapt_ (double tend);
 
    private:
        double *solprev_;
};
 
} // namespace ode
 
#endif