GCC Code Coverage Report

Directory: ../../../builds/dumux-repositories/
File: /builds/dumux-repositories/dumux/dumux/discretization/cellcentered/tpfa/elementvolumevariables.hh
Date: 2024-05-04 19:09:25
Exec Total Coverage
Lines: 102 111 91.9%
Functions: 500 771 64.9%
Branches: 464 620 74.8%
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 CCTpfaDiscretization
10 * \brief The local (stencil) volume variables class for cell centered tpfa models
11 */
12 #ifndef DUMUX_DISCRETIZATION_CCTPFA_ELEMENT_VOLUMEVARIABLES_HH
13 #define DUMUX_DISCRETIZATION_CCTPFA_ELEMENT_VOLUMEVARIABLES_HH
14
15 #include <algorithm>
16 #include <type_traits>
17 #include <vector>
18 #include <utility>
19
20 #include <dumux/discretization/cellcentered/elementsolution.hh>
21
22 namespace Dumux {
23
24 /*!
25 * \ingroup CCTpfaDiscretization
26 * \brief The local (stencil) volume variables class for cell centered tpfa models
27 * \note The class is specialized for versions with and without caching
28 * \tparam GVV the grid volume variables type
29 * \tparam cachingEnabled if the cache is enabled
30 */
31 template<class GVV, bool cachingEnabled>
32 class CCTpfaElementVolumeVariables
33 {};
34
35 /*!
36 * \ingroup CCTpfaDiscretization
37 * \brief The local (stencil) volume variables class for cell centered tpfa models with caching
38 * \note the volume variables are stored for the whole grid view in the corresponding GridVolumeVariables class
39 */
40 template<class GVV>
41
6/12
✓ Branch 1 taken 55 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 55 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 55 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 55 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 55 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 55 times.
✗ Branch 17 not taken.
 17056240 class CCTpfaElementVolumeVariables<GVV, /*cachingEnabled*/true>
42 {
43 public:
44 //! export type of the grid volume variables
45 using GridVolumeVariables = GVV;
46
47 //! export type of the volume variables
48 using VolumeVariables = typename GridVolumeVariables::VolumeVariables;
49
50 //! Constructor
51 CCTpfaElementVolumeVariables(const GridVolumeVariables& gridVolVars)
52 : gridVolVarsPtr_(&gridVolVars)
53
125/157
✓ Branch 1 taken 180486 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 180486 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 1233890 times.
✓ Branch 7 taken 180486 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 1233890 times.
✓ Branch 10 taken 180486 times.
✗ Branch 11 not taken.
✓ Branch 12 taken 1233890 times.
✓ Branch 13 taken 180486 times.
✗ Branch 14 not taken.
✓ Branch 15 taken 1233890 times.
✓ Branch 16 taken 9258850 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 1233890 times.
✓ Branch 19 taken 9258850 times.
✗ Branch 20 not taken.
✓ Branch 21 taken 4140413 times.
✓ Branch 22 taken 9258850 times.
✗ Branch 23 not taken.
✓ Branch 24 taken 4140413 times.
✓ Branch 25 taken 9258850 times.
✓ Branch 26 taken 3470 times.
✓ Branch 27 taken 4140413 times.
✓ Branch 28 taken 9258850 times.
✓ Branch 29 taken 3470 times.
✓ Branch 30 taken 4140413 times.
✓ Branch 31 taken 9697542 times.
✓ Branch 32 taken 3470 times.
✓ Branch 33 taken 4140413 times.
✓ Branch 34 taken 9697542 times.
✓ Branch 35 taken 3470 times.
✓ Branch 36 taken 2926319 times.
✓ Branch 37 taken 9697542 times.
✓ Branch 38 taken 3470 times.
✓ Branch 39 taken 2926319 times.
✓ Branch 40 taken 9697542 times.
✓ Branch 41 taken 3470 times.
✓ Branch 42 taken 2926319 times.
✓ Branch 43 taken 9697542 times.
✓ Branch 44 taken 3470 times.
✓ Branch 45 taken 2926319 times.
✓ Branch 46 taken 5266650 times.
✓ Branch 47 taken 3470 times.
✓ Branch 48 taken 2926319 times.
✓ Branch 49 taken 5266650 times.
✓ Branch 50 taken 3470 times.
✓ Branch 51 taken 19443991 times.
✓ Branch 52 taken 5266650 times.
✓ Branch 53 taken 3470 times.
✓ Branch 54 taken 19443991 times.
✓ Branch 55 taken 5266650 times.
✓ Branch 56 taken 173500 times.
✓ Branch 57 taken 19443991 times.
✓ Branch 58 taken 5266650 times.
✓ Branch 59 taken 173500 times.
✓ Branch 60 taken 19443991 times.
✓ Branch 61 taken 738259 times.
✓ Branch 62 taken 173500 times.
✓ Branch 63 taken 19443991 times.
✓ Branch 64 taken 738259 times.
✓ Branch 65 taken 173500 times.
✓ Branch 66 taken 19443891 times.
✓ Branch 67 taken 738259 times.
✓ Branch 68 taken 173500 times.
✓ Branch 69 taken 19443891 times.
✓ Branch 70 taken 738259 times.
✓ Branch 71 taken 173500 times.
✓ Branch 72 taken 19443891 times.
✓ Branch 73 taken 738259 times.
✓ Branch 74 taken 173500 times.
✓ Branch 75 taken 19443891 times.
✓ Branch 76 taken 320228 times.
✓ Branch 77 taken 173500 times.
✓ Branch 78 taken 19443891 times.
✓ Branch 79 taken 320228 times.
✓ Branch 80 taken 173500 times.
✓ Branch 81 taken 6717 times.
✓ Branch 82 taken 320228 times.
✓ Branch 83 taken 173500 times.
✓ Branch 84 taken 6717 times.
✓ Branch 85 taken 320228 times.
✓ Branch 86 taken 800000 times.
✓ Branch 87 taken 6717 times.
✓ Branch 88 taken 320228 times.
✓ Branch 89 taken 800000 times.
✓ Branch 90 taken 6717 times.
✓ Branch 91 taken 5000 times.
✓ Branch 92 taken 800000 times.
✓ Branch 93 taken 6717 times.
✓ Branch 94 taken 5000 times.
✓ Branch 95 taken 800000 times.
✓ Branch 96 taken 937786 times.
✓ Branch 97 taken 5000 times.
✓ Branch 98 taken 800000 times.
✓ Branch 99 taken 937786 times.
✓ Branch 100 taken 5000 times.
✓ Branch 101 taken 800000 times.
✓ Branch 102 taken 937786 times.
✓ Branch 103 taken 5000 times.
✓ Branch 104 taken 800000 times.
✓ Branch 105 taken 937786 times.
✗ Branch 106 not taken.
✓ Branch 107 taken 800000 times.
✓ Branch 108 taken 937786 times.
✗ Branch 109 not taken.
✓ Branch 110 taken 800000 times.
✓ Branch 111 taken 15642 times.
✗ Branch 112 not taken.
✓ Branch 113 taken 800000 times.
✓ Branch 114 taken 15642 times.
✗ Branch 115 not taken.
✓ Branch 116 taken 1744 times.
✓ Branch 117 taken 15642 times.
✗ Branch 118 not taken.
✓ Branch 119 taken 1744 times.
✓ Branch 120 taken 15642 times.
✗ Branch 121 not taken.
✓ Branch 122 taken 1744 times.
✓ Branch 123 taken 15642 times.
✗ Branch 124 not taken.
✓ Branch 125 taken 1744 times.
✗ Branch 126 not taken.
✓ Branch 128 taken 1744 times.
✗ Branch 129 not taken.
✓ Branch 131 taken 1735 times.
✗ Branch 132 not taken.
✓ Branch 134 taken 1735 times.
✗ Branch 135 not taken.
✓ Branch 137 taken 1735 times.
✗ Branch 138 not taken.
✓ Branch 140 taken 1735 times.
✗ Branch 141 not taken.
✓ Branch 143 taken 1735 times.
✗ Branch 144 not taken.
✓ Branch 146 taken 8000 times.
✗ Branch 147 not taken.
✓ Branch 149 taken 8000 times.
✗ Branch 150 not taken.
✓ Branch 152 taken 8000 times.
✗ Branch 153 not taken.
✓ Branch 155 taken 8000 times.
✗ Branch 156 not taken.
✓ Branch 158 taken 8000 times.
✗ Branch 159 not taken.
✓ Branch 166 taken 101 times.
✗ Branch 167 not taken.
✓ Branch 169 taken 101 times.
✗ Branch 170 not taken.
✓ Branch 172 taken 101 times.
✗ Branch 173 not taken.
✓ Branch 175 taken 101 times.
✗ Branch 176 not taken.
✓ Branch 178 taken 101 times.
✗ Branch 179 not taken.
 405249565 , numScv_(gridVolVars.problem().gridGeometry().numScv())
54 {}
55
56 //! operator for the access with an scv
57 template<class SubControlVolume, typename std::enable_if_t<!std::is_integral<SubControlVolume>::value, int> = 0>
58 2346454003 const VolumeVariables& operator [](const SubControlVolume& scv) const
59 {
60
2/4
✓ Branch 0 taken 1859875898 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1859875898 times.
✗ Branch 3 not taken.
 4692908006 if (scv.dofIndex() < numScv_)
61 7039362009 return gridVolVars().volVars(scv.dofIndex());
62 else
63 ✗ return boundaryVolumeVariables_[getLocalIdx_(scv.dofIndex())];
64 }
65
66 //! operator for the access with an index
67 3696549278 const VolumeVariables& operator [](const std::size_t scvIdx) const
68 {
69
2/2
✓ Branch 0 taken 2721574611 times.
✓ Branch 1 taken 1693516 times.
 3696549278 if (scvIdx < numScv_)
70 7389249204 return gridVolVars().volVars(scvIdx);
71 else
72 1924676 return boundaryVolumeVariables_[getLocalIdx_(scvIdx)];
73 }
74
75 /*!
76 * \brief bind the local view (r-value overload)
77 * This overload is called when an instance of this class is a temporary in the usage context
78 * This allows a usage like this: `const auto view = localView(...).bind(element);`
79 */
80 template<class FVElementGeometry, class SolutionVector>
81 9278289 CCTpfaElementVolumeVariables bind(const typename FVElementGeometry::GridGeometry::GridView::template Codim<0>::Entity& element,
82 const FVElementGeometry& fvGeometry,
83 const SolutionVector& sol) &&
84 {
85 9278289 this->bind_(element, fvGeometry, sol);
86 18556578 return std::move(*this);
87 }
88
89 template<class FVElementGeometry, class SolutionVector>
90 void bind(const typename FVElementGeometry::GridGeometry::GridView::template Codim<0>::Entity& element,
91 const FVElementGeometry& fvGeometry,
92 const SolutionVector& sol) &
93
4/11
✓ Branch 1 taken 1208068 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 55 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2500 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 13 taken 1735 times.
✗ Branch 14 not taken.
 40736398 { this->bind_(element, fvGeometry, sol); }
94
95 /*!
96 * \brief bind the local view (r-value overload)
97 * This overload is called when an instance of this class is a temporary in the usage context
98 * This allows a usage like this: `const auto view = localView(...).bind(element);`
99 */
100 template<class FVElementGeometry, class SolutionVector>
101 ✗ CCTpfaElementVolumeVariables bindElement(const typename FVElementGeometry::GridGeometry::GridView::template Codim<0>::Entity& element,
102 const FVElementGeometry& fvGeometry,
103 const SolutionVector& sol) &&
104 {
105 966862 this->bindElement_(element, fvGeometry, sol);
106 1933724 return std::move(*this);
107 }
108
109 template<class FVElementGeometry, class SolutionVector>
110 ✗ void bindElement(const typename FVElementGeometry::GridGeometry::GridView::template Codim<0>::Entity& element,
111 const FVElementGeometry& fvGeometry,
112 const SolutionVector& sol) &
113
1/2
✓ Branch 1 taken 1205120 times.
✗ Branch 2 not taken.
 8635723 { this->bindElement_(element, fvGeometry, sol); }
114
115 //! The global volume variables object we are a restriction of
116 ✗ const GridVolumeVariables& gridVolVars() const
117 ✗ { return *gridVolVarsPtr_; }
118
119 private:
120
121 //! Clear all local storage
122 void clear_()
123 {
124 15217146 boundaryVolVarIndices_.clear();
125
5/6
✓ Branch 0 taken 7793 times.
✓ Branch 1 taken 1891023 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 5709253 times.
✓ Branch 4 taken 125 times.
✓ Branch 5 taken 379 times.
 7608573 boundaryVolumeVariables_.clear();
126 }
127
128 //! precompute all boundary volume variables in a stencil of an element, the remaining ones are cached
129 template<class FVElementGeometry, class SolutionVector>
130 74626208 void bind_(const typename FVElementGeometry::GridGeometry::GridView::template Codim<0>::Entity& element,
131 const FVElementGeometry& fvGeometry,
132 const SolutionVector& sol)
133 {
134
4/4
✓ Branch 0 taken 7608573 times.
✓ Branch 1 taken 40176999 times.
✓ Branch 2 taken 7608573 times.
✓ Branch 3 taken 40176999 times.
 149252416 if (!fvGeometry.hasBoundaryScvf())
135 return;
136
137
2/2
✓ Branch 0 taken 7918 times.
✓ Branch 1 taken 7600655 times.
 13630892 clear_();
138 27261784 boundaryVolVarIndices_.reserve(fvGeometry.numScvf());
139 27261784 boundaryVolumeVariables_.reserve(fvGeometry.numScvf());
140
141
6/6
✓ Branch 0 taken 41620629 times.
✓ Branch 1 taken 7608573 times.
✓ Branch 2 taken 41620629 times.
✓ Branch 3 taken 7608573 times.
✓ Branch 4 taken 33330906 times.
✓ Branch 5 taken 8289723 times.
 104243370 for (const auto& scvf : scvfs(fvGeometry))
142 {
143
2/2
✓ Branch 0 taken 33330906 times.
✓ Branch 1 taken 8289723 times.
 76981586 if (!scvf.boundary())
144 62048022 continue;
145
146 // check if boundary is a pure dirichlet boundary
147 14933564 const auto& problem = gridVolVars().problem();
148
2/2
✓ Branch 0 taken 183835 times.
✓ Branch 1 taken 6712958 times.
 14933564 const auto bcTypes = problem.boundaryTypes(element, scvf);
149
4/4
✓ Branch 0 taken 255790 times.
✓ Branch 1 taken 8033933 times.
✓ Branch 2 taken 255790 times.
✓ Branch 3 taken 8033933 times.
 29867128 if (bcTypes.hasOnlyDirichlet())
150 {
151
2/2
✓ Branch 0 taken 62349 times.
✓ Branch 1 taken 1874 times.
 462302 const auto dirichletPriVars = elementSolution<FVElementGeometry>(problem.dirichlet(element, scvf));
152 743346 auto&& scvI = fvGeometry.scv(scvf.insideScvIdx());
153
154 371673 VolumeVariables volVars;
155 371673 volVars.update(dirichletPriVars,
156 problem,
157 element,
158 scvI);
159
160 371673 boundaryVolumeVariables_.emplace_back(std::move(volVars));
161 743346 boundaryVolVarIndices_.push_back(scvf.outsideScvIdx());
162 }
163 }
164 }
165
166 //! precompute the volume variables of an element - do nothing: volVars are cached
167 template<class FVElementGeometry, class SolutionVector>
168 ✗ void bindElement_(const typename FVElementGeometry::GridGeometry::GridView::template Codim<0>::Entity& element,
169 const FVElementGeometry& fvGeometry,
170 const SolutionVector& sol)
171 ✗ {}
172
173 const GridVolumeVariables* gridVolVarsPtr_;
174
175 //! map a global scv index to the local storage index
176 1693516 int getLocalIdx_(const int volVarIdx) const
177 {
178 5080548 auto it = std::find(boundaryVolVarIndices_.begin(), boundaryVolVarIndices_.end(), volVarIdx);
179
3/6
✗ Branch 0 not taken.
✓ Branch 1 taken 1462356 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 1462356 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 1462356 times.
 5080548 assert(it != boundaryVolVarIndices_.end() && "Could not find the current volume variables for volVarIdx!");
180 5080548 return std::distance(boundaryVolVarIndices_.begin(), it);
181 }
182
183 std::vector<std::size_t> boundaryVolVarIndices_;
184 std::vector<VolumeVariables> boundaryVolumeVariables_;
185 const std::size_t numScv_;
186 };
187
188 /*!
189 * \ingroup CCTpfaDiscretization
190 * \brief The local (stencil) volume variables class for cell centered tpfa models with caching
191 */
192 template<class GVV>
193
6/12
✓ Branch 1 taken 4 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 4 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 4 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 4 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 4 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 4 times.
✗ Branch 17 not taken.
 348510 class CCTpfaElementVolumeVariables<GVV, /*cachingEnabled*/false>
194 {
195 public:
196 //! export type of the grid volume variables
197 using GridVolumeVariables = GVV;
198
199 //! export type of the volume variables
200 using VolumeVariables = typename GridVolumeVariables::VolumeVariables;
201
202 //! Constructor
203 CCTpfaElementVolumeVariables(const GridVolumeVariables& gridVolVars)
204
44/88
✓ Branch 1 taken 1789 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 41440 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 41740 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 19829793 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 19834051 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 19834315 times.
✗ Branch 17 not taken.
✓ Branch 19 taken 20211330 times.
✗ Branch 20 not taken.
✓ Branch 22 taken 20171668 times.
✗ Branch 23 not taken.
✓ Branch 25 taken 20221416 times.
✗ Branch 26 not taken.
✓ Branch 28 taken 423678 times.
✗ Branch 29 not taken.
✓ Branch 31 taken 423678 times.
✗ Branch 32 not taken.
✓ Branch 34 taken 423678 times.
✗ Branch 35 not taken.
✓ Branch 37 taken 254084 times.
✗ Branch 38 not taken.
✓ Branch 40 taken 519623 times.
✗ Branch 41 not taken.
✓ Branch 43 taken 469623 times.
✗ Branch 44 not taken.
✓ Branch 46 taken 469500 times.
✗ Branch 47 not taken.
✓ Branch 49 taken 469500 times.
✗ Branch 50 not taken.
✓ Branch 52 taken 475236 times.
✗ Branch 53 not taken.
✓ Branch 55 taken 275884 times.
✗ Branch 56 not taken.
✓ Branch 58 taken 6098 times.
✗ Branch 59 not taken.
✓ Branch 61 taken 6098 times.
✗ Branch 62 not taken.
✓ Branch 64 taken 6048 times.
✗ Branch 65 not taken.
✓ Branch 67 taken 6048 times.
✗ Branch 68 not taken.
✓ Branch 70 taken 55781 times.
✗ Branch 71 not taken.
✓ Branch 73 taken 55781 times.
✗ Branch 74 not taken.
✓ Branch 76 taken 55781 times.
✗ Branch 77 not taken.
✓ Branch 79 taken 55781 times.
✗ Branch 80 not taken.
✓ Branch 82 taken 55781 times.
✗ Branch 83 not taken.
✓ Branch 85 taken 55781 times.
✗ Branch 86 not taken.
✓ Branch 88 taken 50000 times.
✗ Branch 89 not taken.
✓ Branch 91 taken 50000 times.
✗ Branch 92 not taken.
✓ Branch 94 taken 50000 times.
✗ Branch 95 not taken.
✓ Branch 97 taken 964000 times.
✗ Branch 98 not taken.
✓ Branch 100 taken 964000 times.
✗ Branch 101 not taken.
✓ Branch 103 taken 964000 times.
✗ Branch 104 not taken.
✓ Branch 106 taken 964000 times.
✗ Branch 107 not taken.
✓ Branch 109 taken 964000 times.
✗ Branch 110 not taken.
✓ Branch 112 taken 964000 times.
✗ Branch 113 not taken.
✓ Branch 115 taken 100 times.
✗ Branch 116 not taken.
✓ Branch 118 taken 100 times.
✗ Branch 119 not taken.
✓ Branch 121 taken 100 times.
✗ Branch 122 not taken.
✓ Branch 127 taken 1 times.
✗ Branch 128 not taken.
✓ Branch 130 taken 1 times.
✗ Branch 131 not taken.
✓ Branch 133 taken 1 times.
✗ Branch 134 not taken.
 130685601 : gridVolVarsPtr_(&gridVolVars) {}
205
206 /*!
207 * \brief bind the local view (r-value overload)
208 * This overload is called when an instance of this class is a temporary in the usage context
209 * This allows a usage like this: `const auto view = localView(...).bind(element);`
210 */
211 template<class FVElementGeometry, class SolutionVector>
212 105599 CCTpfaElementVolumeVariables bind(const typename FVElementGeometry::GridGeometry::GridView::template Codim<0>::Entity& element,
213 const FVElementGeometry& fvGeometry,
214 const SolutionVector& sol) &&
215 {
216 105599 this->bind_(element, fvGeometry, sol);
217 211198 return std::move(*this);
218 }
219
220 template<class FVElementGeometry, class SolutionVector>
221 void bind(const typename FVElementGeometry::GridGeometry::GridView::template Codim<0>::Entity& element,
222 const FVElementGeometry& fvGeometry,
223 const SolutionVector& sol) &
224
6/14
✓ Branch 1 taken 549086 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 87340 times.
✓ Branch 5 taken 36914 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 76900 times.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✓ Branch 17 taken 100 times.
✗ Branch 18 not taken.
✓ Branch 20 taken 1928 times.
✗ Branch 21 not taken.
 22539684 { this->bind_(element, fvGeometry, sol); }
225
226 /*!
227 * \brief bind the local view (r-value overload)
228 * This overload is called when an instance of this class is a temporary in the usage context
229 * This allows a usage like this: `const auto view = localView(...).bind(element);`
230 */
231 template<class FVElementGeometry, class SolutionVector>
232 6316 CCTpfaElementVolumeVariables bindElement(const typename FVElementGeometry::GridGeometry::GridView::template Codim<0>::Entity& element,
233 const FVElementGeometry& fvGeometry,
234 const SolutionVector& sol) &&
235 {
236 6316 this->bindElement_(element, fvGeometry, sol);
237 12632 return std::move(*this);
238 }
239
240 template<class FVElementGeometry, class SolutionVector>
241 void bindElement(const typename FVElementGeometry::GridGeometry::GridView::template Codim<0>::Entity& element,
242 const FVElementGeometry& fvGeometry,
243 const SolutionVector& sol) &
244
7/10
✓ Branch 1 taken 991402 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 26045 times.
✓ Branch 5 taken 3501319 times.
✓ Branch 6 taken 5 times.
✓ Branch 7 taken 21216 times.
✓ Branch 8 taken 7756 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 1100 times.
✗ Branch 11 not taken.
 25501185 { this->bindElement_(element, fvGeometry, sol); }
245
246 //! access operator with scv
247 template<class SubControlVolume, typename std::enable_if_t<!std::is_integral<SubControlVolume>::value, int> = 0>
248 const VolumeVariables& operator [](const SubControlVolume& scv) const
249
28/36
✓ Branch 2 taken 1 times.
✓ Branch 3 taken 4079 times.
✓ Branch 4 taken 1 times.
✓ Branch 5 taken 4079 times.
✓ Branch 6 taken 1512 times.
✓ Branch 7 taken 10584 times.
✓ Branch 8 taken 1512 times.
✓ Branch 9 taken 10584 times.
✓ Branch 16 taken 18192 times.
✓ Branch 17 taken 98479 times.
✓ Branch 18 taken 18192 times.
✓ Branch 19 taken 98479 times.
✓ Branch 20 taken 3037 times.
✓ Branch 21 taken 557794 times.
✓ Branch 22 taken 84445 times.
✓ Branch 23 taken 557794 times.
✓ Branch 24 taken 8256916 times.
✓ Branch 25 taken 47800 times.
✓ Branch 26 taken 8175508 times.
✓ Branch 27 taken 47800 times.
✓ Branch 28 taken 76800 times.
✗ Branch 29 not taken.
✓ Branch 30 taken 76800 times.
✗ Branch 31 not taken.
✓ Branch 40 taken 150000 times.
✗ Branch 41 not taken.
✓ Branch 42 taken 150000 times.
✗ Branch 43 not taken.
✓ Branch 70 taken 150000 times.
✗ Branch 71 not taken.
✓ Branch 72 taken 150000 times.
✗ Branch 73 not taken.
✓ Branch 80 taken 2892000 times.
✗ Branch 81 not taken.
✓ Branch 82 taken 2892000 times.
✗ Branch 83 not taken.
 709414966 { return volumeVariables_[getLocalIdx_(scv.dofIndex())]; }
250
251 //! access operator with scv
252 template<class SubControlVolume, typename std::enable_if_t<!std::is_integral<SubControlVolume>::value, int> = 0>
253 VolumeVariables& operator [](const SubControlVolume& scv)
254
10/16
✓ Branch 5 taken 19456 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 1674817 times.
✓ Branch 8 taken 19456 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 1674817 times.
✓ Branch 11 taken 33000 times.
✗ Branch 12 not taken.
✓ Branch 13 taken 302697 times.
✓ Branch 14 taken 33000 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 302697 times.
✓ Branch 17 taken 1257844 times.
✗ Branch 18 not taken.
✓ Branch 20 taken 1257844 times.
✗ Branch 21 not taken.
 44607594 { return volumeVariables_[getLocalIdx_(scv.dofIndex())]; }
255
256 //! access operator with scv index
257 const VolumeVariables& operator [](std::size_t scvIdx) const
258
135/147
✓ Branch 1 taken 48 times.
✓ Branch 2 taken 9984 times.
✓ Branch 3 taken 48 times.
✓ Branch 4 taken 11391531 times.
✓ Branch 5 taken 19462028 times.
✓ Branch 6 taken 11381547 times.
✓ Branch 7 taken 19783362 times.
✓ Branch 8 taken 304538 times.
✓ Branch 9 taken 5064740 times.
✓ Branch 10 taken 41469676 times.
✓ Branch 11 taken 8728294 times.
✓ Branch 12 taken 43648826 times.
✓ Branch 13 taken 7819022 times.
✓ Branch 14 taken 3866456 times.
✓ Branch 15 taken 37912890 times.
✓ Branch 16 taken 1395376 times.
✓ Branch 17 taken 34814085 times.
✓ Branch 18 taken 30008195 times.
✓ Branch 19 taken 40522394 times.
✓ Branch 20 taken 33201591 times.
✓ Branch 21 taken 44817333 times.
✓ Branch 22 taken 5680491 times.
✓ Branch 23 taken 23828320 times.
✓ Branch 24 taken 15040349 times.
✓ Branch 25 taken 28282380 times.
✓ Branch 26 taken 36009579 times.
✓ Branch 27 taken 25563528 times.
✓ Branch 28 taken 25457939 times.
✓ Branch 29 taken 17812505 times.
✓ Branch 30 taken 44178229 times.
✓ Branch 31 taken 35700693 times.
✓ Branch 32 taken 80222115 times.
✓ Branch 33 taken 35854360 times.
✓ Branch 34 taken 41218095 times.
✓ Branch 35 taken 6800741 times.
✓ Branch 36 taken 9741568 times.
✓ Branch 37 taken 7390697 times.
✓ Branch 38 taken 6589123 times.
✓ Branch 39 taken 10612604 times.
✓ Branch 40 taken 53253374 times.
✓ Branch 41 taken 52141354 times.
✓ Branch 42 taken 54740798 times.
✓ Branch 43 taken 44268156 times.
✓ Branch 44 taken 27963249 times.
✓ Branch 45 taken 56354947 times.
✓ Branch 46 taken 27207015 times.
✓ Branch 47 taken 61119633 times.
✓ Branch 48 taken 2929803 times.
✓ Branch 49 taken 7807607 times.
✓ Branch 50 taken 3418058 times.
✓ Branch 51 taken 2781432 times.
✓ Branch 52 taken 1228206 times.
✓ Branch 53 taken 1368697 times.
✓ Branch 54 taken 89519778 times.
✓ Branch 55 taken 94032705 times.
✓ Branch 56 taken 89579931 times.
✓ Branch 57 taken 93322493 times.
✓ Branch 58 taken 726896 times.
✓ Branch 59 taken 87243377 times.
✓ Branch 60 taken 3579923 times.
✓ Branch 61 taken 103158074 times.
✓ Branch 62 taken 3020555 times.
✓ Branch 63 taken 17450936 times.
✓ Branch 64 taken 106228 times.
✓ Branch 65 taken 7465449 times.
✓ Branch 66 taken 5596 times.
✓ Branch 67 taken 6706168 times.
✓ Branch 68 taken 6646 times.
✓ Branch 69 taken 119462 times.
✓ Branch 70 taken 1122 times.
✗ Branch 71 not taken.
✗ Branch 72 not taken.
✓ Branch 73 taken 8136 times.
✓ Branch 74 taken 21091 times.
✓ Branch 75 taken 53626 times.
✓ Branch 76 taken 34387 times.
✓ Branch 77 taken 57250 times.
✓ Branch 78 taken 13297 times.
✓ Branch 79 taken 27495 times.
✓ Branch 80 taken 71814 times.
✓ Branch 81 taken 16815 times.
✓ Branch 82 taken 73325 times.
✓ Branch 83 taken 1080 times.
✓ Branch 84 taken 1518 times.
✓ Branch 85 taken 6 times.
✓ Branch 86 taken 27978 times.
✓ Branch 87 taken 27978 times.
✓ Branch 88 taken 27972 times.
✓ Branch 89 taken 34308 times.
✓ Branch 90 taken 121248 times.
✓ Branch 91 taken 1325346 times.
✓ Branch 92 taken 1366224 times.
✓ Branch 93 taken 1319010 times.
✓ Branch 94 taken 1279660 times.
✓ Branch 95 taken 531012 times.
✓ Branch 96 taken 830556 times.
✓ Branch 97 taken 567300 times.
✓ Branch 98 taken 2007200 times.
✓ Branch 99 taken 499720 times.
✓ Branch 100 taken 2020844 times.
✓ Branch 101 taken 491404 times.
✓ Branch 102 taken 809716 times.
✓ Branch 103 taken 538067 times.
✓ Branch 104 taken 706964 times.
✓ Branch 105 taken 510095 times.
✓ Branch 106 taken 706764 times.
✓ Branch 107 taken 2763 times.
✓ Branch 108 taken 39056 times.
✓ Branch 109 taken 1214091 times.
✓ Branch 110 taken 39056 times.
✓ Branch 111 taken 1211328 times.
✓ Branch 113 taken 36288 times.
✓ Branch 114 taken 1211328 times.
✓ Branch 115 taken 36288 times.
✓ Branch 116 taken 1211328 times.
✗ Branch 117 not taken.
✓ Branch 118 taken 502725 times.
✓ Branch 119 taken 708603 times.
✓ Branch 120 taken 502725 times.
✓ Branch 121 taken 708603 times.
✗ Branch 123 not taken.
✓ Branch 124 taken 1211328 times.
✗ Branch 125 not taken.
✓ Branch 126 taken 1211328 times.
✗ Branch 127 not taken.
✗ Branch 128 not taken.
✗ Branch 129 not taken.
✗ Branch 134 not taken.
✓ Branch 135 taken 11232 times.
✗ Branch 136 not taken.
✓ Branch 137 taken 11232 times.
✗ Branch 142 not taken.
✓ Branch 143 taken 11232 times.
✗ Branch 144 not taken.
✓ Branch 145 taken 11232 times.
✓ Branch 152 taken 72576 times.
✓ Branch 153 taken 2422656 times.
✓ Branch 154 taken 72576 times.
✓ Branch 155 taken 2422656 times.
✓ Branch 160 taken 1081704 times.
✓ Branch 161 taken 3836184 times.
✓ Branch 162 taken 1081704 times.
✓ Branch 163 taken 3836184 times.
✓ Branch 165 taken 1005450 times.
✓ Branch 166 taken 1417206 times.
✓ Branch 167 taken 1005450 times.
✓ Branch 168 taken 1417206 times.
 1613175953 { return volumeVariables_[getLocalIdx_(scvIdx)]; }
259
260 //! access operator with scv index
261 VolumeVariables& operator [](std::size_t scvIdx)
262
4/8
✗ Branch 1 not taken.
✓ Branch 2 taken 34460 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 34460 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 1056 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 1056 times.
 1574476 { return volumeVariables_[getLocalIdx_(scvIdx)]; }
263
264 //! The global volume variables object we are a restriction of
265 ✗ const GridVolumeVariables& gridVolVars() const
266 ✗ { return *gridVolVarsPtr_; }
267
268 private:
269 //! Clear all local storage
270 void clear_()
271 {
272 97313448 volVarIndices_.clear();
273
15/16
✓ Branch 0 taken 1051237 times.
✓ Branch 1 taken 19547088 times.
✓ Branch 2 taken 4596280 times.
✓ Branch 3 taken 20097015 times.
✓ Branch 4 taken 157168 times.
✓ Branch 5 taken 432734 times.
✓ Branch 6 taken 76849 times.
✓ Branch 7 taken 680734 times.
✓ Branch 8 taken 25049 times.
✓ Branch 9 taken 964013 times.
✗ Branch 10 not taken.
✓ Branch 11 taken 1020045 times.
✓ Branch 12 taken 99 times.
✓ Branch 13 taken 301 times.
✓ Branch 14 taken 1927 times.
✓ Branch 15 taken 6185 times.
 48656724 volumeVariables_.clear();
274 }
275
276 //! Prepares the volume variables within the element stencil
277 template<class FVElementGeometry, class SolutionVector>
278 24999676 void bind_(const typename FVElementGeometry::GridGeometry::GridView::template Codim<0>::Entity& element,
279 const FVElementGeometry& fvGeometry,
280 const SolutionVector& sol)
281 {
282
2/2
✓ Branch 0 taken 1152091 times.
✓ Branch 1 taken 21812020 times.
 24999676 clear_();
283
284 24999676 const auto& problem = gridVolVars().problem();
285 24999676 const auto& gridGeometry = fvGeometry.gridGeometry();
286 49999352 const auto globalI = gridGeometry.elementMapper().index(element);
287 49999352 const auto& connectivityMapI = gridGeometry.connectivityMap()[globalI];
288 24999676 const auto numDofs = connectivityMapI.size() + 1;
289
290 // resize local containers to the required size (for internal elements)
291 24999676 volumeVariables_.resize(numDofs);
292 24999676 volVarIndices_.resize(numDofs);
293 24999676 int localIdx = 0;
294
295 // update the volume variables of the element at hand
296
1/2
✓ Branch 0 taken 16592721 times.
✗ Branch 1 not taken.
 24999676 auto&& scvI = fvGeometry.scv(globalI);
297 25005724 volumeVariables_[localIdx].update(elementSolution(element, sol, gridGeometry),
298 problem,
299 element,
300 scvI);
301 24999676 volVarIndices_[localIdx] = scvI.dofIndex();
302 24999676 ++localIdx;
303
304 // Update the volume variables of the neighboring elements
305
4/4
✓ Branch 0 taken 69500656 times.
✓ Branch 1 taken 22964111 times.
✓ Branch 2 taken 69500656 times.
✓ Branch 3 taken 22964111 times.
 227550784 for (const auto& dataJ : connectivityMapI)
306 {
307
1/3
✗ Branch 0 not taken.
✓ Branch 1 taken 41901486 times.
✗ Branch 2 not taken.
 91271034 const auto& elementJ = gridGeometry.element(dataJ.globalJ);
308
1/3
✗ Branch 0 not taken.
✓ Branch 1 taken 45367402 times.
✗ Branch 2 not taken.
 76275878 auto&& scvJ = fvGeometry.scv(dataJ.globalJ);
309
4/7
✓ Branch 1 taken 12099118 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 29557446 times.
✓ Branch 4 taken 12099118 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 11414638 times.
✗ Branch 8 not taken.
 152564032 volumeVariables_[localIdx].update(elementSolution(elementJ, sol, gridGeometry),
310 problem,
311 elementJ,
312 scvJ);
313 152551756 volVarIndices_[localIdx] = scvJ.dofIndex();
314 76275878 ++localIdx;
315 }
316
317
4/4
✓ Branch 0 taken 3272311 times.
✓ Branch 1 taken 19691800 times.
✓ Branch 2 taken 744708 times.
✓ Branch 3 taken 5626682 times.
 31371066 if (fvGeometry.hasBoundaryScvf())
318 {
319 // Update boundary volume variables
320
6/6
✓ Branch 0 taken 13550617 times.
✓ Branch 1 taken 3272311 times.
✓ Branch 2 taken 13550617 times.
✓ Branch 3 taken 3272311 times.
✓ Branch 4 taken 2198428 times.
✓ Branch 5 taken 756676 times.
 20908843 for (auto&& scvf : scvfs(fvGeometry))
321 {
322 // if we are not on a boundary, skip to the next scvf
323
2/2
✓ Branch 0 taken 8382632 times.
✓ Branch 1 taken 5167985 times.
 14075453 if (!scvf.boundary())
324 8751516 continue;
325
326 // check if boundary is a pure dirichlet boundary
327
2/2
✓ Branch 0 taken 2937257 times.
✓ Branch 1 taken 75431 times.
 5323937 const auto bcTypes = problem.boundaryTypes(element, scvf);
328
4/4
✓ Branch 0 taken 607260 times.
✓ Branch 1 taken 4560725 times.
✓ Branch 2 taken 607260 times.
✓ Branch 3 taken 4560725 times.
 10647874 if (bcTypes.hasOnlyDirichlet())
329 {
330
2/2
✓ Branch 0 taken 6151 times.
✓ Branch 1 taken 21479 times.
 1170884 const auto dirichletPriVars = elementSolution<FVElementGeometry>(problem.dirichlet(element, scvf));
331
332 659913 volumeVariables_.resize(localIdx+1);
333 659913 volVarIndices_.resize(localIdx+1);
334 1319826 volumeVariables_[localIdx].update(dirichletPriVars,
335 problem,
336 element,
337 scvI);
338 1319826 volVarIndices_[localIdx] = scvf.outsideScvIdx();
339 659913 ++localIdx;
340 }
341 }
342 }
343
344 //! Check if user added additional DOF dependencies, i.e. the residual of DOF globalI depends
345 //! on additional DOFs not included in the discretization schemes' occupation pattern
346 // const auto& additionalDofDependencies = problem.getAdditionalDofDependencies(globalI);
347 // if (!additionalDofDependencies.empty())
348 // {
349 // volumeVariables_.resize(volumeVariables_.size() + additionalDofDependencies.size());
350 // volVarIndices_.resize(volVarIndices_.size() + additionalDofDependencies.size());
351 // for (auto globalJ : additionalDofDependencies)
352 // {
353 // const auto& elementJ = gridGeometry.element(globalJ);
354 // auto&& scvJ = fvGeometry.scv(globalJ);
355
356 // volumeVariables_[localIdx].update(elementSolution(elementJ, sol, gridGeometry),
357 // problem,
358 // elementJ,
359 // scvJ);
360 // volVarIndices_[localIdx] = scvJ.dofIndex();
361 // ++localIdx;
362 // }
363 // }
364 24999676 }
365
366 template<class FVElementGeometry, class SolutionVector>
367 27560377 void bindElement_(const typename FVElementGeometry::GridGeometry::GridView::template Codim<0>::Entity& element,
368 const FVElementGeometry& fvGeometry,
369 const SolutionVector& sol)
370 {
371
2/2
✓ Branch 0 taken 4756518 times.
✓ Branch 1 taken 20936095 times.
 27560377 clear_();
372
373 55120754 const auto eIdx = fvGeometry.gridGeometry().elementMapper().index(element);
374 27560377 volumeVariables_.resize(1);
375 27560377 volVarIndices_.resize(1);
376
377 // update the volume variables of the element
378
1/2
✓ Branch 0 taken 18225533 times.
✗ Branch 1 not taken.
 27560377 auto&& scv = fvGeometry.scv(eIdx);
379 49687957 volumeVariables_[0].update(elementSolution(element, sol, fvGeometry.gridGeometry()),
380 27560377 gridVolVars().problem(),
381 element,
382 scv);
383 55120744 volVarIndices_[0] = scv.dofIndex();
384 27560372 }
385
386 const GridVolumeVariables* gridVolVarsPtr_;
387
388 //! map a global scv index to the local storage index
389 2905531569 int getLocalIdx_(const int volVarIdx) const
390 {
391 8716594707 auto it = std::find(volVarIndices_.begin(), volVarIndices_.end(), volVarIdx);
392
3/6
✗ Branch 0 not taken.
✓ Branch 1 taken 2726729436 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 2726729436 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2726729436 times.
 8716594707 assert(it != volVarIndices_.end() && "Could not find the current volume variables for volVarIdx!");
393 8716594707 return std::distance(volVarIndices_.begin(), it);
394 }
395
396 std::vector<std::size_t> volVarIndices_;
397 std::vector<VolumeVariables> volumeVariables_;
398 };
399
400 } // end namespace Dumux
401
402 #endif
403