CLSODA.h

Go to the documentation of this file.

/* Begin CVS Header
   $Source: /Volumes/Home/Users/shoops/cvs/copasi_dev/copasi/odepack++/CLSODA.h,v $
   $Revision: 1.7 $
   $Name:  $
   $Author: gauges $
   $Date: 2006/10/15 07:38:41 $
   End CVS Header */

// Copyright © 2006 by Pedro Mendes, Virginia Tech Intellectual
// Properties, Inc. and EML Research, gGmbH.
// All rights reserved.
//
// This C++ code is based on an f2c conversion of the Fortran
// library ODEPACK available at: http://www.netlib.org/odepack/

#ifndef ODEPACK_CLSODA
#define ODEPACK_CLSODA

#include "odepack++/CInternalSolver.h"

class PJAC;
class SLVS;

class CLSODA : public CInternalSolver
  {
  public:
    CLSODA();
    ~CLSODA();

    C_INT operator() (evalF f,       //  1. evaluate f
                      C_INT *neq,    //  2. number of equations
                      double *y,     //  3. y
                      double *t,     //  4. time
                      double *tout,  //  5. end time
                      C_INT *itol,   //  6. 1 common, 2 atol array, 3 rtol array, 4 both array
                      double *rtol,  //  7. relative tolerance
                      double *atol,  //  8. absolute tolerance
                      C_INT *itask,  //  9. 1
                      C_INT *istate, // 10. status
                      C_INT *iopt,   // 11. 0
                      double *rwork, // 12. work array 22 + neq * max(16, neq + 9)
                      C_INT *lrw,    // 13. size of rwork
                      C_INT *iwork,  // 14. work array 20 + neq
                      C_INT *liw,    // 15. size of iwork
                      evalJ jac,     // 16. evaluate j
                      C_INT *jt);    // 18. type of j evaluation 2 internal full matrix

  private:
    PJAC * mpPJAC;
    SLVS * mpSLVS;

    static const C_INT mxstp0;
    static const C_INT mxhnl0;
    static const C_INT mord[2];
  };

#endif // ODEPACK_CLSODA