#include <ode.hpp>

Inheritance diagram for mfem::SIASolver:

Collaboration diagram for mfem::SIASolver:

Public Member Functions
	SIASolver ()

virtual void	Init (Operator &P, TimeDependentOperator &F)

virtual void	Step (Vector &q, Vector &p, double &t, double &dt)=0

virtual void	Run (Vector &q, Vector &p, double &t, double &dt, double tf)

virtual	~SIASolver ()

Protected Attributes
TimeDependentOperator *	F_

Operator *	P_

Vector	dp_

Vector	dq_

Detailed Description

The SIASolver class is based on the Symplectic Integration Algorithm described in "A Symplectic Integration Algorithm for Separable Hamiltonian Functions" by J. Candy and W. Rozmus, Journal of Computational Physics, Vol. 92, pages 230-256 (1991). The Symplectic Integration Algorithm (SIA) is designed for systems of first order ODEs derived from a Hamiltonian. H(q,p,t) = T(p) + V(q,t) Which leads to the equations: dq/dt = dT/dp dp/dt = -dV/dq In the integrator the operators P and F are defined to be: P = dT/dp F = -dV/dq

Definition at line 570 of file ode.hpp.

Constructor & Destructor Documentation

◆ SIASolver()

mfem::SIASolver::SIASolver ( )

inline

Definition at line 573 of file ode.hpp.

◆ ~SIASolver()

virtual mfem::SIASolver::~SIASolver ( )

inlinevirtual

Definition at line 584 of file ode.hpp.

Member Function Documentation

◆ Init()

void mfem::SIASolver::Init	(	Operator &	P,
		TimeDependentOperator &	F
	)

virtual

Definition at line 916 of file ode.cpp.

◆ Run()

virtual void mfem::SIASolver::Run	(	Vector &	q,
		Vector &	p,
		double &	t,
		double &	dt,
		double	tf
	)

inlinevirtual

Definition at line 579 of file ode.hpp.

◆ Step()

virtual void mfem::SIASolver::Step	(	Vector &	q,
		Vector &	p,
		double &	t,
		double &	dt
	)

pure virtual

Implemented in mfem::SIAVSolver, mfem::SIA2Solver, and mfem::SIA1Solver.

Member Data Documentation

◆ dp_

Vector mfem::SIASolver::dp_

mutableprotected

Definition at line 590 of file ode.hpp.

◆ dq_

Vector mfem::SIASolver::dq_

mutableprotected

Definition at line 591 of file ode.hpp.

◆ F_

TimeDependentOperator* mfem::SIASolver::F_

protected

Definition at line 587 of file ode.hpp.

◆ P_

Operator* mfem::SIASolver::P_