Public Member Functions | Private Attributes | Static Private Attributes
rodas34 Class Reference

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

Inheritance diagram for rodas34:
Inheritance graph
[legend]
Collaboration diagram for rodas34:
Collaboration graph
[legend]

Public Member Functions

 TypeName ("rodas34")
 Runtime type information. More...
 
 rodas34 (const ODESystem &ode, const dictionary &dict)
 Construct from ODE. More...
 
scalar solve (const scalar x0, const scalarField &y0, const scalarField &dydx0, const scalar dx, scalarField &y) const
 Solve a single step dx and return the error. More...
 
void solve (scalar &x, scalarField &y, scalar &dxTry) const
 Solve the ODE system and the update the state. More...
 
virtual void solve (scalar &x, scalarField &y, scalar &dxTry) const
 Inherit solve from ODESolver. More...
 
virtual void solve (scalar &x, scalarField &y, stepState &step) const
 Inherit solve from ODESolver. More...
 
virtual void solve (const scalar xStart, const scalar xEnd, scalarField &y, scalar &dxEst) const
 Inherit solve from ODESolver. 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, stepState &step) 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...
 
- Public Member Functions inherited from adaptiveSolver
 adaptiveSolver (const ODESystem &ode, const dictionary &dict)
 Construct from ODESystem. More...
 
virtual ~adaptiveSolver ()
 Destructor. More...
 
void solve (const ODESystem &ode, scalar &x, scalarField &y, scalar &dxTry) const
 Solve the ODE system and the update the state. More...
 

Private Attributes

scalarField k1_
 
scalarField k2_
 
scalarField k3_
 
scalarField k4_
 
scalarField k5_
 
scalarField dy_
 
scalarField err_
 
scalarField dydx_
 
scalarField dfdx_
 
scalarSquareMatrix dfdy_
 
scalarSquareMatrix a_
 
labelList pivotIndices_
 

Static Private Attributes

static const scalar c2 = 0.386
 
static const scalar c3 = 0.21
 
static const scalar c4 = 0.63
 
static const scalar d1 = 0.25
 
static const scalar d2 = -0.1043
 
static const scalar d3 = 0.1035
 
static const scalar d4 = -0.3620000000000023e-01
 
static const scalar a21 = 0.1544e1
 
static const scalar a31 = 0.9466785280815826
 
static const scalar a32 = 0.2557011698983284
 
static const scalar a41 = 0.3314825187068521e1
 
static const scalar a42 = 0.2896124015972201e1
 
static const scalar a43 = 0.9986419139977817
 
static const scalar a51 = 0.1221224509226641e1
 
static const scalar a52 = 0.6019134481288629e1
 
static const scalar a53 = 0.1253708332932087e2
 
static const scalar a54 = -0.6878860361058950
 
static const scalar c21 = -0.56688e1
 
static const scalar c31 = -0.2430093356833875e1
 
static const scalar c32 = -0.2063599157091915
 
static const scalar c41 = -0.1073529058151375
 
static const scalar c42 = -0.9594562251023355e1
 
static const scalar c43 = -0.2047028614809616e2
 
static const scalar c51 = 0.7496443313967647e1
 
static const scalar c52 = -0.1024680431464352e2
 
static const scalar c53 = -0.3399990352819905e2
 
static const scalar c54 = 0.1170890893206160e2
 
static const scalar c61 = 0.8083246795921522e1
 
static const scalar c62 = -0.7981132988064893e1
 
static const scalar c63 = -0.3152159432874371e2
 
static const scalar c64 = 0.1631930543123136e2
 
static const scalar c65 = -0.6058818238834054e1
 
static const scalar gamma = 0.25
 

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

Detailed Description

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

References: 

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

Definition at line 59 of file rodas34.H.

Constructor & Destructor Documentation

◆ rodas34()

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

Construct from ODE.

Definition at line 75 of file rodas34.C.

Member Function Documentation

◆ TypeName()

TypeName ( "rodas34"  )

Runtime type information.

◆ solve() [1/5]

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

Solve a single step dx and return the error.

Implements adaptiveSolver.

Definition at line 97 of file rodas34.C.

References Foam::constant::physicoChemical::c2, forAll, 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

Solve the ODE system and the update the state.

Reimplemented from ODESolver.

Definition at line 217 of file rodas34.C.

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

Here is the call graph for this function:

◆ solve() [3/5]

virtual void solve
inline

Inherit solve from ODESolver.

Definition at line 178 of file ODESolver.H.

◆ solve() [4/5]

void solve

Inherit solve from ODESolver.

Definition at line 86 of file ODESolver.C.

◆ solve() [5/5]

void solve

Inherit solve from ODESolver.

Definition at line 99 of file ODESolver.C.

Field Documentation

◆ k1_

scalarField k1_
mutableprivate

Definition at line 66 of file rodas34.H.

◆ k2_

scalarField k2_
mutableprivate

Definition at line 67 of file rodas34.H.

◆ k3_

scalarField k3_
mutableprivate

Definition at line 68 of file rodas34.H.

◆ k4_

scalarField k4_
mutableprivate

Definition at line 69 of file rodas34.H.

◆ k5_

scalarField k5_
mutableprivate

Definition at line 70 of file rodas34.H.

◆ dy_

scalarField dy_
mutableprivate

Definition at line 71 of file rodas34.H.

◆ err_

scalarField err_
mutableprivate

Definition at line 72 of file rodas34.H.

◆ dydx_

scalarField dydx_
mutableprivate

Definition at line 73 of file rodas34.H.

◆ dfdx_

scalarField dfdx_
mutableprivate

Definition at line 74 of file rodas34.H.

◆ dfdy_

scalarSquareMatrix dfdy_
mutableprivate

Definition at line 75 of file rodas34.H.

◆ a_

scalarSquareMatrix a_
mutableprivate

Definition at line 76 of file rodas34.H.

◆ pivotIndices_

labelList pivotIndices_
mutableprivate

Definition at line 77 of file rodas34.H.

◆ c2

const scalar c2 = 0.386
staticprivate

Definition at line 80 of file rodas34.H.

◆ c3

const scalar c3 = 0.21
staticprivate

Definition at line 80 of file rodas34.H.

◆ c4

const scalar c4 = 0.63
staticprivate

Definition at line 80 of file rodas34.H.

◆ d1

const scalar d1 = 0.25
staticprivate

Definition at line 81 of file rodas34.H.

◆ d2

const scalar d2 = -0.1043
staticprivate

Definition at line 81 of file rodas34.H.

◆ d3

const scalar d3 = 0.1035
staticprivate

Definition at line 81 of file rodas34.H.

◆ d4

const scalar d4 = -0.3620000000000023e-01
staticprivate

Definition at line 81 of file rodas34.H.

◆ a21

const scalar a21 = 0.1544e1
staticprivate

Definition at line 82 of file rodas34.H.

◆ a31

const scalar a31 = 0.9466785280815826
staticprivate

Definition at line 82 of file rodas34.H.

◆ a32

const scalar a32 = 0.2557011698983284
staticprivate

Definition at line 82 of file rodas34.H.

◆ a41

const scalar a41 = 0.3314825187068521e1
staticprivate

Definition at line 83 of file rodas34.H.

◆ a42

const scalar a42 = 0.2896124015972201e1
staticprivate

Definition at line 83 of file rodas34.H.

◆ a43

const scalar a43 = 0.9986419139977817
staticprivate

Definition at line 83 of file rodas34.H.

◆ a51

const scalar a51 = 0.1221224509226641e1
staticprivate

Definition at line 84 of file rodas34.H.

◆ a52

const scalar a52 = 0.6019134481288629e1
staticprivate

Definition at line 84 of file rodas34.H.

◆ a53

const scalar a53 = 0.1253708332932087e2
staticprivate

Definition at line 84 of file rodas34.H.

◆ a54

const scalar a54 = -0.6878860361058950
staticprivate

Definition at line 84 of file rodas34.H.

◆ c21

const scalar c21 = -0.56688e1
staticprivate

Definition at line 85 of file rodas34.H.

◆ c31

const scalar c31 = -0.2430093356833875e1
staticprivate

Definition at line 85 of file rodas34.H.

◆ c32

const scalar c32 = -0.2063599157091915
staticprivate

Definition at line 85 of file rodas34.H.

◆ c41

const scalar c41 = -0.1073529058151375
staticprivate

Definition at line 86 of file rodas34.H.

◆ c42

const scalar c42 = -0.9594562251023355e1
staticprivate

Definition at line 86 of file rodas34.H.

◆ c43

const scalar c43 = -0.2047028614809616e2
staticprivate

Definition at line 86 of file rodas34.H.

◆ c51

const scalar c51 = 0.7496443313967647e1
staticprivate

Definition at line 87 of file rodas34.H.

◆ c52

const scalar c52 = -0.1024680431464352e2
staticprivate

Definition at line 87 of file rodas34.H.

◆ c53

const scalar c53 = -0.3399990352819905e2
staticprivate

Definition at line 87 of file rodas34.H.

◆ c54

const scalar c54 = 0.1170890893206160e2
staticprivate

Definition at line 87 of file rodas34.H.

◆ c61

const scalar c61 = 0.8083246795921522e1
staticprivate

Definition at line 88 of file rodas34.H.

◆ c62

const scalar c62 = -0.7981132988064893e1
staticprivate

Definition at line 88 of file rodas34.H.

◆ c63

const scalar c63 = -0.3152159432874371e2
staticprivate

Definition at line 88 of file rodas34.H.

◆ c64

const scalar c64 = 0.1631930543123136e2
staticprivate

Definition at line 88 of file rodas34.H.

◆ c65

const scalar c65 = -0.6058818238834054e1
staticprivate

Definition at line 88 of file rodas34.H.

◆ gamma

const scalar gamma = 0.25
staticprivate

Definition at line 89 of file rodas34.H.

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

Copyright © 2011-2018  OpenFOAM | OPENFOAM® is a registered trademark of OpenCFD Ltd.

Generated by   doxygen 1.8.17