BatemanSolver.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 _BatemanSolver_
#define _BatemanSolver_
#include <iosfwd>
using namespace std;
 
//-----------------------------------------------------------------------------//
 
// matrix exponential method
//________________________________________________________________________
 
class BatemanSolver
{
  public:
    BatemanSolver();                            
    BatemanSolver(const BatemanSolver &BS);     
    virtual ~BatemanSolver();                   
 
 
    virtual void Solve(double *N0, double **Matrix, int N, double t1, double t2) = 0;
 
    int GetNumberOfEquationSize()
    {
        return fNVar;   
    }
    void SetNumberOfEquations(int N)
    {
        fNVar = N;
    }
 
    void SetPrecision(double eps = 1e-5)
    {
        fPrecision = eps;   
    }
 
 
    void SetForbidNegativeValue()
    {
        fIsNegativeValueAllowed = false;   
    }
 
 
    void SetXe135EquilibriumIndex(int Te135, int I135, int Xe135)
    {
        fIdxTe135 = Te135;
        fIdxI135 = I135;
        fIdxXe135 = Xe135;
    }
 
  protected:
    int     fNVar;      
    double  fPrecision; 
    double  fTmin;      
    bool    fIsNegativeValueAllowed;
    int     fIdxTe135;  
    int     fIdxI135;   
    int     fIdxXe135;  
};
#endif