Rosenbrock34 Class Reference

L-stable embedded Rosenbrock ODE solver of order (3)4. More...

Inheritance diagram for Rosenbrock34:

Collaboration diagram for Rosenbrock34:

Public Member Functions
	TypeName ("Rosenbrock34")
	Rosenbrock34 (const ODESystem &ode, const dictionary &dict)
virtual	~Rosenbrock34 ()=default
virtual bool	resize ()
virtual scalar	solve (const scalar x0, const scalarField &y0, const scalarField &dydx0, const scalar dx, scalarField &y) const
virtual void	solve (scalar &x, scalarField &y, scalar &dxTry) const
virtual void	solve (scalar &x, scalarField &y, scalar &dxTry) const
virtual void	solve (scalar &x, scalarField &y, stepState &step) const
virtual void	solve (const scalar xStart, const scalar xEnd, scalarField &y, scalar &dxEst) const
Public Member Functions inherited from ODESolver
	TypeName ("ODESolver")
	declareRunTimeSelectionTable (autoPtr, ODESolver, dictionary,(const ODESystem &ode, const dictionary &dict),(ode, dict))
	ODESolver (const ODESystem &ode, const dictionary &dict)
	ODESolver (const ODESystem &ode, const scalarField &absTol, const scalarField &relTol)
virtual	~ODESolver ()=default
label	nEqns () const
scalarField &	absTol ()
scalarField &	relTol ()
template<class Type >
void	resizeField (UList< Type > &f) const
void	resizeMatrix (scalarSquareMatrix &m) const
virtual void	solve (scalar &x, scalarField &y, stepState &step) const
virtual void	solve (const scalar xStart, const scalar xEnd, scalarField &y, scalar &dxEst) const
Public Member Functions inherited from adaptiveSolver
	adaptiveSolver (const ODESystem &ode, const dictionary &dict)
virtual	~adaptiveSolver ()=default
bool	resize (const label n)
void	solve (const ODESystem &ode, scalar &x, scalarField &y, scalar &dxTry) const

Additional Inherited Members
Static Public Member Functions inherited from ODESolver
static autoPtr< ODESolver >	New (const ODESystem &ode, const dictionary &dict)
template<class Type >
static void	resizeField (UList< Type > &f, const label n)
Protected Member Functions inherited from ODESolver
scalar	normalizeError (const scalarField &y0, const scalarField &y, const scalarField &err) const
	ODESolver (const ODESolver &)=delete
void	operator= (const ODESolver &)=delete
Protected Attributes inherited from ODESolver
const ODESystem &	odes_
const label	maxN_
label	n_
scalarField	absTol_
scalarField	relTol_
label	maxSteps_

Detailed Description

L-stable embedded Rosenbrock ODE solver of order (3)4.

    Hairer, E., Nørsett, S. P., & Wanner, G. (1996).
    Solving Ordinary Differential Equations II:
    Stiff and Differential-Algebraic Problems, second edition",
    Springer-Verlag, Berlin.

The default constants are from:

    Shampine, L. F. (1982).
    Implementation of Rosenbrock Methods.
    ACM Transactions on Mathematical Software, vol. 8, pp. 93–113.

with which the scheme is more accurate than with the L-Stable coefficients for small step-size but less stable for large step-size.

The L-Stable scheme constants are provided commented-out in Rosenbrock34.C

Source files

• Rosenbrock34.H
• Rosenbrock34.C

Definition at line 69 of file Rosenbrock34.H.

Constructor & Destructor Documentation

Rosenbrock34()

Rosenbrock34	(	const ODESystem &	ode,
		const dictionary &	dict
	)

Definition at line 104 of file Rosenbrock34.C.

~Rosenbrock34()

virtual ~Rosenbrock34 ( )

virtualdefault

Member Function Documentation

TypeName()

TypeName

(

"Rosenbrock34"

)

resize()

bool resize ( )

virtual

Implements ODESolver.

Definition at line 123 of file Rosenbrock34.C.

References adaptiveSolver::resize(), and ODESolver::resize().

Here is the call graph for this function:

solve() [1/5]

Foam::scalar solve ( const scalar  x0, const scalarField &  y0, const scalarField &  dydx0, const scalar  dx, scalarField &  y  ) const

virtual

Implements adaptiveSolver.

Definition at line 148 of file Rosenbrock34.C.

References Foam::constant::physicoChemical::c2, forAll, gamma, Foam::LUBacksubstitute(), Foam::LUDecompose(), y, and Foam::y0().

Here is the call graph for this function:

solve() [2/5]

void solve ( scalar &  x, scalarField &  y, scalar &  dxTry  ) const

virtual

Reimplemented from ODESolver.

Definition at line 229 of file Rosenbrock34.C.

References adaptiveSolver::solve(), x, and y.

Here is the call graph for this function:

solve() [3/5]

void solve

Definition at line 110 of file ODESolver.C.

solve() [4/5]

void solve

Definition at line 123 of file ODESolver.C.

solve() [5/5]

void solve

Definition at line 136 of file ODESolver.C.

---

The documentation for this class was generated from the following files:

• src/ODE/ODESolvers/Rosenbrock34/Rosenbrock34.H
• src/ODE/ODESolvers/Rosenbrock34/Rosenbrock34.C