RKSolver.hxx

Go to the documentation of this file.

/*
  This file is part of MURE,
  Copyright (C) 2007-2021 MURE developers.
 
  MURE is free software: you can redistribute it and/or modify
  it under the terms of the GNU Lesser General Public License as published by
  the Free Software Foundation, either version 3 of the License, or
  (at your option) any later version.
 
  MURE is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  GNU Lesser General Public License for more details.
 
  You should have received a copy of the GNU Lesser General Public License
  along with MURE.  If not, see <http://www.gnu.org/licenses/>.
*/
 
#ifndef _RKSolver_
#define _RKSolver_
#include <cmath>
#include <vector>
#include <iosfwd>
 
using namespace std;
#include "MUREGlobal.hxx"
#include "MiscFunction.hxx"
#include "BatemanSolver.hxx"
//-----------------------------------------------------------------------------//
 
 
//________________________________________________________________________
 
class RKSolver : public BatemanSolver
{
  public:
    RKSolver();   
    RKSolver(const RKSolver &BS);  
    ~RKSolver() override;                   
 
 
    void Solve(double *N0, double **Matrix, int N, double t1, double t2) override;
 
 
 
    void RungeKutta(double *YStart, double t1, double t2);
 
 
    void BuildEqns(double t, double *N, double *dNdt);
 
    void SetTheMatrix(double **f)
    {
        fTheMatrix = f;
    }
 
 
 
 
  protected:
 
    void RK4(double *y, double *dydx, double x, double h, double *yout);
 
    void AdaptStepSize(double *y, double *dydx, double *x, double htry, double eps, double *yscal, double *hdid, double *hnext);
 
    double  fHestimate; 
    double  fHmin;      
    double  fMaxHdid;   
    double  fMinHdid;   
 
  private:
 
    double      **fTheMatrix;       
 
};
#endif