36#ifndef _RUNGEKUTTAFEHLBERGIVPODESOLVER_HPP_
37#define _RUNGEKUTTAFEHLBERGIVPODESOLVER_HPP_
40#include <boost/serialization/base_object.hpp>
42#include "AbstractOneStepIvpOdeSolver.hpp"
58friend class TestRungeKuttaFehlbergIvpOdeSolver;
69 template<
class Archive>
70 void serialize(Archive & archive,
const unsigned int version)
73 archive & boost::serialization::base_object<AbstractIvpOdeSolver>(*
this);
131 std::vector<double>& rCurrentYValues,
132 std::vector<double>& rWorkingMemory,
138 bool outputSolution);
153 std::vector<double>& rCurrentYValues,
154 std::vector<double>& rNextYValues);
167 const double& rError,
168 const double& rTolerance,
169 const double& rMaxTimeStep,
170 const double& rMinTimeStep);
200 std::vector<double>& rYValues,
204 double ignoredSamplingTime);
222 std::vector<double>& rYValues,
gcov doesn't like this file...
#define CHASTE_CLASS_EXPORT(T)
std::vector< double > mk5
std::vector< double > myk3
std::vector< double > mk2
std::vector< double > mError
void AdjustStepSize(double &rCurrentStepSize, const double &rError, const double &rTolerance, const double &rMaxTimeStep, const double &rMinTimeStep)
void InternalSolve(OdeSolution &rSolution, AbstractOdeSystem *pAbstractOdeSystem, std::vector< double > &rCurrentYValues, std::vector< double > &rWorkingMemory, double startTime, double endTime, double maxTimeStep, double minTimeStep, double tolerance, bool outputSolution)
std::vector< double > myk6
std::vector< double > myk4
std::vector< double > mk1
std::vector< double > mk6
OdeSolution Solve(AbstractOdeSystem *pAbstractOdeSystem, std::vector< double > &rYValues, double startTime, double endTime, double timeStep, double ignoredSamplingTime)
std::vector< double > myk5
std::vector< double > mk3
std::vector< double > myk2
std::vector< double > mk4
void serialize(Archive &archive, const unsigned int version)
RungeKuttaFehlbergIvpOdeSolver()
friend class boost::serialization::access
void CalculateNextYValue(AbstractOdeSystem *pAbstractOdeSystem, double timeStep, double time, std::vector< double > &rCurrentYValues, std::vector< double > &rNextYValues)