ErrorEstimator< ConfiguratorType > Class Template Reference

Class for local error estimation. More...

#include <errorEstimator.h>

Public Types
typedef ConfiguratorType::InitType	MeshType
typedef ConfiguratorType::RealType	RealType
typedef Eigen::VectorXd	VectorType
typedef Eigen::Matrix< RealType, 3, 1 >	VecType

Public Member Functions
	ErrorEstimator (MeshType &grid, RealType gamma=0.4)
RealType	markAndRefineTriangles ()
RealType	markTriangles ()

Protected Member Functions
RealType	f (const typename ConfiguratorType::Point3DType &x) const
void	computeLocalErrorEstimates (VectorType &errEstimates) const
RealType	jumpResidual (int t) const
RealType	getMidpointValue (int t) const

Protected Attributes
MeshType &	_mesh
RealType	_gamma

Detailed Description

template<typename ConfiguratorType>
class ErrorEstimator< ConfiguratorType >

Class for local error estimation.

Definition at line 20 of file errorEstimator.h.

Member Typedef Documentation

template<typename ConfiguratorType>

typedef ConfiguratorType::InitType ErrorEstimator< ConfiguratorType >::MeshType

Definition at line 23 of file errorEstimator.h.

template<typename ConfiguratorType>

typedef ConfiguratorType::RealType ErrorEstimator< ConfiguratorType >::RealType

Definition at line 24 of file errorEstimator.h.

template<typename ConfiguratorType>

typedef Eigen::VectorXd ErrorEstimator< ConfiguratorType >::VectorType

Definition at line 27 of file errorEstimator.h.

template<typename ConfiguratorType>

typedef Eigen::Matrix<RealType, 3, 1> ErrorEstimator< ConfiguratorType >::VecType

Definition at line 28 of file errorEstimator.h.

Constructor & Destructor Documentation

template<typename ConfiguratorType>

ErrorEstimator< ConfiguratorType >::ErrorEstimator ( MeshType &  grid, RealType  gamma = 0.4  )

inline

Definition at line 30 of file errorEstimator.h.

    : _mesh(grid), _gamma(gamma)
    {}

Member Function Documentation

template<typename ConfiguratorType>

void ErrorEstimator< ConfiguratorType >::computeLocalErrorEstimates ( VectorType &  errEstimates ) const

inlineprotected

Definition at line 69 of file errorEstimator.h.

    {
        for (int t = 0; t < _mesh.getNumTriangs(); ++t) {
            RealType h = _mesh.getTriang(t).getEdge(_mesh.getLongestEdgeIndex(t)).norm();
            RealType J_t = jumpResidual(t);
            errEstimates[t] = sqrt(h) * J_t;
        }
    }

template<typename ConfiguratorType>

RealType ErrorEstimator< ConfiguratorType >::f ( const typename ConfiguratorType::Point3DType &  x ) const

inlineprotected

Definition at line 55 of file errorEstimator.h.

    {
        RealType a = 0.8;
        RealType b = 0.0;
        int r = 0;

        for(int i = 0; i < 4; i++) {
            r = r || ( (x[1] > b) && (fabs(x[0]) + fabs(x[1]-b) < a) );
            b += a *= .75;
        }

        return 2.0 * static_cast<RealType>(r);
    }

template<typename ConfiguratorType>

RealType ErrorEstimator< ConfiguratorType >::getMidpointValue ( int  t ) const

inlineprotected

Definition at line 96 of file errorEstimator.h.

                                           {
        typename ConfiguratorType::ElementType &el = _mesh.getTriang(t);

        RealType res = 0.0;
        for (int i = 0; i < 3; ++i) {
            res += f(el.getNode(i));
        }
        res /= 3.0;

        return res;
    }

template<typename ConfiguratorType>

RealType ErrorEstimator< ConfiguratorType >::jumpResidual ( int  t ) const

inlineprotected

Definition at line 78 of file errorEstimator.h.

    {
        RealType Jt = 0.0;

        RealType val = getMidpointValue(t);
        for (int e = 0; e < 3; ++e) {
            int neighbour = _mesh.getNeighbour(t, e);
            // If no neighbour
            if (neighbour == -1)
                continue;

            RealType nVal = getMidpointValue(neighbour);
            Jt += fabs(val - nVal);
        }

        return Jt;
    }

template<typename ConfiguratorType>

RealType ErrorEstimator< ConfiguratorType >::markAndRefineTriangles ( )

inline

Definition at line 34 of file errorEstimator.h.

                                     {
        RealType err = markTriangles();
        _mesh.refineMarkedTriangles();
        return err;
    }

template<typename ConfiguratorType>

RealType ErrorEstimator< ConfiguratorType >::markTriangles ( )

inline

Definition at line 40 of file errorEstimator.h.

    {
        VectorType errEstimates(_mesh.getNumTriangs());

        computeLocalErrorEstimates(errEstimates);
        RealType refineAbove = errEstimates.lpNorm<Eigen::Infinity>() * _gamma;

        for (int t = 0; t < _mesh.getNumTriangs(); ++t)
            if (errEstimates[t] > refineAbove)
                _mesh.mark(t);

        return errEstimates.norm();
    }

Member Data Documentation

template<typename ConfiguratorType>

RealType ErrorEstimator< ConfiguratorType >::_gamma

protected

Definition at line 109 of file errorEstimator.h.

template<typename ConfiguratorType>

MeshType& ErrorEstimator< ConfiguratorType >::_mesh

protected

Definition at line 108 of file errorEstimator.h.

---

The documentation for this class was generated from the following file:

• exercise_2/errorEstimator.h