meshpenalty.hh

#ifndef __DUNE_ACFEM_FUNCTIONS_MODULES_MESHPENALTY_HH__
#define __DUNE_ACFEM_FUNCTIONS_MODULES_MESHPENALTY_HH__
 
#include "../expressions.hh"
#include "elementmeasure.hh"
 
namespace Dune {
 
  namespace ACFem {
 
    namespace GridFunction {
 
      template<class GridPart>
      auto meshPenalty(const GridPart& gridPart)
      {
        auto exponent = -intFraction<1, GridPart::dimension>();
        return pow(elementMeasure(gridPart), exponent);
      }
 
 
    } // NS GridFunction
 
    using GridFunction::meshPenalty;
 
  } // NS ACFem
 
} // NS Dune
 
#endif //  __DUNE_ACFEM_FUNCTIONS_MODULES_MESHPENALTY_HH__