GCC Code Coverage Report

Directory: Bembel/src/
Coverage: low: ≥ 0% medium: ≥ 75.0% high: ≥ 90.0%
Coverage Exec / Excl / Total
Lines: 20.0% 10 / 0 / 50
Functions: 20.0% 1 / 0 / 5
Branches: 11.4% 16 / 0 / 140
DuffyTrick/integrate1.hpp
Line Branch Exec Source
1 // This file is part of Bembel, the higher order C++ boundary element library.
2 //
3 // Copyright (C) 2022 see <http://www.bembel.eu>
4 //
5 // It was written as part of a cooperation of J. Doelz, H. Harbrecht, S. Kurz,
6 // M. Multerer, S. Schoeps, and F. Wolf at Technische Universitaet Darmstadt,
7 // Universitaet Basel, and Universita della Svizzera italiana, Lugano. This
8 // source code is subject to the GNU General Public License version 3 and
9 // provided WITHOUT ANY WARRANTY, see <http://www.bembel.eu> for further
10 // information.
11 #ifndef BEMBEL_SRC_DUFFYTRICK_INTEGRATE1_HPP_
12 #define BEMBEL_SRC_DUFFYTRICK_INTEGRATE1_HPP_
13
14 namespace Bembel {
15 namespace DuffyTrick {
16 /**
17 * \ingroup DuffyTrick
18 * \brief no-problem quadrature routine, elements are sufficiently far away from
19 *each other, but not far-field yet (this is just an isotropic tensor product
20 *quadrature!).
21 *
22 * \todo be sure that map2element computes the weight h*Q.w(i) such that the
23 *integrand may then be scaled by qp1.weight * qp2.weight
24 **/
25 template <typename Derived, class T>
26 256 void integrate1(const LinearOperatorBase<Derived> &LinOp, const T &super_space,
27 const ElementTreeNode &e1, int rot1, const ElementTreeNode &e2,
28 int rot2, const ElementSurfacePoints &ffield_qnodes1,
29 const ElementSurfacePoints &ffield_qnodes2, const Cubature &Q,
30 Eigen::Matrix<typename LinearOperatorTraits<Derived>::Scalar,
31 Eigen::Dynamic, Eigen::Dynamic> *intval) {
32 256 double h = e1.get_h();
33
1/10
void Bembel::DuffyTrick::integrate1<Bembel::DummyOperator, Bembel::SuperSpace<Bembel::DummyOperator> >(Bembel::LinearOperatorBase<Bembel::DummyOperator> const&, Bembel::SuperSpace<Bembel::DummyOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::DummyOperator>::Scalar, -1, -1, 0, -1, -1>*):
✓ Branch 1 taken 256 times.
✗ Branch 2 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::HelmholtzSingleLayerOperator, Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::HelmholtzSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::HomogenisedLaplaceSingleLayerOperator, Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::HomogenisedLaplaceSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::LaplaceSingleLayerOperator, Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::LaplaceSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::LaplaceSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::MaxwellSingleLayerOperator, Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::MaxwellSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::MaxwellSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
 256 intval->setZero();
34
2/20
void Bembel::DuffyTrick::integrate1<Bembel::DummyOperator, Bembel::SuperSpace<Bembel::DummyOperator> >(Bembel::LinearOperatorBase<Bembel::DummyOperator> const&, Bembel::SuperSpace<Bembel::DummyOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::DummyOperator>::Scalar, -1, -1, 0, -1, -1>*):
✓ Branch 1 taken 256 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 256 times.
✗ Branch 5 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::HelmholtzSingleLayerOperator, Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::HelmholtzSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::HomogenisedLaplaceSingleLayerOperator, Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::HomogenisedLaplaceSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::LaplaceSingleLayerOperator, Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::LaplaceSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::LaplaceSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::MaxwellSingleLayerOperator, Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::MaxwellSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::MaxwellSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
 256 SurfacePoint qp1, qp2;
35
2/10
void Bembel::DuffyTrick::integrate1<Bembel::DummyOperator, Bembel::SuperSpace<Bembel::DummyOperator> >(Bembel::LinearOperatorBase<Bembel::DummyOperator> const&, Bembel::SuperSpace<Bembel::DummyOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::DummyOperator>::Scalar, -1, -1, 0, -1, -1>*):
✓ Branch 1 taken 73984 times.
✓ Branch 2 taken 256 times.
void Bembel::DuffyTrick::integrate1<Bembel::HelmholtzSingleLayerOperator, Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::HelmholtzSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::HomogenisedLaplaceSingleLayerOperator, Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::HomogenisedLaplaceSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::LaplaceSingleLayerOperator, Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::LaplaceSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::LaplaceSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::MaxwellSingleLayerOperator, Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::MaxwellSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::MaxwellSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
 74240 for (auto i = 0; i < Q.w_.size(); ++i) {
36
4/40
void Bembel::DuffyTrick::integrate1<Bembel::DummyOperator, Bembel::SuperSpace<Bembel::DummyOperator> >(Bembel::LinearOperatorBase<Bembel::DummyOperator> const&, Bembel::SuperSpace<Bembel::DummyOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::DummyOperator>::Scalar, -1, -1, 0, -1, -1>*):
✓ Branch 1 taken 73984 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 73984 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 73984 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 73984 times.
✗ Branch 11 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::HelmholtzSingleLayerOperator, Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::HelmholtzSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::HomogenisedLaplaceSingleLayerOperator, Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::HomogenisedLaplaceSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::LaplaceSingleLayerOperator, Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::LaplaceSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::LaplaceSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::MaxwellSingleLayerOperator, Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::MaxwellSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::MaxwellSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
 73984 super_space.map2surface(e1, Q.xi_.col(i), h * Q.w_(i), &qp1);
37
2/10
void Bembel::DuffyTrick::integrate1<Bembel::DummyOperator, Bembel::SuperSpace<Bembel::DummyOperator> >(Bembel::LinearOperatorBase<Bembel::DummyOperator> const&, Bembel::SuperSpace<Bembel::DummyOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::DummyOperator>::Scalar, -1, -1, 0, -1, -1>*):
✓ Branch 1 taken 21381376 times.
✓ Branch 2 taken 73984 times.
void Bembel::DuffyTrick::integrate1<Bembel::HelmholtzSingleLayerOperator, Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::HelmholtzSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::HomogenisedLaplaceSingleLayerOperator, Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::HomogenisedLaplaceSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::LaplaceSingleLayerOperator, Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::LaplaceSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::LaplaceSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::MaxwellSingleLayerOperator, Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::MaxwellSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::MaxwellSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
 21455360 for (auto j = 0; j < Q.w_.size(); ++j) {
38
4/40
void Bembel::DuffyTrick::integrate1<Bembel::DummyOperator, Bembel::SuperSpace<Bembel::DummyOperator> >(Bembel::LinearOperatorBase<Bembel::DummyOperator> const&, Bembel::SuperSpace<Bembel::DummyOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::DummyOperator>::Scalar, -1, -1, 0, -1, -1>*):
✓ Branch 1 taken 21381376 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 21381376 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 21381376 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 21381376 times.
✗ Branch 11 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::HelmholtzSingleLayerOperator, Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::HelmholtzSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::HomogenisedLaplaceSingleLayerOperator, Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::HomogenisedLaplaceSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::LaplaceSingleLayerOperator, Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::LaplaceSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::LaplaceSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::MaxwellSingleLayerOperator, Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::MaxwellSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::MaxwellSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
 21381376 super_space.map2surface(e2, Q.xi_.col(j), h * Q.w_(j), &qp2);
39
1/10
void Bembel::DuffyTrick::integrate1<Bembel::DummyOperator, Bembel::SuperSpace<Bembel::DummyOperator> >(Bembel::LinearOperatorBase<Bembel::DummyOperator> const&, Bembel::SuperSpace<Bembel::DummyOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::DummyOperator>::Scalar, -1, -1, 0, -1, -1>*):
✓ Branch 1 taken 21381376 times.
✗ Branch 2 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::HelmholtzSingleLayerOperator, Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::HelmholtzSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::HomogenisedLaplaceSingleLayerOperator, Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::HomogenisedLaplaceSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::LaplaceSingleLayerOperator, Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::LaplaceSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::LaplaceSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
void Bembel::DuffyTrick::integrate1<Bembel::MaxwellSingleLayerOperator, Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> >(Bembel::LinearOperatorBase<Bembel::MaxwellSingleLayerOperator> const&, Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> const&, Bembel::ElementTreeNode const&, int, Bembel::ElementTreeNode const&, int, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > const&, Bembel::Quadrature<2u> const&, Eigen::Matrix<Bembel::LinearOperatorTraits<Bembel::MaxwellSingleLayerOperator>::Scalar, -1, -1, 0, -1, -1>*):
✗ Branch 1 not taken.
✗ Branch 2 not taken.
 21381376 LinOp.evaluateIntegrand(super_space, qp1, qp2, intval);
40 }
41 }
42
43 BEMBEL_UNUSED_(rot1);
44 BEMBEL_UNUSED_(rot2);
45 BEMBEL_UNUSED_(ffield_qnodes1);
46 BEMBEL_UNUSED_(ffield_qnodes2);
47 512 return;
48 }
49 } // namespace DuffyTrick
50 } // namespace Bembel
51 #endif // BEMBEL_SRC_DUFFYTRICK_INTEGRATE1_HPP_
52