GCC Code Coverage Report

Directory: ../../../builds/dumux-repositories/
File: /builds/dumux-repositories/dumux/test/freeflow/navierstokes/analyticalsolutionvectors.hh
Date: 2024-09-21 20:52:54
Exec Total Coverage
Lines: 50 50 100.0%
Functions: 18 18 100.0%
Branches: 86 153 56.2%
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 NavierStokesTests
10 * \copydoc Dumux::NavierStokesAnalyticalSolutionVectors
11 */
12 #ifndef DUMUX_TEST_FREEFLOW_NAVIERSTOKES_ANALYTICALSOLVECTORS_HH
13 #define DUMUX_TEST_FREEFLOW_NAVIERSTOKES_ANALYTICALSOLVECTORS_HH
14
15 #include <dune/common/fvector.hh>
16
17 #include <dumux/discretization/method.hh>
18
19 namespace Dumux {
20 /*!
21 * \ingroup NavierStokesTests
22 * \brief Creates and provides the analytical solution in a form that can be written into the vtk output files
23 */
24 template<class Problem, class Scalar = double>
25 class NavierStokesAnalyticalSolutionVectors
26 {
27 using GridGeometry = std::decay_t<decltype(std::declval<Problem>().gridGeometry())>;
28 static constexpr int dimWorld = GridGeometry::GridView::dimensionworld;
29 using VelocityVector = Dune::FieldVector<Scalar, dimWorld>;
30 using Indices = typename Problem::Indices;
31
32 public:
33 9 NavierStokesAnalyticalSolutionVectors(std::shared_ptr<const Problem> problem, Scalar tInitial = 0.0)
34
3/14
✓ Branch 2 taken 9 times.
✗ Branch 3 not taken.
✓ Branch 5 taken 9 times.
✗ Branch 6 not taken.
✓ Branch 8 taken 9 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.
✗ Branch 16 not taken.
✗ Branch 17 not taken.
 9 : problem_(problem)
35
1/2
✓ Branch 1 taken 9 times.
✗ Branch 2 not taken.
 9 { update(tInitial); }
36
37 /*!
38 * \brief Creates the analytical solution in a form that can be written into the vtk output files
39 */
40 9 void update(Scalar time = 0.0)
41 {
42 27 analyticalPressure_.resize(problem_->gridGeometry().numCellCenterDofs());
43 27 analyticalVelocity_.resize(problem_->gridGeometry().numCellCenterDofs());
44 27 analyticalVelocityAtDofs_.resize(problem_->gridGeometry().numFaceDofs());
45
46 27 auto fvGeometry = localView(problem_->gridGeometry());
47
48
1/2
✓ Branch 4 taken 74209 times.
✗ Branch 5 not taken.
 148436 for (const auto& element : elements(problem_->gridGeometry().gridView()))
49 {
50 74200 fvGeometry.bindElement(element);
51
4/4
✓ Branch 2 taken 74200 times.
✓ Branch 3 taken 74200 times.
✓ Branch 4 taken 74200 times.
✓ Branch 5 taken 74200 times.
 222600 for (const auto& scv : scvs(fvGeometry))
52 {
53 // velocities on faces
54
4/4
✓ Branch 0 taken 296800 times.
✓ Branch 1 taken 74200 times.
✓ Branch 2 taken 296800 times.
✓ Branch 3 taken 74200 times.
 742000 for (const auto& scvf : scvfs(fvGeometry))
55 {
56 890400 analyticalVelocityAtDofs_[scvf.dofIndex()][scvf.directionIndex()] = problem_->analyticalSolution(scvf.center(), time)[Indices::velocity(scvf.directionIndex())];
57 }
58
59 222600 analyticalPressure_[scv.dofIndex()] = problem_->analyticalSolution(scv.dofPosition(), time)[Indices::pressureIdx];
60
61
2/2
✓ Branch 0 taken 148400 times.
✓ Branch 1 taken 74200 times.
 222600 for (int dirIdx = 0; dirIdx < dimWorld; ++dirIdx)
62 296800 analyticalVelocity_[scv.dofIndex()][dirIdx] = problem_->analyticalSolution(scv.center(), time)[Indices::velocity(dirIdx)];
63 }
64 }
65 9 }
66
67 /*!
68 * \brief Returns the analytical solution for the pressure
69 */
70 const std::vector<Scalar>& getAnalyticalPressureSolution() const
71 {
72 return analyticalPressure_;
73 }
74
75 /*!
76 * \brief Returns the analytical solution for the velocity
77 */
78 const std::vector<VelocityVector>& getAnalyticalVelocitySolution() const
79 {
80 return analyticalVelocity_;
81 }
82
83 /*!
84 * \brief Returns the analytical solution for the velocity at the faces
85 */
86 const std::vector<VelocityVector>& getAnalyticalVelocitySolutionAtDofs() const
87 {
88 return analyticalVelocityAtDofs_;
89 }
90
91 private:
92 std::shared_ptr<const Problem> problem_;
93
94 std::vector<Scalar> analyticalPressure_;
95 std::vector<VelocityVector> analyticalVelocity_;
96 std::vector<VelocityVector> analyticalVelocityAtDofs_;
97 };
98
99 template<class Problem>
100 NavierStokesAnalyticalSolutionVectors(std::shared_ptr<Problem> p)
101 -> NavierStokesAnalyticalSolutionVectors<Problem>;
102
103 template<class Problem, class Scalar>
104 NavierStokesAnalyticalSolutionVectors(std::shared_ptr<Problem> p, Scalar t)
105 -> NavierStokesAnalyticalSolutionVectors<Problem, Scalar>;
106
107
108 namespace NavierStokesTest {
109
110 /*!
111 * \ingroup NavierStokesTests
112 * \brief Creates and provides the analytical solution in a form that can be written into the vtk output files
113 * \TODO make this the new NavierStokesAnalyticalSolutionVectors once all tests are ported to the new staggered
114 */
115 template<class MomentumProblem, class MassProblem, class Scalar = double>
116 class AnalyticalSolutionVectors
117 {
118 using MassGridGeometry = std::decay_t<decltype(std::declval<MassProblem>().gridGeometry())>;
119 using MomentumGridGeometry = std::decay_t<decltype(std::declval<MomentumProblem>().gridGeometry())>;
120 static constexpr int dimWorld = MassGridGeometry::GridView::dimensionworld;
121 using VelocityVector = Dune::FieldVector<Scalar, dimWorld>;
122
123 using MassIndices = typename MassProblem::Indices;
124 using MomIndices = typename MomentumProblem::Indices;
125 public:
126 27 AnalyticalSolutionVectors(std::shared_ptr<const MomentumProblem> momentumProblem,
127 std::shared_ptr<const MassProblem> massProblem,
128 Scalar tInitial = 0.0)
129 : momentumProblem_(momentumProblem)
130
3/16
✓ Branch 3 taken 27 times.
✗ Branch 4 not taken.
✓ Branch 6 taken 27 times.
✗ Branch 7 not taken.
✓ Branch 9 taken 27 times.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.
✗ Branch 16 not taken.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
 27 , massProblem_(massProblem)
131
1/2
✓ Branch 1 taken 27 times.
✗ Branch 2 not taken.
 27 { update(tInitial); }
132
133 /*!
134 * \brief Creates the analytical solution in a form that can be written into the vtk output files
135 */
136 83 void update(Scalar time = 0.0)
137 {
138 250 analyticalPressure_.resize(massProblem_->gridGeometry().numDofs());
139 332 analyticalVelocity_.resize(massProblem_->gridGeometry().gridView().size(0));
140 250 analyticalVelocityAtDofs_.resize(momentumProblem_->gridGeometry().numDofs());
141
142 // cell-centers (pressure + velocity)
143 {
144
3/8
✓ Branch 1 taken 7 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 7 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 7 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
 249 auto fvGeometry = localView(massProblem_->gridGeometry());
145
4/8
✓ Branch 1 taken 7 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 7 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 7 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 7 times.
✗ Branch 11 not taken.
 332 const auto gridView = massProblem_->gridGeometry().gridView();
146
10/16
✓ Branch 1 taken 505765 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 7 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 7 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 6736 times.
✓ Branch 10 taken 7 times.
✓ Branch 11 taken 6736 times.
✓ Branch 12 taken 7 times.
✓ Branch 14 taken 6736 times.
✗ Branch 15 not taken.
✓ Branch 17 taken 6736 times.
✗ Branch 18 not taken.
✓ Branch 22 taken 4 times.
✗ Branch 23 not taken.
 1031676 for (const auto& element : elements(gridView))
147 {
148 // output velocity always on elements
149
4/8
✓ Branch 1 taken 6736 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 6736 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 6736 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 6736 times.
✗ Branch 11 not taken.
 2049672 const auto eIdx = massProblem_->gridGeometry().elementMapper().index(element);
150
2/4
✓ Branch 1 taken 6736 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 6736 times.
✗ Branch 5 not taken.
 512418 const auto center = element.geometry().center();
151
2/2
✓ Branch 0 taken 1024804 times.
✓ Branch 1 taken 512418 times.
 1537222 for (int dirIdx = 0; dirIdx < dimWorld; ++dirIdx)
152 1024804 analyticalVelocity_[eIdx][dirIdx]
153 2057140 = momentumProblem_->analyticalSolution(center, time)[MomIndices::velocity(dirIdx)];
154
155 // this works for box and cc methods
156
1/2
✓ Branch 1 taken 6736 times.
✗ Branch 2 not taken.
 512418 fvGeometry.bindElement(element);
157
6/6
✓ Branch 0 taken 13136 times.
✓ Branch 1 taken 6736 times.
✓ Branch 2 taken 517050 times.
✓ Branch 3 taken 510650 times.
✓ Branch 4 taken 505682 times.
✓ Branch 5 taken 505682 times.
 1742136 for (const auto& scv : scvs(fvGeometry))
158 971786 analyticalPressure_[scv.dofIndex()]
159 1884622 = massProblem_->analyticalSolution(scv.dofPosition(), time)[MassIndices::pressureIdx];
160 }
161 }
162
163 // dof positions (velocity)
164 {
165
3/8
✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 1 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
 249 auto fvGeometry = localView(momentumProblem_->gridGeometry());
166
4/8
✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 1 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
 332 const auto gridView = momentumProblem_->gridGeometry().gridView();
167
14/20
✓ Branch 1 taken 505759 times.
✓ Branch 2 taken 6 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 1 times.
✓ Branch 5 taken 6 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 3537 times.
✓ Branch 8 taken 6 times.
✓ Branch 9 taken 6736 times.
✓ Branch 10 taken 7 times.
✓ Branch 11 taken 3200 times.
✓ Branch 12 taken 3537 times.
✗ Branch 13 not taken.
✓ Branch 14 taken 3200 times.
✓ Branch 15 taken 3536 times.
✗ Branch 16 not taken.
✓ Branch 17 taken 3200 times.
✗ Branch 18 not taken.
✓ Branch 22 taken 1 times.
✗ Branch 23 not taken.
 1031670 for (const auto& element : elements(gridView))
168 {
169
1/2
✓ Branch 1 taken 6736 times.
✗ Branch 2 not taken.
 512418 fvGeometry.bindElement(element);
170
171 if constexpr (MomentumGridGeometry::discMethod == DiscretizationMethods::fcstaggered)
172
5/5
✓ Branch 0 taken 7072 times.
✓ Branch 1 taken 2031504 times.
✓ Branch 2 taken 507450 times.
✓ Branch 3 taken 2029736 times.
✓ Branch 4 taken 507450 times.
 4575762 for (const auto& scv : scvs(fvGeometry))
173 2036808 analyticalVelocityAtDofs_[scv.dofIndex()][scv.dofAxis()]
174 6125488 = momentumProblem_->analyticalSolution(scv.center(), time)[MomIndices::velocity(scv.dofAxis())];
175
176 else if constexpr (DiscretizationMethods::isCVFE<typename MomentumGridGeometry::DiscretizationMethod>)
177
4/4
✓ Branch 0 taken 9600 times.
✓ Branch 1 taken 3200 times.
✓ Branch 2 taken 9600 times.
✓ Branch 3 taken 3200 times.
 16000 for (const auto& scv : scvs(fvGeometry))
178
2/2
✓ Branch 0 taken 19200 times.
✓ Branch 1 taken 9600 times.
 28800 for (int dirIdx = 0; dirIdx < dimWorld; ++dirIdx)
179 19200 analyticalVelocityAtDofs_[scv.dofIndex()][dirIdx]
180 57600 = momentumProblem_->analyticalSolution(scv.dofPosition(), time)[MomIndices::velocity(dirIdx)];
181
182 else
183 DUNE_THROW(Dune::Exception, "Unknown discretization method: " << MomentumGridGeometry::discMethod);
184 }
185 }
186 83 }
187
188 /*!
189 * \brief Returns the analytical solution for the pressure
190 */
191 const std::vector<Scalar>& analyticalPressureSolution() const
192 {
193
1/2
✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
 1 return analyticalPressure_;
194 }
195
196 /*!
197 * \brief Returns the analytical solution for the velocity
198 */
199 const std::vector<VelocityVector>& analyticalVelocitySolution() const
200 {
201
1/2
✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
 1 return analyticalVelocity_;
202 }
203
204 /*!
205 * \brief Returns the analytical solution for the velocity at the faces
206 */
207 const std::vector<VelocityVector>& analyticalVelocitySolutionAtDofs() const
208 {
209 return analyticalVelocityAtDofs_;
210 }
211
212 private:
213 std::shared_ptr<const MomentumProblem> momentumProblem_;
214 std::shared_ptr<const MassProblem> massProblem_;
215
216 std::vector<Scalar> analyticalPressure_;
217 std::vector<VelocityVector> analyticalVelocity_;
218 std::vector<VelocityVector> analyticalVelocityAtDofs_;
219 };
220
221 template<class MomentumProblem, class MassProblem>
222 AnalyticalSolutionVectors(std::shared_ptr<MomentumProblem>, std::shared_ptr<MassProblem>)
223 -> AnalyticalSolutionVectors<MomentumProblem, MassProblem>;
224
225 template<class MomentumProblem, class MassProblem, class Scalar>
226 AnalyticalSolutionVectors(std::shared_ptr<MomentumProblem>, std::shared_ptr<MassProblem>, Scalar)
227 -> AnalyticalSolutionVectors<MomentumProblem, MassProblem, Scalar>;
228
229
230 } // end namespace NavierStokesTest
231 } // end namespace Dumux
232
233 #endif
234