TensorProductQuadratureVector.hpp

// This file is part of Bembel, the higher order C++ boundary element library.
//
// Copyright (C) 2022 see <http://www.bembel.eu>
//
// It was written as part of a cooperation of J. Doelz, H. Harbrecht, S. Kurz,
// M. Multerer, S. Schoeps, and F. Wolf at Technische Universitaet Darmstadt,
// Universitaet Basel, and Universita della Svizzera italiana, Lugano. This
// source code is subject to the GNU General Public License version 3 and
// provided WITHOUT ANY WARRANTY, see <http://www.bembel.eu> for further
// information.
 
#ifndef BEMBEL_SRC_QUADRATURE_TENSORPRODUCTQUADRATUREVECTOR_HPP_
#define BEMBEL_SRC_QUADRATURE_TENSORPRODUCTQUADRATUREVECTOR_HPP_
 
namespace Bembel {
 
template <template <unsigned int qrOrder> class QuadratureRule,
          unsigned int Order>
struct TensorProductQuadratureVector {
  TensorProductQuadratureVector() {
    QuadratureRule<Order + 1> GL;
    Q_.xi_.resize(2, GL.xi_.size() * GL.xi_.size());
    Q_.w_.resize(GL.w_.size() * GL.w_.size());
    for (auto k = 0; k < Q_.xi_.cols(); ++k) {
      Q_.xi_.col(k) << GL.xi_[k / GL.xi_.size()], GL.xi_[k % GL.xi_.size()];
      Q_.w_(k) = GL.w_[k / GL.w_.size()] * GL.w_[k % GL.w_.size()];
    }
  }
  Cubature Q_;
  TensorProductQuadratureVector<QuadratureRule, Order - 1>
      remainingQuadratures_;
  const Cubature &operator[](unsigned int i) const {
    return (i == Order) ? Q_ : remainingQuadratures_[i];
  }
};
 
template <template <unsigned int qrOrder> class QuadratureRule>
struct TensorProductQuadratureVector<QuadratureRule, 0> {
  TensorProductQuadratureVector() {
    QuadratureRule<1> GL;
    Q_.xi_.resize(2, GL.xi_.size() * GL.xi_.size());
    Q_.w_.resize(GL.w_.size() * GL.w_.size());
    for (auto k = 0; k < Q_.xi_.cols(); ++k) {
      Q_.xi_.col(k) << GL.xi_[k / GL.xi_.size()], GL.xi_[k % GL.xi_.size()];
      Q_.w_(k) = GL.w_[k / GL.w_.size()] * GL.w_[k % GL.w_.size()];
    }
  }
  Cubature Q_;
  const Cubature &operator[](unsigned int i) const { return Q_; }
};
 
}  // namespace Bembel
#endif  // BEMBEL_SRC_QUADRATURE_TENSORPRODUCTQUADRATUREVECTOR_HPP_