DuffyTrick/farFieldQuadratureNodes.hpp - GCC Code Coverage Report

Directory:	Bembel/src/
Coverage:	low: ≥ 0% medium: ≥ 75.0% high: ≥ 90.0%

	Coverage	Exec / Excl / Total
Lines:	93.8%	75 / 0 / 80
Functions:	100.0%	5 / 0 / 5
Branches:	61.1%	55 / 0 / 90

    DuffyTrick/farFieldQuadratureNodes.hpp
    
        Line
        Branch
        Exec
        Source
      
        // 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_DUFFYTRICK_FARFIELDQUADRATURENODES_HPP_
      
        #define BEMBEL_SRC_DUFFYTRICK_FARFIELDQUADRATURENODES_HPP_
      
        namespace Bembel {
      
        namespace DuffyTrick {
      
        /**
      
         * \ingroup DuffyTrick
      
         * \brief evaluates a given quadrature on all surface panels storage format is
      
         *qNodes.col(k) = [xi, w, Chi(xi); dsChi(xi); dtChi(xi)]
      
         */
      
        template <class T>
      
        9
        std::vector<ElementSurfacePoints> computeFfieldQnodes(const T &super_space,
      
                                                              const Cubature &Q) {
      
        9
          std::vector<ElementSurfacePoints> ffield_qnodes;
      
        9
          int next = 0;
      
          // assume isotropic mesh width h!
      
        9
          double h = (super_space.get_mesh().get_element_tree().cpbegin())->get_h();
      
        9
          auto nE = super_space.get_mesh().get_number_of_elements();
      
        9
          auto pbegin = super_space.get_mesh().get_element_tree().cpbegin();
      
        9
          auto pend = super_space.get_mesh().get_element_tree().cpend();
      
          5/10std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::DummyOperator> >(Bembel::SuperSpace<Bembel::DummyOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 1 taken 3 times.
✗ Branch 2 not taken.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> >(Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.

        9
          ffield_qnodes.reserve(nE);
      
          5/10std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::DummyOperator> >(Bembel::SuperSpace<Bembel::DummyOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 1 taken 3 times.
✗ Branch 2 not taken.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> >(Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.

        9
          SurfacePoint surfpt;
      
          10/10std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::DummyOperator> >(Bembel::SuperSpace<Bembel::DummyOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 2 taken 16 times.
✓ Branch 3 taken 1 time.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 2 taken 288 times.
✓ Branch 3 taken 3 times.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 2 taken 120 times.
✓ Branch 3 taken 2 times.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 2 taken 96 times.
✓ Branch 3 taken 1 time.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> >(Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 2 taken 102 times.
✓ Branch 3 taken 2 times.

        631
          for (auto it = pbegin; it != pend; ++it) {
      
          5/10std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::DummyOperator> >(Bembel::SuperSpace<Bembel::DummyOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 2 taken 16 times.
✗ Branch 3 not taken.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 2 taken 288 times.
✗ Branch 3 not taken.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 2 taken 120 times.
✗ Branch 3 not taken.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 2 taken 96 times.
✗ Branch 3 not taken.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> >(Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 2 taken 102 times.
✗ Branch 3 not taken.

        622
            ffield_qnodes.emplace_back(Q.xi_.cols());
      
          10/10std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::DummyOperator> >(Bembel::SuperSpace<Bembel::DummyOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 1 taken 4624 times.
✓ Branch 2 taken 16 times.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 1 taken 3936 times.
✓ Branch 2 taken 288 times.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 1 taken 1080 times.
✓ Branch 2 taken 120 times.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 1 taken 864 times.
✓ Branch 2 taken 96 times.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> >(Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 1 taken 1632 times.
✓ Branch 2 taken 102 times.

        12758
            for (auto k = 0; k < Q.xi_.cols(); ++k) {
      
              // the quadrature weight is scaled by mesh width
      
              // this corresponds to a scaling of the basis functions
      
              // with respect to the L^2 norm!
      
          20/40std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::DummyOperator> >(Bembel::SuperSpace<Bembel::DummyOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 1 taken 4624 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 4624 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 4624 times.
✗ Branch 8 not taken.
✓ Branch 11 taken 4624 times.
✗ Branch 12 not taken.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 1 taken 3936 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 3936 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 3936 times.
✗ Branch 8 not taken.
✓ Branch 11 taken 3936 times.
✗ Branch 12 not taken.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 1 taken 1080 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 1080 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 1080 times.
✗ Branch 8 not taken.
✓ Branch 11 taken 1080 times.
✗ Branch 12 not taken.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 1 taken 864 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 864 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 864 times.
✗ Branch 8 not taken.
✓ Branch 11 taken 864 times.
✗ Branch 12 not taken.
std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> >(Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> const&, Bembel::Quadrature<2u> const&):
✓ Branch 1 taken 1632 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 1632 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 1632 times.
✗ Branch 8 not taken.
✓ Branch 11 taken 1632 times.
✗ Branch 12 not taken.

        12136
              super_space.map2surface(*it, Q.xi_.col(k), h * Q.w_(k),
      
        12136
                                      &ffield_qnodes[it->id_][k]);
      
            }
      
          }
      
        18
          return ffield_qnodes;
      
        ✗
        }
      
        }  // namespace DuffyTrick
      
        }  // namespace Bembel
      
        #endif  // BEMBEL_SRC_DUFFYTRICK_FARFIELDQUADRATURENODES_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
12			#ifndef BEMBEL_SRC_DUFFYTRICK_FARFIELDQUADRATURENODES_HPP_
13			#define BEMBEL_SRC_DUFFYTRICK_FARFIELDQUADRATURENODES_HPP_
14
15			namespace Bembel {
16			namespace DuffyTrick {
17			/**
18			* \ingroup DuffyTrick
19			* \brief evaluates a given quadrature on all surface panels storage format is
20			*qNodes.col(k) = [xi, w, Chi(xi); dsChi(xi); dtChi(xi)]
21			*/
22			template <class T>
23		9	std::vector<ElementSurfacePoints> computeFfieldQnodes(const T &super_space,
24			const Cubature &Q) {
25		9	std::vector<ElementSurfacePoints> ffield_qnodes;
26		9	int next = 0;
27			// assume isotropic mesh width h!
28		9	double h = (super_space.get_mesh().get_element_tree().cpbegin())->get_h();
29		9	auto nE = super_space.get_mesh().get_number_of_elements();
30		9	auto pbegin = super_space.get_mesh().get_element_tree().cpbegin();
31		9	auto pend = super_space.get_mesh().get_element_tree().cpend();
32	5/10 std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::DummyOperator> >(Bembel::SuperSpace<Bembel::DummyOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 1 taken 1 time. ✗ Branch 2 not taken. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 1 taken 3 times. ✗ Branch 2 not taken. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 1 taken 2 times. ✗ Branch 2 not taken. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 1 taken 1 time. ✗ Branch 2 not taken. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> >(Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 1 taken 2 times. ✗ Branch 2 not taken.	9	ffield_qnodes.reserve(nE);
33	5/10 std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::DummyOperator> >(Bembel::SuperSpace<Bembel::DummyOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 1 taken 1 time. ✗ Branch 2 not taken. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 1 taken 3 times. ✗ Branch 2 not taken. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 1 taken 2 times. ✗ Branch 2 not taken. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 1 taken 1 time. ✗ Branch 2 not taken. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> >(Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 1 taken 2 times. ✗ Branch 2 not taken.	9	SurfacePoint surfpt;
34
35	10/10 std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::DummyOperator> >(Bembel::SuperSpace<Bembel::DummyOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 2 taken 16 times. ✓ Branch 3 taken 1 time. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 2 taken 288 times. ✓ Branch 3 taken 3 times. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 2 taken 120 times. ✓ Branch 3 taken 2 times. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 2 taken 96 times. ✓ Branch 3 taken 1 time. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> >(Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 2 taken 102 times. ✓ Branch 3 taken 2 times.	631	for (auto it = pbegin; it != pend; ++it) {
36	5/10 std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::DummyOperator> >(Bembel::SuperSpace<Bembel::DummyOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 2 taken 16 times. ✗ Branch 3 not taken. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 2 taken 288 times. ✗ Branch 3 not taken. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 2 taken 120 times. ✗ Branch 3 not taken. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 2 taken 96 times. ✗ Branch 3 not taken. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> >(Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 2 taken 102 times. ✗ Branch 3 not taken.	622	ffield_qnodes.emplace_back(Q.xi_.cols());
37	10/10 std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::DummyOperator> >(Bembel::SuperSpace<Bembel::DummyOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 1 taken 4624 times. ✓ Branch 2 taken 16 times. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 1 taken 3936 times. ✓ Branch 2 taken 288 times. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 1 taken 1080 times. ✓ Branch 2 taken 120 times. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 1 taken 864 times. ✓ Branch 2 taken 96 times. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> >(Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 1 taken 1632 times. ✓ Branch 2 taken 102 times.	12758	for (auto k = 0; k < Q.xi_.cols(); ++k) {
38			// the quadrature weight is scaled by mesh width
39			// this corresponds to a scaling of the basis functions
40			// with respect to the L^2 norm!
41	20/40 std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::DummyOperator> >(Bembel::SuperSpace<Bembel::DummyOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 1 taken 4624 times. ✗ Branch 2 not taken. ✓ Branch 4 taken 4624 times. ✗ Branch 5 not taken. ✓ Branch 7 taken 4624 times. ✗ Branch 8 not taken. ✓ Branch 11 taken 4624 times. ✗ Branch 12 not taken. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HelmholtzSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 1 taken 3936 times. ✗ Branch 2 not taken. ✓ Branch 4 taken 3936 times. ✗ Branch 5 not taken. ✓ Branch 7 taken 3936 times. ✗ Branch 8 not taken. ✓ Branch 11 taken 3936 times. ✗ Branch 12 not taken. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 1 taken 1080 times. ✗ Branch 2 not taken. ✓ Branch 4 taken 1080 times. ✗ Branch 5 not taken. ✓ Branch 7 taken 1080 times. ✗ Branch 8 not taken. ✓ Branch 11 taken 1080 times. ✗ Branch 12 not taken. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> >(Bembel::SuperSpace<Bembel::LaplaceSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 1 taken 864 times. ✗ Branch 2 not taken. ✓ Branch 4 taken 864 times. ✗ Branch 5 not taken. ✓ Branch 7 taken 864 times. ✗ Branch 8 not taken. ✓ Branch 11 taken 864 times. ✗ Branch 12 not taken. std::vector<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > >, std::allocator<std::vector<Eigen::Matrix<double, 12, 1, 0, 12, 1>, Eigen::aligned_allocator<Eigen::Matrix<double, 12, 1, 0, 12, 1> > > > > Bembel::DuffyTrick::computeFfieldQnodes<Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> >(Bembel::SuperSpace<Bembel::MaxwellSingleLayerOperator> const&, Bembel::Quadrature<2u> const&): ✓ Branch 1 taken 1632 times. ✗ Branch 2 not taken. ✓ Branch 4 taken 1632 times. ✗ Branch 5 not taken. ✓ Branch 7 taken 1632 times. ✗ Branch 8 not taken. ✓ Branch 11 taken 1632 times. ✗ Branch 12 not taken.	12136	super_space.map2surface(it, Q.xi_.col(k), h Q.w_(k),
42		12136	&ffield_qnodes[it->id_][k]);
43			}
44			}
45		18	return ffield_qnodes;
46		✗	}
47			} // namespace DuffyTrick
48			} // namespace Bembel
49			#endif // BEMBEL_SRC_DUFFYTRICK_FARFIELDQUADRATURENODES_HPP_
50