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 StaggeredDiscretization | ||
10 | * \copydoc Dumux::FreeFlowStaggeredSubControlVolumeFace | ||
11 | */ | ||
12 | #ifndef DUMUX_DISCRETIZATION_STAGGERED_FREE_FLOW_SUBCONTROLVOLUMEFACE_HH | ||
13 | #define DUMUX_DISCRETIZATION_STAGGERED_FREE_FLOW_SUBCONTROLVOLUMEFACE_HH | ||
14 | |||
15 | #include <array> | ||
16 | #include <utility> | ||
17 | #include <dune/common/fvector.hh> | ||
18 | #include <dune/geometry/type.hh> | ||
19 | #include <dune/geometry/multilineargeometry.hh> | ||
20 | |||
21 | #include <dumux/common/indextraits.hh> | ||
22 | #include <dumux/discretization/subcontrolvolumefacebase.hh> | ||
23 | #include <dumux/discretization/staggered/subcontrolvolumeface.hh> | ||
24 | #include <dumux/discretization/staggered/freeflow/staggeredgeometryhelper.hh> | ||
25 | |||
26 | #include <typeinfo> | ||
27 | |||
28 | namespace Dumux { | ||
29 | |||
30 | /*! | ||
31 | * \ingroup StaggeredDiscretization | ||
32 | * \brief Default traits class to be used for the sub-control volume faces | ||
33 | * for the free-flow staggered finite volume scheme | ||
34 | * \tparam GridView the type of the grid view | ||
35 | * \tparam upwindSchemeOrder the order of the upwind scheme | ||
36 | */ | ||
37 | template<class GridView, int upwindSchemeOrder> | ||
38 | struct FreeFlowStaggeredDefaultScvfGeometryTraits | ||
39 | { | ||
40 | using GridIndexType = typename IndexTraits<GridView>::GridIndex; | ||
41 | using LocalIndexType = typename IndexTraits<GridView>::LocalIndex; | ||
42 | using Scalar = typename GridView::ctype; | ||
43 | using GeometryHelper = FreeFlowStaggeredGeometryHelper<GridView, upwindSchemeOrder>; | ||
44 | using PairData = typename GeometryHelper::PairData; | ||
45 | using AxisData = typename GeometryHelper::AxisData; | ||
46 | |||
47 | using Grid = typename GridView::Grid; | ||
48 | static constexpr int dim = Grid::dimension; | ||
49 | static constexpr int dimWorld = Grid::dimensionworld; | ||
50 | |||
51 | using Element = typename GridView::template Codim<0>::Entity; | ||
52 | using GlobalPosition = typename Element::Geometry::GlobalCoordinate; | ||
53 | using Geometry = Dune::AxisAlignedCubeGeometry<Scalar, dim-1, dimWorld>; | ||
54 | }; | ||
55 | |||
56 | /*! | ||
57 | * \ingroup StaggeredDiscretization | ||
58 | * \brief Helper function to turn a given cell scvface into a fake boundary face | ||
59 | * \note This function is considered internal to staggered freeflow and may change or be deleted at any time without deprecation warning | ||
60 | */ | ||
61 | template<class SubControlVolumeFace> | ||
62 | 17557066 | SubControlVolumeFace makeStaggeredBoundaryFace(const SubControlVolumeFace& scvf, | |
63 | const typename SubControlVolumeFace::GlobalPosition& newCenter) | ||
64 | { | ||
65 |
0/9✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✗ Branch 8 not taken.
✗ Branch 11 not taken.
✗ Branch 14 not taken.
|
17557066 | auto bf = scvf; |
66 | 17557066 | bf.setCenter(newCenter); | |
67 | 17557066 | bf.setBoundary(true); | |
68 | 17557066 | bf.setIsGhostFace(true); | |
69 | return bf; | ||
70 | } | ||
71 | |||
72 | /*! | ||
73 | * \ingroup StaggeredDiscretization | ||
74 | * \brief Class for a sub control volume face in the staggered method, i.e a part of the boundary | ||
75 | * of a sub control volume we compute fluxes on. This is a specialization for free flow models. | ||
76 | */ | ||
77 | template<class GV, | ||
78 | int upwindSchemeOrder, | ||
79 | class T = FreeFlowStaggeredDefaultScvfGeometryTraits<GV, upwindSchemeOrder>> | ||
80 | 17939804 | class FreeFlowStaggeredSubControlVolumeFace | |
81 | : public SubControlVolumeFaceBase<FreeFlowStaggeredSubControlVolumeFace<GV, upwindSchemeOrder, T>, T> | ||
82 | { | ||
83 | using ThisType = FreeFlowStaggeredSubControlVolumeFace<GV, upwindSchemeOrder, T>; | ||
84 | using ParentType = SubControlVolumeFaceBase<ThisType, T>; | ||
85 | using Geometry = typename T::Geometry; | ||
86 | using GridIndexType = typename IndexTraits<GV>::GridIndex; | ||
87 | using LocalIndexType = typename IndexTraits<GV>::LocalIndex; | ||
88 | |||
89 | using PairData = typename T::PairData; | ||
90 | using AxisData = typename T::AxisData; | ||
91 | |||
92 | using Scalar = typename T::Scalar; | ||
93 | static const int dim = GV::dimension; | ||
94 | |||
95 | static constexpr int numPairs = 2 * (dim - 1); | ||
96 | |||
97 | static constexpr bool useHigherOrder = upwindSchemeOrder > 1; | ||
98 | |||
99 | public: | ||
100 | using GlobalPosition = typename T::GlobalPosition; | ||
101 | |||
102 | static constexpr int numCornersPerFace = 2 * (dim - 1); | ||
103 | |||
104 | //! State the traits public and thus export all types | ||
105 | using Traits = T; | ||
106 | |||
107 | // The default constructor | ||
108 | FreeFlowStaggeredSubControlVolumeFace() = default; | ||
109 | |||
110 | //! Constructor with intersection | ||
111 | template <class Intersection> | ||
112 | 382738 | FreeFlowStaggeredSubControlVolumeFace(const Intersection& is, | |
113 | const typename Intersection::Geometry& isGeometry, | ||
114 | GridIndexType scvfIndex, | ||
115 | const std::vector<GridIndexType>& scvIndices, | ||
116 | const typename T::GeometryHelper& geometryHelper) | ||
117 | : ParentType(), | ||
118 | 382738 | area_(isGeometry.volume()), | |
119 |
2/2✓ Branch 0 taken 177369 times.
✓ Branch 1 taken 177369 times.
|
382738 | center_(isGeometry.center()), |
120 | 382738 | unitOuterNormal_(is.centerUnitOuterNormal()), | |
121 | 382738 | scvfIndex_(scvfIndex), | |
122 |
2/2✓ Branch 1 taken 28375 times.
✓ Branch 2 taken 375 times.
|
382738 | scvIndices_(scvIndices), |
123 |
2/2✓ Branch 0 taken 375 times.
✓ Branch 1 taken 375 times.
|
382738 | boundary_(is.boundary()), |
124 |
2/2✓ Branch 0 taken 375 times.
✓ Branch 1 taken 375 times.
|
382738 | axisData_(geometryHelper.axisData()), |
125 |
2/2✓ Branch 0 taken 375 times.
✓ Branch 1 taken 375 times.
|
382738 | pairData_(std::move(geometryHelper.pairData())), |
126 |
2/2✓ Branch 0 taken 375 times.
✓ Branch 1 taken 375 times.
|
382738 | localFaceIdx_(geometryHelper.localFaceIndex()), |
127 | 382738 | dirIdx_(geometryHelper.directionIndex()), | |
128 | 382738 | outerNormalSign_(sign(unitOuterNormal_[directionIndex()])), | |
129 | 382738 | isGhostFace_(false) | |
130 | { | ||
131 | 766226 | dimensions[0] = (isGeometry.corner(1) - isGeometry.corner(0)).two_norm(); | |
132 | if constexpr (dim == 3) | ||
133 | 1500 | dimensions[1] = (isGeometry.corner(2) - isGeometry.corner(0)).two_norm(); | |
134 | 382738 | } | |
135 | |||
136 | //! The center of the sub control volume face | ||
137 |
1/2✓ Branch 1 taken 268800 times.
✗ Branch 2 not taken.
|
88557996 | const GlobalPosition& center() const |
138 | { | ||
139 |
28/92✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✓ Branch 7 taken 10 times.
✗ Branch 8 not taken.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✓ Branch 13 taken 9970 times.
✓ Branch 14 taken 115050 times.
✗ Branch 16 not taken.
✗ Branch 17 not taken.
✓ Branch 19 taken 251634 times.
✓ Branch 20 taken 1710 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
✓ Branch 25 taken 10 times.
✗ Branch 26 not taken.
✗ Branch 28 not taken.
✗ Branch 29 not taken.
✓ Branch 31 taken 10 times.
✗ Branch 32 not taken.
✗ Branch 34 not taken.
✗ Branch 35 not taken.
✓ Branch 37 taken 10 times.
✗ Branch 38 not taken.
✗ Branch 40 not taken.
✗ Branch 41 not taken.
✓ Branch 43 taken 10 times.
✗ Branch 44 not taken.
✗ Branch 46 not taken.
✗ Branch 47 not taken.
✓ Branch 49 taken 10 times.
✗ Branch 50 not taken.
✗ Branch 52 not taken.
✗ Branch 53 not taken.
✓ Branch 55 taken 10 times.
✗ Branch 56 not taken.
✗ Branch 58 not taken.
✗ Branch 59 not taken.
✓ Branch 61 taken 10 times.
✗ Branch 62 not taken.
✗ Branch 64 not taken.
✗ Branch 65 not taken.
✓ Branch 67 taken 10 times.
✗ Branch 68 not taken.
✗ Branch 70 not taken.
✗ Branch 71 not taken.
✓ Branch 73 taken 10 times.
✗ Branch 74 not taken.
✗ Branch 76 not taken.
✗ Branch 77 not taken.
✓ Branch 79 taken 10 times.
✗ Branch 80 not taken.
✗ Branch 82 not taken.
✗ Branch 83 not taken.
✓ Branch 85 taken 10 times.
✗ Branch 86 not taken.
✗ Branch 88 not taken.
✗ Branch 89 not taken.
✓ Branch 91 taken 10 times.
✗ Branch 92 not taken.
✗ Branch 94 not taken.
✗ Branch 95 not taken.
✓ Branch 97 taken 10 times.
✗ Branch 98 not taken.
✗ Branch 100 not taken.
✗ Branch 101 not taken.
✓ Branch 103 taken 10 times.
✗ Branch 104 not taken.
✗ Branch 106 not taken.
✗ Branch 107 not taken.
✓ Branch 109 taken 10 times.
✗ Branch 110 not taken.
✗ Branch 112 not taken.
✗ Branch 113 not taken.
✓ Branch 115 taken 10 times.
✗ Branch 116 not taken.
✗ Branch 118 not taken.
✗ Branch 119 not taken.
✓ Branch 121 taken 5 times.
✗ Branch 122 not taken.
✓ Branch 124 taken 5 times.
✗ Branch 125 not taken.
✓ Branch 127 taken 5 times.
✗ Branch 128 not taken.
✓ Branch 130 taken 5 times.
✗ Branch 131 not taken.
✗ Branch 15 not taken.
✓ Branch 18 taken 1680168 times.
✓ Branch 21 taken 20 times.
✗ Branch 12 not taken.
|
859602041 | return center_; |
140 | } | ||
141 | |||
142 | //! The position of the dof living on the face | ||
143 | const GlobalPosition& dofPosition() const | ||
144 | { | ||
145 |
2/2✓ Branch 0 taken 864004 times.
✓ Branch 1 taken 1318948 times.
|
2377675 | return center_; |
146 | } | ||
147 | |||
148 | //! The integration point for flux evaluations in global coordinates | ||
149 |
2/10✓ Branch 1 taken 120800 times.
✗ Branch 2 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✓ Branch 7 taken 475916 times.
✗ Branch 8 not taken.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
|
421176812 | const GlobalPosition& ipGlobal() const |
150 | { | ||
151 | // Return center for now | ||
152 |
120/184✓ Branch 16 taken 37420 times.
✓ Branch 17 taken 6517 times.
✓ Branch 18 taken 10707 times.
✓ Branch 19 taken 2693 times.
✓ Branch 20 taken 968 times.
✓ Branch 21 taken 981 times.
✓ Branch 22 taken 5 times.
✓ Branch 23 taken 15 times.
✓ Branch 24 taken 5 times.
✓ Branch 25 taken 15 times.
✓ Branch 26 taken 5 times.
✓ Branch 27 taken 15 times.
✓ Branch 28 taken 5 times.
✓ Branch 29 taken 15 times.
✓ Branch 30 taken 5 times.
✓ Branch 31 taken 15 times.
✓ Branch 32 taken 5 times.
✓ Branch 33 taken 15 times.
✓ Branch 34 taken 5 times.
✓ Branch 35 taken 15 times.
✓ Branch 36 taken 5 times.
✓ Branch 37 taken 15 times.
✓ Branch 38 taken 5 times.
✓ Branch 39 taken 15 times.
✓ Branch 40 taken 5 times.
✓ Branch 41 taken 15 times.
✓ Branch 42 taken 5 times.
✓ Branch 43 taken 15 times.
✓ Branch 44 taken 5 times.
✓ Branch 45 taken 15 times.
✓ Branch 46 taken 5 times.
✓ Branch 47 taken 15 times.
✗ Branch 48 not taken.
✗ Branch 49 not taken.
✗ Branch 50 not taken.
✗ Branch 51 not taken.
✗ Branch 52 not taken.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✗ Branch 55 not taken.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
✗ Branch 58 not taken.
✗ Branch 59 not taken.
✗ Branch 60 not taken.
✗ Branch 61 not taken.
✗ Branch 62 not taken.
✗ Branch 63 not taken.
✗ Branch 64 not taken.
✗ Branch 65 not taken.
✗ Branch 66 not taken.
✗ Branch 67 not taken.
✗ Branch 68 not taken.
✗ Branch 69 not taken.
✗ Branch 70 not taken.
✗ Branch 71 not taken.
✓ Branch 72 taken 10 times.
✓ Branch 73 taken 30 times.
✓ Branch 74 taken 40 times.
✗ Branch 75 not taken.
✓ Branch 76 taken 1 times.
✗ Branch 77 not taken.
✓ Branch 78 taken 10 times.
✓ Branch 79 taken 30 times.
✓ Branch 80 taken 40 times.
✗ Branch 81 not taken.
✓ Branch 82 taken 1 times.
✗ Branch 83 not taken.
✓ Branch 84 taken 10 times.
✓ Branch 85 taken 30 times.
✓ Branch 86 taken 40 times.
✗ Branch 87 not taken.
✓ Branch 88 taken 1 times.
✗ Branch 89 not taken.
✓ Branch 90 taken 10 times.
✓ Branch 91 taken 30 times.
✓ Branch 92 taken 40 times.
✗ Branch 93 not taken.
✓ Branch 94 taken 1 times.
✗ Branch 95 not taken.
✓ Branch 96 taken 10 times.
✓ Branch 97 taken 30 times.
✓ Branch 98 taken 40 times.
✗ Branch 99 not taken.
✓ Branch 100 taken 1 times.
✗ Branch 101 not taken.
✓ Branch 102 taken 10 times.
✓ Branch 103 taken 30 times.
✓ Branch 104 taken 40 times.
✗ Branch 105 not taken.
✓ Branch 106 taken 1 times.
✗ Branch 107 not taken.
✓ Branch 108 taken 10 times.
✓ Branch 109 taken 30 times.
✓ Branch 110 taken 40 times.
✗ Branch 111 not taken.
✓ Branch 112 taken 1 times.
✗ Branch 113 not taken.
✓ Branch 114 taken 10 times.
✓ Branch 115 taken 30 times.
✓ Branch 116 taken 40 times.
✗ Branch 117 not taken.
✓ Branch 118 taken 1 times.
✗ Branch 119 not taken.
✓ Branch 120 taken 10 times.
✓ Branch 121 taken 30 times.
✓ Branch 122 taken 40 times.
✗ Branch 123 not taken.
✓ Branch 124 taken 1 times.
✗ Branch 125 not taken.
✗ Branch 126 not taken.
✗ Branch 127 not taken.
✓ Branch 128 taken 10 times.
✓ Branch 129 taken 30 times.
✓ Branch 130 taken 40 times.
✗ Branch 131 not taken.
✓ Branch 132 taken 1 times.
✗ Branch 133 not taken.
✗ Branch 134 not taken.
✗ Branch 135 not taken.
✓ Branch 136 taken 10 times.
✓ Branch 137 taken 30 times.
✓ Branch 138 taken 40 times.
✗ Branch 139 not taken.
✓ Branch 140 taken 1 times.
✗ Branch 141 not taken.
✗ Branch 142 not taken.
✗ Branch 143 not taken.
✓ Branch 144 taken 10 times.
✓ Branch 145 taken 30 times.
✓ Branch 146 taken 40 times.
✗ Branch 147 not taken.
✓ Branch 148 taken 1 times.
✗ Branch 149 not taken.
✗ Branch 150 not taken.
✗ Branch 151 not taken.
✓ Branch 152 taken 10 times.
✓ Branch 153 taken 30 times.
✓ Branch 154 taken 40 times.
✗ Branch 155 not taken.
✓ Branch 156 taken 1 times.
✗ Branch 157 not taken.
✓ Branch 158 taken 10 times.
✓ Branch 159 taken 30 times.
✓ Branch 160 taken 40 times.
✗ Branch 161 not taken.
✓ Branch 162 taken 1 times.
✗ Branch 163 not taken.
✓ Branch 164 taken 10 times.
✓ Branch 165 taken 30 times.
✓ Branch 166 taken 40 times.
✗ Branch 167 not taken.
✓ Branch 168 taken 1 times.
✗ Branch 169 not taken.
✓ Branch 170 taken 10 times.
✓ Branch 171 taken 30 times.
✓ Branch 172 taken 40 times.
✗ Branch 173 not taken.
✓ Branch 174 taken 1 times.
✗ Branch 175 not taken.
✓ Branch 216 taken 15 times.
✓ Branch 217 taken 5 times.
✓ Branch 218 taken 15 times.
✓ Branch 219 taken 5 times.
✓ Branch 220 taken 15 times.
✓ Branch 221 taken 5 times.
✓ Branch 222 taken 15 times.
✓ Branch 223 taken 5 times.
✓ Branch 0 taken 4219847 times.
✓ Branch 1 taken 1658998 times.
✓ Branch 2 taken 7973215 times.
✓ Branch 3 taken 35827622 times.
✓ Branch 4 taken 7183327 times.
✓ Branch 5 taken 67929796 times.
✓ Branch 6 taken 7628094 times.
✓ Branch 7 taken 89322116 times.
✓ Branch 8 taken 3418824 times.
✓ Branch 9 taken 101091358 times.
✓ Branch 10 taken 916325 times.
✓ Branch 11 taken 60179200 times.
✓ Branch 12 taken 70642 times.
✓ Branch 13 taken 20370 times.
✓ Branch 14 taken 68895 times.
✓ Branch 15 taken 23102 times.
|
989066090 | return center_; |
153 | } | ||
154 | |||
155 | //! The area of the sub control volume face | ||
156 | 1743375240 | Scalar area() const | |
157 | { | ||
158 |
0/13✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✗ Branch 5 not taken.
✗ Branch 8 not taken.
✗ Branch 11 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 12 not taken.
✗ Branch 15 not taken.
|
1740816670 | return area_; |
159 | } | ||
160 | |||
161 | //! Returns boolean if the sub control volume face is on the boundary | ||
162 | 3054004175 | bool boundary() const | |
163 | { | ||
164 |
182/222✓ Branch 0 taken 9482 times.
✓ Branch 1 taken 305512 times.
✓ Branch 2 taken 276913480 times.
✓ Branch 3 taken 8919454 times.
✓ Branch 4 taken 284552144 times.
✓ Branch 5 taken 346942 times.
✓ Branch 6 taken 22800 times.
✓ Branch 7 taken 406696 times.
✓ Branch 8 taken 12206567 times.
✓ Branch 9 taken 249716310 times.
✓ Branch 10 taken 6410072 times.
✓ Branch 11 taken 268713820 times.
✓ Branch 12 taken 12351308 times.
✓ Branch 13 taken 277463948 times.
✓ Branch 14 taken 193992513 times.
✓ Branch 15 taken 76117775 times.
✓ Branch 16 taken 43115378 times.
✓ Branch 17 taken 148121999 times.
✓ Branch 18 taken 45619644 times.
✓ Branch 19 taken 218280188 times.
✓ Branch 20 taken 9687280 times.
✓ Branch 21 taken 211476840 times.
✓ Branch 22 taken 9894198 times.
✓ Branch 23 taken 100298370 times.
✓ Branch 24 taken 16172880 times.
✓ Branch 25 taken 70506550 times.
✓ Branch 26 taken 79240854 times.
✓ Branch 27 taken 25511588 times.
✓ Branch 28 taken 100101350 times.
✓ Branch 29 taken 12427942 times.
✓ Branch 30 taken 62251260 times.
✓ Branch 31 taken 68097496 times.
✓ Branch 32 taken 50364916 times.
✓ Branch 33 taken 31118796 times.
✓ Branch 34 taken 19489924 times.
✓ Branch 35 taken 27623132 times.
✓ Branch 36 taken 1878794 times.
✓ Branch 37 taken 19458244 times.
✓ Branch 38 taken 374329 times.
✓ Branch 39 taken 2467226 times.
✓ Branch 40 taken 1228806 times.
✓ Branch 41 taken 3044018 times.
✓ Branch 42 taken 839268 times.
✓ Branch 43 taken 3190406 times.
✓ Branch 44 taken 383270 times.
✓ Branch 45 taken 760498 times.
✓ Branch 46 taken 394894 times.
✓ Branch 47 taken 871306 times.
✓ Branch 48 taken 45090 times.
✓ Branch 49 taken 520510 times.
✓ Branch 50 taken 33936 times.
✓ Branch 51 taken 388392 times.
✓ Branch 52 taken 18260 times.
✓ Branch 53 taken 224240 times.
✓ Branch 54 taken 1940 times.
✓ Branch 55 taken 23724 times.
✓ Branch 56 taken 20 times.
✓ Branch 57 taken 44 times.
✓ Branch 58 taken 20 times.
✓ Branch 59 taken 44 times.
✓ Branch 60 taken 20 times.
✓ Branch 61 taken 44 times.
✓ Branch 62 taken 20 times.
✓ Branch 63 taken 44 times.
✓ Branch 64 taken 20 times.
✓ Branch 65 taken 44 times.
✓ Branch 66 taken 20 times.
✓ Branch 67 taken 44 times.
✓ Branch 68 taken 20 times.
✓ Branch 69 taken 44 times.
✓ Branch 70 taken 20 times.
✓ Branch 71 taken 44 times.
✓ Branch 72 taken 20 times.
✓ Branch 73 taken 44 times.
✓ Branch 74 taken 20 times.
✓ Branch 75 taken 44 times.
✓ Branch 76 taken 20 times.
✓ Branch 77 taken 44 times.
✓ Branch 78 taken 20 times.
✓ Branch 79 taken 44 times.
✓ Branch 80 taken 20 times.
✓ Branch 81 taken 44 times.
✓ Branch 82 taken 20 times.
✓ Branch 83 taken 44 times.
✓ Branch 84 taken 20 times.
✓ Branch 85 taken 44 times.
✓ Branch 86 taken 20 times.
✓ Branch 87 taken 44 times.
✓ Branch 88 taken 20 times.
✓ Branch 89 taken 44 times.
✓ Branch 90 taken 20 times.
✓ Branch 91 taken 44 times.
✓ Branch 92 taken 20 times.
✓ Branch 93 taken 44 times.
✗ Branch 94 not taken.
✗ Branch 95 not taken.
✗ Branch 96 not taken.
✗ Branch 97 not taken.
✗ Branch 98 not taken.
✗ Branch 99 not taken.
✗ Branch 100 not taken.
✗ Branch 101 not taken.
✗ Branch 102 not taken.
✗ Branch 103 not taken.
✗ Branch 104 not taken.
✗ Branch 105 not taken.
✗ Branch 106 not taken.
✗ Branch 107 not taken.
✗ Branch 108 not taken.
✗ Branch 109 not taken.
✗ Branch 110 not taken.
✗ Branch 111 not taken.
✗ Branch 112 not taken.
✗ Branch 113 not taken.
✗ Branch 114 not taken.
✗ Branch 115 not taken.
✗ Branch 116 not taken.
✗ Branch 117 not taken.
✓ Branch 118 taken 20 times.
✓ Branch 119 taken 80 times.
✓ Branch 120 taken 40 times.
✓ Branch 121 taken 160 times.
✓ Branch 122 taken 20 times.
✓ Branch 123 taken 80 times.
✓ Branch 124 taken 40 times.
✓ Branch 125 taken 160 times.
✓ Branch 126 taken 20 times.
✓ Branch 127 taken 80 times.
✓ Branch 128 taken 40 times.
✓ Branch 129 taken 160 times.
✓ Branch 130 taken 20 times.
✓ Branch 131 taken 80 times.
✓ Branch 132 taken 40 times.
✓ Branch 133 taken 160 times.
✓ Branch 134 taken 20 times.
✓ Branch 135 taken 80 times.
✓ Branch 136 taken 40 times.
✓ Branch 137 taken 160 times.
✓ Branch 138 taken 20 times.
✓ Branch 139 taken 80 times.
✓ Branch 140 taken 40 times.
✓ Branch 141 taken 160 times.
✓ Branch 142 taken 20 times.
✓ Branch 143 taken 80 times.
✓ Branch 144 taken 40 times.
✓ Branch 145 taken 160 times.
✓ Branch 146 taken 20 times.
✓ Branch 147 taken 80 times.
✓ Branch 148 taken 40 times.
✓ Branch 149 taken 160 times.
✓ Branch 150 taken 20 times.
✓ Branch 151 taken 80 times.
✓ Branch 152 taken 40 times.
✓ Branch 153 taken 160 times.
✗ Branch 154 not taken.
✗ Branch 155 not taken.
✓ Branch 156 taken 20 times.
✓ Branch 157 taken 80 times.
✓ Branch 158 taken 40 times.
✓ Branch 159 taken 160 times.
✗ Branch 160 not taken.
✗ Branch 161 not taken.
✓ Branch 162 taken 20 times.
✓ Branch 163 taken 80 times.
✓ Branch 164 taken 40 times.
✓ Branch 165 taken 160 times.
✗ Branch 166 not taken.
✗ Branch 167 not taken.
✓ Branch 168 taken 20 times.
✓ Branch 169 taken 80 times.
✓ Branch 170 taken 40 times.
✓ Branch 171 taken 160 times.
✗ Branch 172 not taken.
✗ Branch 173 not taken.
✓ Branch 174 taken 20 times.
✓ Branch 175 taken 80 times.
✓ Branch 176 taken 40 times.
✓ Branch 177 taken 160 times.
✓ Branch 178 taken 20 times.
✓ Branch 179 taken 80 times.
✓ Branch 180 taken 40 times.
✓ Branch 181 taken 160 times.
✓ Branch 182 taken 20 times.
✓ Branch 183 taken 80 times.
✓ Branch 184 taken 40 times.
✓ Branch 185 taken 160 times.
✓ Branch 186 taken 20 times.
✓ Branch 187 taken 80 times.
✓ Branch 188 taken 40 times.
✓ Branch 189 taken 160 times.
✓ Branch 190 taken 20 times.
✓ Branch 191 taken 80 times.
✓ Branch 192 taken 40 times.
✓ Branch 193 taken 160 times.
✗ Branch 194 not taken.
✗ Branch 195 not taken.
✓ Branch 196 taken 20 times.
✓ Branch 197 taken 80 times.
✓ Branch 198 taken 40 times.
✓ Branch 199 taken 160 times.
✗ Branch 200 not taken.
✗ Branch 201 not taken.
✓ Branch 202 taken 20 times.
✓ Branch 203 taken 80 times.
✓ Branch 204 taken 40 times.
✓ Branch 205 taken 160 times.
✗ Branch 206 not taken.
✗ Branch 207 not taken.
✓ Branch 208 taken 20 times.
✓ Branch 209 taken 80 times.
✓ Branch 210 taken 40 times.
✓ Branch 211 taken 160 times.
✗ Branch 212 not taken.
✗ Branch 213 not taken.
✓ Branch 214 taken 20 times.
✓ Branch 215 taken 80 times.
✓ Branch 216 taken 20 times.
✓ Branch 217 taken 80 times.
✓ Branch 218 taken 20 times.
✓ Branch 219 taken 80 times.
✓ Branch 220 taken 20 times.
✓ Branch 221 taken 80 times.
|
3054004175 | return boundary_; |
165 | } | ||
166 | |||
167 | //! The unit outer normal vector | ||
168 | 100 | const GlobalPosition& unitOuterNormal() const | |
169 | { | ||
170 | 3156 | return unitOuterNormal_; | |
171 | } | ||
172 | |||
173 | //! Index of the inside sub control volume for spatial param evaluation | ||
174 | 4786421125 | GridIndexType insideScvIdx() const | |
175 | { | ||
176 |
37/56✗ Branch 1 not taken.
✓ Branch 2 taken 100 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 10 taken 238649564 times.
✓ Branch 11 taken 46248708 times.
✓ Branch 13 taken 5 times.
✓ Branch 14 taken 16260349 times.
✓ Branch 16 taken 63837422 times.
✓ Branch 17 taken 40796604 times.
✓ Branch 19 taken 1396466 times.
✓ Branch 20 taken 7777287 times.
✓ Branch 22 taken 17611028 times.
✓ Branch 23 taken 59379936 times.
✓ Branch 25 taken 34067004 times.
✓ Branch 26 taken 54457822 times.
✓ Branch 28 taken 13432822 times.
✓ Branch 29 taken 21045086 times.
✓ Branch 31 taken 3300 times.
✓ Branch 32 taken 3000 times.
✗ Branch 34 not taken.
✓ Branch 35 taken 428996 times.
✓ Branch 37 taken 93296 times.
✓ Branch 38 taken 166600 times.
✓ Branch 40 taken 527932 times.
✓ Branch 41 taken 462924 times.
✓ Branch 43 taken 12360 times.
✓ Branch 44 taken 152440 times.
✓ Branch 46 taken 156340 times.
✗ Branch 47 not taken.
✗ Branch 49 not taken.
✗ Branch 50 not taken.
✗ Branch 52 not taken.
✓ Branch 53 taken 5340 times.
✗ Branch 55 not taken.
✗ Branch 56 not taken.
✗ Branch 58 not taken.
✗ Branch 59 not taken.
✗ Branch 61 not taken.
✗ Branch 62 not taken.
✗ Branch 0 not taken.
✗ Branch 6 not taken.
✗ Branch 9 not taken.
✓ Branch 18 taken 15186246 times.
✓ Branch 21 taken 10466604 times.
✓ Branch 33 taken 6800 times.
✓ Branch 15 taken 140306149 times.
✓ Branch 24 taken 79280022 times.
✓ Branch 39 taken 350220 times.
✓ Branch 36 taken 6192882 times.
✓ Branch 45 taken 13980 times.
✓ Branch 42 taken 230288 times.
✓ Branch 54 taken 55860 times.
✓ Branch 27 taken 508168 times.
|
4641246758 | return scvIndices_[0]; |
177 | } | ||
178 | |||
179 | //! index of the outside sub control volume for spatial param evaluation | ||
180 | 2593914610 | GridIndexType outsideScvIdx() const | |
181 | { | ||
182 |
30/50✓ Branch 0 taken 20 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 20 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 20 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 20 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 20 times.
✓ Branch 9 taken 59452 times.
✓ Branch 10 taken 24660 times.
✓ Branch 11 taken 85684 times.
✓ Branch 12 taken 21140 times.
✓ Branch 13 taken 203956 times.
✓ Branch 14 taken 52676 times.
✓ Branch 15 taken 182196 times.
✓ Branch 16 taken 20 times.
✓ Branch 17 taken 11940 times.
✓ Branch 18 taken 63440 times.
✗ Branch 19 not taken.
✓ Branch 20 taken 20 times.
✗ Branch 21 not taken.
✓ Branch 22 taken 20 times.
✗ Branch 23 not taken.
✓ Branch 24 taken 20 times.
✗ Branch 25 not taken.
✓ Branch 26 taken 20 times.
✗ Branch 27 not taken.
✓ Branch 28 taken 20 times.
✗ Branch 29 not taken.
✓ Branch 30 taken 20 times.
✗ Branch 31 not taken.
✓ Branch 32 taken 20 times.
✗ Branch 33 not taken.
✓ Branch 34 taken 20 times.
✗ Branch 35 not taken.
✓ Branch 36 taken 20 times.
✗ Branch 37 not taken.
✓ Branch 38 taken 20 times.
✗ Branch 39 not taken.
✓ Branch 40 taken 20 times.
✗ Branch 41 not taken.
✓ Branch 42 taken 20 times.
✗ Branch 43 not taken.
✓ Branch 44 taken 20 times.
✗ Branch 45 not taken.
✓ Branch 46 taken 20 times.
✗ Branch 47 not taken.
✓ Branch 50 taken 80 times.
✗ Branch 51 not taken.
|
2652143540 | return scvIndices_[1]; |
183 | } | ||
184 | |||
185 | //! The global index of this sub control volume face | ||
186 | 1751608471 | GridIndexType index() const | |
187 | { | ||
188 |
38/45✓ Branch 1 taken 126520781 times.
✓ Branch 2 taken 76172184 times.
✓ Branch 3 taken 161011852 times.
✓ Branch 4 taken 88929601 times.
✓ Branch 5 taken 54363236 times.
✓ Branch 6 taken 95084685 times.
✓ Branch 7 taken 47270202 times.
✓ Branch 8 taken 97269581 times.
✓ Branch 9 taken 57471149 times.
✓ Branch 10 taken 92932704 times.
✓ Branch 11 taken 95507446 times.
✓ Branch 12 taken 57058141 times.
✓ Branch 13 taken 37102358 times.
✓ Branch 14 taken 12893240 times.
✓ Branch 15 taken 2234880 times.
✓ Branch 16 taken 3160622 times.
✓ Branch 17 taken 594920 times.
✓ Branch 18 taken 411159 times.
✓ Branch 19 taken 636750 times.
✓ Branch 20 taken 810083 times.
✓ Branch 21 taken 1062430 times.
✓ Branch 22 taken 5254845 times.
✓ Branch 23 taken 3403524 times.
✓ Branch 24 taken 645081 times.
✓ Branch 25 taken 660 times.
✓ Branch 26 taken 33 times.
✓ Branch 27 taken 68 times.
✓ Branch 28 taken 665 times.
✓ Branch 29 taken 8 times.
✓ Branch 30 taken 5 times.
✗ Branch 31 not taken.
✓ Branch 32 taken 5 times.
✗ Branch 33 not taken.
✓ Branch 34 taken 5 times.
✗ Branch 35 not taken.
✓ Branch 36 taken 5 times.
✗ Branch 37 not taken.
✓ Branch 38 taken 5 times.
✗ Branch 39 not taken.
✓ Branch 40 taken 5 times.
✗ Branch 41 not taken.
✓ Branch 42 taken 5 times.
✓ Branch 44 taken 4800 times.
✗ Branch 45 not taken.
✓ Branch 0 taken 19 times.
|
1769212855 | return scvfIndex_; |
189 | } | ||
190 | |||
191 | //! The local index of this sub control volume face | ||
192 | 95503390 | LocalIndexType localFaceIdx() const | |
193 | { | ||
194 |
5/9✗ Branch 0 not taken.
✓ Branch 1 taken 12 times.
✓ Branch 4 taken 41368532 times.
✓ Branch 5 taken 372400 times.
✓ Branch 7 taken 41339992 times.
✓ Branch 8 taken 1484000 times.
✗ Branch 2 not taken.
✗ Branch 6 not taken.
✗ Branch 9 not taken.
|
95503390 | return localFaceIdx_; |
195 | } | ||
196 | |||
197 | //! Returns the direction index of the facet (0 = x, 1 = y, 2 = z) | ||
198 | 382738 | unsigned int directionIndex() const | |
199 | { | ||
200 |
138/177✓ Branch 0 taken 8287818 times.
✓ Branch 1 taken 111913 times.
✓ Branch 2 taken 15 times.
✓ Branch 3 taken 8287808 times.
✓ Branch 4 taken 15 times.
✓ Branch 5 taken 5 times.
✓ Branch 6 taken 15 times.
✓ Branch 7 taken 5353 times.
✓ Branch 8 taken 15 times.
✓ Branch 9 taken 5 times.
✓ Branch 10 taken 15 times.
✓ Branch 11 taken 8214059 times.
✓ Branch 12 taken 704795 times.
✓ Branch 13 taken 5 times.
✓ Branch 14 taken 346787 times.
✓ Branch 15 taken 5 times.
✓ Branch 16 taken 15 times.
✓ Branch 17 taken 13506 times.
✓ Branch 18 taken 4989463 times.
✓ Branch 19 taken 261501 times.
✓ Branch 20 taken 5304127 times.
✓ Branch 21 taken 5 times.
✓ Branch 22 taken 243031 times.
✓ Branch 23 taken 3768845 times.
✓ Branch 24 taken 261511 times.
✓ Branch 25 taken 1587601 times.
✓ Branch 26 taken 15 times.
✓ Branch 27 taken 10597 times.
✓ Branch 28 taken 7691 times.
✓ Branch 29 taken 9146402 times.
✓ Branch 30 taken 486587 times.
✓ Branch 31 taken 4186708 times.
✓ Branch 32 taken 3146935 times.
✓ Branch 33 taken 6211879 times.
✓ Branch 34 taken 7376594 times.
✓ Branch 35 taken 2593382 times.
✓ Branch 36 taken 4815421 times.
✓ Branch 37 taken 3641171 times.
✓ Branch 38 taken 1036247 times.
✓ Branch 39 taken 2012639 times.
✓ Branch 40 taken 270879 times.
✓ Branch 41 taken 321049 times.
✓ Branch 42 taken 70693 times.
✓ Branch 43 taken 516271 times.
✓ Branch 44 taken 597079 times.
✓ Branch 45 taken 836873 times.
✓ Branch 46 taken 1124727 times.
✓ Branch 47 taken 5 times.
✓ Branch 49 taken 122612 times.
✗ Branch 50 not taken.
✗ Branch 52 not taken.
✗ Branch 53 not taken.
✗ Branch 55 not taken.
✗ Branch 56 not taken.
✗ Branch 58 not taken.
✗ Branch 59 not taken.
✗ Branch 61 not taken.
✗ Branch 62 not taken.
✗ Branch 64 not taken.
✗ Branch 65 not taken.
✗ Branch 67 not taken.
✗ Branch 68 not taken.
✗ Branch 70 not taken.
✗ Branch 71 not taken.
✗ Branch 73 not taken.
✗ Branch 74 not taken.
✗ Branch 76 not taken.
✗ Branch 77 not taken.
✗ Branch 79 not taken.
✗ Branch 80 not taken.
✗ Branch 82 not taken.
✗ Branch 83 not taken.
✓ Branch 84 taken 20 times.
✓ Branch 85 taken 80 times.
✓ Branch 86 taken 40 times.
✓ Branch 87 taken 160 times.
✓ Branch 89 taken 20 times.
✓ Branch 90 taken 80 times.
✓ Branch 91 taken 40 times.
✓ Branch 92 taken 160 times.
✓ Branch 94 taken 20 times.
✓ Branch 95 taken 80 times.
✓ Branch 96 taken 40 times.
✓ Branch 97 taken 160 times.
✓ Branch 99 taken 20 times.
✓ Branch 100 taken 80 times.
✓ Branch 101 taken 40 times.
✓ Branch 102 taken 160 times.
✓ Branch 104 taken 20 times.
✓ Branch 105 taken 80 times.
✓ Branch 106 taken 40 times.
✓ Branch 107 taken 160 times.
✓ Branch 109 taken 20 times.
✓ Branch 110 taken 80 times.
✓ Branch 111 taken 40 times.
✓ Branch 112 taken 160 times.
✓ Branch 114 taken 20 times.
✓ Branch 115 taken 80 times.
✓ Branch 116 taken 40 times.
✓ Branch 117 taken 160 times.
✓ Branch 119 taken 20 times.
✓ Branch 120 taken 80 times.
✓ Branch 121 taken 40 times.
✓ Branch 122 taken 160 times.
✓ Branch 124 taken 20 times.
✓ Branch 125 taken 80 times.
✓ Branch 126 taken 40 times.
✓ Branch 127 taken 160 times.
✗ Branch 130 not taken.
✗ Branch 131 not taken.
✓ Branch 132 taken 20 times.
✓ Branch 133 taken 80 times.
✓ Branch 134 taken 40 times.
✓ Branch 135 taken 160 times.
✗ Branch 138 not taken.
✗ Branch 139 not taken.
✓ Branch 140 taken 20 times.
✓ Branch 141 taken 80 times.
✓ Branch 142 taken 40 times.
✓ Branch 143 taken 160 times.
✗ Branch 146 not taken.
✗ Branch 147 not taken.
✓ Branch 148 taken 20 times.
✓ Branch 149 taken 80 times.
✓ Branch 150 taken 40 times.
✓ Branch 151 taken 160 times.
✗ Branch 154 not taken.
✗ Branch 155 not taken.
✓ Branch 156 taken 20 times.
✓ Branch 157 taken 80 times.
✓ Branch 158 taken 40 times.
✓ Branch 159 taken 160 times.
✓ Branch 161 taken 20 times.
✓ Branch 162 taken 80 times.
✓ Branch 163 taken 40 times.
✓ Branch 164 taken 160 times.
✓ Branch 166 taken 20 times.
✓ Branch 167 taken 80 times.
✓ Branch 168 taken 40 times.
✓ Branch 169 taken 160 times.
✓ Branch 171 taken 20 times.
✓ Branch 172 taken 80 times.
✓ Branch 173 taken 40 times.
✓ Branch 174 taken 160 times.
✓ Branch 176 taken 20 times.
✓ Branch 177 taken 80 times.
✓ Branch 178 taken 40 times.
✓ Branch 179 taken 160 times.
✗ Branch 182 not taken.
✗ Branch 183 not taken.
✓ Branch 184 taken 20 times.
✓ Branch 185 taken 80 times.
✓ Branch 186 taken 40 times.
✓ Branch 187 taken 160 times.
✗ Branch 190 not taken.
✗ Branch 191 not taken.
✓ Branch 192 taken 20 times.
✓ Branch 193 taken 80 times.
✓ Branch 194 taken 40 times.
✓ Branch 195 taken 160 times.
✗ Branch 198 not taken.
✗ Branch 199 not taken.
✓ Branch 200 taken 20 times.
✓ Branch 201 taken 80 times.
✓ Branch 202 taken 40 times.
✓ Branch 203 taken 160 times.
✗ Branch 206 not taken.
✗ Branch 207 not taken.
✓ Branch 208 taken 20 times.
✓ Branch 209 taken 80 times.
✓ Branch 210 taken 20 times.
✓ Branch 211 taken 80 times.
✓ Branch 212 taken 20 times.
✓ Branch 213 taken 80 times.
✓ Branch 214 taken 20 times.
✓ Branch 215 taken 80 times.
✓ Branch 48 taken 9516 times.
|
269115601 | return dirIdx_; |
201 | } | ||
202 | |||
203 | //! Returns whether the unitNormal of the face points in positive coordinate direction | ||
204 | 284898272 | bool normalInPosCoordDir() const | |
205 | { | ||
206 | 284898272 | return directionSign() > 0; | |
207 | } | ||
208 | |||
209 | //! Returns the sign of the unit outer normal's vector | ||
210 | 3298670525 | int directionSign() const | |
211 | { | ||
212 |
19/19✓ Branch 0 taken 142275652 times.
✓ Branch 1 taken 142622620 times.
✓ Branch 2 taken 27 times.
✓ Branch 3 taken 159857531 times.
✓ Branch 7 taken 71778 times.
✓ Branch 8 taken 119741642 times.
✓ Branch 10 taken 4340649 times.
✓ Branch 11 taken 127294218 times.
✓ Branch 16 taken 15200 times.
✓ Branch 17 taken 1152762 times.
✓ Branch 9 taken 3 times.
✓ Branch 4 taken 21 times.
✓ Branch 5 taken 35919094 times.
✓ Branch 13 taken 3019686 times.
✓ Branch 14 taken 47531653 times.
✓ Branch 19 taken 1835 times.
✓ Branch 20 taken 1024457 times.
✓ Branch 6 taken 44945749 times.
✓ Branch 12 taken 46705 times.
|
2592299151 | return outerNormalSign_; |
213 | } | ||
214 | |||
215 | //! Returns the data for one sub face | ||
216 | 3572575547 | const PairData& pairData(const int idx) const | |
217 | { | ||
218 |
14/31✗ Branch 0 not taken.
✓ Branch 1 taken 48 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 11 taken 346772 times.
✓ Branch 12 taken 280656308 times.
✓ Branch 15 taken 49 times.
✓ Branch 16 taken 285919953 times.
✓ Branch 19 taken 10882948 times.
✓ Branch 20 taken 203778872 times.
✓ Branch 22 taken 68962968 times.
✗ Branch 23 not taken.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✗ Branch 28 not taken.
✗ Branch 29 not taken.
✗ Branch 31 not taken.
✗ Branch 32 not taken.
✗ Branch 2 not taken.
✓ Branch 4 taken 48 times.
✓ Branch 7 taken 96 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 96 times.
✓ Branch 13 taken 3895192 times.
✓ Branch 17 taken 3895190 times.
✓ Branch 21 taken 7531900 times.
✗ Branch 24 not taken.
✗ Branch 27 not taken.
✗ Branch 30 not taken.
✗ Branch 33 not taken.
✗ Branch 34 not taken.
|
2153950111 | return pairData_[idx]; |
219 | } | ||
220 | |||
221 | //! Return an array of all pair data | ||
222 | const std::array<PairData, numPairs>& pairData() const | ||
223 | { | ||
224 |
56/56✓ Branch 20 taken 400 times.
✓ Branch 21 taken 200 times.
✓ Branch 22 taken 400 times.
✓ Branch 23 taken 200 times.
✓ Branch 24 taken 400 times.
✓ Branch 25 taken 200 times.
✓ Branch 26 taken 400 times.
✓ Branch 27 taken 200 times.
✓ Branch 28 taken 400 times.
✓ Branch 29 taken 200 times.
✓ Branch 30 taken 400 times.
✓ Branch 31 taken 200 times.
✓ Branch 32 taken 400 times.
✓ Branch 33 taken 200 times.
✓ Branch 34 taken 400 times.
✓ Branch 35 taken 200 times.
✓ Branch 36 taken 400 times.
✓ Branch 37 taken 200 times.
✓ Branch 38 taken 400 times.
✓ Branch 39 taken 200 times.
✓ Branch 40 taken 400 times.
✓ Branch 41 taken 200 times.
✓ Branch 42 taken 400 times.
✓ Branch 43 taken 200 times.
✓ Branch 44 taken 400 times.
✓ Branch 45 taken 200 times.
✓ Branch 46 taken 400 times.
✓ Branch 47 taken 200 times.
✓ Branch 48 taken 400 times.
✓ Branch 49 taken 200 times.
✓ Branch 50 taken 400 times.
✓ Branch 51 taken 200 times.
✓ Branch 52 taken 400 times.
✓ Branch 53 taken 200 times.
✓ Branch 54 taken 400 times.
✓ Branch 55 taken 200 times.
✓ Branch 56 taken 400 times.
✓ Branch 57 taken 200 times.
✓ Branch 58 taken 400 times.
✓ Branch 59 taken 200 times.
✓ Branch 60 taken 400 times.
✓ Branch 61 taken 200 times.
✓ Branch 62 taken 400 times.
✓ Branch 63 taken 200 times.
✓ Branch 64 taken 400 times.
✓ Branch 65 taken 200 times.
✓ Branch 66 taken 400 times.
✓ Branch 67 taken 200 times.
✓ Branch 0 taken 2809014 times.
✓ Branch 1 taken 142769390 times.
✓ Branch 2 taken 154610512 times.
✓ Branch 3 taken 77305244 times.
✓ Branch 4 taken 137428 times.
✓ Branch 5 taken 1822376 times.
✓ Branch 6 taken 48 times.
✓ Branch 7 taken 12 times.
|
381172720 | return pairData_; |
225 | } | ||
226 | |||
227 | //! Return an array of all pair data | ||
228 | const AxisData& axisData() const | ||
229 | { | ||
230 | return axisData_; | ||
231 | } | ||
232 | |||
233 | //! Returns true if the face is a ghost face | ||
234 | bool isGhostFace() const | ||
235 | { | ||
236 | return isGhostFace_; | ||
237 | } | ||
238 | |||
239 | //! Returns the length of the face in a certain direction (adaptation of area() for 3d) | ||
240 | 289679875 | Scalar faceLength(const int localSubFaceIdx) const | |
241 | { | ||
242 |
2/2✓ Branch 0 taken 24 times.
✓ Branch 1 taken 24 times.
|
48 | if (dim == 3 && localSubFaceIdx > 1) |
243 | 24 | return dimensions[1]; | |
244 | else | ||
245 | 289679851 | return dimensions[0]; | |
246 | } | ||
247 | |||
248 | /*! | ||
249 | * \brief Check if the face has a parallel neighbor | ||
250 | * | ||
251 | * \param localSubFaceIdx The local index of the subface | ||
252 | * \param parallelDegreeIdx The index describing how many faces away from the self face | ||
253 | */ | ||
254 |
48/48✓ Branch 0 taken 12564803 times.
✓ Branch 1 taken 253695864 times.
✓ Branch 2 taken 146930936 times.
✓ Branch 3 taken 8237885 times.
✓ Branch 4 taken 21417 times.
✓ Branch 5 taken 80 times.
✓ Branch 6 taken 1616424 times.
✓ Branch 7 taken 37857 times.
✓ Branch 8 taken 45555 times.
✓ Branch 9 taken 2171605 times.
✓ Branch 10 taken 6393656 times.
✓ Branch 11 taken 190744 times.
✓ Branch 12 taken 320 times.
✓ Branch 13 taken 80 times.
✓ Branch 14 taken 320 times.
✓ Branch 15 taken 80 times.
✓ Branch 16 taken 320 times.
✓ Branch 17 taken 80 times.
✓ Branch 18 taken 320 times.
✓ Branch 19 taken 80 times.
✓ Branch 20 taken 320 times.
✓ Branch 21 taken 80 times.
✓ Branch 22 taken 320 times.
✓ Branch 23 taken 80 times.
✓ Branch 24 taken 320 times.
✓ Branch 25 taken 80 times.
✓ Branch 26 taken 320 times.
✓ Branch 27 taken 80 times.
✓ Branch 28 taken 320 times.
✓ Branch 29 taken 80 times.
✓ Branch 30 taken 320 times.
✓ Branch 31 taken 80 times.
✓ Branch 32 taken 320 times.
✓ Branch 33 taken 80 times.
✓ Branch 34 taken 320 times.
✓ Branch 35 taken 80 times.
✓ Branch 36 taken 320 times.
✓ Branch 37 taken 80 times.
✓ Branch 38 taken 320 times.
✓ Branch 39 taken 80 times.
✓ Branch 40 taken 320 times.
✓ Branch 41 taken 80 times.
✓ Branch 42 taken 320 times.
✓ Branch 43 taken 80 times.
✓ Branch 44 taken 320 times.
✓ Branch 45 taken 80 times.
✓ Branch 46 taken 320 times.
✓ Branch 47 taken 80 times.
|
431914026 | bool hasParallelNeighbor(const int localSubFaceIdx, const int parallelDegreeIdx) const |
255 | { | ||
256 |
48/48✓ Branch 0 taken 12564803 times.
✓ Branch 1 taken 253695864 times.
✓ Branch 2 taken 143123173 times.
✓ Branch 3 taken 9874123 times.
✓ Branch 4 taken 3829180 times.
✓ Branch 5 taken 21177 times.
✓ Branch 6 taken 21417 times.
✓ Branch 7 taken 80 times.
✓ Branch 8 taken 1616424 times.
✓ Branch 9 taken 37857 times.
✓ Branch 10 taken 2116073 times.
✓ Branch 11 taken 53958 times.
✓ Branch 12 taken 45555 times.
✓ Branch 13 taken 2171605 times.
✓ Branch 14 taken 6393656 times.
✓ Branch 15 taken 190744 times.
✓ Branch 16 taken 320 times.
✓ Branch 17 taken 80 times.
✓ Branch 18 taken 320 times.
✓ Branch 19 taken 80 times.
✓ Branch 20 taken 320 times.
✓ Branch 21 taken 80 times.
✓ Branch 22 taken 320 times.
✓ Branch 23 taken 80 times.
✓ Branch 24 taken 320 times.
✓ Branch 25 taken 80 times.
✓ Branch 26 taken 320 times.
✓ Branch 27 taken 80 times.
✓ Branch 28 taken 320 times.
✓ Branch 29 taken 80 times.
✓ Branch 30 taken 320 times.
✓ Branch 31 taken 80 times.
✓ Branch 32 taken 320 times.
✓ Branch 33 taken 80 times.
✓ Branch 34 taken 320 times.
✓ Branch 35 taken 80 times.
✓ Branch 36 taken 320 times.
✓ Branch 37 taken 80 times.
✓ Branch 38 taken 320 times.
✓ Branch 39 taken 80 times.
✓ Branch 40 taken 320 times.
✓ Branch 41 taken 80 times.
✓ Branch 42 taken 320 times.
✓ Branch 43 taken 80 times.
✓ Branch 44 taken 320 times.
✓ Branch 45 taken 80 times.
✓ Branch 46 taken 320 times.
✓ Branch 47 taken 80 times.
|
435762089 | return pairData(localSubFaceIdx).hasParallelNeighbor[parallelDegreeIdx]; |
257 | } | ||
258 | |||
259 | /*! | ||
260 | * \brief Check if the face has a half parallel neighbor | ||
261 | * | ||
262 | * \param localSubFaceIdx The local index of the subface | ||
263 | * | ||
264 | * ------------ | ||
265 | * | | | ||
266 | * | | | ||
267 | * | | | ||
268 | * ----------------------- | ||
269 | * | yyyyyyyy s | | ||
270 | * | yyyyyyyy s | | ||
271 | * | yyyyyyyy s | | ||
272 | * ----------------------- | ||
273 | * In this corner geometry, hasParallelNeighbor will return true for subcontrolvolumeface s belonging to the | ||
274 | * element filled by 'y's, but hasParallelNeighbor will return false for the subcontrolvolumeface that has the | ||
275 | * same dofIndex. We name this situation hasHalfParallelNeighbor. | ||
276 | */ | ||
277 |
7/10✓ Branch 0 taken 10696 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 10696 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 266260267 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 3823512 times.
✓ Branch 7 taken 5348 times.
✓ Branch 8 taken 21097 times.
✓ Branch 9 taken 5348 times.
|
270131616 | bool hasHalfParallelNeighbor(const int localSubFaceIdx) const |
278 | { | ||
279 |
7/10✓ Branch 0 taken 10696 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 10696 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 266260267 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 3823512 times.
✓ Branch 7 taken 5348 times.
✓ Branch 8 taken 21097 times.
✓ Branch 9 taken 5348 times.
|
270131616 | return pairData(localSubFaceIdx).hasHalfParallelNeighbor; |
280 | } | ||
281 | |||
282 | /*! | ||
283 | * \brief Check if the face has a corner parallel neighbor | ||
284 | * | ||
285 | * \param localSubFaceIdx The local index of the subface | ||
286 | * | ||
287 | * ------------ | ||
288 | * | yyyyyyyy s | ||
289 | * | yyyyyyyy s | ||
290 | * | yyyyyyyy s | ||
291 | * ----------------------- | ||
292 | * | | | | ||
293 | * | | | | ||
294 | * | | | | ||
295 | * ----------------------- | ||
296 | * In this corner geometry, hasParallelNeighbor will return true for subcontrolvolumeface s belonging to the | ||
297 | * element filled by 'y's. However, as there also might be a boundary velocity value known at the corner, which | ||
298 | * can be used instead of the standard parallel velocity in some cases, we want to identify this situation. We | ||
299 | * name it cornerParallelNeighbor. | ||
300 | */ | ||
301 |
3/10✗ Branch 0 not taken.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 266260267 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 3823512 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 21097 times.
|
270104876 | bool hasCornerParallelNeighbor(const int localSubFaceIdx) const |
302 | { | ||
303 |
3/10✗ Branch 0 not taken.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 266260267 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 3823512 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 21097 times.
|
270104876 | return pairData(localSubFaceIdx).hasCornerParallelNeighbor; |
304 | } | ||
305 | |||
306 | /*! | ||
307 | * \brief Check if the face has an outer normal neighbor | ||
308 | * | ||
309 | * \param localSubFaceIdx The local index of the subface | ||
310 | */ | ||
311 |
48/48✓ Branch 0 taken 147424376 times.
✓ Branch 1 taken 7186488 times.
✓ Branch 2 taken 320 times.
✓ Branch 3 taken 80 times.
✓ Branch 4 taken 320 times.
✓ Branch 5 taken 80 times.
✓ Branch 6 taken 320 times.
✓ Branch 7 taken 80 times.
✓ Branch 8 taken 320 times.
✓ Branch 9 taken 80 times.
✓ Branch 10 taken 320 times.
✓ Branch 11 taken 80 times.
✓ Branch 12 taken 320 times.
✓ Branch 13 taken 80 times.
✓ Branch 14 taken 320 times.
✓ Branch 15 taken 80 times.
✓ Branch 16 taken 320 times.
✓ Branch 17 taken 80 times.
✓ Branch 18 taken 320 times.
✓ Branch 19 taken 80 times.
✓ Branch 20 taken 320 times.
✓ Branch 21 taken 80 times.
✓ Branch 22 taken 320 times.
✓ Branch 23 taken 80 times.
✓ Branch 24 taken 320 times.
✓ Branch 25 taken 80 times.
✓ Branch 26 taken 320 times.
✓ Branch 27 taken 80 times.
✓ Branch 28 taken 320 times.
✓ Branch 29 taken 80 times.
✓ Branch 30 taken 320 times.
✓ Branch 31 taken 80 times.
✓ Branch 32 taken 320 times.
✓ Branch 33 taken 80 times.
✓ Branch 34 taken 320 times.
✓ Branch 35 taken 80 times.
✓ Branch 36 taken 320 times.
✓ Branch 37 taken 80 times.
✓ Branch 38 taken 320 times.
✓ Branch 39 taken 80 times.
✓ Branch 40 taken 320 times.
✓ Branch 41 taken 80 times.
✓ Branch 42 taken 320 times.
✓ Branch 43 taken 80 times.
✓ Branch 44 taken 320 times.
✓ Branch 45 taken 80 times.
✓ Branch 46 taken 320 times.
✓ Branch 47 taken 80 times.
|
154620064 | bool hasOuterLateral(const int localSubFaceIdx) const |
312 | { | ||
313 |
48/48✓ Branch 0 taken 147424376 times.
✓ Branch 1 taken 7186488 times.
✓ Branch 2 taken 320 times.
✓ Branch 3 taken 80 times.
✓ Branch 4 taken 320 times.
✓ Branch 5 taken 80 times.
✓ Branch 6 taken 320 times.
✓ Branch 7 taken 80 times.
✓ Branch 8 taken 320 times.
✓ Branch 9 taken 80 times.
✓ Branch 10 taken 320 times.
✓ Branch 11 taken 80 times.
✓ Branch 12 taken 320 times.
✓ Branch 13 taken 80 times.
✓ Branch 14 taken 320 times.
✓ Branch 15 taken 80 times.
✓ Branch 16 taken 320 times.
✓ Branch 17 taken 80 times.
✓ Branch 18 taken 320 times.
✓ Branch 19 taken 80 times.
✓ Branch 20 taken 320 times.
✓ Branch 21 taken 80 times.
✓ Branch 22 taken 320 times.
✓ Branch 23 taken 80 times.
✓ Branch 24 taken 320 times.
✓ Branch 25 taken 80 times.
✓ Branch 26 taken 320 times.
✓ Branch 27 taken 80 times.
✓ Branch 28 taken 320 times.
✓ Branch 29 taken 80 times.
✓ Branch 30 taken 320 times.
✓ Branch 31 taken 80 times.
✓ Branch 32 taken 320 times.
✓ Branch 33 taken 80 times.
✓ Branch 34 taken 320 times.
✓ Branch 35 taken 80 times.
✓ Branch 36 taken 320 times.
✓ Branch 37 taken 80 times.
✓ Branch 38 taken 320 times.
✓ Branch 39 taken 80 times.
✓ Branch 40 taken 320 times.
✓ Branch 41 taken 80 times.
✓ Branch 42 taken 320 times.
✓ Branch 43 taken 80 times.
✓ Branch 44 taken 320 times.
✓ Branch 45 taken 80 times.
✓ Branch 46 taken 320 times.
✓ Branch 47 taken 80 times.
|
154620064 | return pairData(localSubFaceIdx).hasOuterLateral; |
314 | } | ||
315 | |||
316 | /*! | ||
317 | * \brief Check if the face has a backward neighbor | ||
318 | * | ||
319 | * \param backwardIdx The index describing how many faces backward this dof is from the opposite face | ||
320 | */ | ||
321 | template<bool enable = useHigherOrder, std::enable_if_t<enable, int> = 0> | ||
322 | 2521770 | bool hasBackwardNeighbor(const int backwardIdx) const | |
323 | { | ||
324 |
4/4✓ Branch 0 taken 18200 times.
✓ Branch 1 taken 550 times.
✓ Branch 2 taken 1612734 times.
✓ Branch 3 taken 33266 times.
|
2521770 | return axisData().hasBackwardNeighbor[backwardIdx]; |
325 | } | ||
326 | |||
327 | /*! | ||
328 | * \brief Check if the face has a forward neighbor | ||
329 | * | ||
330 | * \param forwardIdx The index describing how many faces forward this dof is of the self face | ||
331 | */ | ||
332 | template<bool enable = useHigherOrder, std::enable_if_t<enable, int> = 0> | ||
333 | 2755326 | bool hasForwardNeighbor(const int forwardIdx) const | |
334 | { | ||
335 |
4/4✓ Branch 0 taken 550 times.
✓ Branch 1 taken 18200 times.
✓ Branch 2 taken 1618334 times.
✓ Branch 3 taken 27666 times.
|
2755326 | return axisData().hasForwardNeighbor[forwardIdx]; |
336 | } | ||
337 | |||
338 | //! Returns the dof of the face | ||
339 | 141724051 | GridIndexType dofIndex() const | |
340 | { | ||
341 |
72/81✓ Branch 1 taken 3822668 times.
✓ Branch 2 taken 53412 times.
✓ Branch 28 taken 4996 times.
✗ Branch 29 not taken.
✓ Branch 31 taken 4996 times.
✗ Branch 32 not taken.
✓ Branch 33 taken 18928 times.
✓ Branch 34 taken 22180 times.
✓ Branch 35 taken 18368 times.
✓ Branch 36 taken 22180 times.
✓ Branch 37 taken 40 times.
✓ Branch 38 taken 60 times.
✓ Branch 39 taken 40 times.
✓ Branch 40 taken 60 times.
✓ Branch 41 taken 40 times.
✓ Branch 42 taken 60 times.
✓ Branch 43 taken 40 times.
✓ Branch 44 taken 60 times.
✓ Branch 45 taken 40 times.
✓ Branch 46 taken 60 times.
✓ Branch 47 taken 40 times.
✓ Branch 48 taken 60 times.
✓ Branch 49 taken 40 times.
✓ Branch 50 taken 60 times.
✓ Branch 51 taken 40 times.
✓ Branch 52 taken 60 times.
✓ Branch 53 taken 40 times.
✓ Branch 54 taken 60 times.
✓ Branch 55 taken 40 times.
✓ Branch 56 taken 60 times.
✓ Branch 57 taken 40 times.
✓ Branch 58 taken 60 times.
✓ Branch 59 taken 40 times.
✓ Branch 60 taken 60 times.
✓ Branch 61 taken 40 times.
✓ Branch 62 taken 60 times.
✓ Branch 63 taken 40 times.
✓ Branch 64 taken 60 times.
✓ Branch 65 taken 40 times.
✓ Branch 66 taken 60 times.
✓ Branch 67 taken 40 times.
✓ Branch 68 taken 60 times.
✓ Branch 69 taken 40 times.
✓ Branch 70 taken 60 times.
✓ Branch 71 taken 40 times.
✓ Branch 72 taken 60 times.
✓ Branch 73 taken 40 times.
✓ Branch 74 taken 60 times.
✓ Branch 75 taken 40 times.
✓ Branch 76 taken 60 times.
✓ Branch 77 taken 40 times.
✓ Branch 78 taken 60 times.
✓ Branch 79 taken 40 times.
✓ Branch 80 taken 60 times.
✓ Branch 4 taken 3790128 times.
✓ Branch 5 taken 156160 times.
✓ Branch 9 taken 1570144 times.
✗ Branch 10 not taken.
✓ Branch 13 taken 8678688 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 8961888 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 53635 times.
✓ Branch 19 taken 58765 times.
✓ Branch 20 taken 11362 times.
✓ Branch 21 taken 277438 times.
✓ Branch 22 taken 1453568 times.
✗ Branch 23 not taken.
✗ Branch 3 not taken.
✓ Branch 6 taken 1496000 times.
✗ Branch 7 not taken.
✓ Branch 12 taken 283200 times.
✓ Branch 25 taken 1453568 times.
✗ Branch 26 not taken.
✓ Branch 0 taken 4880 times.
✓ Branch 8 taken 156160 times.
✓ Branch 11 taken 1024 times.
✓ Branch 15 taken 1024 times.
✓ Branch 24 taken 2419200 times.
✓ Branch 27 taken 2419200 times.
✓ Branch 30 taken 560 times.
|
123725507 | return axisData().selfDof; |
342 | } | ||
343 | |||
344 | //! Returns the dof of the opposing face | ||
345 | 77310044 | GridIndexType dofIndexOpposingFace() const | |
346 | { | ||
347 |
1/2✓ Branch 1 taken 12 times.
✗ Branch 2 not taken.
|
77310044 | return axisData().oppositeDof; |
348 | } | ||
349 | |||
350 | //! Returns the dof the first forward face | ||
351 | GridIndexType dofIndexForwardFace() const | ||
352 | { | ||
353 | return axisData().inAxisForwardDofs[0]; | ||
354 | } | ||
355 | |||
356 | //! Returns the dof of the first backward face | ||
357 | GridIndexType dofIndexBackwardFace() const | ||
358 | { | ||
359 | return axisData().inAxisBackwardDofs[0]; | ||
360 | } | ||
361 | |||
362 | //! Returns the distance between the face and the opposite one | ||
363 | 147503500 | Scalar selfToOppositeDistance() const | |
364 | { | ||
365 |
2/3✓ Branch 0 taken 51 times.
✓ Branch 1 taken 145578305 times.
✗ Branch 2 not taken.
|
147503500 | return axisData().selfToOppositeDistance; |
366 | } | ||
367 | |||
368 | /*! | ||
369 | * \brief Returns the distance between the parallel dofs | ||
370 | * | ||
371 | * \param localSubFaceIdx The local index of the subface | ||
372 | * \param parallelDegreeIdx The index describing how many faces away from the self | ||
373 | */ | ||
374 | 289679923 | Scalar parallelDofsDistance(const int localSubFaceIdx, const int parallelDegreeIdx) const | |
375 | { | ||
376 |
2/2✓ Branch 0 taken 48 times.
✓ Branch 1 taken 48 times.
|
96 | if (parallelDegreeIdx == 0) |
377 |
2/2✓ Branch 0 taken 21097 times.
✓ Branch 1 taken 1657335 times.
|
285784683 | return (faceLength(localSubFaceIdx) + pairData(localSubFaceIdx).parallelCellWidths[0]) * 0.5; |
378 | // pairData(localSubFaceIdx).parallelCellWidths[0]) will return 0.0 if the subface perpendicular the scvf lies on a boundary | ||
379 | else | ||
380 | { | ||
381 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 48 times.
|
48 | assert((parallelDegreeIdx == 1) && "Only the width of the first two parallel cells (indices 0 and 1) is stored for each scvf."); |
382 | 2171573 | return (pairData(localSubFaceIdx).parallelCellWidths[0] + pairData(localSubFaceIdx).parallelCellWidths[1]) * 0.5; | |
383 | } | ||
384 | } | ||
385 | |||
386 | //! set the center to a different position | ||
387 | 17557066 | void setCenter(const GlobalPosition& center) | |
388 | 17557066 | { center_ = center; } | |
389 | |||
390 | //! set the boundary flag | ||
391 | 17557066 | void setBoundary(bool boundaryFlag) | |
392 | 17557066 | { boundary_ = boundaryFlag; } | |
393 | |||
394 | //! set the ghost face flag | ||
395 | 17557066 | void setIsGhostFace(bool isGhostFaceFlag) | |
396 |
5/21✓ Branch 1 taken 120800 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 11928 times.
✓ Branch 5 taken 2940 times.
✓ Branch 7 taken 475916 times.
✗ Branch 8 not taken.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 16 not taken.
✗ Branch 17 not taken.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
✗ Branch 0 not taken.
✗ Branch 3 not taken.
✓ Branch 6 taken 508848 times.
✗ Branch 9 not taken.
✗ Branch 12 not taken.
|
17557066 | { isGhostFace_ = isGhostFaceFlag; } |
397 | |||
398 | private: | ||
399 | std::array<Scalar, dim-1> dimensions; | ||
400 | Scalar area_; | ||
401 | GlobalPosition center_; | ||
402 | GlobalPosition unitOuterNormal_; | ||
403 | GridIndexType scvfIndex_; | ||
404 | std::vector<GridIndexType> scvIndices_; | ||
405 | bool boundary_; | ||
406 | |||
407 | Scalar selfToOppositeDistance_; | ||
408 | AxisData axisData_; | ||
409 | std::array<PairData, numPairs> pairData_; | ||
410 | |||
411 | int localFaceIdx_; | ||
412 | unsigned int dirIdx_; | ||
413 | int outerNormalSign_; | ||
414 | bool isGhostFace_; | ||
415 | }; | ||
416 | |||
417 | } // end namespace Dumux | ||
418 | |||
419 | #endif | ||
420 |