#include <ode.hpp>

Inheritance diagram for mfem::TrapezoidalRuleSolver:

Collaboration diagram for mfem::TrapezoidalRuleSolver:

Public Member Functions
virtual void	Init (TimeDependentOperator &f_)
	Associate a TimeDependentOperator with the ODE solver. More...

virtual void	Step (Vector &x, double &t, double &dt)
	Perform a time step from time t [in] to time t [out] based on the requested step size dt [in]. More...

Public Member Functions inherited from mfem::ODESolver
	ODESolver ()

virtual void	Run (Vector &x, double &t, double &dt, double tf)
	Perform time integration from time t [in] to time tf [in]. More...

virtual int	GetMaxStateSize ()
	Function for getting and setting the state vectors. More...

virtual int	GetStateSize ()

virtual const Vector &	GetStateVector (int i)

virtual void	GetStateVector (int i, Vector &state)

virtual void	SetStateVector (int i, Vector &state)

virtual	~ODESolver ()

Protected Attributes
Vector	k

Vector	y

Protected Attributes inherited from mfem::ODESolver
TimeDependentOperator *	f
	Pointer to the associated TimeDependentOperator. More...

MemoryType	mem_type

Detailed Description

Two stage, explicit singly diagonal implicit Runge-Kutta (ESDIRK) method of order 2. A-stable.

Definition at line 467 of file ode.hpp.

Member Function Documentation

void mfem::TrapezoidalRuleSolver::Init ( TimeDependentOperator & f_ )

virtual

Associate a TimeDependentOperator with the ODE solver.

This method has to be called:

Reimplemented from mfem::ODESolver.

Definition at line 688 of file ode.cpp.

void mfem::TrapezoidalRuleSolver::Step	(	Vector &	x,
		double &	t,
		double &	dt
	)

virtual

Perform a time step from time t [in] to time t [out] based on the requested step size dt [in].

Parameters

[in,out]	x	Approximate solution.
[in,out]	t	Time associated with the approximate solution x.
[in,out]	dt	Time step size.

The following rules describe the common behavior of the method:

The input x [in] is the approximate solution for the input time t [in].
The input dt [in] is the desired time step size, defining the desired target time: t [target] = t [in] + dt [in].
The output x [out] is the approximate solution for the output time t [out].
The output dt [out] is the last time step taken by the method which may be smaller or larger than the input dt [in] value, e.g. because of time step control.
The method may perform more than one time step internally; in this case dt [out] is the last internal time step size.
The output value of t [out] may be smaller or larger than t [target], however, it is not smaller than t [in] + dt [out], if at least one internal time step was performed.
The value x [out] may be obtained by interpolation using internally stored data.
In some cases, the contents of x [in] may not be used, e.g. when x [out] from a previous Step() call was obtained by interpolation.
In consecutive calls to this method, the output t [out] of one Step() call has to be the same as the input t [in] to the next Step() call.
If the previous rule has to be broken, e.g. to restart a time stepping sequence, then the ODE solver must be re-initialized by calling Init() between the two Step() calls.