GCC Code Coverage Report


Directory: ../../../builds/dumux-repositories/
File: dumux/dumux/discretization/porenetwork/subcontrolvolumeface.hh
Date: 2025-06-14 19:21:29
Exec Total Coverage
Lines: 13 13 100.0%
Functions: 0 0 -%
Branches: 80 109 73.4%

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-FileCopyrightText: 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 PoreNetworkDiscretization
10 * \brief Base class for a sub control volume face
11 */
12 #ifndef DUMUX_DISCRETIZATION_PNM_SUBCONTROLVOLUMEFACE_HH
13 #define DUMUX_DISCRETIZATION_PNM_SUBCONTROLVOLUMEFACE_HH
14
15 #include <dune/geometry/axisalignedcubegeometry.hh>
16 #include <dumux/common/indextraits.hh>
17
18 namespace Dumux::PoreNetwork {
19
20 /*!
21 * \ingroup PoreNetworkDiscretization
22 * \brief Default traits class
23 * \tparam GV the type of the grid view
24 */
25 template<class GridView>
26 struct PNMDefaultScvfGeometryTraits
27 {
28 using Grid = typename GridView::Grid;
29 static constexpr int dim = Grid::dimension;
30 static constexpr int dimWorld = Grid::dimensionworld;
31
32 using GridIndexType = typename IndexTraits<GridView>::GridIndex;
33 using LocalIndexType = typename IndexTraits<GridView>::LocalIndex;
34
35 using Scalar = typename Grid::ctype;
36 using GlobalPosition = Dune::FieldVector<Scalar, dimWorld>;
37 using CornerStorage = std::array<Dune::FieldVector<Scalar, dimWorld>, 1>;
38 using Geometry = Dune::AxisAlignedCubeGeometry<Scalar, dim-1, dimWorld>;
39
40 };
41
42 /*!
43 * \ingroup PoreNetworkDiscretization
44 * \brief Class for a sub control volume face for porenetworks
45 * \tparam GV the type of the grid view
46 * \tparam T the scvf geometry traits
47 */
48 template<class GV,
49 class T = PNMDefaultScvfGeometryTraits<GV> >
50 class PNMSubControlVolumeFace
51 {
52 using ThisType = PNMSubControlVolumeFace<GV, T>;
53 using GridIndexType = typename T::GridIndexType;
54 using LocalIndexType = typename T::LocalIndexType;
55 using Scalar = typename T::Scalar;
56
57 public:
58 //! export the type used for global coordinates
59 using GlobalPosition = typename T::GlobalPosition;
60 //! state the traits public and thus export all types
61 using Traits = T;
62
63 //! The default constructor
64
15/33
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✓ Branch 15 taken 1 times.
✗ Branch 16 not taken.
✓ Branch 18 taken 1 times.
✗ Branch 19 not taken.
✓ Branch 21 taken 2 times.
✗ Branch 22 not taken.
✓ Branch 24 taken 1 times.
✗ Branch 25 not taken.
✓ Branch 29 taken 338742 times.
✗ Branch 30 not taken.
✓ Branch 32 taken 338742 times.
✗ Branch 33 not taken.
✓ Branch 35 taken 182830 times.
✗ Branch 36 not taken.
✓ Branch 38 taken 182830 times.
✗ Branch 39 not taken.
✓ Branch 43 taken 2 times.
✗ Branch 44 not taken.
✓ Branch 46 taken 2 times.
✗ Branch 47 not taken.
✓ Branch 49 taken 1 times.
✗ Branch 50 not taken.
✓ Branch 52 taken 1 times.
✗ Branch 53 not taken.
✓ Branch 6 taken 952 times.
2807832 PNMSubControlVolumeFace() = default;
65
66 //! Constructor for inner scvfs
67 3955347 PNMSubControlVolumeFace(const GlobalPosition& center,
68 const GlobalPosition& unitOuterNormal,
69 const Scalar& area,
70 GridIndexType scvfIndex,
71 std::array<LocalIndexType, 2>&& scvIndices)
72 3955347 : center_(center),
73 3955347 unitOuterNormal_(unitOuterNormal),
74 area_(area),
75 scvfIndex_(scvfIndex),
76 scvIndices_(std::move(scvIndices))
77 {}
78
79 //! The center of the sub control volume face
80 const GlobalPosition& center() const
81 1872 { return center_; }
82
83 //! The integration point for flux evaluations in global coordinates
84 const GlobalPosition& ipGlobal() const
85 1242046 { return center_; }
86
87 //! The area of the sub control volume face
88 Scalar area() const
89 { return area_; }
90
91 //! We assume to always have a pore body and not a pore throat at the boundary
92 bool boundary() const
93 { return false; }
94
95 //! The unit outer normal of the sub control volume face
96 const GlobalPosition& unitOuterNormal() const
97 1872 { return unitOuterNormal_; }
98
99 //! Index of the inside sub control volume for spatial param evaluation
100 89464026 LocalIndexType insideScvIdx() const
101
26/28
✓ Branch 0 taken 31854 times.
✓ Branch 1 taken 22413 times.
✓ Branch 2 taken 1228437 times.
✓ Branch 3 taken 896284 times.
✓ Branch 4 taken 1776342 times.
✓ Branch 5 taken 2223717 times.
✓ Branch 6 taken 369536 times.
✓ Branch 7 taken 953223 times.
✓ Branch 8 taken 2532950 times.
✓ Branch 9 taken 823435 times.
✓ Branch 10 taken 7330 times.
✓ Branch 11 taken 1174600 times.
✓ Branch 12 taken 1092264 times.
✗ Branch 13 not taken.
✓ Branch 14 taken 4385 times.
✓ Branch 15 taken 10999 times.
✓ Branch 16 taken 396 times.
✓ Branch 17 taken 2524594 times.
✓ Branch 18 taken 22032 times.
✓ Branch 19 taken 2539784 times.
✓ Branch 20 taken 6681 times.
✓ Branch 21 taken 1129442 times.
✓ Branch 22 taken 10470328 times.
✓ Branch 23 taken 21478346 times.
✗ Branch 24 not taken.
✓ Branch 25 taken 32058562 times.
✓ Branch 26 taken 18206 times.
✓ Branch 27 taken 25653 times.
85813058 { return scvIndices_[0]; }
102
103 //! Index of the outside sub control volume for spatial param evaluation
104 // This results in undefined behaviour if boundary is true
105 89381325 LocalIndexType outsideScvIdx() const
106
22/24
✓ Branch 0 taken 31854 times.
✓ Branch 1 taken 22413 times.
✓ Branch 2 taken 1228437 times.
✓ Branch 3 taken 896284 times.
✓ Branch 4 taken 1776342 times.
✓ Branch 5 taken 2223717 times.
✓ Branch 6 taken 369536 times.
✓ Branch 7 taken 953223 times.
✓ Branch 8 taken 2532950 times.
✓ Branch 9 taken 825855 times.
✓ Branch 10 taken 4910 times.
✓ Branch 11 taken 1174600 times.
✓ Branch 12 taken 1092264 times.
✗ Branch 13 not taken.
✓ Branch 14 taken 4385 times.
✓ Branch 15 taken 10999 times.
✓ Branch 16 taken 396 times.
✓ Branch 17 taken 2524594 times.
✓ Branch 18 taken 6681 times.
✓ Branch 19 taken 2560154 times.
✓ Branch 20 taken 10470328 times.
✓ Branch 21 taken 22570610 times.
✗ Branch 22 not taken.
✓ Branch 23 taken 32058562 times.
85730357 { return scvIndices_[1]; }
107
108 //! The local index of this sub control volume face
109 140633325 GridIndexType index() const
110
17/24
✓ Branch 0 taken 20364 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 24066409 times.
✓ Branch 4 taken 302005 times.
✓ Branch 5 taken 10231662 times.
✓ Branch 7 taken 132634 times.
✓ Branch 8 taken 2546135 times.
✓ Branch 10 taken 48084 times.
✓ Branch 11 taken 1174589 times.
✓ Branch 13 taken 61904 times.
✓ Branch 14 taken 32058562 times.
✓ Branch 17 taken 1 times.
✓ Branch 18 taken 31948673 times.
✓ Branch 19 taken 1 times.
✓ Branch 20 taken 31948673 times.
✗ Branch 21 not taken.
✓ Branch 22 taken 8736 times.
✗ Branch 23 not taken.
✓ Branch 24 taken 8056 times.
✗ Branch 12 not taken.
✓ Branch 15 taken 61536 times.
✗ Branch 6 not taken.
✗ Branch 9 not taken.
108684652 { return scvfIndex_; }
111
112 private:
113 GlobalPosition center_;
114 GlobalPosition unitOuterNormal_;
115 Scalar area_;
116 GridIndexType scvfIndex_;
117 std::array<LocalIndexType, 2> scvIndices_;
118 };
119
120 } // end namespace Dumux::PoreNetwork
121
122 #endif
123