LehrFEM++ 1.0.0
A simple Finite Element Library for teaching
Loading...
Searching...
No Matches
gmsh_file_v4.cc
1
9#include "gmsh_file_v4.h"
10
11namespace lf::io {
12
14std::ostream& operator<<(std::ostream& stream, GMshFileV4::ElementType et) {
15 switch (et) {
16 default:
17 break;
18 case GMshFileV4::ElementType::EDGE2:
19 stream << "EDGE2";
20 break;
21 case GMshFileV4::ElementType::TRIA3:
22 stream << "TRIA3";
23 break;
24 case GMshFileV4::ElementType::QUAD4:
25 stream << "QUAD4";
26 break;
27 case GMshFileV4::ElementType::TET4:
28 stream << "TET4";
29 break;
30 case GMshFileV4::ElementType::HEX8:
31 stream << "HEX8";
32 break;
33 case GMshFileV4::ElementType::PRISM6:
34 stream << "PRISM6";
35 break;
36 case GMshFileV4::ElementType::PYRAMID5:
37 stream << "PYRAMID5";
38 break;
39 case GMshFileV4::ElementType::EDGE3:
40 stream << "EDGE3";
41 break;
42 case GMshFileV4::ElementType::TRIA6:
43 stream << "TRIA6";
44 break;
45 case GMshFileV4::ElementType::QUAD9:
46 stream << "QUAD9";
47 break;
48 case GMshFileV4::ElementType::TET10:
49 stream << "TET10";
50 break;
51 case GMshFileV4::ElementType::HEX27:
52 stream << "HEX27";
53 break;
54 case GMshFileV4::ElementType::PRISM18:
55 stream << "PRISM18";
56 break;
57 case GMshFileV4::ElementType::PYRAMID14:
58 stream << "PYRAMID14";
59 break;
60 case GMshFileV4::ElementType::POINT:
61 stream << "POINT";
62 break;
63 case GMshFileV4::ElementType::QUAD8:
64 stream << "QUAD8";
65 break;
66 case GMshFileV4::ElementType::HEX20:
67 stream << "HEX20";
68 break;
69 case GMshFileV4::ElementType::PRISM15:
70 stream << "PRISM15";
71 break;
72 case GMshFileV4::ElementType::PYRAMID13:
73 stream << "PYRAMID13";
74 break;
75 case GMshFileV4::ElementType::TRIA9:
76 stream << "TRIA9";
77 break;
78 case GMshFileV4::ElementType::TRIA10:
79 stream << "TRIA10";
80 break;
81 case GMshFileV4::ElementType::TRIA12:
82 stream << "TRIA12";
83 break;
84 case GMshFileV4::ElementType::TRIA15:
85 stream << "TRIA15";
86 break;
87 case GMshFileV4::ElementType::TRIA15_5:
88 stream << "TRIA15_5";
89 break;
90 case GMshFileV4::ElementType::TRIA21:
91 stream << "TRIA21";
92 break;
93 case GMshFileV4::ElementType::EDGE4:
94 stream << "EDGE4";
95 break;
96 case GMshFileV4::ElementType::EDGE5:
97 stream << "EDGE5";
98 break;
99 case GMshFileV4::ElementType::EDGE6:
100 stream << "EDGE6";
101 break;
102 case GMshFileV4::ElementType::TET20:
103 stream << "TET20";
104 break;
105 case GMshFileV4::ElementType::TET35:
106 stream << "TET35";
107 break;
108 case GMshFileV4::ElementType::TET56:
109 stream << "TET56";
110 break;
111 case GMshFileV4::ElementType::HEX64:
112 stream << "HEX64";
113 break;
114 case GMshFileV4::ElementType::HEX125:
115 stream << "HEX125";
116 break;
117 }
118 return stream;
119}
120
122int NumNodes(GMshFileV4::ElementType et) {
123 switch (et) {
124 default:
125 break;
126 case GMshFileV4::ElementType::EDGE2:
127 return 2;
128 case GMshFileV4::ElementType::TRIA3:
129 return 3;
130 case GMshFileV4::ElementType::QUAD4:
131 case GMshFileV4::ElementType::TET4:
132 return 4;
133 case GMshFileV4::ElementType::HEX8:
134 return 8;
135 case GMshFileV4::ElementType::PRISM6:
136 return 6;
137 case GMshFileV4::ElementType::PYRAMID5:
138 return 5;
139 case GMshFileV4::ElementType::EDGE3:
140 return 3;
141 case GMshFileV4::ElementType::TRIA6:
142 return 6;
143 case GMshFileV4::ElementType::QUAD9:
144 return 9;
145 case GMshFileV4::ElementType::TET10:
146 return 10;
147 case GMshFileV4::ElementType::HEX27:
148 return 27;
149 case GMshFileV4::ElementType::PRISM18:
150 return 18;
151 case GMshFileV4::ElementType::PYRAMID14:
152 return 14;
153 case GMshFileV4::ElementType::POINT:
154 return 1;
155 case GMshFileV4::ElementType::QUAD8:
156 return 8;
157 case GMshFileV4::ElementType::HEX20:
158 return 20;
159 case GMshFileV4::ElementType::PRISM15:
160 return 15;
161 case GMshFileV4::ElementType::PYRAMID13:
162 return 13;
163 case GMshFileV4::ElementType::TRIA9:
164 return 9;
165 case GMshFileV4::ElementType::TRIA10:
166 return 10;
167 case GMshFileV4::ElementType::TRIA12:
168 return 12;
169 case GMshFileV4::ElementType::TRIA15:
170 case GMshFileV4::ElementType::TRIA15_5:
171 return 15;
172 case GMshFileV4::ElementType::TRIA21:
173 return 21;
174 case GMshFileV4::ElementType::EDGE4:
175 return 4;
176 case GMshFileV4::ElementType::EDGE5:
177 return 5;
178 case GMshFileV4::ElementType::EDGE6:
179 return 6;
180 case GMshFileV4::ElementType::TET20:
181 return 20;
182 case GMshFileV4::ElementType::TET35:
183 return 35;
184 case GMshFileV4::ElementType::TET56:
185 return 56;
186 case GMshFileV4::ElementType::HEX64:
187 return 64;
188 case GMshFileV4::ElementType::HEX125:
189 return 125;
190 }
191 LF_VERIFY_MSG(false, "unknown Gmsh element type");
192 // Make compiler happy:
193 return 0;
194}
195
196base::RefEl RefElOf(GMshFileV4::ElementType et) {
197 switch (et) {
198 case GMshFileV4::ElementType::POINT:
199 return base::RefEl::kPoint();
200
201 case GMshFileV4::ElementType::EDGE2:
202 case GMshFileV4::ElementType::EDGE3:
203 case GMshFileV4::ElementType::EDGE4:
204 case GMshFileV4::ElementType::EDGE5:
205 case GMshFileV4::ElementType::EDGE6:
206 return base::RefEl::kSegment();
207
208 case GMshFileV4::ElementType::TRIA3:
209 case GMshFileV4::ElementType::TRIA6:
210 case GMshFileV4::ElementType::TRIA9:
211 case GMshFileV4::ElementType::TRIA10:
212 case GMshFileV4::ElementType::TRIA12:
213 case GMshFileV4::ElementType::TRIA15:
214 case GMshFileV4::ElementType::TRIA15_5:
215 case GMshFileV4::ElementType::TRIA21:
216 return base::RefEl::kTria();
217
218 case GMshFileV4::ElementType::QUAD4:
219 case GMshFileV4::ElementType::QUAD8:
220 case GMshFileV4::ElementType::QUAD9:
221 return base::RefEl::kQuad();
222
223 case GMshFileV4::ElementType::TET4:
224 case GMshFileV4::ElementType::HEX8:
225 case GMshFileV4::ElementType::PRISM6:
226 case GMshFileV4::ElementType::PYRAMID5:
227 case GMshFileV4::ElementType::TET10:
228 case GMshFileV4::ElementType::HEX27:
229 case GMshFileV4::ElementType::PRISM18:
230 case GMshFileV4::ElementType::PYRAMID14:
231 case GMshFileV4::ElementType::HEX20:
232 case GMshFileV4::ElementType::PRISM15:
233 case GMshFileV4::ElementType::PYRAMID13:
234 case GMshFileV4::ElementType::TET20:
235 case GMshFileV4::ElementType::TET35:
236 case GMshFileV4::ElementType::TET56:
237 case GMshFileV4::ElementType::HEX64:
238 case GMshFileV4::ElementType::HEX125:
239 default:
240 LF_VERIFY_MSG(
241 false, "Reference element not supported for GmshElement type " << et);
242 }
243}
244
246int DimOf(GMshFileV4::ElementType et) {
247 switch (et) {
248 case GMshFileV4::ElementType::POINT:
249 return 0;
250 case GMshFileV4::ElementType::EDGE2:
251 case GMshFileV4::ElementType::EDGE3:
252 case GMshFileV4::ElementType::EDGE4:
253 case GMshFileV4::ElementType::EDGE5:
254 case GMshFileV4::ElementType::EDGE6:
255 return 1;
256 case GMshFileV4::ElementType::TRIA3:
257 case GMshFileV4::ElementType::QUAD4:
258 case GMshFileV4::ElementType::TRIA6:
259 case GMshFileV4::ElementType::QUAD9:
260 case GMshFileV4::ElementType::QUAD8:
261 case GMshFileV4::ElementType::TRIA9:
262 case GMshFileV4::ElementType::TRIA10:
263 case GMshFileV4::ElementType::TRIA12:
264 case GMshFileV4::ElementType::TRIA15:
265 case GMshFileV4::ElementType::TRIA15_5:
266 case GMshFileV4::ElementType::TRIA21:
267 return 2;
268 case GMshFileV4::ElementType::TET4:
269 case GMshFileV4::ElementType::HEX8:
270 case GMshFileV4::ElementType::PRISM6:
271 case GMshFileV4::ElementType::PYRAMID5:
272 case GMshFileV4::ElementType::TET10:
273 case GMshFileV4::ElementType::HEX27:
274 case GMshFileV4::ElementType::PRISM18:
275 case GMshFileV4::ElementType::PYRAMID14:
276 case GMshFileV4::ElementType::HEX20:
277 case GMshFileV4::ElementType::PRISM15:
278 case GMshFileV4::ElementType::PYRAMID13:
279 case GMshFileV4::ElementType::TET20:
280 case GMshFileV4::ElementType::TET35:
281 case GMshFileV4::ElementType::TET56:
282 case GMshFileV4::ElementType::HEX64:
283 case GMshFileV4::ElementType::HEX125:
284 return 3;
285 default:
286 LF_VERIFY_MSG(false, "Unknown GmshElement Type.");
287 }
288 // Make compiler happy:
289 return -1;
290}
291
292namespace detail {
293
294// defined in gmsh_file_v4_text.cc
295bool ParseGmshFileV4Text(std::string::const_iterator begin,
296 std::string::const_iterator end, GMshFileV4* result);
297
298bool ParseGmshFileV4Binary(std::string::const_iterator begin,
299 std::string::const_iterator end, int one,
300 GMshFileV4* result);
301
302} // namespace detail
303
304GMshFileV4 ReadGmshFileV4(std::string::const_iterator begin,
305 std::string::const_iterator end,
306 const std::string& version, bool is_binary,
307 int size_t_size, int one,
308 const std::string& filename) {
309 LF_VERIFY_MSG(version == "4.1", "Only version 4.1 is supported so far");
310 LF_VERIFY_MSG(size_t_size == sizeof(std::size_t),
311 "size of size_t must be " << sizeof(std::size_t));
312
313 GMshFileV4 result;
314 result.version_number = version;
315 result.is_binary = is_binary;
316 result.size_t_size = size_t_size;
317
318 bool succesful = false;
319 if (!is_binary) {
320 succesful = detail::ParseGmshFileV4Text(begin, end, &result);
321 } else {
322 succesful = detail::ParseGmshFileV4Binary(begin, end, one, &result);
323 }
324
325 LF_VERIFY_MSG(succesful, "Could not parse file " << filename);
326 // LF_VERIFY_MSG(iter == end, "Could not parse all of file " << filename);
327
328 // transpose all periodic matrices because they are read in column-first mode
329 // but gmsh writes them in row-first mode
330 for (auto& p : result.periodic_links) {
331 if (p.affine_transform) {
332 p.affine_transform->transposeInPlace();
333 }
334 }
335
336 return result;
337}
338
339} // namespace lf::io
lf::base::RefEl
Represents a reference element with all its properties.
Definition ref_el.h:109
lf::base::RefEl::kSegment
static constexpr RefEl kSegment()
Returns the (1-dimensional) reference segment.
Definition ref_el.h:153
lf::base::RefEl::kPoint
static constexpr RefEl kPoint()
Returns the (0-dimensional) reference point.
Definition ref_el.h:144
lf::base::RefEl::kTria
static constexpr RefEl kTria()
Returns the reference triangle.
Definition ref_el.h:161
lf::base::RefEl::kQuad
static constexpr RefEl kQuad()
Returns the reference quadrilateral.
Definition ref_el.h:169
lf::io
Mesh input (from file) and output (in various formats) facilities.
Definition gmsh_file_v2.cc:35
lf::io::ReadGmshFileV4
GMshFileV4 ReadGmshFileV4(std::string::const_iterator begin, std::string::const_iterator end, const std::string &version, bool is_binary, int size_t_size, int one, const std::string &filename)
Read a GmshFile with format 4 and return it as an in-memory struct.
Definition gmsh_file_v4.cc:304
lf::io::DimOf
int DimOf(GMshFileV2::ElementType et)
Dimension of the GmshElement type.
Definition gmsh_file_v2.cc:313
lf::io::RefElOf
base::RefEl RefElOf(GMshFileV2::ElementType et)
Reference element type of a GmshElementType.
Definition gmsh_file_v2.cc:263
lf::io::operator<<
std::ostream & operator<<(std::ostream &stream, GMshFileV2::ElementType et)
Output the element type onto the console:
Definition gmsh_file_v2.cc:38
lf::io::NumNodes
size_type NumNodes(GMshFileV2::ElementType et)
Number of nodes that this element type has.
Definition gmsh_file_v2.cc:189
lf::io::GMshFileV4
A representation of a .msh file (V4) in a c++ data structure.
Definition gmsh_file_v4.h:20
lf::io::GMshFileV4::periodic_links
std::vector< PeriodicLink > periodic_links
Definition gmsh_file_v4.h:481
lf::io::GMshFileV4::ElementType
ElementType
All possible element types (see GMSH documentation)
Definition gmsh_file_v4.h:317
lf::io::GMshFileV4::version_number
std::string version_number
The version of GMSH of the msh file, equals usually 4.1.
Definition gmsh_file_v4.h:22
lf::io::GMshFileV4::size_t_size
int size_t_size
how many bytes is a size_t?
Definition gmsh_file_v4.h:26
lf::io::GMshFileV4::is_binary
bool is_binary
Is it a binary file?
Definition gmsh_file_v4.h:24
Generated by doxygen 1.10.0