geometry.hh

#ifndef __DUNE_ACFEM_GEOMETRY_HH__
#define __DUNE_ACFEM_GEOMETRY_HH__
 
#include <cmath>
 
namespace Dune {
 
  namespace ACFem {
 
    template<class Geometry>
    typename Geometry::ctype h2Estimate(const Geometry& geometry)
    {
      const int dimension = Geometry::mydimension;
      typename Geometry::ctype volume = geometry.volume();
 
      switch(dimension) {
      case 1:
        return volume*volume;
      case 2:
        return volume;
      case 3:
        return std::pow(volume, 2.0/(double)dimension);
      default:
        return -1.0;
      }
    }
 
    template<class Geometry>
    typename Geometry::ctype hEstimate(const Geometry& geometry)
    {
      typedef typename Geometry::ctype ctype;
      const int dimension = Geometry::mydimension;
      const ctype volume = geometry.volume();
 
      switch(dimension) {
      case 1:
        return volume;
      case 2:
        return std::sqrt(volume);
      case 3:
        return std::pow(volume, 1.0/(double)dimension);
      default:
        return -1.0;
      }
    }
 
  } // namespace ACFem
 
} // namespace Dune
 
#endif // __DUNE_ACFEM_GEOMETRY_HH__