AMRTGA< T > Class Template Reference

#include <AMRTGA.H>

Detailed Description

template<class T>
class AMRTGA< T >

Template implementation of the TGA algorithm to solve the heat equation.


Public Member Functions
TGAHelmOp< T > *	newOp (const ProblemDomain &a_indexSpace, const AMRLevelOpFactory< T > &a_opFact)
	~AMRTGA ()
	AMRTGA (const RefCountedPtr< AMRMultiGrid< T > > &a_solver, const AMRLevelOpFactory< T > &a_factory, const ProblemDomain &a_level0Domain, const Vector< int > &a_refRat, int a_numLevels=-1, int a_verbosity=0)
void	oneStep (Vector< T * > &a_phiNew, Vector< T * > &a_phiOld, Vector< T * > &a_source, const Real &a_dt, int a_lbase, int a_lmax, Real a_timeOld=0)
void	resetAlphaAndBeta (const Real &a_alpha, const Real &a_beta)
void	setTime (Real time)
Protected Member Functions
void	solveHelm (Vector< T * > &a_ans, Vector< T * > &a_rhs, int a_lbase, int a_lmax, Real a_mu, Real a_dt)
void	applyHelm (Vector< T * > &a_ans, Vector< T * > &a_source, int a_lbase, int a_lmax, Real a_mu, Real a_dt, bool a_homogeneousBC)
void	setMu ()
void	createData (Vector< T * > &a_source, int a_lbase, int a_lmax)
Private Member Functions
	AMRTGA (const AMRTGA< T > &a_opin)
void	operator= (const AMRTGA< T > &a_opin)
	AMRTGA ()
	weak construction is bad. Ref Counted pointers are your friends.
Private Attributes
Vector< TGAHelmOp< T > * >	m_ops
ProblemDomain	m_level0Domain
Vector< int >	m_refRat
RefCountedPtr< AMRMultiGrid< T > >	m_solver
Real	m_mu1
Real	m_mu2
Real	m_mu3
Real	m_mu4
int	m_verbosity
int	m_numLevels
bool	m_dataCreated
Vector< T * >	m_rhst
Vector< T * >	m_srct

Constructor & Destructor Documentation

template<class T>

AMRTGA< T >::~AMRTGA ( ) [inline]

References AMRTGA< T >::m_rhst, AMRTGA< T >::m_srct, and Vector< T >::size().

template<class T>

AMRTGA< T >::AMRTGA	(	const RefCountedPtr< AMRMultiGrid< T > > &	a_solver,
		const AMRLevelOpFactory< T > &	a_factory,
		const ProblemDomain &	a_level0Domain,
		const Vector< int > &	a_refRat,
		int	a_numLevels = `-1`,
		int	a_verbosity = `0`
	)			`[inline]`

References MayDay::Error(), AMRTGA< T >::m_dataCreated, AMRTGA< T >::m_level0Domain, AMRTGA< T >::m_numLevels, AMRTGA< T >::m_ops, AMRTGA< T >::m_refRat, AMRTGA< T >::m_solver, AMRTGA< T >::m_verbosity, Vector< T >::resize(), AMRTGA< T >::setMu(), and Vector< T >::size().

template<class T>

AMRTGA< T >::AMRTGA ( const AMRTGA< T > & a_opin ) [inline, private]

References MayDay::Error().

template<class T>

AMRTGA< T >::AMRTGA ( ) [inline, private]

weak construction is bad. Ref Counted pointers are your friends.

References MayDay::Error().

Member Function Documentation

template<class T>

TGAHelmOp<T>* AMRTGA< T >::newOp	(	const ProblemDomain &	a_indexSpace,
		const AMRLevelOpFactory< T > &	a_opFact
	)			`[inline]`

References AMRLevelOpFactory< T >::AMRnewOp().

template<class T>

void AMRTGA< T >::oneStep	(	Vector< T * > &	a_phiNew,
		Vector< T * > &	a_phiOld,
		Vector< T * > &	a_source,
		const Real &	a_dt,
		int	a_lbase,
		int	a_lmax,
		Real	a_timeOld = `0`
	)			`[inline]`

This advances a parabolic pde from a_phiOld to a_phiNew using TGA on a non-moving domain with source term a_source If your operator is time-dependent, be sure to send the time.

References AMRTGA< T >::applyHelm(), AMRTGA< T >::createData(), AMRTGA< T >::m_dataCreated, AMRTGA< T >::m_mu1, AMRTGA< T >::m_mu2, AMRTGA< T >::m_mu3, AMRTGA< T >::m_mu4, AMRTGA< T >::m_ops, AMRTGA< T >::m_rhst, AMRTGA< T >::m_srct, AMRTGA< T >::m_verbosity, pout(), AMRTGA< T >::setTime(), and AMRTGA< T >::solveHelm().

template<class T>

void AMRTGA< T >::resetAlphaAndBeta	(	const Real &	a_alpha,
		const Real &	a_beta
	)			`[inline]`

References AMRTGA< T >::m_solver, TGAHelmOp< T >::setAlphaAndBeta(), and Vector< T >::size().

Referenced by AMRTGA< T >::applyHelm(), and AMRTGA< T >::solveHelm().

template<class T>

void AMRTGA< T >::setTime ( Real time ) [inline]

References AMRTGA< T >::m_solver, TGAHelmOp< T >::setTime(), and Vector< T >::size().

Referenced by AMRTGA< T >::oneStep().

template<class T>

void AMRTGA< T >::solveHelm	(	Vector< T * > &	a_ans,
		Vector< T * > &	a_rhs,
		int	a_lbase,
		int	a_lmax,
		Real	a_mu,
		Real	a_dt
	)			`[inline, protected]`

References AMRTGA< T >::m_solver, AMRTGA< T >::m_verbosity, pout(), and AMRTGA< T >::resetAlphaAndBeta().

Referenced by AMRTGA< T >::oneStep().

template<class T>

void AMRTGA< T >::applyHelm	(	Vector< T * > &	a_ans,
		Vector< T * > &	a_source,
		int	a_lbase,
		int	a_lmax,
		Real	a_mu,
		Real	a_dt,
		bool	a_homogeneousBC
	)			`[inline, protected]`