Inheritance diagram for estimator::ImplicitResidual< F >:

Public Member Functions
virtual ImplicitResidual< F > *	clone () const =0

void	addBoundaryData (const concepts::Set< uint > &nSet, const concepts::ElementFormula< F > &g)

void	addBoundaryData (const concepts::Set< uint > &rSet, const concepts::ElementFormula< F > &h1, const concepts::ElementFormula< F > &h2)

Protected Member Functions
	ImplicitResidual (const concepts::SpaceOnCells< F > &spc, const concepts::Vector< F > &sol)

virtual std::ostream &	info (std::ostream &os) const

Protected Attributes
concepts::Sequence< const concepts::ElementFormula< F > * >	g_

concepts::Sequence< concepts::Set< uint > >	gAttrbs_

concepts::Set< uint >	eAttrbs_

template<class F>
class estimator::ImplicitResidual< F >

Abstract class representing the implicit residual a posteriori error estimator.

At the moment just one underlying space is supported, i.e. with that robin conditions cannot be handeled for the residuals, atm

This uses the flux approximation proposed by Bank and Weise, that is

n_k \cdot [\mathcal{A} Grad(u)] = 0.5 \cdot n_K ( \mathcal{A}Grad[u_h]_{|K} + \mathcal{A}Grad[u_h]_{K'} )

where n_k is the outer normal on a cell K, and (K, K') are neighboured Elements.

Definition at line 37 of file implicitResidual.hh.

template<class F >

estimator::ImplicitResidual< F >::ImplicitResidual	(	const concepts::SpaceOnCells< F > &	spc,
		const concepts::Vector< F > &	sol
	)

inlineprotected

Constructur of the explicit residual a posteriori error estimator.

The input residual is builded outside the class. For given elliptic PDE Lu = f, res = f - L(u_h), where u_h is the FEM approximation.

Example: L = -Laplace + c(x)*ID, use hp2D::Functions as Laplacian, Value and ParsedFormulas

Parameters

spc	Underlying 2d Space, i.e. a AdaptiveSpace
sol	solution vector of the approximated FEM solution
res	strong Residual.

Definition at line 125 of file implicitResidual.hh.

template<class F >

void estimator::ImplicitResidual< F >::addBoundaryData	(	const concepts::Set< uint > &	nSet,
		const concepts::ElementFormula< F > &	g
	)

inline

Add boundary informations (i.e. on point/edge/face) to build boundary residuals on a given requested set of edges through edge attributes.

This applicates i.e for Neumann b.c. :

Let the boundary condition be of the form :

 n^T * [A*Grad(u)] = g

Then the L2 edge residual representation has the form :

$|| g - n_e^T * [A*Grad(u_{hp})] ||_{L^2(e)}$

where n_e is the outward pointing normal vector, that comes from the underlying element that has edge e.

Definition at line 64 of file implicitResidual.hh.

template<class F >

void estimator::ImplicitResidual< F >::addBoundaryData	(	const concepts::Set< uint > &	rSet,
		const concepts::ElementFormula< F > &	h1,
		const concepts::ElementFormula< F > &	h2
	)