80 y[i] =
y0[i] + dx*dydx0[i];
84 odes_.derivatives(x0 + dx,
y, err_);
90 y[i] =
y0[i] + 0.5*dx*(dydx0[i] + err_[i]);
91 err_[i] = 0.5*dx*(err_[i] - dydx0[i]);
94 return normalizeError(
y0,
y, err_);
virtual scalar solve(const scalar x0, const scalarField &y0, const scalarField &dydx0, const scalar dx, scalarField &y) const =0
Solve a single step dx and return the error.
Abstract base class for the systems of ordinary differential equations.
virtual bool resize()=0
Resize the ODE solver.
virtual scalar solve(const scalar x0, const scalarField &y0, const scalarField &dydx0, const scalar dx, scalarField &y) const
Solve a single step dx and return the error.
A list of keyword definitions, which are a keyword followed by a number of values (eg...
dimensionedScalar y0(const dimensionedScalar &ds)
An ODE solver for chemistry.
Macros for easy insertion into run-time selection tables.
#define forAll(list, i)
Loop across all elements in list.
bool resize(const label n)
Resize the ODE solver.
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
defineTypeNameAndDebug(combustionModel, 0)
Trapezoidal ODE solver of order (1)2.
Abstract base-class for ODE system solvers.
static void resizeField(UList< Type > &f, const label n)
Trapezoid(const ODESystem &ode, const dictionary &dict)
Construct from ODESystem.
label n_
Size of the ODESystem (adjustable)
virtual bool resize()
Resize the ODE solver.
addToRunTimeSelectionTable(functionObject, pointHistory, dictionary)