ode_embedded.h

Go to the documentation of this file.

#ifndef ODE_EMBEDDED_H_
#define ODE_EMBEDDED_H_
 
 
#include <cmath>
 
#include "ode_io.h"
#include "ode_util.h"
#include "ode_adaptive.h"
 
namespace ode {
 
class OdeEmbedded : public OdeAdaptive {
 
    public:
 
        OdeEmbedded (unsigned long neq, bool need_jac, int lowerord);
        virtual ~OdeEmbedded ();
 
        //-------------------
        //getters and setters
 
        double get_facsafe () { return(facsafe_); }
        double get_facmin () { return(facmin_); }
        double get_facmax () { return(facmax_); }
 
        void set_facsafe (double facsafe) { facsafe_ = facsafe; }
        void set_facmin (double facmin) {
            if ( facmin >= 1.0 )
                ode_print_exit("facmin must be a number less than 1 in set_facmin()");
            facmin_ = facmin;
        }
        void set_facmax (double facmax) {
            if ( facmax <= 1.0 )
                ode_print_exit("facmax must be a number greater than 1 in set_facmax()");
            facmax_ = facmax;
        }
 
    protected:
 
        double facsafe_;
        double facmin_;
        double facmax_;
 
        double *solemb_;
        double error (double abstol, double reltol);
        double facopt (double err);
 
        //------------------------------------------------
        //implementation of inherited methods for adapting
 
        virtual void adapt (double abstol, double reltol);
        virtual bool is_rejected ();
        virtual double dt_adapt ();
 
    private:
        int lowerord_;
        bool isrej_;
        double dtopt_;
};
 
} // namespace ode
 
#endif