A semi-implicit mid-point solver for stiff systems of ordinary differential equations. More...

Inheritance diagram for SIBS:

Collaboration diagram for SIBS:

Public Member Functions
	TypeName ("SIBS")
	Runtime type information. More...

	SIBS (const ODESystem &ode, const dictionary &dict)
	Construct from ODE system. More...

virtual	~SIBS ()
	Destructor. More...

virtual bool	resize ()
	Resize the ODE solver. More...

virtual void	solve (scalar &x, scalarField &y, scalar &dxTry) const
	Solve the ODE system as far as possible up to dxTry. 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...

label	nEqns () const
	Return the number of equations to solve. More...

scalarField &	absTol ()
	Return access to the absolute tolerance field. More...

scalarField &	relTol ()
	Return access to the relative tolerance field. More...

template<class Type >
void	resizeField (UList< Type > &f) const

void	resizeMatrix (scalarSquareMatrix &m) const

virtual void	solve (scalar &x, scalarField &y, stepState &step) const
	Solve the ODE system as far as possible up to 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...

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
	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...

const label	maxN_
	Maximum size of the ODESystem. More...

label	n_
	Size of the ODESystem (adjustable) 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...

Detailed Description

A semi-implicit mid-point solver for stiff systems of ordinary differential equations.

Reference:

    Bader, G., & Deuflhard, P. (1983).
    A semi-implicit mid-point rule for stiff systems
    of ordinary differential equations.
    Numerische Mathematik, 41(3), 373-398.

Source files

Definition at line 58 of file SIBS.H.

Constructor & Destructor Documentation

◆ SIBS()

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

Construct from ODE system.

Definition at line 49 of file SIBS.C.

◆ ~SIBS()

virtual ~SIBS ( )

inlinevirtual

Destructor.

Definition at line 125 of file SIBS.H.

References SIBS::resize(), SIBS::solve(), and x.

Here is the call graph for this function:

Member Function Documentation

◆ TypeName()

TypeName ( "SIBS" )

Runtime type information.

◆ resize()

bool resize ( )

virtual

Resize the ODE solver.

Implements ODESolver.

Definition at line 71 of file SIBS.C.

References ODESolver::resize(), ODESolver::resizeField(), ODESolver::resizeMatrix(), and SIBS::solve().

Referenced by SIBS::~SIBS().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ solve()

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

virtual

Solve the ODE system as far as possible up to dxTry.

adjusting the step as necessary to provide a solution within the specified tolerance. Update the state and return an estimate for the next step in dxTry

Reimplemented from ODESolver.

Definition at line 92 of file SIBS.C.

References ODESolver::absTol_, ODESystem::derivatives(), Foam::exit(), Foam::FatalError, FatalErrorInFunction, h, ODESystem::jacobian(), k, Foam::mag(), Foam::max(), Foam::min(), ODESolver::n_, ODESolver::odes_, Foam::pow(), ODESolver::relTol_, Foam::sqr(), and y.