GCC Code Coverage Report

Directory:	../../../builds/dumux-repositories/
File:	/builds/dumux-repositories/dumux/dumux/multidomain/embedded/integrationpointsource.hh
Date:	2024-09-21 20:52:54

	Exec	Total	Coverage
Lines:	30	49	61.2%
Functions:	9	16	56.2%
Branches:	44	124	35.5%

  
      Line
      Branch
      Exec
      Source
    
      // -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
    
      // vi: set et ts=4 sw=4 sts=4:
    
      //
    
      // SPDX-FileCopyrightInfo: Copyright © DuMux Project contributors, see AUTHORS.md in root folder
    
      // SPDX-License-Identifier: GPL-3.0-or-later
    
      //
    
      /*!
    
       * \file
    
       * \ingroup EmbeddedCoupling
    
       * \brief An integration point source class,
    
       *        i.e. sources located at a single point in space
    
       *        associated with a quadrature point
    
       */
    
      #ifndef DUMUX_INTEGRATION_POINTSOURCE_HH
    
      #define DUMUX_INTEGRATION_POINTSOURCE_HH
    
      #include <type_traits>
    
      #include <dune/common/reservedvector.hh>
    
      #include <dumux/common/pointsource.hh>
    
      #include <dumux/common/parameters.hh>
    
      #include <dumux/discretization/method.hh>
    
      namespace Dumux {
    
      /*!
    
       * \ingroup EmbeddedCoupling
    
       * \brief An integration point source class with an identifier to attach data
    
       *        and a quadrature weight and integration element
    
       */
    
      template<class GlobalPosition, class SourceValues, typename IdType = std::size_t>
    
      class IntegrationPointSource : public IdPointSource<GlobalPosition, SourceValues, IdType>
    
      {
    
          using ParentType = IdPointSource<GlobalPosition, SourceValues, IdType>;
    
          using Scalar = std::decay_t<decltype(std::declval<SourceValues>()[0])>;
    
      public:
    
          //! Constructor for integration point sources
    
          IntegrationPointSource(GlobalPosition pos, SourceValues values, IdType id,
    
                                 Scalar qpweight, Scalar integrationElement,
    
                                 std::size_t elementIndex)
    
          : ParentType(pos, values, id)
    
          , qpweight_(qpweight)
    
          , integrationElement_(integrationElement)
    
          , elementIndex_(elementIndex)
    
          {}
    
          //! Constructor for integration point sources, when there is no
    
          //! value known at the time of initialization
    
      169609
          IntegrationPointSource(GlobalPosition pos, IdType id,
    
                                 Scalar qpweight, Scalar integrationElement,
    
                                 std::size_t elementIndex)
    
          : ParentType(pos, id)
    
          , qpweight_(qpweight)
    
          , integrationElement_(integrationElement)
    
      339218
          , elementIndex_(elementIndex)
    
          {}
    
      ✗
          Scalar quadratureWeight() const
    
      ✗
          { return qpweight_; }
    
      ✗
          Scalar integrationElement() const
    
      ✗
          { return integrationElement_; }
    
          void setQuadratureWeight(const Scalar qpWeight)
    
          { return qpweight_ = qpWeight; }
    
      ✗
          void setIntegrationElement(const Scalar ie)
    
      ✗
          { integrationElement_ = ie; }
    
          //! The index of the element this integration point source is associated with
    
      ✗
          std::size_t elementIndex() const
    
      ✗
          { return elementIndex_; }
    
          //! Convenience = operator overload modifying only the values
    
          IntegrationPointSource& operator= (const SourceValues& values)
    
          {
    
      7048064
              ParentType::operator=(values);
    
              return *this;
    
          }
    
          //! Convenience = operator overload modifying only the values
    
          IntegrationPointSource& operator= (Scalar s)
    
          {
    
      61177792
              ParentType::operator=(s);
    
              return *this;
    
          }
    
      private:
    
          Scalar qpweight_;
    
          Scalar integrationElement_;
    
          std::size_t elementIndex_;
    
      };
    
      /*!
    
       * \ingroup EmbeddedCoupling
    
       * \brief A helper class calculating a DOF-index to point source map
    
       */
    
      class IntegrationPointSourceHelper
    
      {
    
      public:
    
          //! calculate a DOF index to point source map from given vector of point sources
    
          template<class GridGeometry, class PointSource, class PointSourceMap>
    
      96
          static void computePointSourceMap(const GridGeometry& gridGeometry,
    
                                            const std::vector<PointSource>& sources,
    
                                            PointSourceMap& pointSourceMap,
    
                                            const std::string& paramGroup = "")
    
          {
    
        4/4✓ Branch 0 taken 169609 times.
✓ Branch 1 taken 48 times.
✓ Branch 2 taken 169609 times.
✓ Branch 3 taken 48 times.

      344010
              for (const auto& source : sources)
    
              {
    
                  // get the index of the element in which the point source falls
    
      339218
                  const auto eIdx = source.elementIndex();
    
                  if constexpr (GridGeometry::discMethod == DiscretizationMethods::box
    
                      || GridGeometry::discMethod == DiscretizationMethods::fcdiamond)
    
                  {
    
        1/4✓ Branch 1 taken 4494 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.

      13964
                      auto fvGeometry = localView(gridGeometry);
    
                      // check in which subcontrolvolume(s) we are
    
        4/8✓ Branch 1 taken 4494 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 4494 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 4494 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 2252 times.
✗ Branch 11 not taken.

      32884
                      const auto element = gridGeometry.boundingBoxTree().entitySet().entity(eIdx);
    
        1/2✓ Branch 1 taken 2252 times.
✗ Branch 2 not taken.

      9460
                      fvGeometry.bindElement(element);
    
      9460
                      const auto globalPos = source.position();
    
        4/6✓ Branch 0 taken 6 times.
✓ Branch 1 taken 4724 times.
✓ Branch 3 taken 6 times.
✗ Branch 4 not taken.
✓ Branch 6 taken 6 times.
✗ Branch 7 not taken.

      9460
                      static const bool boxPointSourceLumping = getParamFromGroup<bool>(paramGroup, "PointSource.EnableBoxLumping", true);
    
        1/2✓ Branch 0 taken 4730 times.
✗ Branch 1 not taken.

      9460
                      if (boxPointSourceLumping)
    
                      {
    
                          // loop over all sub control volumes and check if the point source is inside
    
      8988
                          constexpr int dim = GridGeometry::GridView::dimension;
    
        2/2✓ Branch 0 taken 236 times.
✓ Branch 1 taken 236 times.

      9932
                          Dune::ReservedVector<std::size_t, 1<<dim> scvIndices;
    
        4/4✓ Branch 0 taken 27416 times.
✓ Branch 1 taken 4730 times.
✓ Branch 2 taken 27416 times.
✓ Branch 3 taken 4730 times.

      73752
                          for (const auto& scv : scvs(fvGeometry))
    
        3/4✓ Branch 1 taken 27416 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 6979 times.
✓ Branch 5 taken 20437 times.

      54832
                              if (intersectsPointGeometry(globalPos, fvGeometry.geometry(scv)))
    
      27916
                                  scvIndices.push_back(scv.indexInElement());
    
                          // for all scvs that where tested positive add the point sources
    
                          // to the element/scv to point source map
    
        2/2✓ Branch 0 taken 6979 times.
✓ Branch 1 taken 4730 times.

      42338
                          for (auto scvIdx : scvIndices)
    
                          {
    
      13958
                              const auto key = std::make_pair(eIdx, scvIdx);
    
      13958
                              if (pointSourceMap.count(key))
    
        2/4✓ Branch 1 taken 2063 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2063 times.
✗ Branch 5 not taken.

      4126
                                  pointSourceMap.at(key).push_back(source);
    
                              else
    
        6/16✓ Branch 1 taken 4916 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 4916 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 4916 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 4916 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 4916 times.
✓ Branch 14 taken 4916 times.
✗ Branch 15 not taken.
✗ Branch 16 not taken.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✗ Branch 19 not taken.

      39328
                                  pointSourceMap.insert({key, {source}});
    
                              // split equally on the number of matched scvs
    
        1/2✓ Branch 1 taken 6979 times.
✗ Branch 2 not taken.

      13958
                              auto& s = pointSourceMap.at(key).back();
    
      27916
                              s.setEmbeddings(scvIndices.size()*s.embeddings());
    
                          }
    
                      }
    
                      // distribute the sources using the basis function weights
    
                      else
    
                      {
    
      ✗
                          const auto& localBasis = fvGeometry.feLocalBasis();
    
      ✗
                          const auto ipLocal = element.geometry().local(globalPos);
    
                          using Scalar = std::decay_t<decltype(source.values()[0])>;
    
      ✗
                          std::vector<typename Dune::FieldVector<Scalar, 1>> shapeValues;
    
      ✗
                          localBasis.evaluateFunction(ipLocal, shapeValues);
    
      ✗
                          for (const auto& scv : scvs(fvGeometry))
    
                          {
    
      ✗
                              const auto key = std::make_pair(eIdx, scv.indexInElement());
    
      ✗
                              if (pointSourceMap.count(key))
    
      ✗
                                  pointSourceMap.at(key).push_back(source);
    
                              else
    
      ✗
                                  pointSourceMap.insert({key, {source}});
    
                              // adjust the integration element
    
      ✗
                              auto& s = pointSourceMap.at(key).back();
    
      ✗
                              s.setIntegrationElement(shapeValues[scv.indexInElement()]*s.integrationElement());
    
                          }
    
                      }
    
                  }
    
                  else
    
                  {
    
                      // add the pointsource to the DOF map
    
      329758
                      const auto key = std::make_pair(eIdx, /*scvIdx=*/ 0);
    
      659516
                      if (pointSourceMap.count(key))
    
        3/6✓ Branch 1 taken 148858 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 148858 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 148858 times.
✗ Branch 8 not taken.

      595432
                          pointSourceMap.at(key).push_back(source);
    
                      else
    
        6/16✓ Branch 1 taken 16021 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 16021 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 16021 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 16021 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 16021 times.
✓ Branch 14 taken 16021 times.
✗ Branch 15 not taken.
✗ Branch 16 not taken.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✗ Branch 19 not taken.

      128168
                          pointSourceMap.insert({key, {source}});
    
                  }
    
              }
    
      96
          }
    
      };
    
      } // end namespace Dumux
    
      #endif

Line	Branch	Exec	Source
1			// -- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 --
2			// vi: set et ts=4 sw=4 sts=4:
3			//
4			// SPDX-FileCopyrightInfo: Copyright © DuMux Project contributors, see AUTHORS.md in root folder
5			// SPDX-License-Identifier: GPL-3.0-or-later
6			//
7			/*!
8			* \file
9			* \ingroup EmbeddedCoupling
10			* \brief An integration point source class,
11			* i.e. sources located at a single point in space
12			* associated with a quadrature point
13			*/
14
15			#ifndef DUMUX_INTEGRATION_POINTSOURCE_HH
16			#define DUMUX_INTEGRATION_POINTSOURCE_HH
17
18			#include <type_traits>
19			#include <dune/common/reservedvector.hh>
20			#include <dumux/common/pointsource.hh>
21			#include <dumux/common/parameters.hh>
22			#include <dumux/discretization/method.hh>
23
24			namespace Dumux {
25
26			/*!
27			* \ingroup EmbeddedCoupling
28			* \brief An integration point source class with an identifier to attach data
29			* and a quadrature weight and integration element
30			*/
31			template<class GlobalPosition, class SourceValues, typename IdType = std::size_t>
32			class IntegrationPointSource : public IdPointSource<GlobalPosition, SourceValues, IdType>
33			{
34			using ParentType = IdPointSource<GlobalPosition, SourceValues, IdType>;
35			using Scalar = std::decay_t<decltype(std::declval<SourceValues>()[0])>;
36
37			public:
38			//! Constructor for integration point sources
39			IntegrationPointSource(GlobalPosition pos, SourceValues values, IdType id,
40			Scalar qpweight, Scalar integrationElement,
41			std::size_t elementIndex)
42			: ParentType(pos, values, id)
43			, qpweight_(qpweight)
44			, integrationElement_(integrationElement)
45			, elementIndex_(elementIndex)
46			{}
47
48			//! Constructor for integration point sources, when there is no
49			//! value known at the time of initialization
50		169609	IntegrationPointSource(GlobalPosition pos, IdType id,
51			Scalar qpweight, Scalar integrationElement,
52			std::size_t elementIndex)
53			: ParentType(pos, id)
54			, qpweight_(qpweight)
55			, integrationElement_(integrationElement)
56		339218	, elementIndex_(elementIndex)
57			{}
58
59		✗	Scalar quadratureWeight() const
60		✗	{ return qpweight_; }
61
62		✗	Scalar integrationElement() const
63		✗	{ return integrationElement_; }
64
65			void setQuadratureWeight(const Scalar qpWeight)
66			{ return qpweight_ = qpWeight; }
67
68		✗	void setIntegrationElement(const Scalar ie)
69		✗	{ integrationElement_ = ie; }
70
71			//! The index of the element this integration point source is associated with
72		✗	std::size_t elementIndex() const
73		✗	{ return elementIndex_; }
74
75			//! Convenience = operator overload modifying only the values
76			IntegrationPointSource& operator= (const SourceValues& values)
77			{
78		7048064	ParentType::operator=(values);
79			return *this;
80			}
81
82			//! Convenience = operator overload modifying only the values
83			IntegrationPointSource& operator= (Scalar s)
84			{
85		61177792	ParentType::operator=(s);
86			return *this;
87			}
88
89			private:
90			Scalar qpweight_;
91			Scalar integrationElement_;
92			std::size_t elementIndex_;
93			};
94
95			/*!
96			* \ingroup EmbeddedCoupling
97			* \brief A helper class calculating a DOF-index to point source map
98			*/
99			class IntegrationPointSourceHelper
100			{
101
102			public:
103			//! calculate a DOF index to point source map from given vector of point sources
104			template<class GridGeometry, class PointSource, class PointSourceMap>
105		96	static void computePointSourceMap(const GridGeometry& gridGeometry,
106			const std::vector<PointSource>& sources,
107			PointSourceMap& pointSourceMap,
108			const std::string& paramGroup = "")
109			{
110	4/4 ✓ Branch 0 taken 169609 times. ✓ Branch 1 taken 48 times. ✓ Branch 2 taken 169609 times. ✓ Branch 3 taken 48 times.	344010	for (const auto& source : sources)
111			{
112			// get the index of the element in which the point source falls
113		339218	const auto eIdx = source.elementIndex();
114			if constexpr (GridGeometry::discMethod == DiscretizationMethods::box
115			\|\| GridGeometry::discMethod == DiscretizationMethods::fcdiamond)
116			{
117	1/4 ✓ Branch 1 taken 4494 times. ✗ Branch 2 not taken. ✗ Branch 3 not taken. ✗ Branch 4 not taken.	13964	auto fvGeometry = localView(gridGeometry);
118			// check in which subcontrolvolume(s) we are
119	4/8 ✓ Branch 1 taken 4494 times. ✗ Branch 2 not taken. ✓ Branch 4 taken 4494 times. ✗ Branch 5 not taken. ✓ Branch 7 taken 4494 times. ✗ Branch 8 not taken. ✓ Branch 10 taken 2252 times. ✗ Branch 11 not taken.	32884	const auto element = gridGeometry.boundingBoxTree().entitySet().entity(eIdx);
120	1/2 ✓ Branch 1 taken 2252 times. ✗ Branch 2 not taken.	9460	fvGeometry.bindElement(element);
121		9460	const auto globalPos = source.position();
122
123	4/6 ✓ Branch 0 taken 6 times. ✓ Branch 1 taken 4724 times. ✓ Branch 3 taken 6 times. ✗ Branch 4 not taken. ✓ Branch 6 taken 6 times. ✗ Branch 7 not taken.	9460	static const bool boxPointSourceLumping = getParamFromGroup<bool>(paramGroup, "PointSource.EnableBoxLumping", true);
124	1/2 ✓ Branch 0 taken 4730 times. ✗ Branch 1 not taken.	9460	if (boxPointSourceLumping)
125			{
126			// loop over all sub control volumes and check if the point source is inside
127		8988	constexpr int dim = GridGeometry::GridView::dimension;
128	2/2 ✓ Branch 0 taken 236 times. ✓ Branch 1 taken 236 times.	9932	Dune::ReservedVector<std::size_t, 1<<dim> scvIndices;
129	4/4 ✓ Branch 0 taken 27416 times. ✓ Branch 1 taken 4730 times. ✓ Branch 2 taken 27416 times. ✓ Branch 3 taken 4730 times.	73752	for (const auto& scv : scvs(fvGeometry))
130	3/4 ✓ Branch 1 taken 27416 times. ✗ Branch 2 not taken. ✓ Branch 4 taken 6979 times. ✓ Branch 5 taken 20437 times.	54832	if (intersectsPointGeometry(globalPos, fvGeometry.geometry(scv)))
131		27916	scvIndices.push_back(scv.indexInElement());
132
133			// for all scvs that where tested positive add the point sources
134			// to the element/scv to point source map
135	2/2 ✓ Branch 0 taken 6979 times. ✓ Branch 1 taken 4730 times.	42338	for (auto scvIdx : scvIndices)
136			{
137		13958	const auto key = std::make_pair(eIdx, scvIdx);
138		13958	if (pointSourceMap.count(key))
139	2/4 ✓ Branch 1 taken 2063 times. ✗ Branch 2 not taken. ✓ Branch 4 taken 2063 times. ✗ Branch 5 not taken.	4126	pointSourceMap.at(key).push_back(source);
140			else
141	6/16 ✓ Branch 1 taken 4916 times. ✗ Branch 2 not taken. ✓ Branch 4 taken 4916 times. ✗ Branch 5 not taken. ✓ Branch 7 taken 4916 times. ✗ Branch 8 not taken. ✓ Branch 10 taken 4916 times. ✗ Branch 11 not taken. ✗ Branch 12 not taken. ✓ Branch 13 taken 4916 times. ✓ Branch 14 taken 4916 times. ✗ Branch 15 not taken. ✗ Branch 16 not taken. ✗ Branch 17 not taken. ✗ Branch 18 not taken. ✗ Branch 19 not taken.	39328	pointSourceMap.insert({key, {source}});
142			// split equally on the number of matched scvs
143	1/2 ✓ Branch 1 taken 6979 times. ✗ Branch 2 not taken.	13958	auto& s = pointSourceMap.at(key).back();
144		27916	s.setEmbeddings(scvIndices.size()*s.embeddings());
145			}
146			}
147
148			// distribute the sources using the basis function weights
149			else
150			{
151		✗	const auto& localBasis = fvGeometry.feLocalBasis();
152		✗	const auto ipLocal = element.geometry().local(globalPos);
153			using Scalar = std::decay_t<decltype(source.values()[0])>;
154		✗	std::vector<typename Dune::FieldVector<Scalar, 1>> shapeValues;
155		✗	localBasis.evaluateFunction(ipLocal, shapeValues);
156		✗	for (const auto& scv : scvs(fvGeometry))
157			{
158		✗	const auto key = std::make_pair(eIdx, scv.indexInElement());
159		✗	if (pointSourceMap.count(key))
160		✗	pointSourceMap.at(key).push_back(source);
161			else
162		✗	pointSourceMap.insert({key, {source}});
163
164			// adjust the integration element
165		✗	auto& s = pointSourceMap.at(key).back();
166		✗	s.setIntegrationElement(shapeValues[scv.indexInElement()]*s.integrationElement());
167			}
168			}
169			}
170			else
171			{
172			// add the pointsource to the DOF map
173		329758	const auto key = std::make_pair(eIdx, /scvIdx=/ 0);
174		659516	if (pointSourceMap.count(key))
175	3/6 ✓ Branch 1 taken 148858 times. ✗ Branch 2 not taken. ✓ Branch 4 taken 148858 times. ✗ Branch 5 not taken. ✓ Branch 7 taken 148858 times. ✗ Branch 8 not taken.	595432	pointSourceMap.at(key).push_back(source);
176			else
177	6/16 ✓ Branch 1 taken 16021 times. ✗ Branch 2 not taken. ✓ Branch 4 taken 16021 times. ✗ Branch 5 not taken. ✓ Branch 7 taken 16021 times. ✗ Branch 8 not taken. ✓ Branch 10 taken 16021 times. ✗ Branch 11 not taken. ✗ Branch 12 not taken. ✓ Branch 13 taken 16021 times. ✓ Branch 14 taken 16021 times. ✗ Branch 15 not taken. ✗ Branch 16 not taken. ✗ Branch 17 not taken. ✗ Branch 18 not taken. ✗ Branch 19 not taken.	128168	pointSourceMap.insert({key, {source}});
178			}
179			}
180		96	}
181			};
182
183			} // end namespace Dumux
184
185			#endif
186