Adaptive limiter. More...

#include <advection.hh>

Inheritance diagram for FvdgAdaptiveLimiter:

Collaboration diagram for FvdgAdaptiveLimiter:

Public Types
typedef Cmplxtype< Real >::type	c_type
	Real type of data type.
typedef Realtype< Real >::type	r_type
	Real type of data type.
typedef Real	type
	Type of data, e.g. matrix entries.
Public Member Functions
virtual const uint	dimX () const
	Returns the size of the image space of the operator (number of rows of the corresponding matrix)
virtual const uint	dimY () const
	Returns the size of the source space of the operator (number of columns of the corresponding matrix)
	FvdgAdaptiveLimiter (const linDG3D::FvdgSpace &spc, const Real alpha, const Real lambda, concepts::BoundaryConditions *bc)
	Constructor.
void	getIndicator (concepts::Vector< Real > &indicator)
	Gets values of element indicators.
virtual void	operator() (const Function< c_type > &fncY, Function< c_type > &fncX)
	Application operator for complex function `fncY`.
virtual void	operator() (const concepts::Function< Real > &fncY, concepts::Function< Real > &fncX)
	Application operator.
virtual void	operator() (const Function< r_type > &fncY, Function< Real > &fncX)
	Application operator for real function `fncY`.
virtual const concepts::Space < Real > &	spaceX () const
	Returns the image space of the operator.
virtual const concepts::Space < Real > &	spaceY () const
	Returns the source space of the operator.
Protected Member Functions
void	bndryContrib_ (const concepts::Vector< Real > &vecY, concepts::Vector< Real > &vecX)
	Calculates boundary face contributions.
Real	elmVolume_ (const linDG3D::FvdgElement &elm) const
	Calculates volume of an element.
virtual std::ostream &	info (std::ostream &os) const
	Returns information in an output stream.
void	innerContrib_ (const concepts::Vector< Real > &vecY, concepts::Vector< Real > &vecX)
	Calculates inner face contributions.
void	limit_ (concepts::Vector< Real > &vecX, const linDG3D::FvdgElement &elm) const
void	reset_ ()
	Resets element indicator to zero.
Protected Attributes
const Real	alpha_
concepts::BoundaryConditions *	bc_
uint	dimX_
	Dimension of image space and the source space.
uint	dimY_
std::map< const linDG3D::FvdgElement *, Real >	elmIndicator_
Real	lambda_
const linDG3D::FvdgSpace &	spc_

Detailed Description

Adaptive limiter.

Author:: Adrian Burri, 2002

Definition at line 308 of file advection.hh.

Member Typedef Documentation

typedef Cmplxtype<Real >::type concepts::Operator< Real >::c_type [inherited]

Real type of data type.

Reimplemented in concepts::VecOperator< Real >, concepts::DenseMatrix< Real >, concepts::DiagonalMatrix< Real >, concepts::Matrix< Real >, concepts::Permutation< Real >, and concepts::SparseMatrix< Real >.

Definition at line 47 of file compositions.hh.

typedef Realtype<Real >::type concepts::Operator< Real >::r_type [inherited]

Real type of data type.

Reimplemented in concepts::VecOperator< Real >, concepts::DenseMatrix< Real >, concepts::DiagonalMatrix< Real >, concepts::Matrix< Real >, concepts::Permutation< Real >, and concepts::SparseMatrix< Real >.

Definition at line 45 of file compositions.hh.

typedef Real concepts::Operator< Real >::type [inherited]

Type of data, e.g. matrix entries.

Definition at line 43 of file compositions.hh.

Constructor & Destructor Documentation

FvdgAdaptiveLimiter::FvdgAdaptiveLimiter	(	const linDG3D::FvdgSpace &	spc,
		const Real	alpha,
		const Real	lambda,
		concepts::BoundaryConditions *	bc
	)

Constructor.

Parameters:

spc	FvdgSpace
alpha	Weight factor for indicator
lambda	Weight factor for gradient damping
bc	Boundary conditions

Member Function Documentation

void FvdgAdaptiveLimiter::bndryContrib_	(	const concepts::Vector< Real > &	vecY,
		concepts::Vector< Real > &	vecX
	)		`[protected]`

Calculates boundary face contributions.

virtual const uint concepts::Operator< Real >::dimX ( ) const [inline, virtual, inherited]

Returns the size of the image space of the operator (number of rows of the corresponding matrix)

Examples:: hpFEM2d-simple.cc, and hpFEM2d.cc.

Definition at line 87 of file compositions.hh.

virtual const uint concepts::Operator< Real >::dimY ( ) const [inline, virtual, inherited]

Returns the size of the source space of the operator (number of columns of the corresponding matrix)

Definition at line 92 of file compositions.hh.

Real FvdgAdaptiveLimiter::elmVolume_ ( const linDG3D::FvdgElement & elm ) const [protected]

Calculates volume of an element.

Parameters:

elm	Element of which the volume is calculated

void FvdgAdaptiveLimiter::getIndicator ( concepts::Vector< Real > & indicator )

Gets values of element indicators.

Parameters:

inidicator Vector to store element indicator values

virtual std::ostream& concepts::Operator< Real >::info ( std::ostream & os ) const [protected, virtual, inherited]

Returns information in an output stream.

Reimplemented from concepts::OutputOperator.

Reimplemented in EquivalentStress2d, Advection< NumFlux >, hp3D::SpaceTransition, concepts::Constrained, concepts::LapackChol, concepts::NRLSolver, concepts::Pardiso, concepts::PETSc, concepts::PETScMat, concepts::Spooles, concepts::Umfpack, sparseqr::Driver, vectorial::BramblePasciak, vectorial::CGUzawa, concepts::VecOperator< Real >, concepts::DenseMatrix< Real >, concepts::DiagonalMatrix< Real >, concepts::Permutation< Real >, concepts::TrivExtendRestrict< Real >, concepts::SparseMatrix< Real >, and sparseqr::GivensRotations< Real >.

void FvdgAdaptiveLimiter::innerContrib_	(	const concepts::Vector< Real > &	vecY,
		concepts::Vector< Real > &	vecX
	)		`[protected]`

Calculates inner face contributions.

void FvdgAdaptiveLimiter::limit_	(	concepts::Vector< Real > &	vecX,
		const linDG3D::FvdgElement &	elm
	)		const `[protected]`

virtual void concepts::Operator< Real >::operator()	(	const Function< c_type > &	fncY,
		Function< c_type > &	fncX
	)		`[virtual, inherited]`

Application operator for complex function fncY.

Computes fncX = A(fncY) where A is this operator. fncX becomes complex.

In derived classes its enough to implement the operator() for complex Operator's. If a real counterpart is not implemented, the function fncY is splitted into real and imaginary part and the application operator for real functions is called for each. Then the result is combined.

If in a derived class the operator() for complex Operator's is not implemented, a exception is thrown from here.

Reimplemented in concepts::VecOperator< Real >, concepts::DenseMatrix< Real >, concepts::DiagonalMatrix< Real >, concepts::Matrix< Real >, concepts::Permutation< Real >, and concepts::SparseMatrix< Real >.

virtual void FvdgAdaptiveLimiter::operator()	(	const concepts::Function< Real > &	fncY,
		concepts::Function< Real > &	fncX
	)		`[virtual]`

Application operator.

Limits according to an indicator function.

Parameters:

fncY	Function containing unrestricted solution
fncX	Function for the limited solution

virtual void concepts::Operator< Real >::operator()	(	const Function< r_type > &	fncY,
		Function< Real > &	fncX
	)		`[virtual, inherited]`

Application operator for real function fncY.

Computes fncX = A(fncY) where A is this operator.

fncX becomes the type of the operator, for real data it becomes real, for complex data it becomes complex.

In derived classes its enough to implement the operator() for real Operator's. If a complex counterpart is not implemented, the function fncY is transformed to a complex function and then the application operator for complex functions is called.

If in a derived class the operator() for real Operator's is not implemented, a exception is thrown from here.

Reimplemented in concepts::VecOperator< Real >, concepts::DenseMatrix< Real >, concepts::DiagonalMatrix< Real >, concepts::Matrix< Real >, concepts::Permutation< Real >, and concepts::SparseMatrix< Real >.

void FvdgAdaptiveLimiter::reset_ ( ) [protected]

Resets element indicator to zero.

virtual const concepts::Space<Real>& FvdgAdaptiveLimiter::spaceX ( ) const [inline, virtual]

Returns the image space of the operator.

Definition at line 321 of file advection.hh.

virtual const concepts::Space<Real>& FvdgAdaptiveLimiter::spaceY ( ) const [inline, virtual]

Returns the source space of the operator.

Definition at line 325 of file advection.hh.