GCC Code Coverage Report

Directory: ../../../builds/dumux-repositories/
File: /builds/dumux-repositories/dumux/dumux/io/vtkoutputmodule.hh
Date: 2024-05-04 19:09:25
Exec Total Coverage
Lines: 232 265 87.5%
Functions: 1527 2920 52.3%
Branches: 1044 2352 44.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-FileCopyrightInfo: 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 InputOutput
10 * \brief A VTK output module to simplify writing dumux simulation data to VTK format
11 */
12 #ifndef DUMUX_IO_VTK_OUTPUT_MODULE_HH
13 #define DUMUX_IO_VTK_OUTPUT_MODULE_HH
14
15 #include <functional>
16 #include <memory>
17 #include <string>
18
19 #include <dune/common/timer.hh>
20 #include <dune/common/fvector.hh>
21 #include <dune/common/typetraits.hh>
22
23 #include <dune/geometry/type.hh>
24 #include <dune/geometry/multilineargeometry.hh>
25
26 #include <dune/grid/common/mcmgmapper.hh>
27 #include <dune/grid/common/partitionset.hh>
28 #include <dune/grid/io/file/vtk/vtkwriter.hh>
29 #include <dune/grid/io/file/vtk/vtksequencewriter.hh>
30 #include <dune/grid/common/partitionset.hh>
31
32 #include <dumux/common/parameters.hh>
33 #include <dumux/io/format.hh>
34 #include <dumux/discretization/method.hh>
35
36 #include <dumux/io/vtk/function.hh>
37 #include <dumux/io/vtk/fieldtype.hh>
38 #include "velocityoutput.hh"
39
40 namespace Dumux {
41
42 /*!
43 * \ingroup InputOutput
44 * \brief A VTK output module to simplify writing dumux simulation data to VTK format
45 * \note This is a base class providing only rudimentary features
46 */
47 template<class GridGeometry>
48 class VtkOutputModuleBase
49 {
50 using GridView = typename GridGeometry::GridView;
51 static constexpr int dim = GridView::dimension;
52
53 public:
54 //! the type of Field that can be added to this writer
55 using Field = Vtk::template Field<GridView>;
56
57 956 VtkOutputModuleBase(const GridGeometry& gridGeometry,
58 const std::string& name,
59 const std::string& paramGroup = "",
60 Dune::VTK::DataMode dm = Dune::VTK::conforming,
61 bool verbose = true)
62 : gridGeometry_(gridGeometry)
63 , name_(name)
64 , paramGroup_(paramGroup)
65 , dm_(dm)
66
17/30
✓ Branch 2 taken 708 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 6 times.
✓ Branch 5 taken 607 times.
✓ Branch 6 taken 95 times.
✓ Branch 7 taken 613 times.
✓ Branch 8 taken 547 times.
✓ Branch 9 taken 161 times.
✓ Branch 10 taken 458 times.
✓ Branch 11 taken 244 times.
✗ Branch 12 not taken.
✓ Branch 13 taken 553 times.
✓ Branch 14 taken 141 times.
✓ Branch 15 taken 472 times.
✓ Branch 16 taken 95 times.
✓ Branch 17 taken 141 times.
✓ Branch 18 taken 472 times.
✓ Branch 19 taken 95 times.
✗ Branch 20 not taken.
✓ Branch 21 taken 472 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✗ Branch 31 not taken.
 1016 , verbose_(gridGeometry.gridView().comm().rank() == 0 && verbose)
67 {
68
1/2
✓ Branch 1 taken 708 times.
✗ Branch 2 not taken.
 956 const auto precisionString = getParamFromGroup<std::string>(paramGroup, "Vtk.Precision", "Float32");
69
2/4
✓ Branch 1 taken 708 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 708 times.
✗ Branch 5 not taken.
 1912 precision_ = Dumux::Vtk::stringToPrecision(precisionString);
70
4/10
✓ Branch 1 taken 708 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 708 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 708 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✓ Branch 10 taken 708 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
 956 const auto coordPrecision = Dumux::Vtk::stringToPrecision(getParamFromGroup<std::string>(paramGroup, "Vtk.CoordPrecision", precisionString));
71
4/8
✓ Branch 1 taken 708 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 708 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 708 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 708 times.
 1912 writer_ = std::make_shared<Dune::VTKWriter<GridView>>(gridGeometry.gridView(), dm, coordPrecision);
72
3/6
✓ Branch 1 taken 708 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 708 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 708 times.
 956 sequenceWriter_ = std::make_unique<Dune::VTKSequenceWriter<GridView>>(writer_, name);
73
2/6
✓ Branch 1 taken 708 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 708 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
 956 addProcessRank_ = getParamFromGroup<bool>(this->paramGroup(), "Vtk.AddProcessRank", true);
74 956 }
75
76
5/8
✓ Branch 1 taken 708 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 708 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✓ Branch 6 taken 708 times.
✓ Branch 7 taken 630 times.
✓ Branch 8 taken 78 times.
 4248 virtual ~VtkOutputModuleBase() = default;
77
78 //! the parameter group for getting parameter from the parameter tree
79 const std::string& paramGroup() const
80
43/86
✓ Branch 1 taken 50 times.
✓ Branch 2 taken 30 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 30 times.
✗ Branch 6 not taken.
✓ Branch 14 taken 1 times.
✗ Branch 15 not taken.
✓ Branch 17 taken 1 times.
✗ Branch 18 not taken.
✓ Branch 21 taken 1 times.
✗ Branch 22 not taken.
✓ Branch 24 taken 1 times.
✗ Branch 25 not taken.
✓ Branch 28 taken 1 times.
✗ Branch 29 not taken.
✓ Branch 31 taken 1 times.
✗ Branch 32 not taken.
✓ Branch 35 taken 1 times.
✗ Branch 36 not taken.
✓ Branch 38 taken 1 times.
✗ Branch 39 not taken.
✓ Branch 45 taken 1 times.
✗ Branch 46 not taken.
✓ Branch 48 taken 1 times.
✗ Branch 49 not taken.
✓ Branch 51 taken 1 times.
✗ Branch 52 not taken.
✓ Branch 54 taken 1 times.
✗ Branch 55 not taken.
✓ Branch 57 taken 1 times.
✗ Branch 58 not taken.
✓ Branch 60 taken 1 times.
✗ Branch 61 not taken.
✓ Branch 63 taken 1 times.
✗ Branch 64 not taken.
✓ Branch 66 taken 1 times.
✗ Branch 67 not taken.
✓ Branch 69 taken 1 times.
✗ Branch 70 not taken.
✓ Branch 72 taken 1 times.
✗ Branch 73 not taken.
✓ Branch 75 taken 1 times.
✗ Branch 76 not taken.
✓ Branch 78 taken 1 times.
✗ Branch 79 not taken.
✓ Branch 81 taken 1 times.
✗ Branch 82 not taken.
✓ Branch 84 taken 1 times.
✗ Branch 85 not taken.
✓ Branch 87 taken 1 times.
✗ Branch 88 not taken.
✓ Branch 90 taken 1 times.
✗ Branch 91 not taken.
✓ Branch 93 taken 1 times.
✗ Branch 94 not taken.
✓ Branch 96 taken 1 times.
✗ Branch 97 not taken.
✓ Branch 102 taken 1 times.
✗ Branch 103 not taken.
✓ Branch 105 taken 1 times.
✗ Branch 106 not taken.
✓ Branch 108 taken 1 times.
✗ Branch 109 not taken.
✓ Branch 111 taken 1 times.
✗ Branch 112 not taken.
✓ Branch 114 taken 1 times.
✗ Branch 115 not taken.
✓ Branch 117 taken 1 times.
✗ Branch 118 not taken.
✓ Branch 120 taken 1 times.
✗ Branch 121 not taken.
✓ Branch 123 taken 1 times.
✗ Branch 124 not taken.
✓ Branch 126 taken 1 times.
✗ Branch 127 not taken.
✓ Branch 129 taken 1 times.
✗ Branch 130 not taken.
✓ Branch 132 taken 1 times.
✗ Branch 133 not taken.
✓ Branch 135 taken 1 times.
✗ Branch 136 not taken.
✓ Branch 139 taken 1 times.
✗ Branch 140 not taken.
✓ Branch 142 taken 1 times.
✗ Branch 143 not taken.
 780 { return paramGroup_; }
81
82 /*!
83 * \brief Add a scalar or vector valued vtk field
84 *
85 * \param v The field to be added. Can be any indexable container. Its value type can be a number or itself an indexable container.
86 * \param name The name of the field
87 * \param fieldType The type of the field.
88 * This determines whether the values are associated with vertices or elements.
89 * By default, the method automatically deduces the correct type for the given input.
90 */
91 template<typename Vector>
92 void addField(const Vector& v,
93 const std::string& name,
94 Vtk::FieldType fieldType = Vtk::FieldType::automatic)
95
26/40
✓ Branch 1 taken 171 times.
✓ Branch 2 taken 17 times.
✓ Branch 3 taken 1 times.
✓ Branch 4 taken 115 times.
✓ Branch 5 taken 17 times.
✓ Branch 6 taken 1 times.
✓ Branch 7 taken 56 times.
✓ Branch 8 taken 17 times.
✓ Branch 9 taken 1 times.
✓ Branch 10 taken 37 times.
✗ Branch 11 not taken.
✓ Branch 12 taken 2 times.
✓ Branch 13 taken 6 times.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 2 times.
✓ Branch 19 taken 2 times.
✗ Branch 20 not taken.
✓ Branch 21 taken 2 times.
✗ Branch 22 not taken.
✓ Branch 24 taken 2 times.
✗ Branch 25 not taken.
✓ Branch 27 taken 2 times.
✗ Branch 28 not taken.
✓ Branch 30 taken 2 times.
✗ Branch 31 not taken.
✓ Branch 33 taken 2 times.
✗ Branch 34 not taken.
✓ Branch 36 taken 2 times.
✗ Branch 37 not taken.
✓ Branch 39 taken 1 times.
✗ Branch 40 not taken.
✓ Branch 42 taken 1 times.
✗ Branch 43 not taken.
✓ Branch 45 taken 1 times.
✗ Branch 46 not taken.
✓ Branch 48 taken 1 times.
✗ Branch 49 not taken.
 551 { addField(v, name, this->precision(), fieldType); }
96
97 /*!
98 * \brief Add a scalar or vector valued vtk field
99 *
100 * \param v The field to be added. Can be any indexable container. Its value type can be a number or itself an indexable container.
101 * \param name The name of the field
102 * \param fieldType The type of the field.
103 * This determines whether the values are associated with vertices or elements.
104 * By default, the method automatically deduces the correct type for the given input.
105 * \param precision The output precision of this field (see Dune::VTK::Precision)
106 */
107 template<typename Vector>
108 921 void addField(const Vector& v,
109 const std::string& name,
110 Dumux::Vtk::Precision precision,
111 Vtk::FieldType fieldType = Vtk::FieldType::automatic)
112 {
113 // Deduce the number of components from the given vector type
114 1056 const auto nComp = getNumberOfComponents_(v);
115
116
4/6
✓ Branch 0 taken 387 times.
✓ Branch 1 taken 166 times.
✓ Branch 2 taken 268 times.
✓ Branch 3 taken 165 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
 1657 const auto numElemDofs = gridGeometry().elementMapper().size();
117
4/4
✓ Branch 0 taken 397 times.
✓ Branch 1 taken 166 times.
✓ Branch 2 taken 328 times.
✓ Branch 3 taken 166 times.
 1756 const auto numVertexDofs = gridGeometry().vertexMapper().size();
118
119 // Automatically deduce the field type ...
120
2/2
✓ Branch 0 taken 397 times.
✓ Branch 1 taken 166 times.
 921 if(fieldType == Vtk::FieldType::automatic)
121 {
122
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 397 times.
 589 if(numElemDofs == numVertexDofs)
123 ✗ DUNE_THROW(Dune::InvalidStateException, "Automatic deduction of FieldType failed. Please explicitly specify FieldType::element or FieldType::vertex.");
124
125
4/4
✓ Branch 0 taken 61 times.
✓ Branch 1 taken 336 times.
✓ Branch 2 taken 52 times.
✓ Branch 3 taken 336 times.
 1165 if(v.size() == numElemDofs)
126 fieldType = Vtk::FieldType::element;
127
1/4
✗ Branch 0 not taken.
✓ Branch 1 taken 61 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
 72 else if(v.size() == numVertexDofs)
128 fieldType = Vtk::FieldType::vertex;
129 else
130 ✗ DUNE_THROW(Dune::RangeError, "Size mismatch of added field!");
131 }
132 // ... or check if the user-specified type matches the size of v
133 else
134 {
135
2/2
✓ Branch 0 taken 115 times.
✓ Branch 1 taken 51 times.
 332 if(fieldType == Vtk::FieldType::element)
136
2/4
✗ Branch 0 not taken.
✓ Branch 1 taken 115 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 115 times.
 460 if(v.size() != numElemDofs)
137 ✗ DUNE_THROW(Dune::RangeError, "Size mismatch of added field!");
138
139
2/2
✓ Branch 0 taken 51 times.
✓ Branch 1 taken 115 times.
 332 if(fieldType == Vtk::FieldType::vertex)
140
2/4
✗ Branch 0 not taken.
✓ Branch 1 taken 51 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 51 times.
 204 if(v.size() != numVertexDofs)
141 ✗ DUNE_THROW(Dune::RangeError, "Size mismatch of added field!");
142 }
143
144 // add the appropriate field
145
2/2
✓ Branch 0 taken 115 times.
✓ Branch 1 taken 51 times.
 332 if (fieldType == Vtk::FieldType::element)
146 2241 fields_.emplace_back(gridGeometry_.gridView(), gridGeometry_.elementMapper(), v, name, nComp, 0, dm_, precision);
147 else
148 522 fields_.emplace_back(gridGeometry_.gridView(), gridGeometry_.vertexMapper(), v, name, nComp, dim, dm_, precision);
149 921 }
150
151 /*!
152 * \brief Add a scalar or vector valued vtk field
153 * \param field The parameters passed on to the field constructor
154 */
155 void addField(Field&& field)
156 {
157 fields_.push_back(std::move(field));
158 }
159
160 //! Write the data for this timestep to file in four steps
161 //! (1) We assemble all registered variable fields
162 //! (2) We register them with the vtk writer
163 //! (3) The writer writes the output for us
164 //! (4) Clear the writer for the next time step
165 7737 void write(double time, Dune::VTK::OutputType type = Dune::VTK::ascii)
166 {
167 7737 Dune::Timer timer;
168
169 // write to file depending on data mode
170
2/2
✓ Branch 0 taken 6621 times.
✓ Branch 1 taken 52 times.
 7737 if (dm_ == Dune::VTK::conforming)
171 7665 writeConforming_(time, type);
172
1/2
✓ Branch 0 taken 52 times.
✗ Branch 1 not taken.
 72 else if (dm_ == Dune::VTK::nonconforming)
173 72 writeNonConforming_(time, type);
174 else
175 ✗ DUNE_THROW(Dune::NotImplemented, "Output for provided VTK data mode");
176
177 //! output
178 7737 timer.stop();
179
2/2
✓ Branch 0 taken 6300 times.
✓ Branch 1 taken 373 times.
 7737 if (verbose_)
180
2/6
✓ Branch 4 taken 6300 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 6300 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
 15832 std::cout << Fmt::format("Writing output for problem \"{}\". Took {:.2g} seconds.\n", name_, timer.elapsed());
181 7737 }
182
183 protected:
184 ✗ const GridGeometry& gridGeometry() const { return gridGeometry_; }
185
186 ✗ bool verbose() const { return verbose_; }
187
1/2
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
 212 const std::string& name() const { return name_; }
188 ✗ Dune::VTK::DataMode dataMode() const { return dm_; }
189 ✗ Dumux::Vtk::Precision precision() const { return precision_; }
190
191 13758 Dune::VTKWriter<GridView>& writer() { return *writer_; }
192
6/40
✓ Branch 1 taken 2060 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2060 times.
✗ Branch 5 not taken.
✗ Branch 7 not taken.
✗ 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 taken 206 times.
✗ Branch 20 not taken.
✓ Branch 22 taken 206 times.
✗ Branch 23 not taken.
✓ Branch 29 taken 206 times.
✗ Branch 30 not taken.
✓ Branch 32 taken 206 times.
✗ Branch 33 not taken.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 41 not taken.
✗ Branch 42 not taken.
✗ Branch 44 not taken.
✗ Branch 45 not taken.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✗ Branch 50 not taken.
✗ Branch 51 not taken.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
✗ Branch 63 not taken.
✗ Branch 64 not taken.
✗ Branch 66 not taken.
✗ Branch 67 not taken.
 206828 Dune::VTKSequenceWriter<GridView>& sequenceWriter() { return *sequenceWriter_; }
193
194
0/4
✗ Branch 0 not taken.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
 6878 const std::vector<Field>& fields() const { return fields_; }
195
196 private:
197 //! Assembles the fields and adds them to the writer (conforming output)
198 5 virtual void writeConforming_(double time, Dune::VTK::OutputType type)
199 {
200 //////////////////////////////////////////////////////////////
201 //! (1) Assemble all variable fields and add to writer
202 //////////////////////////////////////////////////////////////
203
204 // process rank
205
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 5 times.
 5 std::vector<int> rank;
206
207 //! Abort if no data was registered
208
2/6
✗ Branch 0 not taken.
✓ Branch 1 taken 5 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 5 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
 10 if (!fields_.empty() || addProcessRank_)
209 {
210
4/6
✓ Branch 0 taken 1 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 4 times.
✗ Branch 5 not taken.
 10 const auto numCells = gridGeometry_.gridView().size(0);
211
212 // maybe allocate space for the process rank
213
1/2
✓ Branch 0 taken 5 times.
✗ Branch 1 not taken.
 5 if (addProcessRank_)
214 {
215
1/2
✓ Branch 1 taken 5 times.
✗ Branch 2 not taken.
 5 rank.resize(numCells);
216
217
5/17
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 10 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 13 times.
✓ Branch 5 taken 1 times.
✓ Branch 6 taken 20 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
 67 for (const auto& element : elements(gridGeometry_.gridView(), Dune::Partitions::interior))
218 {
219
0/6
✗ Branch 0 not taken.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
 50 const auto eIdxGlobal = gridGeometry_.elementMapper().index(element);
220
2/13
✗ Branch 0 not taken.
✓ Branch 1 taken 16 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 16 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
 77 rank[eIdxGlobal] = gridGeometry_.gridView().comm().rank();
221 }
222 }
223
224 //////////////////////////////////////////////////////////////
225 //! (2) Register data fields with the vtk writer
226 //////////////////////////////////////////////////////////////
227
228 // the process rank
229
1/2
✓ Branch 0 taken 5 times.
✗ Branch 1 not taken.
 5 if (addProcessRank_)
230
10/26
✓ Branch 1 taken 5 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 5 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 5 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 5 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 5 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 5 times.
✗ Branch 17 not taken.
✓ Branch 20 taken 5 times.
✗ Branch 21 not taken.
✓ Branch 22 taken 5 times.
✗ Branch 23 not taken.
✓ Branch 24 taken 5 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✓ Branch 27 taken 5 times.
✗ Branch 28 not taken.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✗ Branch 31 not taken.
✗ Branch 32 not taken.
✗ Branch 33 not taken.
 25 this->sequenceWriter().addCellData(Field(gridGeometry_.gridView(), gridGeometry_.elementMapper(), rank, "process rank", 1, 0).get());
231
232 // also register additional (non-standardized) user fields if any
233
4/4
✓ Branch 0 taken 18 times.
✓ Branch 1 taken 5 times.
✓ Branch 2 taken 18 times.
✓ Branch 3 taken 5 times.
 33 for (auto&& field : fields_)
234 {
235
1/2
✓ Branch 0 taken 18 times.
✗ Branch 1 not taken.
 18 if (field.codim() == 0)
236
2/6
✓ Branch 3 taken 18 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 18 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
 54 this->sequenceWriter().addCellData(field.get());
237 ✗ else if (field.codim() == dim)
238 ✗ this->sequenceWriter().addVertexData(field.get());
239 else
240 ✗ DUNE_THROW(Dune::RangeError, "Cannot add wrongly sized vtk scalar field!");
241 }
242 }
243
244 //////////////////////////////////////////////////////////////
245 //! (2) The writer writes the output for us
246 //////////////////////////////////////////////////////////////
247
2/4
✓ Branch 1 taken 5 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 5 times.
✗ Branch 5 not taken.
 10 this->sequenceWriter().write(time, type);
248
249 //////////////////////////////////////////////////////////////
250 //! (3) Clear the writer
251 //////////////////////////////////////////////////////////////
252
1/2
✓ Branch 2 taken 5 times.
✗ Branch 3 not taken.
 10 this->writer().clear();
253 5 }
254
255 //! Assembles the fields and adds them to the writer (nonconforming output)
256 ✗ virtual void writeNonConforming_(double time, Dune::VTK::OutputType type)
257 {
258 ✗ DUNE_THROW(Dune::NotImplemented, "Non-conforming VTK output");
259 }
260
261 //! Deduces the number of components of the value type of a vector of values
262 template<class Vector>
263 ✗ std::size_t getNumberOfComponents_(const Vector& v)
264 {
265 if constexpr (Dune::IsIndexable<decltype(std::declval<Vector>()[0])>::value)
266
3/4
✓ Branch 0 taken 48 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 19 times.
✗ Branch 3 not taken.
 71 return v[0].size();
267 else
268 ✗ return 1;
269 }
270
271 const GridGeometry& gridGeometry_;
272 std::string name_;
273 const std::string paramGroup_;
274 Dune::VTK::DataMode dm_;
275 bool verbose_;
276 Dumux::Vtk::Precision precision_;
277
278 std::shared_ptr<Dune::VTKWriter<GridView>> writer_;
279 std::unique_ptr<Dune::VTKSequenceWriter<GridView>> sequenceWriter_;
280
281 std::vector<Field> fields_; //!< Registered scalar and vector fields
282
283 bool addProcessRank_ = true;
284 };
285
286 /*!
287 * \ingroup InputOutput
288 * \brief A VTK output module to simplify writing dumux simulation data to VTK format
289 *
290 * \tparam GridVariables The grid variables
291 * \tparam SolutionVector The solution vector
292 *
293 * Handles the output of scalar and vector fields to VTK formatted file for multiple
294 * variables and timesteps. Certain predefined fields can be registered on
295 * initialization and/or be turned on/off using the designated properties. Additionally
296 * non-standardized scalar and vector fields can be added to the writer manually.
297 */
298 template<class GridVariables, class SolutionVector>
299 class VtkOutputModule : public VtkOutputModuleBase<typename GridVariables::GridGeometry>
300 {
301 using ParentType = VtkOutputModuleBase<typename GridVariables::GridGeometry>;
302 using GridGeometry = typename GridVariables::GridGeometry;
303
304 using VV = typename GridVariables::VolumeVariables;
305 using Scalar = typename GridVariables::Scalar;
306
307 using GridView = typename GridGeometry::GridView;
308
309 enum {
310 dim = GridView::dimension,
311 dimWorld = GridView::dimensionworld
312 };
313
314 using Element = typename GridView::template Codim<0>::Entity;
315 using VolVarsVector = Dune::FieldVector<Scalar, dimWorld>;
316
317 static constexpr bool isBox = GridGeometry::discMethod == DiscretizationMethods::box;
318 static constexpr bool isDiamond = GridGeometry::discMethod == DiscretizationMethods::fcdiamond;
319 static constexpr bool isPQ1Bubble = GridGeometry::discMethod == DiscretizationMethods::pq1bubble;
320
321 struct VolVarScalarDataInfo { std::function<Scalar(const VV&)> get; std::string name; Dumux::Vtk::Precision precision_; };
322 struct VolVarVectorDataInfo { std::function<VolVarsVector(const VV&)> get; std::string name; Dumux::Vtk::Precision precision_; };
323
324 using VelocityOutputType = Dumux::VelocityOutput<GridVariables>;
325
326 public:
327 //! the type of Field that can be added to this writer
328 using Field = Vtk::template Field<GridView>;
329 //! export type of the volume variables for the outputfields
330 using VolumeVariables = VV;
331
332 987 VtkOutputModule(const GridVariables& gridVariables,
333 const SolutionVector& sol,
334 const std::string& name,
335 const std::string& paramGroup = "",
336 Dune::VTK::DataMode dm = Dune::VTK::conforming,
337 bool verbose = true)
338 : ParentType(gridVariables.gridGeometry(), name, paramGroup, dm, verbose)
339 , gridVariables_(gridVariables)
340 , sol_(sol)
341
3/6
✓ Branch 3 taken 703 times.
✗ Branch 4 not taken.
✓ Branch 6 taken 703 times.
✗ Branch 7 not taken.
✓ Branch 9 taken 703 times.
✗ Branch 10 not taken.
 1974 , velocityOutput_(std::make_shared<VelocityOutputType>())
342 {
343
1/2
✓ Branch 1 taken 703 times.
✗ Branch 2 not taken.
 987 enableVelocityOutput_ = getParamFromGroup<bool>(this->paramGroup(), "Vtk.AddVelocity", false);
344
1/4
✓ Branch 1 taken 703 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
 987 addProcessRank_ = getParamFromGroup<bool>(this->paramGroup(), "Vtk.AddProcessRank", true);
345 987 }
346
347 //////////////////////////////////////////////////////////////////////////////////////////////
348 //! Methods to conveniently add primary and secondary variables upon initialization
349 //! Do not call these methods after initialization i.e. _not_ within the time loop
350 //////////////////////////////////////////////////////////////////////////////////////////////
351
352 /*!
353 * \brief Add a velocity output policy
354 *
355 * \param velocityOutput the output policy
356 * \note the default policy does not add any velocity output
357 */
358 void addVelocityOutput(std::shared_ptr<VelocityOutputType> velocityOutput)
359
24/48
✓ Branch 1 taken 366 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 30 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 1 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 1 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 1 times.
✗ Branch 17 not taken.
✓ Branch 19 taken 1 times.
✗ Branch 20 not taken.
✓ Branch 22 taken 1 times.
✗ Branch 23 not taken.
✓ Branch 25 taken 1 times.
✗ Branch 26 not taken.
✓ Branch 28 taken 1 times.
✗ Branch 29 not taken.
✓ Branch 31 taken 1 times.
✗ Branch 32 not taken.
✓ Branch 34 taken 1 times.
✗ Branch 35 not taken.
✓ Branch 37 taken 1 times.
✗ Branch 38 not taken.
✓ Branch 40 taken 1 times.
✗ Branch 41 not taken.
✓ Branch 43 taken 1 times.
✗ Branch 44 not taken.
✓ Branch 46 taken 1 times.
✗ Branch 47 not taken.
✓ Branch 49 taken 1 times.
✗ Branch 50 not taken.
✓ Branch 52 taken 1 times.
✗ Branch 53 not taken.
✓ Branch 55 taken 1 times.
✗ Branch 56 not taken.
✓ Branch 58 taken 1 times.
✗ Branch 59 not taken.
✓ Branch 61 taken 1 times.
✗ Branch 62 not taken.
✓ Branch 64 taken 1 times.
✗ Branch 65 not taken.
✓ Branch 67 taken 1 times.
✗ Branch 68 not taken.
✓ Branch 70 taken 1 times.
✗ Branch 71 not taken.
 418 { velocityOutput_ = velocityOutput; }
360
361 //! Output a scalar volume variable
362 //! \param name The name of the vtk field
363 //! \param f A function taking a VolumeVariables object and returning the desired scalar
364 5412 void addVolumeVariable(std::function<Scalar(const VolumeVariables&)>&& f,
365 const std::string& name)
366 {
367
1/2
✓ Branch 3 taken 4372 times.
✗ Branch 4 not taken.
 5412 volVarScalarDataInfo_.push_back(VolVarScalarDataInfo{f, name, this->precision()});
368 5412 }
369
370 //! Add a vector-valued variable
371 //! \param f A function taking a VolumeVariables object and returning the desired vector
372 //! \param name The name of the vtk field
373 //! \note This method is only available for dimWorld > 1. For 1-D problems, the overload for volVar methods returning a Scalar will be used.
374 template<class VVV = VolVarsVector, typename std::enable_if_t<(VVV::dimension > 1), int> = 0>
375 150 void addVolumeVariable(std::function<VolVarsVector(const VolumeVariables&)>&& f,
376 const std::string& name)
377 {
378
1/2
✓ Branch 3 taken 110 times.
✗ Branch 4 not taken.
 150 volVarVectorDataInfo_.push_back(VolVarVectorDataInfo{f, name, this->precision()});
379 150 }
380
381 protected:
382 // some return functions for differing implementations to use
383
22/39
✓ Branch 3 taken 17 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
✓ Branch 6 taken 17 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 1 times.
✓ Branch 9 taken 17 times.
✗ Branch 10 not taken.
✓ Branch 11 taken 1 times.
✓ Branch 12 taken 17 times.
✗ Branch 13 not taken.
✓ Branch 14 taken 1 times.
✓ Branch 15 taken 17 times.
✗ Branch 16 not taken.
✓ Branch 17 taken 1 times.
✓ Branch 18 taken 17 times.
✗ Branch 19 not taken.
✓ Branch 20 taken 1 times.
✗ Branch 21 not taken.
✓ Branch 23 taken 2 times.
✗ Branch 24 not taken.
✓ Branch 26 taken 2 times.
✗ Branch 27 not taken.
✓ Branch 29 taken 2 times.
✗ Branch 30 not taken.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✓ Branch 35 taken 2 times.
✗ Branch 36 not taken.
✓ Branch 38 taken 2 times.
✗ Branch 39 not taken.
✓ Branch 43 taken 4368 times.
✗ Branch 44 not taken.
✓ Branch 46 taken 4368 times.
✗ Branch 47 not taken.
✓ Branch 51 taken 2014 times.
✗ Branch 52 not taken.
✓ Branch 54 taken 2014 times.
✗ Branch 55 not taken.
 455268 const auto& problem() const { return gridVariables_.curGridVolVars().problem(); }
384 ✗ const GridVariables& gridVariables() const { return gridVariables_; }
385
152/315
✗ Branch 0 not taken.
✓ Branch 1 taken 206 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 206 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 375739 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 48675 times.
✓ Branch 8 taken 327064 times.
✗ Branch 9 not taken.
✓ Branch 11 taken 2123 times.
✗ Branch 12 not taken.
✓ Branch 13 taken 60 times.
✓ Branch 14 taken 2123 times.
✓ Branch 15 taken 60 times.
✓ Branch 17 taken 63 times.
✗ Branch 18 not taken.
✓ Branch 20 taken 63 times.
✓ Branch 21 taken 60 times.
✗ Branch 22 not taken.
✓ Branch 23 taken 60 times.
✗ Branch 24 not taken.
✗ Branch 26 not taken.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✓ Branch 29 taken 266 times.
✗ Branch 30 not taken.
✓ Branch 31 taken 60 times.
✓ Branch 32 taken 206 times.
✗ Branch 33 not taken.
✗ Branch 34 not taken.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✓ Branch 37 taken 60 times.
✗ Branch 38 not taken.
✓ Branch 39 taken 60 times.
✗ Branch 40 not taken.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✓ Branch 45 taken 60 times.
✗ Branch 46 not taken.
✓ Branch 47 taken 60 times.
✗ Branch 48 not taken.
✗ Branch 49 not taken.
✗ Branch 51 not taken.
✗ Branch 52 not taken.
✓ Branch 53 taken 60 times.
✗ Branch 54 not taken.
✓ Branch 55 taken 60 times.
✗ Branch 57 not taken.
✗ Branch 58 not taken.
✗ Branch 60 not taken.
✓ Branch 61 taken 60 times.
✗ Branch 62 not taken.
✓ Branch 63 taken 60 times.
✗ Branch 64 not taken.
✗ Branch 66 not taken.
✗ Branch 67 not taken.
✗ Branch 68 not taken.
✓ Branch 69 taken 60 times.
✗ Branch 70 not taken.
✓ Branch 71 taken 60 times.
✗ Branch 76 not taken.
✓ Branch 77 taken 60 times.
✗ Branch 78 not taken.
✓ Branch 79 taken 60 times.
✗ Branch 84 not taken.
✓ Branch 85 taken 60 times.
✗ Branch 86 not taken.
✓ Branch 87 taken 60 times.
✗ Branch 92 not taken.
✓ Branch 93 taken 60 times.
✗ Branch 94 not taken.
✓ Branch 95 taken 60 times.
✗ Branch 100 not taken.
✓ Branch 101 taken 60 times.
✗ Branch 102 not taken.
✓ Branch 103 taken 60 times.
✗ Branch 108 not taken.
✓ Branch 109 taken 60 times.
✗ Branch 110 not taken.
✓ Branch 111 taken 60 times.
✗ Branch 116 not taken.
✓ Branch 117 taken 60 times.
✗ Branch 118 not taken.
✓ Branch 119 taken 60 times.
✗ Branch 124 not taken.
✓ Branch 125 taken 60 times.
✗ Branch 126 not taken.
✓ Branch 127 taken 60 times.
✗ Branch 132 not taken.
✓ Branch 133 taken 60 times.
✗ Branch 134 not taken.
✓ Branch 135 taken 60 times.
✗ Branch 140 not taken.
✓ Branch 141 taken 60 times.
✗ Branch 142 not taken.
✓ Branch 143 taken 60 times.
✗ Branch 148 not taken.
✓ Branch 149 taken 60 times.
✗ Branch 150 not taken.
✓ Branch 151 taken 60 times.
✗ Branch 156 not taken.
✓ Branch 157 taken 60 times.
✗ Branch 158 not taken.
✓ Branch 159 taken 60 times.
✗ Branch 164 not taken.
✓ Branch 165 taken 60 times.
✗ Branch 166 not taken.
✓ Branch 167 taken 60 times.
✗ Branch 172 not taken.
✓ Branch 173 taken 60 times.
✗ Branch 174 not taken.
✓ Branch 175 taken 60 times.
✗ Branch 180 not taken.
✓ Branch 181 taken 60 times.
✗ Branch 182 not taken.
✓ Branch 183 taken 60 times.
✗ Branch 188 not taken.
✓ Branch 189 taken 60 times.
✗ Branch 190 not taken.
✓ Branch 191 taken 60 times.
✓ Branch 195 taken 1 times.
✗ Branch 196 not taken.
✓ Branch 198 taken 1 times.
✗ Branch 199 not taken.
✓ Branch 201 taken 1 times.
✗ Branch 202 not taken.
✓ Branch 204 taken 1 times.
✗ Branch 205 not taken.
✓ Branch 209 taken 1 times.
✗ Branch 210 not taken.
✓ Branch 212 taken 1 times.
✗ Branch 213 not taken.
✓ Branch 215 taken 1 times.
✗ Branch 216 not taken.
✓ Branch 218 taken 1 times.
✗ Branch 219 not taken.
✓ Branch 223 taken 1 times.
✗ Branch 224 not taken.
✓ Branch 226 taken 1 times.
✗ Branch 227 not taken.
✓ Branch 229 taken 1 times.
✗ Branch 230 not taken.
✓ Branch 232 taken 1 times.
✗ Branch 233 not taken.
✓ Branch 237 taken 1 times.
✗ Branch 238 not taken.
✓ Branch 240 taken 1 times.
✗ Branch 241 not taken.
✓ Branch 243 taken 1 times.
✗ Branch 244 not taken.
✓ Branch 246 taken 1 times.
✗ Branch 247 not taken.
✓ Branch 251 taken 1 times.
✗ Branch 252 not taken.
✓ Branch 254 taken 1 times.
✗ Branch 255 not taken.
✓ Branch 257 taken 1 times.
✗ Branch 258 not taken.
✓ Branch 260 taken 1 times.
✗ Branch 261 not taken.
✓ Branch 265 taken 1 times.
✗ Branch 266 not taken.
✓ Branch 268 taken 1 times.
✗ Branch 269 not taken.
✓ Branch 271 taken 1 times.
✗ Branch 272 not taken.
✓ Branch 274 taken 1 times.
✗ Branch 275 not taken.
✓ Branch 279 taken 1 times.
✗ Branch 280 not taken.
✓ Branch 282 taken 1 times.
✗ Branch 283 not taken.
✓ Branch 285 taken 1 times.
✗ Branch 286 not taken.
✓ Branch 288 taken 1 times.
✗ Branch 289 not taken.
✓ Branch 293 taken 1 times.
✗ Branch 294 not taken.
✓ Branch 296 taken 1 times.
✗ Branch 297 not taken.
✓ Branch 299 taken 1 times.
✗ Branch 300 not taken.
✓ Branch 302 taken 1 times.
✗ Branch 303 not taken.
✓ Branch 307 taken 1 times.
✗ Branch 308 not taken.
✓ Branch 310 taken 1 times.
✗ Branch 311 not taken.
✓ Branch 313 taken 1 times.
✗ Branch 314 not taken.
✓ Branch 316 taken 1 times.
✗ Branch 317 not taken.
✓ Branch 321 taken 1 times.
✗ Branch 322 not taken.
✓ Branch 324 taken 1 times.
✗ Branch 325 not taken.
✓ Branch 327 taken 1 times.
✗ Branch 328 not taken.
✓ Branch 330 taken 1 times.
✗ Branch 331 not taken.
✓ Branch 335 taken 1 times.
✗ Branch 336 not taken.
✓ Branch 338 taken 1 times.
✗ Branch 339 not taken.
✓ Branch 341 taken 1 times.
✗ Branch 342 not taken.
✓ Branch 344 taken 1 times.
✗ Branch 345 not taken.
✓ Branch 349 taken 1 times.
✗ Branch 350 not taken.
✓ Branch 352 taken 1 times.
✗ Branch 353 not taken.
✓ Branch 355 taken 1 times.
✗ Branch 356 not taken.
✓ Branch 358 taken 1 times.
✗ Branch 359 not taken.
✓ Branch 363 taken 1 times.
✗ Branch 364 not taken.
✓ Branch 366 taken 1 times.
✗ Branch 367 not taken.
✓ Branch 369 taken 1 times.
✗ Branch 370 not taken.
✓ Branch 372 taken 1 times.
✗ Branch 373 not taken.
✓ Branch 377 taken 1 times.
✗ Branch 378 not taken.
✓ Branch 380 taken 1 times.
✗ Branch 381 not taken.
✓ Branch 383 taken 1 times.
✗ Branch 384 not taken.
✓ Branch 386 taken 1 times.
✗ Branch 387 not taken.
✓ Branch 391 taken 1 times.
✗ Branch 392 not taken.
✓ Branch 394 taken 1 times.
✗ Branch 395 not taken.
✓ Branch 397 taken 1 times.
✗ Branch 398 not taken.
✓ Branch 400 taken 1 times.
✗ Branch 401 not taken.
✓ Branch 405 taken 1 times.
✗ Branch 406 not taken.
✓ Branch 408 taken 1 times.
✗ Branch 409 not taken.
✓ Branch 411 taken 1 times.
✗ Branch 412 not taken.
✓ Branch 414 taken 1 times.
✗ Branch 415 not taken.
✓ Branch 419 taken 1 times.
✗ Branch 420 not taken.
✓ Branch 422 taken 1 times.
✗ Branch 423 not taken.
✓ Branch 425 taken 1 times.
✗ Branch 426 not taken.
✓ Branch 428 taken 1 times.
✗ Branch 429 not taken.
✓ Branch 433 taken 1 times.
✗ Branch 434 not taken.
✓ Branch 436 taken 1 times.
✗ Branch 437 not taken.
✓ Branch 439 taken 1 times.
✗ Branch 440 not taken.
✓ Branch 442 taken 1 times.
✗ Branch 443 not taken.
✓ Branch 447 taken 1 times.
✗ Branch 448 not taken.
✓ Branch 450 taken 1 times.
✗ Branch 451 not taken.
✓ Branch 453 taken 1 times.
✗ Branch 454 not taken.
✓ Branch 456 taken 1 times.
✗ Branch 457 not taken.
✓ Branch 461 taken 1 times.
✗ Branch 462 not taken.
✓ Branch 464 taken 1 times.
✗ Branch 465 not taken.
✓ Branch 467 taken 1 times.
✗ Branch 468 not taken.
✓ Branch 470 taken 1 times.
✗ Branch 471 not taken.
✓ Branch 475 taken 1 times.
✗ Branch 476 not taken.
✓ Branch 478 taken 1 times.
✗ Branch 479 not taken.
✓ Branch 481 taken 1 times.
✗ Branch 482 not taken.
✓ Branch 484 taken 1 times.
✗ Branch 485 not taken.
✓ Branch 489 taken 1 times.
✗ Branch 490 not taken.
✓ Branch 492 taken 1 times.
✗ Branch 493 not taken.
✓ Branch 495 taken 1 times.
✗ Branch 496 not taken.
✓ Branch 498 taken 1 times.
✗ Branch 499 not taken.
✓ Branch 503 taken 1 times.
✗ Branch 504 not taken.
✓ Branch 506 taken 1 times.
✗ Branch 507 not taken.
✓ Branch 509 taken 1 times.
✗ Branch 510 not taken.
✓ Branch 512 taken 1 times.
✗ Branch 513 not taken.
✓ Branch 517 taken 1 times.
✗ Branch 518 not taken.
✓ Branch 520 taken 1 times.
✗ Branch 521 not taken.
✓ Branch 523 taken 1 times.
✗ Branch 524 not taken.
✓ Branch 526 taken 1 times.
✗ Branch 527 not taken.
 25527542 const GridGeometry& gridGeometry() const { return gridVariables_.gridGeometry(); }
386 ✗ const SolutionVector& sol() const { return sol_; }
387
388
2/8
✗ Branch 0 not taken.
✓ Branch 1 taken 206 times.
✓ Branch 2 taken 206 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
 412 const std::vector<VolVarScalarDataInfo>& volVarScalarDataInfo() const { return volVarScalarDataInfo_; }
389
2/8
✗ Branch 0 not taken.
✓ Branch 1 taken 206 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 206 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
 412 const std::vector<VolVarVectorDataInfo>& volVarVectorDataInfo() const { return volVarVectorDataInfo_; }
390
391 using VelocityOutput = VelocityOutputType;
392
8/80
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✓ Branch 7 taken 206 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 206 times.
✗ Branch 11 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 16 not taken.
✗ Branch 17 not taken.
✓ Branch 19 taken 327064 times.
✗ Branch 20 not taken.
✓ Branch 22 taken 327064 times.
✗ Branch 23 not taken.
✓ Branch 25 taken 327064 times.
✗ Branch 26 not taken.
✓ Branch 28 taken 327064 times.
✗ Branch 29 not taken.
✗ Branch 31 not taken.
✗ Branch 32 not taken.
✗ Branch 34 not taken.
✗ Branch 35 not taken.
✗ Branch 37 not taken.
✗ Branch 38 not taken.
✗ Branch 40 not taken.
✗ Branch 41 not taken.
✓ Branch 43 taken 206 times.
✗ Branch 44 not taken.
✓ Branch 46 taken 206 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.
 1309080 const VelocityOutput& velocityOutput() const { return *velocityOutput_; }
393
394 private:
395
396 //! Assembles the fields and adds them to the writer (conforming output)
397 7836 void writeConforming_(double time, Dune::VTK::OutputType type) override
398 {
399 7836 const Dune::VTK::DataMode dm = Dune::VTK::conforming;
400 //////////////////////////////////////////////////////////////
401 //! (1) Assemble all variable fields and add to writer
402 //////////////////////////////////////////////////////////////
403
404 // instantiate the velocity output
405 using VelocityVector = typename VelocityOutput::VelocityVector;
406
4/8
✓ Branch 1 taken 6616 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 6616 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 6616 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 6616 times.
✗ Branch 11 not taken.
 23508 std::vector<VelocityVector> velocity(velocityOutput_->numFluidPhases());
407
408 // process rank
409
2/4
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 6606 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
 15672 std::vector<double> rank;
410
411 // volume variable data
412
4/4
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 6606 times.
✓ Branch 3 taken 6576 times.
✓ Branch 4 taken 40 times.
 15672 std::vector<std::vector<Scalar>> volVarScalarData;
413
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 6606 times.
 15672 std::vector<std::vector<VolVarsVector>> volVarVectorData;
414
415 //! Abort if no data was registered
416
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 6606 times.
 7836 if (!volVarScalarDataInfo_.empty()
417
4/4
✓ Branch 0 taken 4 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 4 times.
✓ Branch 3 taken 6 times.
 20 || !volVarVectorDataInfo_.empty()
418
2/4
✗ Branch 0 not taken.
✓ Branch 1 taken 4 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 4 times.
 8 || !this->fields().empty()
419 ✗ || velocityOutput_->enableOutput()
420
2/4
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 6606 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
 7836 || addProcessRank_)
421 {
422
5/8
✓ Branch 0 taken 11 times.
✓ Branch 1 taken 5319 times.
✓ Branch 2 taken 11 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 5330 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 5319 times.
✗ Branch 8 not taken.
 23508 const auto numCells = gridGeometry().gridView().size(0);
423 15672 const auto numDofs = numDofs_();
424
425 // get fields for all volume variables
426
4/4
✓ Branch 0 taken 6606 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 5311 times.
✓ Branch 3 taken 10 times.
 14097 if (!volVarScalarDataInfo_.empty())
427
5/12
✓ Branch 1 taken 6606 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 6606 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 6606 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 6606 times.
✗ Branch 11 not taken.
✓ Branch 12 taken 6606 times.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.
 31304 volVarScalarData.resize(volVarScalarDataInfo_.size(), std::vector<Scalar>(numDofs));
428
4/4
✓ Branch 0 taken 1139 times.
✓ Branch 1 taken 5477 times.
✓ Branch 2 taken 1139 times.
✓ Branch 3 taken 5477 times.
 15672 if (!volVarVectorDataInfo_.empty())
429
5/12
✓ Branch 1 taken 1139 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 1139 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 1139 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 1139 times.
✗ Branch 11 not taken.
✓ Branch 12 taken 1139 times.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.
 4732 volVarVectorData.resize(volVarVectorDataInfo_.size(), std::vector<VolVarsVector>(numDofs));
430
431
4/6
✓ Branch 1 taken 6616 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 6616 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 3078 times.
✓ Branch 7 taken 3538 times.
 15672 if (velocityOutput_->enableOutput())
432 {
433
4/6
✓ Branch 1 taken 6693 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 6693 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 3615 times.
✓ Branch 7 taken 3078 times.
 11200 for (int phaseIdx = 0; phaseIdx < velocityOutput_->numFluidPhases(); ++phaseIdx)
434 {
435
4/6
✓ Branch 1 taken 3615 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 3615 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 266 times.
✓ Branch 7 taken 3349 times.
 8216 if (velocityOutput_->fieldType() == VelocityOutput::FieldType::element)
436
2/4
✓ Branch 1 taken 266 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 266 times.
✗ Branch 5 not taken.
 592 velocity[phaseIdx].resize(numCells);
437
3/6
✓ Branch 1 taken 3349 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 3349 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 3349 times.
 7624 else if (velocityOutput_->fieldType() == VelocityOutput::FieldType::vertex)
438 ✗ velocity[phaseIdx].resize(numDofs);
439 else
440 {
441 if(isBox && dim == 1)
442
2/4
✓ Branch 1 taken 667 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 667 times.
✗ Branch 5 not taken.
 1342 velocity[phaseIdx].resize(numCells);
443 else
444
2/4
✓ Branch 1 taken 2682 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2682 times.
✗ Branch 5 not taken.
 6282 velocity[phaseIdx].resize(numDofs);
445 }
446 }
447 }
448
449 // maybe allocate space for the process rank
450
3/4
✓ Branch 0 taken 6576 times.
✓ Branch 1 taken 40 times.
✓ Branch 3 taken 6576 times.
✗ Branch 4 not taken.
 7836 if (addProcessRank_) rank.resize(numCells);
451
452
3/6
✓ Branch 1 taken 5815 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 5815 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 98 times.
✗ Branch 7 not taken.
 18418 auto fvGeometry = localView(gridGeometry());
453
5/8
✓ Branch 1 taken 5819 times.
✓ Branch 2 taken 548 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 5819 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 1291 times.
✓ Branch 7 taken 151 times.
✗ Branch 8 not taken.
 24549 auto elemVolVars = localView(gridVariables_.curGridVolVars());
454
19/22
✓ Branch 1 taken 5819 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 443532 times.
✓ Branch 4 taken 6605 times.
✓ Branch 5 taken 443521 times.
✓ Branch 6 taken 797 times.
✓ Branch 7 taken 5819 times.
✓ Branch 8 taken 179267 times.
✓ Branch 9 taken 4809624 times.
✓ Branch 10 taken 1011 times.
✓ Branch 11 taken 480730 times.
✓ Branch 12 taken 446 times.
✓ Branch 13 taken 52317 times.
✓ Branch 14 taken 480730 times.
✓ Branch 15 taken 1375366 times.
✓ Branch 16 taken 565 times.
✓ Branch 17 taken 1856096 times.
✓ Branch 18 taken 565 times.
✓ Branch 20 taken 1375366 times.
✗ Branch 21 not taken.
✓ Branch 23 taken 1375366 times.
✗ Branch 24 not taken.
 15679557 for (const auto& element : elements(gridGeometry().gridView(), Dune::Partitions::interior))
455 {
456
3/6
✓ Branch 1 taken 2087115 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2087115 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 2087115 times.
✗ Branch 8 not taken.
 27503670 const auto eIdxGlobal = gridGeometry().elementMapper().index(element);
457 // If velocity output is enabled we need to bind to the whole stencil
458 // otherwise element-local data is sufficient
459
4/6
✓ Branch 1 taken 6623703 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 6623703 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2747294 times.
✓ Branch 7 taken 3876409 times.
 18335780 if (velocityOutput_->enableOutput())
460 {
461
1/3
✗ Branch 0 not taken.
✓ Branch 1 taken 1087076 times.
✗ Branch 2 not taken.
 3738630 fvGeometry.bind(element);
462
1/2
✓ Branch 1 taken 2637294 times.
✗ Branch 2 not taken.
 3738630 elemVolVars.bind(element, fvGeometry, sol_);
463 }
464 else
465 {
466
2/3
✓ Branch 0 taken 175134 times.
✓ Branch 1 taken 2332419 times.
✗ Branch 2 not taken.
 5429260 fvGeometry.bindElement(element);
467
1/2
✓ Branch 1 taken 1946654 times.
✗ Branch 2 not taken.
 4931824 elemVolVars.bindElement(element, fvGeometry, sol_);
468 }
469
470
8/8
✓ Branch 0 taken 54488 times.
✓ Branch 1 taken 6569215 times.
✓ Branch 2 taken 54488 times.
✓ Branch 3 taken 6569215 times.
✓ Branch 4 taken 42000 times.
✓ Branch 5 taken 12488 times.
✓ Branch 6 taken 42000 times.
✓ Branch 7 taken 12488 times.
 18335780 if (!volVarScalarDataInfo_.empty() || !volVarVectorDataInfo_.empty())
471 {
472
4/4
✓ Branch 0 taken 10800651 times.
✓ Branch 1 taken 6611215 times.
✓ Branch 2 taken 8937835 times.
✓ Branch 3 taken 4861639 times.
 29511406 for (const auto& scv : scvs(fvGeometry))
473 {
474 13803684 const auto dofIdxGlobal = scv.dofIndex();
475 13803684 const auto& volVars = elemVolVars[scv];
476
477 // get the scalar-valued data
478
5/5
✓ Branch 0 taken 83088838 times.
✓ Branch 1 taken 31540762 times.
✓ Branch 2 taken 86501077 times.
✓ Branch 3 taken 31540762 times.
✓ Branch 4 taken 3412239 times.
 252021472 for (std::size_t i = 0; i < volVarScalarDataInfo_.size(); ++i)
479
2/4
✗ Branch 0 not taken.
✓ Branch 1 taken 107241188 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 107241188 times.
 595544470 volVarScalarData[i][dofIdxGlobal] = volVarScalarDataInfo_[i].get(volVars);
480
481 // get the vector-valued data
482
4/4
✓ Branch 0 taken 713800 times.
✓ Branch 1 taken 10800651 times.
✓ Branch 2 taken 713800 times.
✓ Branch 3 taken 10800651 times.
 29061096 for (std::size_t i = 0; i < volVarVectorDataInfo_.size(); ++i)
483
2/4
✗ Branch 0 not taken.
✓ Branch 1 taken 713800 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 713800 times.
 3634320 volVarVectorData[i][dofIdxGlobal] = volVarVectorDataInfo_[i].get(volVars);
484 }
485 }
486
487 // velocity output
488
4/6
✓ Branch 1 taken 6623703 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 6623703 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2747294 times.
✓ Branch 7 taken 3876409 times.
 18335780 if (velocityOutput_->enableOutput())
489 {
490
4/10
✓ Branch 1 taken 958290 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 926938 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 926938 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 20294 times.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
 12576102 const auto elemFluxVarsCache = localView(gridVariables_.gridFluxVarsCache()).bind(element, fvGeometry, elemVolVars);
491
492
4/6
✓ Branch 1 taken 5568440 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 5568440 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2821146 times.
✓ Branch 7 taken 2747294 times.
 7561112 for (int phaseIdx = 0; phaseIdx < velocityOutput_->numFluidPhases(); ++phaseIdx)
493
3/6
✓ Branch 1 taken 2821146 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2821146 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 2821146 times.
✗ Branch 8 not taken.
 11467446 velocityOutput_->calculateVelocity(velocity[phaseIdx], element, fvGeometry, elemVolVars, elemFluxVarsCache, phaseIdx);
494 }
495
496 //! the rank
497
2/2
✓ Branch 0 taken 6308703 times.
✓ Branch 1 taken 315000 times.
 9167890 if (addProcessRank_)
498
3/6
✗ Branch 0 not taken.
✓ Branch 1 taken 4853237 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 4853237 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 4853237 times.
 26798970 rank[eIdxGlobal] = static_cast<double>(gridGeometry().gridView().comm().rank());
499 }
500
501 //////////////////////////////////////////////////////////////
502 //! (2) Register data fields with the vtk writer
503 //////////////////////////////////////////////////////////////
504
505 // volume variables if any
506 if constexpr (isBox || isPQ1Bubble)
507 {
508
4/4
✓ Branch 0 taken 28673 times.
✓ Branch 1 taken 2699 times.
✓ Branch 2 taken 28673 times.
✓ Branch 3 taken 2699 times.
 60843 for (std::size_t i = 0; i < volVarScalarDataInfo_.size(); ++i)
509
9/22
✓ Branch 1 taken 28673 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 28673 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 28673 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 28673 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 28673 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 28673 times.
✗ Branch 17 not taken.
✓ Branch 20 taken 28673 times.
✗ Branch 21 not taken.
✓ Branch 22 taken 28673 times.
✗ Branch 23 not taken.
✓ Branch 24 taken 28673 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✗ Branch 29 not taken.
 232184 this->sequenceWriter().addVertexData( Field(gridGeometry().gridView(), gridGeometry().dofMapper(), volVarScalarData[i],
510
2/4
✓ Branch 1 taken 28673 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 28673 times.
✗ Branch 5 not taken.
 58046 volVarScalarDataInfo_[i].name, /*numComp*/1, /*codim*/dim, dm, this->precision()).get() );
511
4/4
✓ Branch 0 taken 32 times.
✓ Branch 1 taken 2699 times.
✓ Branch 2 taken 32 times.
✓ Branch 3 taken 2699 times.
 5706 for (std::size_t i = 0; i < volVarVectorDataInfo_.size(); ++i)
512
9/22
✓ Branch 1 taken 32 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 32 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 32 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 32 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 32 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 32 times.
✗ Branch 17 not taken.
✓ Branch 20 taken 32 times.
✗ Branch 21 not taken.
✓ Branch 22 taken 32 times.
✗ Branch 23 not taken.
✓ Branch 24 taken 32 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✗ Branch 29 not taken.
 448 this->sequenceWriter().addVertexData( Field(gridGeometry().gridView(), gridGeometry().dofMapper(), volVarVectorData[i],
513
2/4
✓ Branch 1 taken 32 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 32 times.
✗ Branch 5 not taken.
 112 volVarVectorDataInfo_[i].name, /*numComp*/dimWorld, /*codim*/dim, dm, this->precision()).get() );
514
515 if constexpr (isPQ1Bubble)
516 {
517
4/4
✓ Branch 0 taken 4 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 4 times.
✓ Branch 3 taken 10 times.
 28 for (std::size_t i = 0; i < volVarScalarDataInfo_.size(); ++i)
518
9/22
✓ Branch 1 taken 4 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 4 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 4 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 4 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 4 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 4 times.
✗ Branch 17 not taken.
✓ Branch 20 taken 4 times.
✗ Branch 21 not taken.
✓ Branch 22 taken 4 times.
✗ Branch 23 not taken.
✓ Branch 24 taken 4 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✗ Branch 29 not taken.
 32 this->sequenceWriter().addCellData( Field(gridGeometry().gridView(), gridGeometry().dofMapper(), volVarScalarData[i],
519
2/4
✓ Branch 1 taken 4 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 4 times.
✗ Branch 5 not taken.
 8 volVarScalarDataInfo_[i].name, /*numComp*/1, /*codim*/0,dm, this->precision()).get() );
520
4/4
✓ Branch 0 taken 6 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 6 times.
✓ Branch 3 taken 10 times.
 32 for (std::size_t i = 0; i < volVarVectorDataInfo_.size(); ++i)
521
9/22
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 6 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 6 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 6 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 6 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 6 times.
✗ Branch 17 not taken.
✓ Branch 20 taken 6 times.
✗ Branch 21 not taken.
✓ Branch 22 taken 6 times.
✗ Branch 23 not taken.
✓ Branch 24 taken 6 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✗ Branch 29 not taken.
 48 this->sequenceWriter().addCellData( Field(gridGeometry().gridView(), gridGeometry().dofMapper(), volVarVectorData[i],
522
2/4
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 6 times.
✗ Branch 5 not taken.
 12 volVarVectorDataInfo_[i].name, /*numComp*/dimWorld, /*codim*/0,dm, this->precision()).get() );
523 }
524
525 }
526 else
527 {
528
4/4
✓ Branch 0 taken 44700 times.
✓ Branch 1 taken 3917 times.
✓ Branch 2 taken 44700 times.
✓ Branch 3 taken 3917 times.
 109757 for (std::size_t i = 0; i < volVarScalarDataInfo_.size(); ++i)
529
9/22
✓ Branch 1 taken 44700 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 44700 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 44700 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 44700 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 44700 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 44700 times.
✗ Branch 17 not taken.
✓ Branch 20 taken 44700 times.
✗ Branch 21 not taken.
✓ Branch 22 taken 44700 times.
✗ Branch 23 not taken.
✓ Branch 24 taken 44700 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✗ Branch 29 not taken.
 418872 this->sequenceWriter().addCellData( Field(gridGeometry().gridView(), gridGeometry().elementMapper(), volVarScalarData[i],
530
2/4
✓ Branch 1 taken 44700 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 44700 times.
✗ Branch 5 not taken.
 104718 volVarScalarDataInfo_[i].name, /*numComp*/1, /*codim*/0,dm, this->precision()).get() );
531
4/4
✓ Branch 0 taken 2214 times.
✓ Branch 1 taken 3917 times.
✓ Branch 2 taken 2214 times.
✓ Branch 3 taken 3917 times.
 14586 for (std::size_t i = 0; i < volVarVectorDataInfo_.size(); ++i)
532
9/22
✓ Branch 1 taken 2214 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2214 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 2214 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 2214 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 2214 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 2214 times.
✗ Branch 17 not taken.
✓ Branch 20 taken 2214 times.
✗ Branch 21 not taken.
✓ Branch 22 taken 2214 times.
✗ Branch 23 not taken.
✓ Branch 24 taken 2214 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✗ Branch 29 not taken.
 18032 this->sequenceWriter().addCellData( Field(gridGeometry().gridView(), gridGeometry().elementMapper(), volVarVectorData[i],
533
2/4
✓ Branch 1 taken 2214 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2214 times.
✗ Branch 5 not taken.
 4508 volVarVectorDataInfo_[i].name, /*numComp*/dimWorld, /*codim*/0,dm, this->precision()).get() );
534 }
535
536 // the velocity field
537
4/6
✓ Branch 1 taken 6616 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 6616 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 3078 times.
✓ Branch 7 taken 3538 times.
 15672 if (velocityOutput_->enableOutput())
538 {
539
2/4
✓ Branch 1 taken 3078 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 3078 times.
✗ Branch 5 not taken.
 7092 if (velocityOutput_->fieldType() == VelocityOutput::FieldType::vertex
540
5/8
✓ Branch 0 taken 3078 times.
✗ Branch 1 not taken.
✓ Branch 3 taken 3078 times.
✗ Branch 4 not taken.
✓ Branch 6 taken 3078 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 508 times.
✓ Branch 9 taken 9 times.
 3546 || ( (velocityOutput_->fieldType() == VelocityOutput::FieldType::automatic) && dim > 1 && isBox ))
541 {
542
4/6
✓ Branch 1 taken 1153 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 1153 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 645 times.
✓ Branch 7 taken 508 times.
 2306 for (int phaseIdx = 0; phaseIdx < velocityOutput_->numFluidPhases(); ++phaseIdx)
543
15/40
✓ Branch 1 taken 645 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 645 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 645 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 645 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 645 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 645 times.
✗ Branch 17 not taken.
✓ Branch 19 taken 645 times.
✗ Branch 20 not taken.
✓ Branch 22 taken 645 times.
✗ Branch 23 not taken.
✓ Branch 25 taken 645 times.
✗ Branch 26 not taken.
✓ Branch 29 taken 645 times.
✗ Branch 30 not taken.
✓ Branch 31 taken 645 times.
✗ Branch 32 not taken.
✓ Branch 33 taken 645 times.
✗ Branch 34 not taken.
✓ Branch 35 taken 645 times.
✗ Branch 36 not taken.
✗ Branch 37 not taken.
✓ Branch 38 taken 645 times.
✗ Branch 39 not taken.
✓ Branch 40 taken 645 times.
✗ Branch 41 not taken.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✗ Branch 45 not taken.
✗ Branch 46 not taken.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✗ Branch 49 not taken.
✗ Branch 50 not taken.
 2580 this->sequenceWriter().addVertexData( Field(gridGeometry().gridView(), gridGeometry().vertexMapper(), velocity[phaseIdx],
544
2/4
✓ Branch 1 taken 645 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 645 times.
✗ Branch 5 not taken.
 1290 "velocity_" + velocityOutput_->phaseName(phaseIdx) + " (m/s)",
545 /*numComp*/dimWorld, /*codim*/dim, dm, this->precision()).get() );
546 }
547 // cell-centered models
548 else
549 {
550
4/6
✓ Branch 1 taken 5540 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 5540 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2970 times.
✓ Branch 7 taken 2570 times.
 13002 for (int phaseIdx = 0; phaseIdx < velocityOutput_->numFluidPhases(); ++phaseIdx)
551
15/40
✓ Branch 1 taken 2970 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2970 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 2970 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 2970 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 2970 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 2970 times.
✗ Branch 17 not taken.
✓ Branch 19 taken 2970 times.
✗ Branch 20 not taken.
✓ Branch 22 taken 2970 times.
✗ Branch 23 not taken.
✓ Branch 25 taken 2970 times.
✗ Branch 26 not taken.
✓ Branch 29 taken 2970 times.
✗ Branch 30 not taken.
✓ Branch 31 taken 2970 times.
✗ Branch 32 not taken.
✓ Branch 33 taken 2970 times.
✗ Branch 34 not taken.
✓ Branch 35 taken 2970 times.
✗ Branch 36 not taken.
✗ Branch 37 not taken.
✓ Branch 38 taken 2970 times.
✗ Branch 39 not taken.
✓ Branch 40 taken 2970 times.
✗ Branch 41 not taken.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✗ Branch 45 not taken.
✗ Branch 46 not taken.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✗ Branch 49 not taken.
✗ Branch 50 not taken.
 13852 this->sequenceWriter().addCellData( Field(gridGeometry().gridView(), gridGeometry().elementMapper(), velocity[phaseIdx],
552
2/4
✓ Branch 1 taken 2970 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 2970 times.
✗ Branch 5 not taken.
 6926 "velocity_" + velocityOutput_->phaseName(phaseIdx) + " (m/s)",
553 /*numComp*/dimWorld, /*codim*/0, dm, this->precision()).get() );
554 }
555 }
556
557 // the process rank
558
2/2
✓ Branch 0 taken 6576 times.
✓ Branch 1 taken 40 times.
 7836 if (addProcessRank_)
559
12/30
✓ Branch 1 taken 6576 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 6576 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 6576 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 6576 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 6576 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 6576 times.
✗ Branch 17 not taken.
✓ Branch 19 taken 6576 times.
✗ Branch 20 not taken.
✓ Branch 22 taken 6576 times.
✗ Branch 23 not taken.
✓ Branch 26 taken 6576 times.
✗ Branch 27 not taken.
✓ Branch 28 taken 6576 times.
✗ Branch 29 not taken.
✓ Branch 30 taken 6576 times.
✗ Branch 31 not taken.
✗ Branch 32 not taken.
✓ Branch 33 taken 6576 times.
✗ Branch 34 not taken.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✗ Branch 37 not taken.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
 38980 this->sequenceWriter().addCellData(Field(gridGeometry().gridView(), gridGeometry().elementMapper(), rank, "process rank", 1, 0).get());
560
561 // also register additional (non-standardized) user fields if any
562
4/4
✓ Branch 0 taken 8124 times.
✓ Branch 1 taken 6616 times.
✓ Branch 2 taken 8124 times.
✓ Branch 3 taken 6616 times.
 32223 for (auto&& field : this->fields())
563 {
564
2/2
✓ Branch 0 taken 5831 times.
✓ Branch 1 taken 2293 times.
 8715 if (field.codim() == 0)
565
2/6
✓ Branch 3 taken 5831 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 5831 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
 19122 this->sequenceWriter().addCellData(field.get());
566
1/2
✓ Branch 0 taken 2293 times.
✗ Branch 1 not taken.
 2341 else if (field.codim() == dim)
567
2/8
✓ Branch 3 taken 2293 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2293 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
 7023 this->sequenceWriter().addVertexData(field.get());
568 else
569 ✗ DUNE_THROW(Dune::RangeError, "Cannot add wrongly sized vtk scalar field!");
570 }
571 }
572
573 //////////////////////////////////////////////////////////////
574 //! (2) The writer writes the output for us
575 //////////////////////////////////////////////////////////////
576
2/4
✓ Branch 1 taken 6616 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 6616 times.
✗ Branch 5 not taken.
 15672 this->sequenceWriter().write(time, type);
577
578 //////////////////////////////////////////////////////////////
579 //! (3) Clear the writer
580 //////////////////////////////////////////////////////////////
581 15672 this->writer().clear();
582 7836 }
583
584 //! Assembles the fields and adds them to the writer (nonconforming output)
585 72 void writeNonConforming_(double time, Dune::VTK::OutputType type) override
586 {
587 72 const Dune::VTK::DataMode dm = Dune::VTK::nonconforming;
588
589 // only supports finite-element-like discretization schemes
590 if(!isBox && !isDiamond)
591 ✗ DUNE_THROW(Dune::NotImplemented,
592 "Non-conforming output for discretization scheme " << GridGeometry::discMethod
593 );
594
595 //////////////////////////////////////////////////////////////
596 //! (1) Assemble all variable fields and add to writer
597 //////////////////////////////////////////////////////////////
598
599 // check the velocity output
600
1/4
✗ Branch 0 not taken.
✓ Branch 1 taken 52 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
 72 if (enableVelocityOutput_ && !velocityOutput_->enableOutput())
601 std::cerr << "Warning! Velocity output was enabled in the input file"
602 << " but no velocity output policy was set for the VTK output module:"
603 << " There will be no velocity output."
604 ✗ << " Use the addVelocityOutput member function of the VTK output module." << std::endl;
605 using VelocityVector = typename VelocityOutput::VelocityVector;
606
4/8
✓ Branch 1 taken 52 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 52 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 52 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 52 times.
✗ Branch 11 not taken.
 216 std::vector<VelocityVector> velocity(velocityOutput_->numFluidPhases());
607
608 // process rank
609
2/4
✓ Branch 0 taken 6 times.
✓ Branch 1 taken 46 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
 144 std::vector<double> rank;
610
611 // volume variable data (indexing: volvardata/element/localindex)
612 using ScalarDataContainer = std::vector< std::vector<Scalar> >;
613 using VectorDataContainer = std::vector< std::vector<VolVarsVector> >;
614
3/4
✓ Branch 0 taken 6 times.
✓ Branch 1 taken 46 times.
✓ Branch 3 taken 52 times.
✗ Branch 4 not taken.
 144 std::vector< ScalarDataContainer > volVarScalarData;
615
2/2
✓ Branch 0 taken 6 times.
✓ Branch 1 taken 46 times.
 144 std::vector< VectorDataContainer > volVarVectorData;
616
617 //! Abort if no data was registered
618
2/2
✓ Branch 0 taken 6 times.
✓ Branch 1 taken 46 times.
 72 if (!volVarScalarDataInfo_.empty()
619
2/4
✗ Branch 0 not taken.
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 6 times.
 12 || !volVarVectorDataInfo_.empty()
620 ✗ || !this->fields().empty()
621 ✗ || velocityOutput_->enableOutput()
622
2/4
✓ Branch 0 taken 6 times.
✓ Branch 1 taken 46 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
 72 || addProcessRank_)
623 {
624
3/8
✗ Branch 0 not taken.
✓ Branch 1 taken 46 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 46 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 46 times.
✗ Branch 8 not taken.
 216 const auto numCells = gridGeometry().gridView().size(0);
625 144 const auto outputSize = numDofs_();
626
627 // get fields for all volume variables
628
3/4
✓ Branch 0 taken 46 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 18 times.
✗ Branch 3 not taken.
 90 if (!volVarScalarDataInfo_.empty())
629
4/8
✓ Branch 1 taken 46 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 46 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 46 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 46 times.
✗ Branch 11 not taken.
 132 volVarScalarData.resize(volVarScalarDataInfo_.size(), ScalarDataContainer(numCells));
630
4/4
✓ Branch 0 taken 6 times.
✓ Branch 1 taken 46 times.
✓ Branch 2 taken 6 times.
✓ Branch 3 taken 46 times.
 144 if (!volVarVectorDataInfo_.empty())
631
4/8
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 6 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 6 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 6 times.
✗ Branch 11 not taken.
 12 volVarVectorData.resize(volVarVectorDataInfo_.size(), VectorDataContainer(numCells));
632
633
4/6
✓ Branch 1 taken 52 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 52 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 4 times.
✓ Branch 7 taken 48 times.
 144 if (velocityOutput_->enableOutput())
634 {
635
4/6
✓ Branch 1 taken 8 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 8 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 4 times.
✓ Branch 7 taken 4 times.
 12 for (int phaseIdx = 0; phaseIdx < velocityOutput_->numFluidPhases(); ++phaseIdx)
636 {
637 if((isBox && dim == 1) || isDiamond)
638
2/4
✓ Branch 1 taken 4 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 4 times.
✗ Branch 5 not taken.
 8 velocity[phaseIdx].resize(numCells);
639 else
640 ✗ velocity[phaseIdx].resize(outputSize);
641 }
642 }
643
644 // maybe allocate space for the process rank
645
2/4
✓ Branch 0 taken 52 times.
✗ Branch 1 not taken.
✓ Branch 3 taken 52 times.
✗ Branch 4 not taken.
 72 if (addProcessRank_) rank.resize(numCells);
646
647 // now we go element-local to extract values at local dof locations
648
2/6
✓ Branch 1 taken 52 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 52 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
 144 auto fvGeometry = localView(gridGeometry());
649
3/7
✓ Branch 1 taken 52 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 52 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 48 times.
✗ Branch 7 not taken.
 280 auto elemVolVars = localView(gridVariables_.curGridVolVars());
650
15/22
✓ Branch 1 taken 52 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 52 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 52 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 125378 times.
✓ Branch 10 taken 26 times.
✓ Branch 11 taken 55504 times.
✓ Branch 12 taken 6 times.
✓ Branch 13 taken 20 times.
✓ Branch 14 taken 55504 times.
✓ Branch 15 taken 36000 times.
✓ Branch 16 taken 20 times.
✓ Branch 17 taken 91504 times.
✓ Branch 18 taken 20 times.
✓ Branch 20 taken 36000 times.
✗ Branch 21 not taken.
✓ Branch 23 taken 36000 times.
✗ Branch 24 not taken.
 267416 for (const auto& element : elements(gridGeometry().gridView(), Dune::Partitions::interior))
651 {
652
3/6
✓ Branch 1 taken 91504 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 91504 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 91504 times.
✗ Branch 8 not taken.
 592056 const auto eIdxGlobal = gridGeometry().elementMapper().index(element);
653 // If velocity output is enabled we need to bind to the whole stencil
654 // otherwise element-local data is sufficient
655
4/6
✓ Branch 1 taken 161352 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 161352 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 55304 times.
✓ Branch 7 taken 106048 times.
 394704 if (velocityOutput_->enableOutput())
656 {
657
1/2
✓ Branch 1 taken 55304 times.
✗ Branch 2 not taken.
 55304 fvGeometry.bind(element);
658
1/2
✓ Branch 1 taken 55304 times.
✗ Branch 2 not taken.
 55304 elemVolVars.bind(element, fvGeometry, sol_);
659 }
660 else
661 {
662
1/2
✓ Branch 1 taken 63848 times.
✗ Branch 2 not taken.
 142048 fvGeometry.bindElement(element);
663
1/2
✓ Branch 1 taken 106048 times.
✗ Branch 2 not taken.
 142048 elemVolVars.bindElement(element, fvGeometry, sol_);
664 }
665
666 197352 const auto numLocalDofs = fvGeometry.numScv();
667 // resize element-local data containers
668
4/4
✓ Branch 0 taken 465904 times.
✓ Branch 1 taken 161352 times.
✓ Branch 2 taken 465904 times.
✓ Branch 3 taken 161352 times.
 741672 for (std::size_t i = 0; i < volVarScalarDataInfo_.size(); ++i)
669
3/6
✓ Branch 1 taken 465904 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 465904 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 465904 times.
✗ Branch 8 not taken.
 1632960 volVarScalarData[i][eIdxGlobal].resize(numLocalDofs);
670
4/4
✓ Branch 0 taken 42000 times.
✓ Branch 1 taken 161352 times.
✓ Branch 2 taken 42000 times.
✓ Branch 3 taken 161352 times.
 436704 for (std::size_t i = 0; i < volVarVectorDataInfo_.size(); ++i)
671
3/6
✓ Branch 1 taken 42000 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 42000 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 42000 times.
✗ Branch 8 not taken.
 126000 volVarVectorData[i][eIdxGlobal].resize(numLocalDofs);
672
673
6/8
✓ Branch 0 taken 42000 times.
✓ Branch 1 taken 119352 times.
✓ Branch 2 taken 42000 times.
✓ Branch 3 taken 119352 times.
✓ Branch 4 taken 42000 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 42000 times.
✗ Branch 7 not taken.
 394704 if (!volVarScalarDataInfo_.empty() || !volVarVectorDataInfo_.empty())
674 {
675
4/4
✓ Branch 0 taken 609408 times.
✓ Branch 1 taken 161352 times.
✓ Branch 2 taken 609408 times.
✓ Branch 3 taken 161352 times.
 1112112 for (const auto& scv : scvs(fvGeometry))
676 {
677 717408 const auto& volVars = elemVolVars[scv];
678
679 // get the scalar-valued data
680
4/4
✓ Branch 0 taken 1785200 times.
✓ Branch 1 taken 609408 times.
✓ Branch 2 taken 1785200 times.
✓ Branch 3 taken 609408 times.
 4758304 for (std::size_t i = 0; i < volVarScalarDataInfo_.size(); ++i)
681
2/4
✗ Branch 0 not taken.
✓ Branch 1 taken 1785200 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 1785200 times.
 12122688 volVarScalarData[i][eIdxGlobal][scv.localDofIndex()] = volVarScalarDataInfo_[i].get(volVars);
682
683 // get the vector-valued data
684
4/4
✓ Branch 0 taken 168000 times.
✓ Branch 1 taken 609408 times.
✓ Branch 2 taken 168000 times.
✓ Branch 3 taken 609408 times.
 1770816 for (std::size_t i = 0; i < volVarVectorDataInfo_.size(); ++i)
685
2/4
✗ Branch 0 not taken.
✓ Branch 1 taken 168000 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 168000 times.
 1008000 volVarVectorData[i][eIdxGlobal][scv.localDofIndex()] = volVarVectorDataInfo_[i].get(volVars);
686 }
687 }
688
689 // velocity output
690
4/6
✓ Branch 1 taken 161352 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 161352 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 55304 times.
✓ Branch 7 taken 106048 times.
 394704 if (velocityOutput_->enableOutput())
691 {
692
3/10
✓ Branch 1 taken 55304 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 55304 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 55304 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
 165912 const auto elemFluxVarsCache = localView(gridVariables_.gridFluxVarsCache()).bind(element, fvGeometry, elemVolVars);
693
4/6
✓ Branch 1 taken 110608 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 110608 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 55304 times.
✓ Branch 7 taken 55304 times.
 110608 for (int phaseIdx = 0; phaseIdx < velocityOutput_->numFluidPhases(); ++phaseIdx)
694
3/6
✓ Branch 1 taken 55304 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 55304 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 55304 times.
✗ Branch 8 not taken.
 165912 velocityOutput_->calculateVelocity(velocity[phaseIdx], element, fvGeometry, elemVolVars, elemFluxVarsCache, phaseIdx);
695 }
696
697 //! the rank
698
1/2
✓ Branch 0 taken 161352 times.
✗ Branch 1 not taken.
 197352 if (addProcessRank_)
699
3/6
✗ Branch 0 not taken.
✓ Branch 1 taken 125352 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 125352 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 125352 times.
 592056 rank[eIdxGlobal] = static_cast<double>(gridGeometry().gridView().comm().rank());
700 }
701
702 //////////////////////////////////////////////////////////////
703 //! Register data fields with the vtk writer
704 //////////////////////////////////////////////////////////////
705
706 // volume variables if any
707 static constexpr int dofLocCodim = isDiamond ? 1 : dim;
708
4/4
✓ Branch 0 taken 234 times.
✓ Branch 1 taken 52 times.
✓ Branch 2 taken 234 times.
✓ Branch 3 taken 52 times.
 624 for (std::size_t i = 0; i < volVarScalarDataInfo_.size(); ++i)
709
7/18
✓ Branch 1 taken 234 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 234 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 234 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 234 times.
✗ Branch 11 not taken.
✓ Branch 14 taken 234 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 234 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 234 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
 1656 this->sequenceWriter().addVertexData(Field(
710
1/2
✓ Branch 1 taken 234 times.
✗ Branch 2 not taken.
 276 gridGeometry().gridView(), gridGeometry().elementMapper(),
711
3/6
✓ Branch 1 taken 234 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 234 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 234 times.
✗ Branch 8 not taken.
 828 volVarScalarData[i], volVarScalarDataInfo_[i].name,
712 /*numComp*/1, /*codim*/dofLocCodim, /*nonconforming*/dm, this->precision()
713 ).get());
714
715
4/4
✓ Branch 0 taken 6 times.
✓ Branch 1 taken 52 times.
✓ Branch 2 taken 6 times.
✓ Branch 3 taken 52 times.
 156 for (std::size_t i = 0; i < volVarVectorDataInfo_.size(); ++i)
716
7/18
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 6 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 6 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 6 times.
✗ Branch 11 not taken.
✓ Branch 14 taken 6 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 6 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 6 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
 36 this->sequenceWriter().addVertexData(Field(
717
1/2
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
 6 gridGeometry().gridView(), gridGeometry().elementMapper(),
718
3/6
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 6 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 6 times.
✗ Branch 8 not taken.
 18 volVarVectorData[i], volVarVectorDataInfo_[i].name,
719 /*numComp*/dimWorld, /*codim*/dofLocCodim, /*nonconforming*/dm, this->precision()
720 ).get());
721
722 // the velocity field
723
4/6
✓ Branch 1 taken 52 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 52 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 4 times.
✓ Branch 7 taken 48 times.
 144 if (velocityOutput_->enableOutput())
724 {
725 // node-wise velocities
726 if (dim > 1 && !isDiamond)
727 ✗ for (int phaseIdx = 0; phaseIdx < velocityOutput_->numFluidPhases(); ++phaseIdx)
728 ✗ this->sequenceWriter().addVertexData(Field(
729 ✗ gridGeometry().gridView(), gridGeometry().vertexMapper(), velocity[phaseIdx],
730 ✗ "velocity_" + velocityOutput_->phaseName(phaseIdx) + " (m/s)",
731 /*numComp*/dimWorld, /*codim*/dofLocCodim, dm, this->precision()
732 ).get());
733
734 // cell-wise velocities
735 else
736
4/6
✓ Branch 1 taken 8 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 8 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 4 times.
✓ Branch 7 taken 4 times.
 16 for (int phaseIdx = 0; phaseIdx < velocityOutput_->numFluidPhases(); ++phaseIdx)
737
12/36
✓ Branch 1 taken 4 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 4 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 4 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 4 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 4 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 4 times.
✗ Branch 17 not taken.
✓ Branch 20 taken 4 times.
✗ Branch 21 not taken.
✓ Branch 22 taken 4 times.
✗ Branch 23 not taken.
✓ Branch 24 taken 4 times.
✗ Branch 25 not taken.
✓ Branch 26 taken 4 times.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✓ Branch 29 taken 4 times.
✗ Branch 30 not taken.
✓ Branch 31 taken 4 times.
✗ Branch 32 not taken.
✗ Branch 33 not taken.
✗ Branch 34 not taken.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✗ Branch 37 not taken.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✗ Branch 41 not taken.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
 16 this->sequenceWriter().addCellData(Field(
738
3/6
✓ Branch 1 taken 4 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 4 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 4 times.
✗ Branch 8 not taken.
 12 gridGeometry().gridView(), gridGeometry().elementMapper(), velocity[phaseIdx],
739
2/4
✓ Branch 1 taken 4 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 4 times.
✗ Branch 5 not taken.
 8 "velocity_" + velocityOutput_->phaseName(phaseIdx) + " (m/s)",
740 /*numComp*/dimWorld, /*codim*/0, dm, this->precision()
741 ).get());
742 }
743 }
744
745 // the process rank
746
1/2
✓ Branch 0 taken 52 times.
✗ Branch 1 not taken.
 72 if (addProcessRank_)
747
10/26
✓ Branch 1 taken 52 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 52 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 52 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 52 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 52 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 52 times.
✗ Branch 17 not taken.
✓ Branch 20 taken 52 times.
✗ Branch 21 not taken.
✓ Branch 22 taken 52 times.
✗ Branch 23 not taken.
✓ Branch 24 taken 52 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✓ Branch 27 taken 52 times.
✗ Branch 28 not taken.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✗ Branch 31 not taken.
✗ Branch 32 not taken.
✗ Branch 33 not taken.
 360 this->sequenceWriter().addCellData(Field(
748
2/4
✓ Branch 1 taken 52 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 52 times.
✗ Branch 5 not taken.
 144 gridGeometry().gridView(), gridGeometry().elementMapper(),
749 rank, "process rank", /*numComp*/1, /*codim*/0
750 ).get());
751
752 // also register additional (non-standardized) user fields if any
753
4/4
✓ Branch 0 taken 6 times.
✓ Branch 1 taken 52 times.
✓ Branch 2 taken 6 times.
✓ Branch 3 taken 52 times.
 222 for (const auto& field : this->fields())
754 {
755
1/2
✓ Branch 0 taken 6 times.
✗ Branch 1 not taken.
 6 if (field.codim() == 0)
756
2/6
✓ Branch 3 taken 6 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 6 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
 18 this->sequenceWriter().addCellData(field.get());
757 ✗ else if (field.codim() == dim || field.codim() == 1)
758 ✗ this->sequenceWriter().addVertexData(field.get());
759 else
760 ✗ DUNE_THROW(Dune::RangeError, "Cannot add wrongly sized vtk scalar field!");
761 }
762
763 //////////////////////////////////////////////////////////////
764 //! (2) The writer writes the output for us
765 //////////////////////////////////////////////////////////////
766
2/4
✓ Branch 1 taken 52 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 52 times.
✗ Branch 5 not taken.
 144 this->sequenceWriter().write(time, type);
767
768 //////////////////////////////////////////////////////////////
769 //! (3) Clear the writer
770 //////////////////////////////////////////////////////////////
771 144 this->writer().clear();
772 72 }
773
774 //! return the number of dofs, we only support vertex and cell data
775 std::size_t numDofs_() const
776 {
777 // TODO this should actually always be dofMapper.size()
778 // maybe some discretizations needs special treatment (?)
779 if constexpr (isBox || isDiamond || isPQ1Bubble)
780
29/52
✓ Branch 0 taken 39 times.
✓ Branch 1 taken 12 times.
✓ Branch 2 taken 39 times.
✓ Branch 3 taken 12 times.
✓ Branch 4 taken 39 times.
✓ Branch 5 taken 12 times.
✓ Branch 6 taken 2628 times.
✓ Branch 7 taken 16 times.
✓ Branch 8 taken 2634 times.
✓ Branch 9 taken 4 times.
✓ Branch 10 taken 2646 times.
✓ Branch 11 taken 4 times.
✓ Branch 12 taken 509 times.
✓ Branch 13 taken 12 times.
✓ Branch 14 taken 6 times.
✗ Branch 15 not taken.
✗ Branch 16 not taken.
✗ Branch 17 not taken.
✓ Branch 18 taken 23 times.
✗ Branch 19 not taken.
✓ Branch 20 taken 26 times.
✗ Branch 21 not taken.
✓ Branch 22 taken 26 times.
✗ Branch 23 not taken.
✓ Branch 24 taken 24 times.
✗ Branch 25 not taken.
✓ Branch 26 taken 14 times.
✗ Branch 27 not taken.
✓ Branch 28 taken 11 times.
✗ Branch 29 not taken.
✓ Branch 30 taken 13 times.
✗ Branch 31 not taken.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✓ Branch 34 taken 2 times.
✗ Branch 35 not taken.
✓ Branch 36 taken 2 times.
✗ Branch 37 not taken.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✗ Branch 41 not taken.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✓ Branch 44 taken 11 times.
✗ Branch 45 not taken.
✓ Branch 46 taken 11 times.
✗ Branch 47 not taken.
✓ Branch 48 taken 11 times.
✗ Branch 49 not taken.
✓ Branch 50 taken 11 times.
✗ Branch 51 not taken.
 8799 return gridGeometry().dofMapper().size();
781 else
782
78/144
✓ Branch 0 taken 3242 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 3230 times.
✓ Branch 3 taken 12 times.
✓ Branch 4 taken 3230 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 455 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 405 times.
✓ Branch 9 taken 101 times.
✓ Branch 10 taken 405 times.
✓ Branch 11 taken 101 times.
✓ Branch 12 taken 90 times.
✓ Branch 13 taken 101 times.
✓ Branch 14 taken 45 times.
✓ Branch 15 taken 101 times.
✓ Branch 16 taken 25 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 126 times.
✗ Branch 19 not taken.
✓ Branch 20 taken 25 times.
✓ Branch 21 taken 101 times.
✓ Branch 22 taken 23 times.
✗ Branch 23 not taken.
✓ Branch 24 taken 23 times.
✗ Branch 25 not taken.
✓ Branch 26 taken 23 times.
✗ Branch 27 not taken.
✓ Branch 28 taken 1 times.
✗ Branch 29 not taken.
✓ Branch 30 taken 1 times.
✗ Branch 31 not taken.
✓ Branch 32 taken 1 times.
✗ Branch 33 not taken.
✓ Branch 34 taken 1 times.
✗ Branch 35 not taken.
✓ Branch 36 taken 1 times.
✗ Branch 37 not taken.
✓ Branch 38 taken 1 times.
✗ Branch 39 not taken.
✓ Branch 40 taken 1 times.
✗ Branch 41 not taken.
✓ Branch 42 taken 1 times.
✗ Branch 43 not taken.
✓ Branch 44 taken 1 times.
✗ Branch 45 not taken.
✓ Branch 46 taken 1 times.
✗ Branch 47 not taken.
✓ Branch 48 taken 1 times.
✗ Branch 49 not taken.
✓ Branch 50 taken 1 times.
✗ Branch 51 not taken.
✓ Branch 52 taken 1 times.
✗ Branch 53 not taken.
✓ Branch 54 taken 1 times.
✗ Branch 55 not taken.
✓ Branch 56 taken 1 times.
✗ Branch 57 not taken.
✓ Branch 58 taken 1 times.
✗ Branch 59 not taken.
✓ Branch 60 taken 1 times.
✗ Branch 61 not taken.
✓ Branch 62 taken 1 times.
✗ Branch 63 not taken.
✓ Branch 64 taken 1 times.
✗ Branch 65 not taken.
✓ Branch 66 taken 1 times.
✗ Branch 67 not taken.
✓ Branch 68 taken 1 times.
✗ Branch 69 not taken.
✓ Branch 70 taken 1 times.
✗ Branch 71 not taken.
✓ Branch 72 taken 1 times.
✗ Branch 73 not taken.
✓ Branch 74 taken 1 times.
✗ Branch 75 not taken.
✓ Branch 76 taken 1 times.
✗ Branch 77 not taken.
✓ Branch 78 taken 1 times.
✗ Branch 79 not taken.
✓ Branch 80 taken 1 times.
✗ Branch 81 not taken.
✓ Branch 82 taken 1 times.
✗ Branch 83 not taken.
✓ Branch 84 taken 1 times.
✗ Branch 85 not taken.
✓ Branch 86 taken 1 times.
✗ Branch 87 not taken.
✓ Branch 88 taken 1 times.
✗ Branch 89 not taken.
✓ Branch 90 taken 1 times.
✗ Branch 91 not taken.
✓ Branch 92 taken 1 times.
✗ Branch 93 not taken.
✓ Branch 94 taken 1 times.
✗ Branch 95 not taken.
✓ Branch 96 taken 1 times.
✗ Branch 97 not taken.
✓ Branch 98 taken 1 times.
✗ Branch 99 not taken.
✓ Branch 100 taken 1 times.
✗ Branch 101 not taken.
✓ Branch 102 taken 1 times.
✗ Branch 103 not taken.
✓ Branch 104 taken 1 times.
✗ Branch 105 not taken.
✓ Branch 106 taken 1 times.
✗ Branch 107 not taken.
✓ Branch 108 taken 1 times.
✗ Branch 109 not taken.
✓ Branch 110 taken 1 times.
✗ Branch 111 not taken.
✓ Branch 112 taken 1 times.
✗ Branch 113 not taken.
✓ Branch 114 taken 1 times.
✗ Branch 115 not taken.
✓ Branch 116 taken 1 times.
✗ Branch 117 not taken.
✓ Branch 118 taken 1 times.
✗ Branch 119 not taken.
✓ Branch 120 taken 1 times.
✗ Branch 121 not taken.
✓ Branch 122 taken 1 times.
✗ Branch 123 not taken.
✓ Branch 124 taken 1 times.
✗ Branch 125 not taken.
✓ Branch 126 taken 1 times.
✗ Branch 127 not taken.
✓ Branch 128 taken 1 times.
✗ Branch 129 not taken.
✓ Branch 130 taken 1 times.
✗ Branch 131 not taken.
✓ Branch 132 taken 1 times.
✗ Branch 133 not taken.
✓ Branch 134 taken 1 times.
✗ Branch 135 not taken.
✓ Branch 136 taken 1 times.
✗ Branch 137 not taken.
✓ Branch 138 taken 1 times.
✗ Branch 139 not taken.
✓ Branch 140 taken 1 times.
✗ Branch 141 not taken.
✓ Branch 142 taken 1 times.
✗ Branch 143 not taken.
 11922 return gridGeometry().elementMapper().size();
783 }
784
785 const GridVariables& gridVariables_;
786 const SolutionVector& sol_;
787
788 std::vector<VolVarScalarDataInfo> volVarScalarDataInfo_; //!< Registered volume variables (scalar)
789 std::vector<VolVarVectorDataInfo> volVarVectorDataInfo_; //!< Registered volume variables (vector)
790
791 std::shared_ptr<VelocityOutput> velocityOutput_; //!< The velocity output policy
792 bool enableVelocityOutput_ = false;
793 bool addProcessRank_ = true;
794 };
795
796 } // end namespace Dumux
797
798 #endif
799