16 std::vector<double> & errors)
const {
17 const unsigned int nvar = s.
variable.size();
22 d1.time = s.
time + dt/2.0;
44 if (h<=0)
throw std::runtime_error (
"SimpleRKStepper: negative stepsize");
45 const unsigned int nvar = s.
variable.size();
48 std::vector<std::vector<double> >k(tableau.
nSteps());
49 for (
unsigned int i=0;i<tableau.
nSteps();i++) {
52 for (
unsigned int v=0;v<nvar;v++) arg[v]=s.
variable[v];
53 for (
unsigned int j=0;j<i;j++) {
54 for (
unsigned int v=0;v<nvar;v++) arg[v] += h*tableau.
A(i,j)*k[j][v];
56 for (
unsigned int v=0;v<nvar;v++) k[i][v]=(*data->
_diffEqn[v])(arg);
62 for (
unsigned int i=0;i<tableau.
nSteps();i++) {
63 for (
unsigned int v=0;v<nvar;v++) d.
firstDerivative[v] += tableau.
b(i)*k[i][v];
77 return tableau.
order();
double & A(unsigned int i, unsigned int j)
unsigned int nSteps() const
unsigned int order() const
double & b(unsigned int i)
std::vector< const AbsFunction * > _diffEqn
virtual void step(const RKIntegrator::RKData *data, const RKIntegrator::RKData::Data &sdata, RKIntegrator::RKData::Data &ddata, std::vector< double > &errors) const
virtual unsigned int order() const
void doStep(const RKIntegrator::RKData *data, const RKIntegrator::RKData::Data &s, RKIntegrator::RKData::Data &d) const
virtual StepDoublingRKStepper * clone() const
StepDoublingRKStepper(const ButcherTableau &tableau)
virtual ~StepDoublingRKStepper()
std::vector< double > firstDerivative
std::vector< double > variable