An extrapolation-algorithm, based on the linearly implicit Euler method with step size control and order selection. More...
Public Member Functions | |
| TypeName ("seulex") | |
| Runtime type information. More... | |
| seulex (const ODESystem &ode, const dictionary &dict) | |
| Construct from ODE. More... | |
| void | solve (scalar &x, scalarField &y, stepState &step) const |
| Solve the ODE system and the update the state. More... | |
 Public Member Functions inherited from ODESolver | |
| TypeName ("ODESolver") | |
| Runtime type information. More... | |
| declareRunTimeSelectionTable (autoPtr, ODESolver, dictionary,(const ODESystem &ode, const dictionary &dict),(ode, dict)) | |
| ODESolver (const ODESystem &ode, const dictionary &dict) | |
| Construct for given ODESystem. More... | |
| ODESolver (const ODESystem &ode, const scalarField &absTol, const scalarField &relTol) | |
| Construct for given ODESystem specifying tolerances. More... | |
| virtual | ~ODESolver () |
| Destructor. More... | |
| scalarField & | absTol () |
| scalarField & | relTol () |
| virtual void | solve (scalar &x, scalarField &y, scalar &dxTry) const |
| Solve the ODE system as far as possible upto dxTry. More... | |
| virtual void | solve (const scalar xStart, const scalar xEnd, scalarField &y, scalar &dxEst) const |
| Solve the ODE system from xStart to xEnd, update the state. More... | |
Additional Inherited Members | |
| Static Public Member Functions inherited from ODESolver | |
| static autoPtr< ODESolver > | New (const ODESystem &ode, const dictionary &dict) |
| Select null constructed. More... | |
| Protected Member Functions inherited from ODESolver | |
| scalar | normalizeError (const scalarField &y0, const scalarField &y, const scalarField &err) const |
| Return the nomalized scalar error. More... | |
| ODESolver (const ODESolver &) | |
| Disallow default bitwise copy construct. More... | |
| void | operator= (const ODESolver &) |
| Disallow default bitwise assignment. More... | |
| Protected Attributes inherited from ODESolver | |
| const ODESystem & | odes_ |
| Reference to ODESystem. More... | |
| label | n_ |
| Size of the ODESystem. More... | |
| scalarField | absTol_ |
| Absolute convergence tolerance per step. More... | |
| scalarField | relTol_ |
| Relative convergence tolerance per step. More... | |
| label | maxSteps_ |
| The maximum number of sub-steps allowed for the integration step. More... | |
An extrapolation-algorithm, based on the linearly implicit Euler method with step size control and order selection.
The implementation is based on the SEULEX algorithm in
"Solving Ordinary Differential Equations II: Stiff
and Differential-Algebraic Problems, second edition",
Hairer, E.,
Nørsett, S.,
Wanner, G.,
Springer-Verlag, Berlin. 1996.
| seulex | ( | const ODESystem & | ode, |
| const dictionary & | dict | ||
| ) |
Construct from ODE.
Definition at line 50 of file seulex.C.
References ODESystem::derivatives(), k, Foam::LUBacksubstitute(), Foam::LUDecompose(), Foam::mag(), Foam::min(), ODESolver::n_, ODESolver::odes_, seulex::solve(), Foam::sqr(), Foam::sqrt(), y, and Foam::y0().
| TypeName | ( | "seulex" | ) |
Runtime type information.
|
virtual |
Solve the ODE system and the update the state.
Reimplemented from ODESolver.
Definition at line 215 of file seulex.C.
References ODESolver::absTol_, ODESolver::stepState::dxDid, ODESolver::stepState::dxTry, Foam::endl(), ODESolver::stepState::first, forAll, ODESolver::stepState::forward, ODESystem::jacobian(), k, ODESolver::stepState::last, Foam::log10(), Foam::mag(), Foam::max(), Foam::min(), ODESolver::n_, ODESolver::odes_, Foam::pow(), ODESolver::stepState::prevReject, ODESolver::stepState::reject, ODESolver::relTol_, Foam::sqr(), Foam::sqrt(), success, WarningInFunction, and y.
Referenced by seulex::seulex().
1.8.11