geometry.hh

#ifndef DUNE_FEM_QUADRATURE_GEOMETRIC_GEOMETRY_HH
#define DUNE_FEM_QUADRATURE_GEOMETRIC_GEOMETRY_HH
 
#include <cstddef>
 
#include <functional>
 
#include <dune/geometry/type.hh>
 
#include "quadrature.hh"
 
namespace Dune
{
 
  namespace Fem
  {
 
    // GeometryQuadrature
    // ------------------
 
    template< class QuadratureRule >
    class GeometryQuadrature
    : public GeometricQuadrature< typename QuadratureRule::CoordType, QuadratureRule::d, QuadratureRule::d, GeometryQuadrature< QuadratureRule > >
    {
      using ThisType = GeometryQuadrature< QuadratureRule >;
      using BaseType = GeometricQuadrature< typename QuadratureRule::CoordType, QuadratureRule::d, QuadratureRule::d, GeometryQuadrature< QuadratureRule > >;
 
    public:
      using QuadratureRuleType = QuadratureRule;
 
      using FieldType = typename BaseType::FieldType;
 
      using CoordinateType = typename BaseType::CoordinateType;
      using LocalCoordinateType = typename BaseType::LocalCoordinateType;
 
      explicit GeometryQuadrature ( const QuadratureRuleType &quadratureRule )
        : quadratureRule_( quadratureRule )
      {}
 
      GeometryQuadrature ( const ThisType & ) = default;
 
      GeometryQuadrature ( ThisType && ) = default;
 
      GeometryQuadrature &operator= ( const ThisType & ) = default;
 
      GeometryQuadrature &operator= ( ThisType && ) = default;
 
      Dune::GeometryType type () const { return quadratureRule().type(); }
 
      int order () const { return quadratureRule().order(); }
 
      std::size_t nop () const { return quadratureRule().size(); }
 
      const CoordinateType &point ( std::size_t i ) const
      {
        return quadratureRule()[ i ].position();
      }
 
      const LocalCoordinateType &localPoint ( std::size_t i ) const
      {
        return point( i );
      }
 
      FieldType weight ( std::size_t i ) const
      {
        return quadratureRule()[ i ].weight();
      }
 
    private:
      const QuadratureRuleType &quadratureRule () const { return quadratureRule_.get(); }
 
      std::reference_wrapper< const QuadratureRuleType > quadratureRule_;
    };
 
  } // namespace Fem
 
} // namespace Dune
 
#endif // #ifndef DUNE_FEM_QUADRATURE_GEOMETRIC_GEOMETRY_HH