MFEM  v4.5.2
Finite element discretization library
Public Member Functions | Protected Attributes | List of all members
mfem::SecondOrderODESolver Class Referenceabstract

Abstract class for solving systems of ODEs: d2x/dt2 = f(x,dx/dt,t) More...

#include <ode.hpp>

Inheritance diagram for mfem::SecondOrderODESolver:
[legend]
Collaboration diagram for mfem::SecondOrderODESolver:
[legend]

Public Member Functions

 SecondOrderODESolver ()
 
virtual void Init (SecondOrderTimeDependentOperator &f)
 Associate a TimeDependentOperator with the ODE solver. More...
 
virtual void Step (Vector &x, Vector &dxdt, double &t, double &dt)=0
 Perform a time step from time t [in] to time t [out] based on the requested step size dt [in]. More...
 
virtual void Run (Vector &x, Vector &dxdt, 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 VectorGetStateVector (int i)
 
virtual void GetStateVector (int i, Vector &state)
 
virtual void SetStateVector (int i, Vector &state)
 
virtual ~SecondOrderODESolver ()
 

Protected Attributes

SecondOrderTimeDependentOperatorf
 Pointer to the associated TimeDependentOperator. More...
 
MemoryType mem_type
 

Detailed Description

Abstract class for solving systems of ODEs: d2x/dt2 = f(x,dx/dt,t)

Definition at line 607 of file ode.hpp.

Constructor & Destructor Documentation

◆ SecondOrderODESolver()

mfem::SecondOrderODESolver::SecondOrderODESolver ( )
inline

Definition at line 615 of file ode.hpp.

◆ ~SecondOrderODESolver()

virtual mfem::SecondOrderODESolver::~SecondOrderODESolver ( )
inlinevirtual

Definition at line 705 of file ode.hpp.

Member Function Documentation

◆ GetMaxStateSize()

virtual int mfem::SecondOrderODESolver::GetMaxStateSize ( )
inlinevirtual

Function for getting and setting the state vectors.

Reimplemented in mfem::GeneralizedAlpha2Solver.

Definition at line 689 of file ode.hpp.

◆ GetStateSize()

virtual int mfem::SecondOrderODESolver::GetStateSize ( )
inlinevirtual

Reimplemented in mfem::GeneralizedAlpha2Solver.

Definition at line 690 of file ode.hpp.

◆ GetStateVector() [1/2]

virtual const Vector& mfem::SecondOrderODESolver::GetStateVector ( int  i)
inlinevirtual

Reimplemented in mfem::GeneralizedAlpha2Solver.

Definition at line 691 of file ode.hpp.

◆ GetStateVector() [2/2]

virtual void mfem::SecondOrderODESolver::GetStateVector ( int  i,
Vector state 
)
inlinevirtual

Reimplemented in mfem::GeneralizedAlpha2Solver.

Definition at line 696 of file ode.hpp.

◆ Init()

void mfem::SecondOrderODESolver::Init ( SecondOrderTimeDependentOperator f)
virtual

Associate a TimeDependentOperator with the ODE solver.

This method has to be called:

Reimplemented in mfem::GeneralizedAlpha2Solver, and mfem::NewmarkSolver.

Definition at line 989 of file ode.cpp.

◆ Run()

virtual void mfem::SecondOrderODESolver::Run ( Vector x,
Vector dxdt,
double &  t,
double &  dt,
double  tf 
)
inlinevirtual

Perform time integration from time t [in] to time tf [in].

Parameters
[in,out]xApproximate solution.
[in,out]dxdtApproximate rate.
[in,out]tTime associated with the approximate solution x.
[in,out]dtTime step size.
[in]tfRequested final time.

The default implementation makes consecutive calls to Step() until reaching tf. 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 dxdt [in] is the approximate rate for the input time t [in].
  • The input dt [in] is the initial time step size.
  • 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 output value of t [out] is not smaller than tf [in].

Definition at line 683 of file ode.hpp.

◆ SetStateVector()

virtual void mfem::SecondOrderODESolver::SetStateVector ( int  i,
Vector state 
)
inlinevirtual

Reimplemented in mfem::GeneralizedAlpha2Solver.

Definition at line 700 of file ode.hpp.

◆ Step()

virtual void mfem::SecondOrderODESolver::Step ( Vector x,
Vector dxdt,
double &  t,
double &  dt 
)
pure virtual

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

Parameters
[in,out]xApproximate solution.
[in,out]dxdtApproximate rate.
[in,out]tTime associated with the approximate solution x and rate @ dxdt
[in,out]dtTime 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 dxdt [in] is the approximate rate 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 dxdt [out] is the approximate rate 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.

Implemented in mfem::GeneralizedAlpha2Solver, and mfem::NewmarkSolver.

Member Data Documentation

◆ f

SecondOrderTimeDependentOperator* mfem::SecondOrderODESolver::f
protected

Pointer to the associated TimeDependentOperator.

Definition at line 611 of file ode.hpp.

◆ mem_type

MemoryType mfem::SecondOrderODESolver::mem_type
protected

Definition at line 612 of file ode.hpp.


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