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 |