CVodeSolver.h

Go to the documentation of this file.

 
#pragma once
 
#include <memory>
 
#include "ODESolverTypes.h"
#include "FunctionHandles.h"
 
namespace BaseLib
{
class ConfigTree;
}
 
namespace MathLib
{
namespace ODE
{
 
class CVodeSolverImpl;
 
class CVodeSolver
{
protected:
    CVodeSolver(BaseLib::ConfigTree const& config,
                unsigned const num_equations);
 
    void setTolerance(double const* const abstol, const double reltol);
    void setTolerance(const double abstol, const double reltol);
 
    void setFunction(std::unique_ptr<detail::FunctionHandles>&& f);
 
    void setIC(const double t0, double const* const y0);
 
    void preSolve();
    bool solve(const double t_end);
 
    double const* getSolution() const;
    double getTime() const;
    void getYDot(const double t,
                 double const* const y,
                 double* const y_dot) const;
 
    ~CVodeSolver();
 
private:
    std::unique_ptr<CVodeSolverImpl> impl_;
};
 
 
}  // namespace ODE
}  // namespace MathLib