quadraturerules.hh

#ifndef DUNE_PDELAB_COMMON_QUADRATURERULES_HH
#define DUNE_PDELAB_COMMON_QUADRATURERULES_HH
 
#include <dune/geometry/quadraturerules.hh>
 
namespace Dune {
  namespace PDELab {
 
 
    template<typename QR>
    class QuadratureRuleWrapper
    {
 
    public:
 
      using CoordType = typename QR::CoordType;
 
      using size_type = typename QR::size_type;
 
      using const_iterator = typename QR::const_iterator;
 
      using iterator = const_iterator;
 
      int order() const
      {
        return _quadrature_rule->order();
      }
 
      GeometryType type() const
      {
        return _quadrature_rule->type();
      }
 
      size_type size() const
      {
        return _quadrature_rule->size();
      }
 
      const_iterator begin() const
      {
        return _quadrature_rule->begin();
      }
 
      const_iterator end() const
      {
        return _quadrature_rule->end();
      }
 
#ifndef DOXYGEN
 
      QuadratureRuleWrapper(const QR& quadrature_rule)
        : _quadrature_rule(&quadrature_rule)
      {}
 
#endif
 
    private:
 
      const QR* _quadrature_rule;
 
    };
 
 
    template<typename Geometry>
    QuadratureRuleWrapper<
      QuadratureRule<
        typename Geometry::ctype,
        Geometry::mydimension
        >
      >
    quadratureRule(const Geometry& geo, std::size_t order, QuadratureType::Enum quadrature_type = QuadratureType::GaussLegendre)
    {
      return { QuadratureRules<typename Geometry::ctype,Geometry::mydimension>::rule(geo.type(),order,quadrature_type) };
    }
 
  } // namespace PDELab
 
  // inject the function into the Dune namespace to enable ADL.
  // TODO: Remove this injection once the core modules gain their own version of this function!
  using Dune::PDELab::quadratureRule;
 
} // namespace Dune
 
#endif // DUNE_PDELAB_COMMON_QUADRATURERULES_HH