LinearForm/DiscreteLinearForm.hpp - GCC Code Coverage Report

Directory:	Bembel/src/
Coverage:	low: ≥ 0% medium: ≥ 75.0% high: ≥ 90.0%
	Coverage	Exec / Excl / Total
Lines:	100.0%	273 / 0 / 273
Functions:	100.0%	35 / 0 / 35
Branches:	53.6%	210 / 0 / 392
    LinearForm/DiscreteLinearForm.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_LINEARFORM_DISCRETELINEARFORM_HPP_
      
        #define BEMBEL_SRC_LINEARFORM_DISCRETELINEARFORM_HPP_
      
        namespace Bembel {
      
        /**
      
         *  \ingroup LinearForm
      
         *  \brief This class, given a LinearForm with defined traits, takes care of the
      
         * assembly of the right hand side.
      
         */
      
        template <typename Derived, typename LinOp>
      
        class DiscreteLinearForm {
      
         public:
      
          //////////////////////////////////////////////////////////////////////////////
      
          //    constructors
      
          //////////////////////////////////////////////////////////////////////////////
      
          DiscreteLinearForm() {}
      
          14/28Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::HomogenisedLaplaceSingleLayerOperator>::DiscreteLinearForm(Bembel::AnsatzSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&):
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✓ Branch 5 taken 2 times.
✗ Branch 6 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceBeltramiOperator>::DiscreteLinearForm(Bembel::AnsatzSpace<Bembel::LaplaceBeltramiOperator> const&):
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✓ Branch 5 taken 2 times.
✗ Branch 6 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceSingleLayerOperator>::DiscreteLinearForm(Bembel::AnsatzSpace<Bembel::LaplaceSingleLayerOperator> const&):
✓ Branch 2 taken 1 time.
✗ Branch 3 not taken.
✓ Branch 5 taken 1 time.
✗ Branch 6 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::MassMatrixScalarCont>::DiscreteLinearForm(Bembel::AnsatzSpace<Bembel::MassMatrixScalarCont> const&):
✓ Branch 2 taken 1 time.
✗ Branch 3 not taken.
✓ Branch 5 taken 1 time.
✗ Branch 6 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<std::complex<double> >, Bembel::HelmholtzSingleLayerOperator>::DiscreteLinearForm(Bembel::AnsatzSpace<Bembel::HelmholtzSingleLayerOperator> const&):
✓ Branch 2 taken 1 time.
✗ Branch 3 not taken.
✓ Branch 5 taken 1 time.
✗ Branch 6 not taken.
Bembel::DiscreteLinearForm<Bembel::RotatedTangentialTrace<std::complex<double> >, Bembel::MaxwellSingleLayerOperator>::DiscreteLinearForm(Bembel::AnsatzSpace<Bembel::MaxwellSingleLayerOperator> const&):
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✓ Branch 5 taken 2 times.
✗ Branch 6 not taken.
Bembel::DiscreteLinearForm<Bembel::TangentialTrace<double>, TestOperatorDivC>::DiscreteLinearForm(Bembel::AnsatzSpace<TestOperatorDivC> const&):
✓ Branch 2 taken 1 time.
✗ Branch 3 not taken.
✓ Branch 5 taken 1 time.
✗ Branch 6 not taken.

        10
          explicit DiscreteLinearForm(const AnsatzSpace<LinOp> &ansatz_space) {
      
          7/14Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::HomogenisedLaplaceSingleLayerOperator>::DiscreteLinearForm(Bembel::AnsatzSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&):
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceBeltramiOperator>::DiscreteLinearForm(Bembel::AnsatzSpace<Bembel::LaplaceBeltramiOperator> const&):
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceSingleLayerOperator>::DiscreteLinearForm(Bembel::AnsatzSpace<Bembel::LaplaceSingleLayerOperator> const&):
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::MassMatrixScalarCont>::DiscreteLinearForm(Bembel::AnsatzSpace<Bembel::MassMatrixScalarCont> const&):
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<std::complex<double> >, Bembel::HelmholtzSingleLayerOperator>::DiscreteLinearForm(Bembel::AnsatzSpace<Bembel::HelmholtzSingleLayerOperator> const&):
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::RotatedTangentialTrace<std::complex<double> >, Bembel::MaxwellSingleLayerOperator>::DiscreteLinearForm(Bembel::AnsatzSpace<Bembel::MaxwellSingleLayerOperator> const&):
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::TangentialTrace<double>, TestOperatorDivC>::DiscreteLinearForm(Bembel::AnsatzSpace<TestOperatorDivC> const&):
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.

        10
            init_DiscreteLinearForm(ansatz_space);
      
        10
          }
      
          //////////////////////////////////////////////////////////////////////////////
      
          //    init_DiscreteLinearForm
      
          //////////////////////////////////////////////////////////////////////////////
      
        10
          void init_DiscreteLinearForm(const AnsatzSpace<LinOp> &ansatz_space) {
      
          14/28Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::HomogenisedLaplaceSingleLayerOperator>::init_DiscreteLinearForm(Bembel::AnsatzSpace<Bembel::HomogenisedLaplaceSingleLayerOperator> const&):
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceBeltramiOperator>::init_DiscreteLinearForm(Bembel::AnsatzSpace<Bembel::LaplaceBeltramiOperator> const&):
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceSingleLayerOperator>::init_DiscreteLinearForm(Bembel::AnsatzSpace<Bembel::LaplaceSingleLayerOperator> const&):
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 time.
✗ Branch 5 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::MassMatrixScalarCont>::init_DiscreteLinearForm(Bembel::AnsatzSpace<Bembel::MassMatrixScalarCont> const&):
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 time.
✗ Branch 5 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<std::complex<double> >, Bembel::HelmholtzSingleLayerOperator>::init_DiscreteLinearForm(Bembel::AnsatzSpace<Bembel::HelmholtzSingleLayerOperator> const&):
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 time.
✗ Branch 5 not taken.
Bembel::DiscreteLinearForm<Bembel::RotatedTangentialTrace<std::complex<double> >, Bembel::MaxwellSingleLayerOperator>::init_DiscreteLinearForm(Bembel::AnsatzSpace<Bembel::MaxwellSingleLayerOperator> const&):
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
Bembel::DiscreteLinearForm<Bembel::TangentialTrace<double>, TestOperatorDivC>::init_DiscreteLinearForm(Bembel::AnsatzSpace<TestOperatorDivC> const&):
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 time.
✗ Branch 5 not taken.

        10
            ansatz_space_ = ansatz_space;
      
        10
            deg_ = ansatz_space_.get_polynomial_degree() + 1;
      
        10
            return;
      
          }
      
          //////////////////////////////////////////////////////////////////////////////
      
          //    compute
      
          //////////////////////////////////////////////////////////////////////////////
      
          /**
      
           *    \todo Add inline commentary
      
           **/
      
        10
          void compute() {
      
          7/14Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::HomogenisedLaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceBeltramiOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::MassMatrixScalarCont>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<std::complex<double> >, Bembel::HelmholtzSingleLayerOperator>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::RotatedTangentialTrace<std::complex<double> >, Bembel::MaxwellSingleLayerOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::TangentialTrace<double>, TestOperatorDivC>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.

        10
            GaussSquare<Constants::maximum_quadrature_degree> GS;
      
          14/28Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::HomogenisedLaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceBeltramiOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 time.
✗ Branch 5 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::MassMatrixScalarCont>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 time.
✗ Branch 5 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<std::complex<double> >, Bembel::HelmholtzSingleLayerOperator>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 time.
✗ Branch 5 not taken.
Bembel::DiscreteLinearForm<Bembel::RotatedTangentialTrace<std::complex<double> >, Bembel::MaxwellSingleLayerOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
Bembel::DiscreteLinearForm<Bembel::TangentialTrace<double>, TestOperatorDivC>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 time.
✗ Branch 5 not taken.

        10
            auto Q = GS[deg_];
      
          7/14Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::HomogenisedLaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceBeltramiOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::MassMatrixScalarCont>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<std::complex<double> >, Bembel::HelmholtzSingleLayerOperator>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::RotatedTangentialTrace<std::complex<double> >, Bembel::MaxwellSingleLayerOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::TangentialTrace<double>, TestOperatorDivC>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.

        10
            SurfacePoint qp;
      
        10
            auto super_space = ansatz_space_.get_superspace();
      
        10
            const ElementTree &element_tree = super_space.get_mesh().get_element_tree();
      
        10
            auto number_of_elements = element_tree.get_number_of_elements();
      
        10
            auto polynomial_degree = super_space.get_polynomial_degree();
      
        10
            auto polynomial_degree_plus_one_squared =
      
        10
                (polynomial_degree + 1) * (polynomial_degree + 1);
      
        10
            const auto function_space_dimension =
      
                getFunctionSpaceVectorDimension<LinearOperatorTraits<LinOp>::Form>();
      
            Eigen::Matrix<typename LinearFormTraits<Derived>::Scalar, Eigen::Dynamic,
      
                          function_space_dimension>
      
          7/14Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::HomogenisedLaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceBeltramiOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::MassMatrixScalarCont>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<std::complex<double> >, Bembel::HelmholtzSingleLayerOperator>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::RotatedTangentialTrace<std::complex<double> >, Bembel::MaxwellSingleLayerOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::TangentialTrace<double>, TestOperatorDivC>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.

        10
                intval(polynomial_degree_plus_one_squared, function_space_dimension);
      
            Eigen::Matrix<typename LinearFormTraits<Derived>::Scalar, Eigen::Dynamic,
      
                          function_space_dimension>
      
          7/14Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::HomogenisedLaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceBeltramiOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::MassMatrixScalarCont>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<std::complex<double> >, Bembel::HelmholtzSingleLayerOperator>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::RotatedTangentialTrace<std::complex<double> >, Bembel::MaxwellSingleLayerOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::TangentialTrace<double>, TestOperatorDivC>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.

        10
                disc_lf_matrix(polynomial_degree_plus_one_squared * number_of_elements,
      
                               function_space_dimension);
      
          7/14Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::HomogenisedLaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceBeltramiOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::MassMatrixScalarCont>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<std::complex<double> >, Bembel::HelmholtzSingleLayerOperator>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::RotatedTangentialTrace<std::complex<double> >, Bembel::MaxwellSingleLayerOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::TangentialTrace<double>, TestOperatorDivC>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.

        10
            disc_lf_matrix.setZero();
      
          14/14Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::HomogenisedLaplaceSingleLayerOperator>::compute():
✓ Branch 3 taken 120 times.
✓ Branch 4 taken 2 times.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceBeltramiOperator>::compute():
✓ Branch 3 taken 12 times.
✓ Branch 4 taken 2 times.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceSingleLayerOperator>::compute():
✓ Branch 3 taken 96 times.
✓ Branch 4 taken 1 time.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::MassMatrixScalarCont>::compute():
✓ Branch 3 taken 6 times.
✓ Branch 4 taken 1 time.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<std::complex<double> >, Bembel::HelmholtzSingleLayerOperator>::compute():
✓ Branch 3 taken 96 times.
✓ Branch 4 taken 1 time.
Bembel::DiscreteLinearForm<Bembel::RotatedTangentialTrace<std::complex<double> >, Bembel::MaxwellSingleLayerOperator>::compute():
✓ Branch 3 taken 102 times.
✓ Branch 4 taken 2 times.
Bembel::DiscreteLinearForm<Bembel::TangentialTrace<double>, TestOperatorDivC>::compute():
✓ Branch 3 taken 5 times.
✓ Branch 4 taken 1 time.

        447
            for (auto element = element_tree.cpbegin(); element != element_tree.cpend();
      
        437
                 ++element) {
      
          7/14Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::HomogenisedLaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 120 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceBeltramiOperator>::compute():
✓ Branch 1 taken 12 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 96 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::MassMatrixScalarCont>::compute():
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<std::complex<double> >, Bembel::HelmholtzSingleLayerOperator>::compute():
✓ Branch 1 taken 96 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::RotatedTangentialTrace<std::complex<double> >, Bembel::MaxwellSingleLayerOperator>::compute():
✓ Branch 1 taken 102 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::TangentialTrace<double>, TestOperatorDivC>::compute():
✓ Branch 1 taken 5 times.
✗ Branch 2 not taken.

        437
              intval.setZero();
      
          14/14Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::HomogenisedLaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 120 times.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceBeltramiOperator>::compute():
✓ Branch 1 taken 108 times.
✓ Branch 2 taken 12 times.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 384 times.
✓ Branch 2 taken 96 times.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::MassMatrixScalarCont>::compute():
✓ Branch 1 taken 54 times.
✓ Branch 2 taken 6 times.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<std::complex<double> >, Bembel::HelmholtzSingleLayerOperator>::compute():
✓ Branch 1 taken 384 times.
✓ Branch 2 taken 96 times.
Bembel::DiscreteLinearForm<Bembel::RotatedTangentialTrace<std::complex<double> >, Bembel::MaxwellSingleLayerOperator>::compute():
✓ Branch 1 taken 918 times.
✓ Branch 2 taken 102 times.
Bembel::DiscreteLinearForm<Bembel::TangentialTrace<double>, TestOperatorDivC>::compute():
✓ Branch 1 taken 45 times.
✓ Branch 2 taken 5 times.

        2810
              for (auto i = 0; i < Q.w_.size(); ++i) {
      
          14/28Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::HomogenisedLaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 480 times.
✗ Branch 2 not taken.
✓ Branch 5 taken 480 times.
✗ Branch 6 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceBeltramiOperator>::compute():
✓ Branch 1 taken 108 times.
✗ Branch 2 not taken.
✓ Branch 5 taken 108 times.
✗ Branch 6 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 384 times.
✗ Branch 2 not taken.
✓ Branch 5 taken 384 times.
✗ Branch 6 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::MassMatrixScalarCont>::compute():
✓ Branch 1 taken 54 times.
✗ Branch 2 not taken.
✓ Branch 5 taken 54 times.
✗ Branch 6 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<std::complex<double> >, Bembel::HelmholtzSingleLayerOperator>::compute():
✓ Branch 1 taken 384 times.
✗ Branch 2 not taken.
✓ Branch 5 taken 384 times.
✗ Branch 6 not taken.
Bembel::DiscreteLinearForm<Bembel::RotatedTangentialTrace<std::complex<double> >, Bembel::MaxwellSingleLayerOperator>::compute():
✓ Branch 1 taken 918 times.
✗ Branch 2 not taken.
✓ Branch 5 taken 918 times.
✗ Branch 6 not taken.
Bembel::DiscreteLinearForm<Bembel::TangentialTrace<double>, TestOperatorDivC>::compute():
✓ Branch 1 taken 45 times.
✗ Branch 2 not taken.
✓ Branch 5 taken 45 times.
✗ Branch 6 not taken.

        2373
                super_space.map2surface(*element, Q.xi_.col(i),
      
          14/28Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::HomogenisedLaplaceSingleLayerOperator>::compute():
✓ Branch 3 taken 480 times.
✗ Branch 4 not taken.
✓ Branch 6 taken 480 times.
✗ Branch 7 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceBeltramiOperator>::compute():
✓ Branch 3 taken 108 times.
✗ Branch 4 not taken.
✓ Branch 6 taken 108 times.
✗ Branch 7 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceSingleLayerOperator>::compute():
✓ Branch 3 taken 384 times.
✗ Branch 4 not taken.
✓ Branch 6 taken 384 times.
✗ Branch 7 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::MassMatrixScalarCont>::compute():
✓ Branch 3 taken 54 times.
✗ Branch 4 not taken.
✓ Branch 6 taken 54 times.
✗ Branch 7 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<std::complex<double> >, Bembel::HelmholtzSingleLayerOperator>::compute():
✓ Branch 3 taken 384 times.
✗ Branch 4 not taken.
✓ Branch 6 taken 384 times.
✗ Branch 7 not taken.
Bembel::DiscreteLinearForm<Bembel::RotatedTangentialTrace<std::complex<double> >, Bembel::MaxwellSingleLayerOperator>::compute():
✓ Branch 3 taken 918 times.
✗ Branch 4 not taken.
✓ Branch 6 taken 918 times.
✗ Branch 7 not taken.
Bembel::DiscreteLinearForm<Bembel::TangentialTrace<double>, TestOperatorDivC>::compute():
✓ Branch 3 taken 45 times.
✗ Branch 4 not taken.
✓ Branch 6 taken 45 times.
✗ Branch 7 not taken.

        2373
                                        element->get_h() * Q.w_(i), &qp);
      
          7/14Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::HomogenisedLaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 480 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceBeltramiOperator>::compute():
✓ Branch 1 taken 108 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 384 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::MassMatrixScalarCont>::compute():
✓ Branch 1 taken 54 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<std::complex<double> >, Bembel::HelmholtzSingleLayerOperator>::compute():
✓ Branch 1 taken 384 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::RotatedTangentialTrace<std::complex<double> >, Bembel::MaxwellSingleLayerOperator>::compute():
✓ Branch 1 taken 918 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::TangentialTrace<double>, TestOperatorDivC>::compute():
✓ Branch 1 taken 45 times.
✗ Branch 2 not taken.

        2373
                lf_.evaluateIntegrand_impl(super_space, qp, &intval);
      
              }
      
          7/14Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::HomogenisedLaplaceSingleLayerOperator>::compute():
✓ Branch 2 taken 120 times.
✗ Branch 3 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceBeltramiOperator>::compute():
✓ Branch 2 taken 12 times.
✗ Branch 3 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceSingleLayerOperator>::compute():
✓ Branch 2 taken 96 times.
✗ Branch 3 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::MassMatrixScalarCont>::compute():
✓ Branch 2 taken 6 times.
✗ Branch 3 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<std::complex<double> >, Bembel::HelmholtzSingleLayerOperator>::compute():
✓ Branch 2 taken 96 times.
✗ Branch 3 not taken.
Bembel::DiscreteLinearForm<Bembel::RotatedTangentialTrace<std::complex<double> >, Bembel::MaxwellSingleLayerOperator>::compute():
✓ Branch 2 taken 102 times.
✗ Branch 3 not taken.
Bembel::DiscreteLinearForm<Bembel::TangentialTrace<double>, TestOperatorDivC>::compute():
✓ Branch 2 taken 5 times.
✗ Branch 3 not taken.

        874
              disc_lf_matrix.block(polynomial_degree_plus_one_squared * element->id_, 0,
      
                                   polynomial_degree_plus_one_squared,
      
          7/14Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::HomogenisedLaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 120 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceBeltramiOperator>::compute():
✓ Branch 1 taken 12 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 96 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::MassMatrixScalarCont>::compute():
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<std::complex<double> >, Bembel::HelmholtzSingleLayerOperator>::compute():
✓ Branch 1 taken 96 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::RotatedTangentialTrace<std::complex<double> >, Bembel::MaxwellSingleLayerOperator>::compute():
✓ Branch 1 taken 102 times.
✗ Branch 2 not taken.
Bembel::DiscreteLinearForm<Bembel::TangentialTrace<double>, TestOperatorDivC>::compute():
✓ Branch 1 taken 5 times.
✗ Branch 2 not taken.

        437
                                   function_space_dimension) = intval;
      
            }
      
          14/28Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::HomogenisedLaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceBeltramiOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 time.
✗ Branch 5 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::MassMatrixScalarCont>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 time.
✗ Branch 5 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<std::complex<double> >, Bembel::HelmholtzSingleLayerOperator>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 time.
✗ Branch 5 not taken.
Bembel::DiscreteLinearForm<Bembel::RotatedTangentialTrace<std::complex<double> >, Bembel::MaxwellSingleLayerOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
Bembel::DiscreteLinearForm<Bembel::TangentialTrace<double>, TestOperatorDivC>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 time.
✗ Branch 5 not taken.

        10
            disc_lf_ =
      
          7/14Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::HomogenisedLaplaceSingleLayerOperator>::compute():
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceBeltramiOperator>::compute():
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceSingleLayerOperator>::compute():
✓ Branch 2 taken 1 time.
✗ Branch 3 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::MassMatrixScalarCont>::compute():
✓ Branch 2 taken 1 time.
✗ Branch 3 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<std::complex<double> >, Bembel::HelmholtzSingleLayerOperator>::compute():
✓ Branch 2 taken 1 time.
✗ Branch 3 not taken.
Bembel::DiscreteLinearForm<Bembel::RotatedTangentialTrace<std::complex<double> >, Bembel::MaxwellSingleLayerOperator>::compute():
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
Bembel::DiscreteLinearForm<Bembel::TangentialTrace<double>, TestOperatorDivC>::compute():
✓ Branch 2 taken 1 time.
✗ Branch 3 not taken.

        10
                ansatz_space_.get_transformation_matrix().transpose() *
      
          21/42Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::HomogenisedLaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 2 times.
✗ Branch 8 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceBeltramiOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 2 times.
✗ Branch 8 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::LaplaceSingleLayerOperator>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 time.
✗ Branch 5 not taken.
✓ Branch 7 taken 1 time.
✗ Branch 8 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<double>, Bembel::MassMatrixScalarCont>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 time.
✗ Branch 5 not taken.
✓ Branch 7 taken 1 time.
✗ Branch 8 not taken.
Bembel::DiscreteLinearForm<Bembel::DirichletTrace<std::complex<double> >, Bembel::HelmholtzSingleLayerOperator>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 time.
✗ Branch 5 not taken.
✓ Branch 7 taken 1 time.
✗ Branch 8 not taken.
Bembel::DiscreteLinearForm<Bembel::RotatedTangentialTrace<std::complex<double> >, Bembel::MaxwellSingleLayerOperator>::compute():
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 2 times.
✗ Branch 8 not taken.
Bembel::DiscreteLinearForm<Bembel::TangentialTrace<double>, TestOperatorDivC>::compute():
✓ Branch 1 taken 1 time.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 time.
✗ Branch 5 not taken.
✓ Branch 7 taken 1 time.
✗ Branch 8 not taken.

        20
                Eigen::Map<Eigen::Matrix<typename LinearOperatorTraits<LinOp>::Scalar,
      
                                         Eigen::Dynamic, 1>>(
      
                    disc_lf_matrix.data(),
      
        10
                    disc_lf_matrix.rows() * disc_lf_matrix.cols());
      
        20
            return;
      
        10
          }
      
          //////////////////////////////////////////////////////////////////////////////
      
          //    setter
      
          //////////////////////////////////////////////////////////////////////////////
      
          void set_degree(const int &deg) { deg_ = deg; }
      
          //////////////////////////////////////////////////////////////////////////////
      
          //    getter
      
          //////////////////////////////////////////////////////////////////////////////
      
        10
          Derived &get_linear_form() { return lf_; }
      
          const Eigen::Matrix<typename LinearFormTraits<Derived>::Scalar,
      
                              Eigen::Dynamic, 1>
      
        11
              &get_discrete_linear_form() const {
      
        11
            return disc_lf_;
      
          }
      
          //////////////////////////////////////////////////////////////////////////////
      
          //    private member variables
      
          //////////////////////////////////////////////////////////////////////////////
      
         private:
      
          int deg_;
      
          Derived lf_;
      
          Eigen::Matrix<typename LinearFormTraits<Derived>::Scalar, Eigen::Dynamic, 1>
      
              disc_lf_;
      
          AnsatzSpace<LinOp> ansatz_space_;
      
        };
      
        }  // namespace Bembel
      
        #endif  // BEMBEL_SRC_LINEARFORM_DISCRETELINEARFORM_HPP_