To calculate an a posteriori estimator of the solution. More...

#include <Estimator.h>

Public Types

enum  EstimatorType {
  ESTIM_ZZ = 0 ,
  ESTIM_ND_JUMP = 1
}
 

Public Member Functions

 Estimator ()
 Default Constructor.
 
 Estimator (Mesh &m)
 Constructor using finite element mesh.
 
 ~Estimator ()
 Destructor.
 
void setType (EstimatorType t=ESTIM_ZZ)
 Select type of a posteriori estimator.
 
void setSolution (const Vect< real_t > &u)
 Provide solution vector in order to determine error index.
 
void getElementWiseIndex (Vect< real_t > &e)
 Get vector containing elementwise error index.
 
void getNodeWiseIndex (Vect< real_t > &e)
 Get vector containing nodewise error index.
 
void getSideWiseIndex (Vect< real_t > &e)
 Get vector containing sidewise error index.
 
real_t getAverage () const
 Return averaged error.
 
MeshgetMesh () const
 Return a reference to the finite element mesh.
 

Detailed Description

To calculate an a posteriori estimator of the solution.

This class enables calculating an estimator of a solution in order to evaluate reliability. Estimation uses the so-called Zienkiewicz-Zhu estimator.

Author
Rachid Touzani

Member Enumeration Documentation

◆ EstimatorType

Enumerate variable that selects an error estimator for mesh adaptation purposes

Enumerator
ESTIM_ZZ 

Zhu-Zienckiewicz elementwise estimator

ESTIM_ND_JUMP 

Normal derivative jump sidewise estimator

Constructor & Destructor Documentation

◆ Estimator()

Estimator ( Mesh m)

Constructor using finite element mesh.

Parameters
[in]mMesh instance

Member Function Documentation

◆ getElementWiseIndex()

void getElementWiseIndex ( Vect< real_t > &  e)

Get vector containing elementwise error index.

Parameters
[in,out]eVector that contains once the member function setError is invoked a posteriori estimator at each element

◆ getNodeWiseIndex()

void getNodeWiseIndex ( Vect< real_t > &  e)

Get vector containing nodewise error index.

Parameters
[in,out]eVector that contains once the member function setError is invoked a posteriori estimator at each node

◆ getSideWiseIndex()

void getSideWiseIndex ( Vect< real_t > &  e)

Get vector containing sidewise error index.

Parameters
[in,out]eVector that contains once the member function setError is invoked a posteriori estimator at each side

◆ setSolution()

void setSolution ( const Vect< real_t > &  u)

Provide solution vector in order to determine error index.

Parameters
[in]uVector containing solution at mesh nodes

◆ setType()

void setType ( EstimatorType  t = ESTIM_ZZ)

Select type of a posteriori estimator.

Parameters
[in]tType of estimator. It has to be chosen among the enumerated values:
  • ESTIM_ZZ: The Zhu-Zienckiewicz estimator (Default value)
  • ESTIM_ND_JUMP: An estimator based on the jump of normal derivatives of the solution across mesh sides