subentityquadraturerule.hh

Go to the documentation of this file.

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
 
// SPDX-FileCopyrightText: Copyright © DUNE-FUFEM Project contributors, see file AUTHORS.md
// SPDX-License-Identifier: LicenseRef-GPL-2.0-only-with-DUNE-exception OR LGPL-3.0-or-later
 
#ifndef DUNE_FUFEM_QUADRATURERULES_SUBENTITYQUADRATURERULE_HH
#define DUNE_FUFEM_QUADRATURERULES_SUBENTITYQUADRATURERULE_HH
 
#include <dune/geometry/type.hh>
#include <dune/geometry/referenceelements.hh>
#include <dune/geometry/quadraturerules.hh>
 
namespace Dune::Fufem {
 
template <class ct, int dim, int codim>
class SubEntityQuadratureRule:
  public Dune::QuadratureRule<ct,dim>
{
  using Base = Dune::QuadratureRule<ct,dim>;
 
public:
 
  using Base::type;
 
  SubEntityQuadratureRule(Dune::GeometryType gt, std::size_t subEntityIndex, const Dune::QuadratureRule<ct, dim-codim>& rule)
    : Base(gt),
      subEntityRule_(&rule)
  {
    [[maybe_unused]] auto subEntityGeometryType = Dune::referenceElement<ct,dim>(type()).type(subEntityIndex, codim);
    assert(subEntityGeometryType == subEntityRule().type());
    auto subEntityGeometry = Dune::referenceElement<ct,dim>(type()).template geometry<codim>(subEntityIndex);
    for(auto qp : subEntityRule())
      this->emplace_back(subEntityGeometry.global(qp.position()), qp.weight());
    this->delivered_order = subEntityRule().order();
  }
 
  const Dune::QuadratureRule<ct, dim-1>& subEntityRule() const
  {
    return *subEntityRule_;
  }
 
protected:
  const Dune::QuadratureRule<ct, dim-1>* subEntityRule_;
};
 
 
} // namespace Dune::Fufem
 
#endif // DUNE_FUFEM_QUADRATURERULES_SUBENTITYQUADRATURERULE_HH