---

rdRKF.h

// rdRKF.h
//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
// Copyright 2000, 2003 Realistic Dynamics, Inc.
// All rights reserved.
//
// Please do not read, copy, distribute, or use without permission.
// Contact: Frank C. Anderson, fca@RealisticDynamics.com
//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
#ifndef __rdRKF_h__
#define __rdRKF_h__

// INCLUDES
#include "rdIntegrator.h"
#include <rdModel.h>
#include <rdControlSet.h>


// EXPORTED CONSTANTS
extern RDSIMULATION_API const int rdRKF_NORMAL;
extern RDSIMULATION_API const int rdRKF_FINE;
extern RDSIMULATION_API const int rdRKF_POOR;
extern RDSIMULATION_API const int rdRKF_NAN;
extern RDSIMULATION_API const int rdRKF_ERROR;
extern RDSIMULATION_API const int rdRKF_TOO_MANY_STEPS;


//=============================================================================
//=============================================================================
class RDSIMULATION_API rdRKF
{
//=============================================================================
// DATA
//=============================================================================
protected:
   double _tol;
   double _tolFine;
   rdModel *_model;
   int _nx;
   double *_x;
   int _ny;
   double *_yv,*_ye,*_dy;
   double *_k1,*_k2,*_k3,*_k4,*_k5,*_k6;

//=============================================================================
// METHODS
//=============================================================================
   //--------------------------------------------------------------------------
   // CONSTRUCTION
   //--------------------------------------------------------------------------
public:
   rdRKF(rdModel *aModel,double aTol=1.0e-4,double aTolFine=-1.0);
   virtual ~rdRKF();

   //--------------------------------------------------------------------------
   // MEMORY
   //--------------------------------------------------------------------------
private:
   int allocateMemory();
   int freeMemory();

   //--------------------------------------------------------------------------
   // GET AND SET
   //--------------------------------------------------------------------------
public:
   // MODEL
   rdModel* getModel();
   // NUMBER OF STATES
   int getNY();
   // TOLERANCE
   void setTolerance(double aTol,double aTolFine=-1.0);
   double getTolerance();
   void setFineTolerance(double aFineTol);
   double getFineTolerance();

   //--------------------------------------------------------------------------
   // INTEGRATION
   //--------------------------------------------------------------------------
   int step(double dt,double t,rdControlSet &x,double *y);
   int step(double dt,double t,double *x,double *y);
   int stepFixed(double dt,double t,rdControlSet &x,double *y);
   int stepFixed(double dt,double t,double *x,double *y);

//=============================================================================
}; // END class rdRKF
//=============================================================================
//=============================================================================

#endif  // __rdRKF_h__

---