Go to the documentation of this file.
16 std::vector<double> & errors)
const {
17 const unsigned int nvar =
s.variable.size();
21 double dt = (d.
time -
s.time);
22 d1.time =
s.time + dt/2.0;
41 double h = (d.
time -
s.time);
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];
77 return tableau.
order();
double & b(unsigned int i)
virtual unsigned int order() const
double & A(unsigned int i, unsigned int j)
virtual void step(const RKIntegrator::RKData *data, const RKIntegrator::RKData::Data &sdata, RKIntegrator::RKData::Data &ddata, std::vector< double > &errors) const
virtual ~StepDoublingRKStepper()
std::vector< double > variable
they are gone ZOOM Features Discontinued The following features of the ZOOM package were felt to be extreme overkill These have been after checking that no existing user code was utilizing as in SpaceVector v
unsigned int order() const
std::vector< const AbsFunction * > _diffEqn
std::vector< double > firstDerivative
Methods applicble to containers of as in std::list< LorentzVector > s
StepDoublingRKStepper(const ButcherTableau &tableau)
virtual StepDoublingRKStepper * clone() const
unsigned int nSteps() const
void doStep(const RKIntegrator::RKData *data, const RKIntegrator::RKData::Data &s, RKIntegrator::RKData::Data &d) const