tensorproductquadrule.hh

Go to the documentation of this file.

// SPDX-FileCopyrightText: 2021-2024 The Ikarus Developers mueller@ibb.uni-stuttgart.de
// SPDX-License-Identifier: LGPL-3.0-or-later
 
#pragma once
 
#include <dune/geometry/quadraturerules.hh>
namespace Ikarus {
 
template <class BaseQuadrature, class Quadrature>
auto tensorProductQuadrature(const BaseQuadrature& baseQuad, const Quadrature& onedQuad) {
  constexpr int baseQuadDim       = BaseQuadrature::d;
  auto rule                       = Dune::QuadratureRule<double, baseQuadDim + 1>();
  const unsigned int baseQuadSize = baseQuad.size();
  for (unsigned int bqi = 0; bqi < baseQuadSize; ++bqi) {
    const typename Dune::QuadraturePoint<double, baseQuadDim>::Vector& basePoint = baseQuad[bqi].position();
    const double& baseWeight                                                     = baseQuad[bqi].weight();
 
    typename Dune::QuadraturePoint<double, baseQuadDim + 1>::Vector point;
    for (unsigned int i = 0; i < baseQuadDim; ++i)
      point[i] = basePoint[i];
 
    const unsigned int onedQuadSize = onedQuad.size();
    for (unsigned int oqi = 0; oqi < onedQuadSize; ++oqi) {
      point[baseQuadDim] = onedQuad[oqi].position()[0];
      rule.emplace_back(Dune::QuadraturePoint(point, baseWeight * onedQuad[oqi].weight()));
    }
  }
  return rule;
}
} // namespace Ikarus