GCC Code Coverage Report


Directory: gridformat/
File: gridformat/vtk/pvtu_writer.hpp
Date: 2026-10-12 17:50:57
Exec Total Coverage
Lines: 51 53 96.2%
Functions: 115 217 53.0%
Branches: 128 428 29.9%

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 VTK
6 * \copydoc GridFormat::PVTUWriter
7 */
8 #ifndef GRIDFORMAT_VTK_PVTU_WRITER_HPP_
9 #define GRIDFORMAT_VTK_PVTU_WRITER_HPP_
10
11 #include <ostream>
12 #include <string>
13 #include <filesystem>
14
15 #include <gridformat/common/exceptions.hpp>
16 #include <gridformat/common/lvalue_reference.hpp>
17 #include <gridformat/common/output_file.hpp>
18 #include <gridformat/parallel/communication.hpp>
19 #include <gridformat/parallel/helpers.hpp>
20
21 #include <gridformat/grid/grid.hpp>
22 #include <gridformat/xml/element.hpp>
23 #include <gridformat/vtk/vtu_writer.hpp>
24 #include <gridformat/vtk/parallel.hpp>
25
26 namespace GridFormat {
27
28 /*!
29 * \ingroup VTK
30 * \brief Writer for parallel .pvtu files
31 */
32 template<Concepts::UnstructuredGrid Grid,
33 Concepts::Communicator Communicator>
34 class PVTUWriter : public VTK::XMLWriterBase<Grid, PVTUWriter<Grid, Communicator>> {
35 using ParentType = VTK::XMLWriterBase<Grid, PVTUWriter<Grid, Communicator>>;
36
37 public:
38 300 explicit PVTUWriter(LValueReferenceOf<const Grid> grid,
39 Communicator comm,
40 VTK::XMLOptions xml_opts = {})
41 300 : ParentType(grid.get(), ".pvtu", false, xml_opts)
42
2/4
✓ Branch 1 taken 213 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 213 times.
✗ Branch 5 not taken.
900 , _comm(comm)
43 300 {}
44
45 180 const Communicator& communicator() const {
46 180 return _comm;
47 }
48
49 private:
50 Communicator _comm;
51
52 70 PVTUWriter _with(VTK::XMLOptions xml_opts) const override {
53
1/2
✓ Branch 4 taken 70 times.
✗ Branch 5 not taken.
70 return PVTUWriter{this->grid(), _comm, std::move(xml_opts)};
54 }
55
56 ✗ void _write(std::ostream&) const override {
57 ✗ throw InvalidState(
58 "PVTUWriter does not support direct export into stream. "
59 "Use overload with filename instead!"
60 );
61 }
62
63 423 virtual void _write(const std::string& filename_with_ext) const override {
64 // the collective invocations also ensure that all pieces are written before the main file,
65 // and that the main file is written before returning
66
2/2
✓ Branch 1 taken 258 times.
✓ Branch 2 taken 6 times.
741 Parallel::invoke_collectively(_comm, [&] () {
67 264 _write_piece(filename_with_ext);
68
2/4
✓ Branch 1 taken 264 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 264 times.
✗ Branch 5 not taken.
858 }, "Writing the pieces of '" + filename_with_ext + "'");
69
2/2
✓ Branch 1 taken 255 times.
✓ Branch 2 taken 3 times.
735 Parallel::invoke_collectively(_comm, [&] () {
70
8/8
✓ Branch 1 taken 96 times.
✓ Branch 2 taken 89 times.
✓ Branch 4 taken 32 times.
✓ Branch 5 taken 37 times.
✓ Branch 7 taken 1 times.
✓ Branch 8 taken 1 times.
✓ Branch 10 taken 1 times.
✓ Branch 11 taken 1 times.
258 if (Parallel::rank(_comm) == 0)
71 130 _write_pvtu_file(filename_with_ext);
72
2/4
✓ Branch 1 taken 258 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 258 times.
✗ Branch 5 not taken.
840 }, "Writing '" + filename_with_ext + "'");
73 414 }
74
75 423 void _write_piece(const std::string& par_filename) const {
76
1/2
✓ Branch 3 taken 264 times.
✗ Branch 4 not taken.
423 VTUWriter writer{this->grid(), this->_xml_opts};
77
1/2
✓ Branch 1 taken 264 times.
✗ Branch 2 not taken.
423 this->copy_fields(writer);
78
6/8
✓ Branch 1 taken 263 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 263 times.
✓ Branch 5 taken 1 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 258 times.
✓ Branch 8 taken 5 times.
428 writer.write(PVTK::piece_basefilename(par_filename, Parallel::rank(_comm)));
79 423 }
80
81 211 void _write_pvtu_file(const std::string& filename_with_ext) const {
82
1/2
✓ Branch 1 taken 130 times.
✗ Branch 2 not taken.
422 XMLElement pvtk_xml("VTKFile");
83
2/4
✓ Branch 1 taken 130 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 130 times.
✗ Branch 5 not taken.
422 pvtk_xml.set_attribute("type", "PUnstructuredGrid");
84
85
2/4
✓ Branch 1 taken 130 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 130 times.
✗ Branch 5 not taken.
422 XMLElement& grid = pvtk_xml.add_child("PUnstructuredGrid");
86
2/4
✓ Branch 1 taken 130 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 130 times.
✗ Branch 5 not taken.
422 XMLElement& ppoint_data = grid.add_child("PPointData");
87
2/4
✓ Branch 1 taken 130 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 130 times.
✗ Branch 5 not taken.
211 XMLElement& pcell_data = grid.add_child("PCellData");
88 780 std::visit([&] (const auto& encoder) {
89 488 std::visit([&] (const auto& data_format) {
90 130 PVTK::PDataArrayHelper pdata_helper{encoder, data_format, ppoint_data};
91 130 PVTK::PDataArrayHelper cdata_helper{encoder, data_format, pcell_data};
92
10/96
✓ Branch 1 taken 4 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 4 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 88 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 88 times.
✗ Branch 11 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 16 not taken.
✗ Branch 17 not taken.
✓ Branch 19 taken 33 times.
✗ Branch 20 not taken.
✓ Branch 22 taken 33 times.
✗ Branch 23 not taken.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✗ Branch 28 not taken.
✗ Branch 29 not taken.
✓ Branch 31 taken 4 times.
✗ Branch 32 not taken.
✓ Branch 34 taken 4 times.
✗ Branch 35 not taken.
✗ Branch 37 not taken.
✗ Branch 38 not taken.
✗ Branch 40 not taken.
✗ Branch 41 not taken.
✓ Branch 43 taken 1 times.
✗ Branch 44 not taken.
✓ Branch 46 taken 1 times.
✗ Branch 47 not taken.
✗ Branch 49 not taken.
✗ 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 85 not taken.
✗ Branch 86 not taken.
✗ Branch 88 not taken.
✗ Branch 89 not taken.
✗ Branch 91 not taken.
✗ Branch 92 not taken.
✗ Branch 94 not taken.
✗ Branch 95 not taken.
✗ Branch 97 not taken.
✗ Branch 98 not taken.
✗ Branch 100 not taken.
✗ Branch 101 not taken.
✗ Branch 103 not taken.
✗ Branch 104 not taken.
✗ Branch 106 not taken.
✗ Branch 107 not taken.
✗ Branch 109 not taken.
✗ Branch 110 not taken.
✗ Branch 112 not taken.
✗ Branch 113 not taken.
✗ Branch 115 not taken.
✗ Branch 116 not taken.
✗ Branch 118 not taken.
✗ Branch 119 not taken.
✗ Branch 121 not taken.
✗ Branch 122 not taken.
✗ Branch 124 not taken.
✗ Branch 125 not taken.
✗ Branch 127 not taken.
✗ Branch 128 not taken.
✗ Branch 130 not taken.
✗ Branch 131 not taken.
✗ Branch 133 not taken.
✗ Branch 134 not taken.
✗ Branch 136 not taken.
✗ Branch 137 not taken.
✗ Branch 139 not taken.
✗ Branch 140 not taken.
✗ Branch 142 not taken.
✗ Branch 143 not taken.
1596 std::ranges::for_each(this->_point_field_names(), [&] (const std::string& name) {
93 492 pdata_helper.add(name, this->_get_point_field(name));
94 });
95
10/96
✓ Branch 1 taken 4 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 4 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 88 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 88 times.
✗ Branch 11 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 16 not taken.
✗ Branch 17 not taken.
✓ Branch 19 taken 33 times.
✗ Branch 20 not taken.
✓ Branch 22 taken 33 times.
✗ Branch 23 not taken.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✗ Branch 28 not taken.
✗ Branch 29 not taken.
✓ Branch 31 taken 4 times.
✗ Branch 32 not taken.
✓ Branch 34 taken 4 times.
✗ Branch 35 not taken.
✗ Branch 37 not taken.
✗ Branch 38 not taken.
✗ Branch 40 not taken.
✗ Branch 41 not taken.
✓ Branch 43 taken 1 times.
✗ Branch 44 not taken.
✓ Branch 46 taken 1 times.
✗ Branch 47 not taken.
✗ Branch 49 not taken.
✗ 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 85 not taken.
✗ Branch 86 not taken.
✗ Branch 88 not taken.
✗ Branch 89 not taken.
✗ Branch 91 not taken.
✗ Branch 92 not taken.
✗ Branch 94 not taken.
✗ Branch 95 not taken.
✗ Branch 97 not taken.
✗ Branch 98 not taken.
✗ Branch 100 not taken.
✗ Branch 101 not taken.
✗ Branch 103 not taken.
✗ Branch 104 not taken.
✗ Branch 106 not taken.
✗ Branch 107 not taken.
✗ Branch 109 not taken.
✗ Branch 110 not taken.
✗ Branch 112 not taken.
✗ Branch 113 not taken.
✗ Branch 115 not taken.
✗ Branch 116 not taken.
✗ Branch 118 not taken.
✗ Branch 119 not taken.
✗ Branch 121 not taken.
✗ Branch 122 not taken.
✗ Branch 124 not taken.
✗ Branch 125 not taken.
✗ Branch 127 not taken.
✗ Branch 128 not taken.
✗ Branch 130 not taken.
✗ Branch 131 not taken.
✗ Branch 133 not taken.
✗ Branch 134 not taken.
✗ Branch 136 not taken.
✗ Branch 137 not taken.
✗ Branch 139 not taken.
✗ Branch 140 not taken.
✗ Branch 142 not taken.
✗ Branch 143 not taken.
1590 std::ranges::for_each(this->_cell_field_names(), [&] (const std::string& name) {
96 489 cdata_helper.add(name, this->_get_cell_field(name));
97 });
98
6/24
✓ Branch 1 taken 4 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 89 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 3 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 32 times.
✗ Branch 11 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✓ Branch 16 taken 1 times.
✗ Branch 17 not taken.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 22 taken 1 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 34 not taken.
✗ Branch 35 not taken.
130 }, this->_xml_settings.data_format);
99
1/2
✓ Branch 1 taken 130 times.
✗ Branch 2 not taken.
211 }, this->_xml_settings.encoder);
100
101
4/8
✓ Branch 1 taken 130 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 130 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 130 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 130 times.
✗ Branch 11 not taken.
633 XMLElement& point_array = grid.add_child("PPoints").add_child("PDataArray");
102
2/4
✓ Branch 1 taken 130 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 130 times.
✗ Branch 5 not taken.
422 point_array.set_attribute("NumberOfComponents", "3");
103 358 std::visit([&] <typename T> (const Precision<T>& prec) {
104
12/48
✓ Branch 2 taken 94 times.
✗ Branch 3 not taken.
✓ Branch 5 taken 94 times.
✗ Branch 6 not taken.
✓ Branch 8 taken 94 times.
✗ Branch 9 not taken.
✓ Branch 16 taken 34 times.
✗ Branch 17 not taken.
✓ Branch 19 taken 34 times.
✗ Branch 20 not taken.
✓ Branch 22 taken 34 times.
✗ Branch 23 not taken.
✓ Branch 30 taken 1 times.
✗ Branch 31 not taken.
✓ Branch 33 taken 1 times.
✗ Branch 34 not taken.
✓ Branch 36 taken 1 times.
✗ Branch 37 not taken.
✓ Branch 44 taken 1 times.
✗ Branch 45 not taken.
✓ Branch 47 taken 1 times.
✗ Branch 48 not taken.
✓ Branch 50 taken 1 times.
✗ Branch 51 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 72 not taken.
✗ Branch 73 not taken.
✗ Branch 75 not taken.
✗ Branch 76 not taken.
✗ Branch 78 not taken.
✗ Branch 79 not taken.
✗ Branch 86 not taken.
✗ Branch 87 not taken.
✗ Branch 89 not taken.
✗ Branch 90 not taken.
✗ Branch 92 not taken.
✗ Branch 93 not taken.
✗ Branch 100 not taken.
✗ Branch 101 not taken.
✗ Branch 103 not taken.
✗ Branch 104 not taken.
✗ Branch 106 not taken.
✗ Branch 107 not taken.
390 point_array.set_attribute("type", VTK::attribute_name(prec));
105
1/2
✓ Branch 1 taken 130 times.
✗ Branch 2 not taken.
211 }, this->_xml_settings.coordinate_precision);
106
107
2/4
✓ Branch 1 taken 130 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 130 times.
✗ Branch 5 not taken.
529 std::ranges::for_each(Parallel::ranks(_comm), [&] (int rank) {
108
20/40
✓ Branch 1 taken 185 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 185 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 185 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 185 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 185 times.
✗ Branch 14 not taken.
✓ Branch 21 taken 69 times.
✗ Branch 22 not taken.
✓ Branch 24 taken 69 times.
✗ Branch 25 not taken.
✓ Branch 27 taken 69 times.
✗ Branch 28 not taken.
✓ Branch 30 taken 69 times.
✗ Branch 31 not taken.
✓ Branch 33 taken 69 times.
✗ Branch 34 not taken.
✓ Branch 41 taken 2 times.
✗ Branch 42 not taken.
✓ Branch 44 taken 2 times.
✗ Branch 45 not taken.
✓ Branch 47 taken 2 times.
✗ Branch 48 not taken.
✓ Branch 50 taken 2 times.
✗ Branch 51 not taken.
✓ Branch 53 taken 2 times.
✗ Branch 54 not taken.
✓ Branch 61 taken 2 times.
✗ Branch 62 not taken.
✓ Branch 64 taken 2 times.
✗ Branch 65 not taken.
✓ Branch 67 taken 2 times.
✗ Branch 68 not taken.
✓ Branch 70 taken 2 times.
✗ Branch 71 not taken.
✓ Branch 73 taken 2 times.
✗ Branch 74 not taken.
1548 grid.add_child("Piece").set_attribute("Source", std::filesystem::path{
109
8/16
✓ Branch 1 taken 185 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 185 times.
✗ Branch 5 not taken.
✓ Branch 11 taken 69 times.
✗ Branch 12 not taken.
✓ Branch 14 taken 69 times.
✗ Branch 15 not taken.
✓ Branch 21 taken 2 times.
✗ Branch 22 not taken.
✓ Branch 24 taken 2 times.
✗ Branch 25 not taken.
✓ Branch 31 taken 2 times.
✗ Branch 32 not taken.
✓ Branch 34 taken 2 times.
✗ Branch 35 not taken.
516 PVTK::piece_basefilename(filename_with_ext, rank) + ".vtu"
110
4/8
✓ Branch 1 taken 185 times.
✗ Branch 2 not taken.
✓ Branch 7 taken 69 times.
✗ Branch 8 not taken.
✓ Branch 13 taken 2 times.
✗ Branch 14 not taken.
✓ Branch 19 taken 2 times.
✗ Branch 20 not taken.
774 }.filename());
111 });
112
113
2/4
✓ Branch 2 taken 130 times.
✗ Branch 3 not taken.
✓ Branch 5 taken 130 times.
✗ Branch 6 not taken.
211 this->_set_default_active_fields(pvtk_xml.get_child("PUnstructuredGrid"));
114
2/2
✓ Branch 1 taken 128 times.
✓ Branch 2 taken 2 times.
373 write_to(filename_with_ext, [&] (std::ostream& s) {
115
8/16
✓ Branch 1 taken 94 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 94 times.
✗ Branch 5 not taken.
✓ Branch 9 taken 32 times.
✗ Branch 10 not taken.
✓ Branch 12 taken 32 times.
✗ Branch 13 not taken.
✓ Branch 17 taken 1 times.
✗ Branch 18 not taken.
✓ Branch 20 taken 1 times.
✗ Branch 21 not taken.
✓ Branch 25 taken 1 times.
✗ Branch 26 not taken.
✓ Branch 28 taken 1 times.
✗ Branch 29 not taken.
128 write_xml_with_version_header(pvtk_xml, s, Indentation{{.width = 2}});
116 });
117 211 }
118 };
119
120 template<typename G, Concepts::Communicator C>
121 PVTUWriter(G&&, const C&, VTK::XMLOptions = {}) -> PVTUWriter<std::remove_cvref_t<G>, C>;
122
123 } // namespace GridFormat
124
125 #endif // GRIDFORMAT_VTK_PVTU_WRITER_HPP_
126