121 forward(dx > 0 ? true : false),
175 inline label
nEqns()
const;
184 virtual bool resize() = 0;
label nEqns() const
Return the number of equations to solve.
virtual void solve(scalar &x, scalarField &y, scalar &dxTry) const
Solve the ODE system as far as possible up to dxTry.
scalarField & relTol()
Return access to the relative tolerance field.
Abstract base class for the systems of ordinary differential equations.
virtual bool resize()=0
Resize the ODE solver.
void operator=(const ODESolver &)=delete
No copy assignment.
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.
stepState(const scalar dx)
void resizeMatrix(scalarSquareMatrix &m) const
const ODESystem & odes_
Reference to ODESystem.
TypeName("ODESolver")
Runtime type information.
label maxSteps_
The maximum number of sub-steps allowed for the integration step.
const label maxN_
Maximum size of the ODESystem.
declareRunTimeSelectionTable(autoPtr, ODESolver, dictionary,(const ODESystem &ode, const dictionary &dict),(ode, dict))
virtual ~ODESolver()=default
Destructor.
scalarField relTol_
Relative convergence tolerance per step.
scalar normalizeError(const scalarField &y0, const scalarField &y, const scalarField &err) const
Return the nomalized scalar error.
Abstract base-class for ODE system solvers.
scalarField absTol_
Absolute convergence tolerance per step.
static autoPtr< ODESolver > New(const ODESystem &ode, const dictionary &dict)
Select null constructed.
scalarField & absTol()
Return access to the absolute tolerance field.
static void resizeField(UList< Type > &f, const label n)
Pointer management similar to std::unique_ptr, with some additional methods and type checking...
ODESolver(const ODESolver &)=delete
No copy construct.
label n_
Size of the ODESystem (adjustable)