GCC Code Coverage Report

Directory: gridformat/
File: gridformat/grid/grid.hpp
Date: 2026-10-18 00:09:19
Exec Total Coverage
Lines: 52 53 98.1%
Functions: 1079 1139 94.7%
Branches: 32 47 68.1%
Line Branch Exec Source
1 // SPDX-FileCopyrightText: 2022-2023 Dennis Gläser <dennis.glaeser@iws.uni-stuttgart.de>
2 // SPDX-License-Identifier: MIT
3 /*!
4 * \file
5 * \ingroup Grid
6 * \brief Defines the grid interface
7 */
8 #ifndef GRIDFORMAT_GRID_GRID_HPP_
9 #define GRIDFORMAT_GRID_GRID_HPP_
10
11 #include <ranges>
12 #include <concepts>
13 #include <cassert>
14 #include <numeric>
15 #include <algorithm>
16 #include <unordered_map>
17
18 #include <gridformat/common/ranges.hpp>
19 #include <gridformat/common/concepts.hpp>
20 #include <gridformat/common/flat_index_mapper.hpp>
21
22 #include <gridformat/grid/cell_type.hpp>
23 #include <gridformat/grid/traits.hpp>
24 #include <gridformat/grid/type_traits.hpp>
25 #include <gridformat/grid/concepts.hpp>
26
27 namespace GridFormat {
28
29 #ifndef DOXYGEN
30 namespace GridDetail {
31
32 template<typename T, std::size_t dim>
33 19280 std::array<std::array<T, dim>, dim> standard_basis() {
34 std::array<std::array<T, dim>, dim> result;
35
2/3
✓ Branch 1 taken 102 times.
✓ Branch 2 taken 922 times.
✗ Branch 3 not taken.
 46329 std::ranges::for_each(result, [&] (auto& b) { std::ranges::fill(b, T{0.0}); });
36
2/2
✓ Branch 0 taken 27049 times.
✓ Branch 1 taken 9864 times.
 72456 for (std::size_t i = 0; i < dim; ++i)
37 53176 result[i][i] = 1.0;
38 19280 return result;
39 }
40
41 } // namespace GridDetail
42 #endif // DOXYGEN
43
44 //! \addtogroup Grid
45 //! \{
46
47 template<GridDetail::ExposesPointRange Grid>
48 1118813 std::ranges::range decltype(auto) points(const Grid& grid) {
49 1118813 return Traits::Points<Grid>::get(grid);
50 }
51
52 template<GridDetail::ExposesCellRange Grid>
53 1079455 std::ranges::range decltype(auto) cells(const Grid& grid) {
54 1079455 return Traits::Cells<Grid>::get(grid);
55 }
56
57 template<GridDetail::ExposesPointId Grid>
58 515212 std::integral auto id(const Grid& grid, const Point<Grid>& point) {
59 515212 return Traits::PointId<Grid, Point<Grid>>::get(grid, point);
60 }
61
62 template<GridDetail::ExposesPointCoordinates Grid>
63 2742584 std::ranges::range decltype(auto) coordinates(const Grid& grid, const Point<Grid>& point) {
64 2742584 return Traits::PointCoordinates<Grid, Point<Grid>>::get(grid, point);
65 }
66
67 template<GridDetail::ExposesCellPoints Grid>
68 651963 std::ranges::range decltype(auto) points(const Grid& grid, const Cell<Grid>& cell) {
69 651963 return Traits::CellPoints<Grid, Cell<Grid>>::get(grid, cell);
70 }
71
72 template<GridDetail::ExposesCellType Grid>
73 310372 CellType type(const Grid& grid, const Cell<Grid>& cell) {
74 310372 return Traits::CellType<Grid, Cell<Grid>>::get(grid, cell);
75 }
76
77 template<typename Grid>
78 requires(GridDetail::ExposesNumberOfCells<Grid> or
79 GridDetail::ExposesCellRange<Grid>)
80 956589 std::size_t number_of_cells(const Grid& grid) {
81 if constexpr (GridDetail::ExposesNumberOfCells<Grid>)
82 6118 return Traits::NumberOfCells<Grid>::get(grid);
83 else
84
2/4
✓ Branch 1 taken 2224 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2224 times.
✗ Branch 5 not taken.
 950471 return Ranges::size(cells(grid));
85 }
86
87 template<typename Grid>
88 requires(GridDetail::ExposesNumberOfPoints<Grid> or
89 GridDetail::ExposesPointRange<Grid>)
90 984886 std::size_t number_of_points(const Grid& grid) {
91 if constexpr (GridDetail::ExposesNumberOfPoints<Grid>)
92 6040 return Traits::NumberOfPoints<Grid>::get(grid);
93 else
94
2/4
✓ Branch 1 taken 3219 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 3219 times.
✗ Branch 5 not taken.
 978846 return Ranges::size(points(grid));
95 }
96
97 template<typename Grid>
98 requires(GridDetail::ExposesNumberOfCellPoints<Grid> or
99 GridDetail::ExposesCellPoints<Grid>)
100 376172 std::size_t number_of_points(const Grid& grid, const Cell<Grid>& cell) {
101 if constexpr (GridDetail::ExposesNumberOfCellPoints<Grid>)
102 74680 return Traits::NumberOfCellPoints<Grid, Cell<Grid>>::get(grid, cell);
103 else
104
2/4
✓ Branch 1 taken 151599 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 20648 times.
✗ Branch 5 not taken.
 301492 return Ranges::size(points(grid, cell));
105 }
106
107 template<GridDetail::ExposesExtents Grid>
108 190043 Concepts::StaticallySizedRange decltype(auto) extents(const Grid& grid) {
109 190043 return Traits::Extents<Grid>::get(grid);
110 }
111
112 template<GridDetail::ExposesExtents Grid>
113 82373 auto point_extents(const Grid& grid) {
114
1/2
✓ Branch 1 taken 152 times.
✗ Branch 2 not taken.
 82373 std::ranges::range auto result = extents(grid);
115
1/2
✓ Branch 1 taken 10349 times.
✗ Branch 2 not taken.
 194462 std::ranges::for_each(result, [] (std::integral auto& ext) { ext++; });
116 82373 return result;
117 }
118
119 template<GridDetail::ExposesOrigin Grid>
120 5385 Concepts::StaticallySizedRange decltype(auto) origin(const Grid& grid) {
121 5385 return Traits::Origin<Grid>::get(grid);
122 }
123
124 template<GridDetail::ExposesOrdinates Grid>
125 24720 Concepts::MDRange<1> decltype(auto) ordinates(const Grid& grid, unsigned dir) {
126 24720 return Traits::Ordinates<Grid>::get(grid, dir);
127 }
128
129 template<GridDetail::ExposesSpacing Grid>
130 10964 Concepts::StaticallySizedRange decltype(auto) spacing(const Grid& grid) {
131 10964 return Traits::Spacing<Grid>::get(grid);
132 }
133
134 template<Concepts::StructuredEntitySet Grid>
135 11164 Concepts::StaticallySizedMDRange<2> decltype(auto) basis(const Grid& grid) {
136 static constexpr std::size_t dim = dimension<Grid>;
137 if constexpr (GridDetail::ExposesBasis<Grid>) {
138 using ResultType = std::remove_cvref_t<decltype(Traits::Basis<Grid>::get(grid))>;
139 static_assert(static_size<ResultType> == dim);
140 static_assert(static_size<std::ranges::range_value_t<ResultType>> == dim);
141 124 return Traits::Basis<Grid>::get(grid);
142 }
143 else {
144 11040 return GridDetail::standard_basis<CoordinateType<Grid>, dim>();
145 }
146 }
147
148 template<GridDetail::ExposesCellLocation Grid>
149 1029526 Concepts::StaticallySizedRange decltype(auto) location(const Grid& grid, const Cell<Grid>& c) {
150 1029526 return Traits::Location<Grid, Cell<Grid>>::get(grid, c);
151 }
152
153 template<GridDetail::ExposesPointLocation Grid>
154 requires(not std::same_as<Cell<Grid>, Point<Grid>>) // avoid ambiguity if cell & point type are the same
155 1638133 Concepts::StaticallySizedRange decltype(auto) location(const Grid& grid, const Point<Grid>& p) {
156 1638133 return Traits::Location<Grid, Point<Grid>>::get(grid, p);
157 }
158
159 template<GridDetail::ExposesPointRange Grid> requires(GridDetail::ExposesPointId<Grid>)
160 1698 std::unordered_map<std::size_t, std::size_t> make_point_id_map(const Grid& grid) {
161 1698 std::size_t i = 0;
162 1698 std::unordered_map<std::size_t, std::size_t> point_id_to_running_idx;
163
13/15
✓ Branch 1 taken 308 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 308 times.
✓ Branch 5 taken 11215 times.
✓ Branch 6 taken 700 times.
✓ Branch 7 taken 271 times.
✓ Branch 8 taken 4380 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 23301 times.
✓ Branch 11 taken 18222 times.
✓ Branch 12 taken 126 times.
✓ Branch 13 taken 23411 times.
✓ Branch 14 taken 72 times.
✓ Branch 15 taken 13051 times.
✓ Branch 16 taken 110 times.
 97710 for (const auto& p : points(grid)) {
164
3/4
✓ Branch 1 taken 18970 times.
✓ Branch 2 taken 3204 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 3593 times.
 95966 assert(id(grid, p) >= 0);
165
4/7
✓ Branch 1 taken 18970 times.
✓ Branch 2 taken 33768 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 18970 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 17431 times.
✗ Branch 8 not taken.
 95966 point_id_to_running_idx[id(grid, p)] = i++;
166 }
167 1698 return point_id_to_running_idx;
168 ✗ }
169
170 //! \} group Grid
171
172 } // namespace GridFormat
173
174 #endif // GRIDFORMAT_GRID_GRID_HPP_
175