Integrator.h

Go to the documentation of this file.

/**
 *  @file Integrator.h
 */

/**
 * @defgroup odeGroup ODE Integrators
 */

// Copyright 2001  California Institute of Technology

#ifndef CT_INTEGRATOR_H
#define CT_INTEGRATOR_H
#include "FuncEval.h"

#include "cantera/base/global.h"

namespace Cantera
{

const int DIAG  = 1;
const int DENSE = 2;
const int NOJAC = 4;
const int JAC   = 8;
const int GMRES =16;
const int BAND  =32;

/**
 * Specifies the method used to integrate the system of equations.
 * Not all methods are supported by all integrators.
 */
enum MethodType {
    BDF_Method,          /**< Backward Differentiation */
    Adams_Method         /**< Adams                    */
};

//! Specifies the method used for iteration.
/*!
 * Not all methods are supported by all integrators.
 */
enum IterType {
    //!  Newton Iteration
    Newton_Iter,
    //! Functional Iteration
    Functional_Iter
};


//!  Abstract base class for ODE system integrators.
/*!
 *  @ingroup odeGroup
 */
class Integrator
{
public:
    //! Default Constructor
    Integrator() {
    }

    //!  Destructor
    virtual ~Integrator() {
    }

    /** Set or reset the number of equations. */
    //virtual void resize(int n)=0;

    //!  Set error tolerances.
    /*!
     * @param reltol scalar relative tolerance
     * @param n      Number of equations
     * @param abstol array of N absolute tolerance values
     */
    virtual void setTolerances(doublereal reltol, size_t n,
                               doublereal* abstol) {
        warn("setTolerances");
    }

    //!  Set error tolerances.
    /*!
     * @param reltol scalar relative tolerance
     * @param abstol scalar absolute tolerance
     */
    virtual void setTolerances(doublereal reltol, doublereal abstol) {
        warn("setTolerances");
    }

    //!  Set the sensitivity error tolerances
    /*!
     * @param reltol scalar relative tolerance
     * @param abstol scalar absolute tolerance
     */
    virtual void setSensitivityTolerances(doublereal reltol, doublereal abstol)
    { }

    //!    Set the problem type.
    /*!
     * @param probtype    Type of the problem
     */
    virtual void setProblemType(int probtype) {
        warn("setProblemType");
    }

    /**
     * Initialize the integrator for a new problem. Call after
     * all options have been set.
     * @param t0   initial time
     * @param func RHS evaluator object for system of equations.
     */
    virtual void initialize(doublereal t0, FuncEval& func) {
        warn("initialize");
    }

    virtual void reinitialize(doublereal t0, FuncEval& func) {
        warn("reinitialize");
    }

    //! Integrate the system of equations.
    /*!
     * @param tout Integrate to this time. Note that this is the
     *             absolute time value, not a time interval.
     */
    virtual void integrate(doublereal tout) {
        warn("integrate");
    }

    /**
     * Integrate the system of equations.
     * @param tout integrate to this time. Note that this is the
     * absolute time value, not a time interval.
     */
    virtual doublereal step(doublereal tout) {
        warn("step");
        return 0.0;
    }

    /** The current value of the solution of equation k. */
    virtual doublereal& solution(size_t k) {
        warn("solution");
        return m_dummy;
    }

    /** The current value of the solution of the system of equations. */
    virtual doublereal* solution() {
        warn("solution");
        return 0;
    }

    /** The number of equations. */
    virtual int nEquations() const {
        warn("nEquations");
        return 0;
    }

    /** The number of function evaluations. */
    virtual int nEvals() const {
        warn("nEvals");
        return 0;
    }

    /** Set the maximum integration order that will be used. **/
    virtual void setMaxOrder(int n) {
        warn("setMaxorder");
    }

    /** Set the solution method */
    virtual void setMethod(MethodType t) {
        warn("setMethodType");
    }

    /** Set the linear iterator. */
    virtual void setIterator(IterType t) {
        warn("setInterator");
    }

    /** Set the maximum step size */
    virtual void setMaxStepSize(double hmax) {
        warn("setMaxStepSize");
    }

    /** Set the minimum step size */
    virtual void setMinStepSize(double hmin) {
        warn("setMinStepSize");
    }

    //! Set the maximum permissible number of error test failures
    virtual void setMaxErrTestFails(int n) {
        warn("setMaxErrTestFails");
    }

    virtual void setMaxSteps(int nmax) {
        warn("setMaxStep");
    }

    virtual void setBandwidth(int N_Upper, int N_Lower) {
        warn("setBandwidth");
    }

    virtual int nSensParams() {
        warn("nSensParams()");
        return 0;
    }

    virtual double sensitivity(size_t k, size_t p) {
        warn("sensitivity");
        return 0.0;
    }

private:

    doublereal m_dummy;
    void warn(const std::string& msg) const {
        writelog(">>>> Warning: method "+msg+" of base class "
                 +"Integrator called. Nothing done.\n");
    }

};

// defined in ODE_integrators.cpp
Integrator* newIntegrator(const std::string& itype);

}    // namespace

#endif