errorEstimator.h

Go to the documentation of this file.

#ifndef EXERCISE_3_ERRORESTIMATOR_H_
#define EXERCISE_3_ERRORESTIMATOR_H_

#include <typeinfo>
#include <cxxabi.h>

#include <Eigen/Dense>

template<typename ConfiguratorType>
class ErrorEstimator
{
public:
    typedef typename ConfiguratorType::InitType       MeshType;
    typedef typename ConfiguratorType::RealType       RealType;
    typedef typename ConfiguratorType::DomVecType     DomVecType;
    typedef typename ConfiguratorType::TangentVecType TangentVecType;
    // Strangely, lpNorm fails if this is set to
    // "typename ConfiguratorType::VectorType". Bug?
    typedef typename Eigen::VectorXd VectorType;
    typedef typename Eigen::Matrix<RealType, 3, 1> VecType;

    ErrorEstimator(MeshType &mesh, const ConfiguratorType &config,
            const VectorType &numericalSolution,
            RealType gamma = 0.4)
    : _mesh(mesh),
      _numericalSolution(config, numericalSolution), _gamma(gamma)
    {}

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

    RealType markTriangles()
    {
        VectorType errEstimates(_mesh.getNumTriangs());

        computeLocalErrorEstimates(errEstimates);


        return errEstimates.norm();
    }

    void computeLocalErrorEstimates(VectorType &errEstimates) const
    {
    }

    void computeH(VectorType &h) const
    {
        for (int t = 0; t < _mesh.getNumTriangs(); ++t) {
            h[t] = _mesh.getTriang(t).getEdge(_mesh.getLongestEdgeIndex(t)).norm();
        }
    }

    MeshType& getMesh()
    {
        return _mesh;
    }

protected:

    MeshType &_mesh;
    const shellFE::DiscreteFunctionDefaultShellFE<ConfiguratorType> _numericalSolution;
    RealType _gamma;
};

#endif /* EXERCISE_3_ERRORESTIMATOR_H_ */