Foam::RK::RK(ODE& ode) : ODESolver(ode), yTemp_(ode.nEqns()), ak2_(ode.nEqns()), ak3_(ode.nEqns()), ak4_(ode.nEqns()), ak5_(ode.nEqns()), ak6_(ode.nEqns()), yErr_(ode.nEqns()), yTemp2_(ode.nEqns()) {}
Foam::ODESolver::ODESolver(const ODE& ode) : n_(ode.nEqns()), yScale_(n_), dydx_(n_) {}
Foam::ODESolver::ODESolver(ODE& ode) : ode_(ode), yScale_(ode.nEqns()), dydx_(ode.nEqns()) {}