1 // Copyright (C) 2007-2024 CEA, EDF, OPEN CASCADE
3 // Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
6 // This library is free software; you can redistribute it and/or
7 // modify it under the terms of the GNU Lesser General Public
8 // License as published by the Free Software Foundation; either
9 // version 2.1 of the License, or (at your option) any later version.
11 // This library is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 // Lesser General Public License for more details.
16 // You should have received a copy of the GNU Lesser General Public
17 // License along with this library; if not, write to the Free Software
18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
23 // SMESH SMDS : implementation of Salome mesh data structure
24 // File : SMDS_Mesh.hxx
27 #ifndef _SMDS_Mesh_HeaderFile
28 #define _SMDS_Mesh_HeaderFile
30 #include "SMESH_SMDS.hxx"
32 #include "SMDS_BallElement.hxx"
33 #include "SMDS_ElemIterator.hxx"
34 #include "SMDS_Mesh0DElement.hxx"
35 #include "SMDS_MeshCell.hxx"
36 #include "SMDS_MeshEdge.hxx"
37 #include "SMDS_MeshFace.hxx"
38 #include "SMDS_MeshInfo.hxx"
39 #include "SMDS_MeshNode.hxx"
40 #include "SMDS_MeshVolume.hxx"
41 #include "SMDS_UnstructuredGrid.hxx"
46 #include <smIdType.hxx>
48 class SMDS_ElementHolder;
49 class SMDS_ElementFactory;
50 class SMDS_NodeFactory;
52 class SMDS_EXPORT SMDS_Mesh : public SMDS_MeshObject
58 //! actual nodes coordinates, cells definition and reverse connectivity are stored in a vtkUnstructuredGrid
59 inline SMDS_UnstructuredGrid* GetGrid() { return myGrid; }
61 virtual SMDS_NodeIteratorPtr nodesIterator () const;
62 virtual SMDS_EdgeIteratorPtr edgesIterator () const;
63 virtual SMDS_FaceIteratorPtr facesIterator () const;
64 virtual SMDS_VolumeIteratorPtr volumesIterator() const;
66 virtual SMDS_ElemIteratorPtr elementsIterator(SMDSAbs_ElementType type=SMDSAbs_All) const;
67 virtual SMDS_ElemIteratorPtr elementGeomIterator(SMDSAbs_GeometryType type) const;
68 virtual SMDS_ElemIteratorPtr elementEntityIterator(SMDSAbs_EntityType type) const;
70 virtual SMDS_NodeIteratorPtr shapeNodesIterator (int shapeID,
71 size_t nbElemsToReturn=-1,
72 const SMDS_MeshNode* sm1stNode=0) const;
73 virtual SMDS_ElemIteratorPtr shapeElementsIterator(int shapeID,
74 size_t nbElemsToReturn=-1,
75 const SMDS_MeshElement* sm1stElem=0) const;
77 SMDSAbs_ElementType GetElementType( const smIdType id, const bool iselem ) const;
79 SMDS_Mesh *AddSubMesh();
81 virtual SMDS_MeshNode* AddNodeWithID(double x, double y, double z, smIdType ID);
82 virtual SMDS_MeshNode* AddNode (double x, double y, double z);
84 virtual SMDS_Mesh0DElement* Add0DElementWithID(smIdType n, smIdType ID);
85 virtual SMDS_Mesh0DElement* Add0DElementWithID(const SMDS_MeshNode * n, smIdType ID);
86 virtual SMDS_Mesh0DElement* Add0DElement (const SMDS_MeshNode * n);
88 virtual SMDS_BallElement* AddBallWithID(smIdType n, double diameter, smIdType ID);
89 virtual SMDS_BallElement* AddBallWithID(const SMDS_MeshNode * n, double diameter, smIdType ID);
90 virtual SMDS_BallElement* AddBall (const SMDS_MeshNode * n, double diameter);
92 virtual SMDS_MeshEdge* AddEdgeWithID(smIdType n1, smIdType n2, smIdType ID);
93 virtual SMDS_MeshEdge* AddEdgeWithID(const SMDS_MeshNode * n1,
94 const SMDS_MeshNode * n2,
96 virtual SMDS_MeshEdge* AddEdge(const SMDS_MeshNode * n1,
97 const SMDS_MeshNode * n2);
99 // 2d order edge with 3 nodes: n12 - node between n1 and n2
100 virtual SMDS_MeshEdge* AddEdgeWithID(smIdType n1, smIdType n2, smIdType n12, smIdType ID);
101 virtual SMDS_MeshEdge* AddEdgeWithID(const SMDS_MeshNode * n1,
102 const SMDS_MeshNode * n2,
103 const SMDS_MeshNode * n12,
105 virtual SMDS_MeshEdge* AddEdge(const SMDS_MeshNode * n1,
106 const SMDS_MeshNode * n2,
107 const SMDS_MeshNode * n12);
109 virtual SMDS_MeshFace* AddFaceWithID(smIdType n1, smIdType n2, smIdType n3, smIdType ID);
110 virtual SMDS_MeshFace* AddFaceWithID(const SMDS_MeshNode * n1,
111 const SMDS_MeshNode * n2,
112 const SMDS_MeshNode * n3,
114 virtual SMDS_MeshFace* AddFace(const SMDS_MeshNode * n1,
115 const SMDS_MeshNode * n2,
116 const SMDS_MeshNode * n3);
118 virtual SMDS_MeshFace* AddFaceWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4, smIdType ID);
119 virtual SMDS_MeshFace* AddFaceWithID(const SMDS_MeshNode * n1,
120 const SMDS_MeshNode * n2,
121 const SMDS_MeshNode * n3,
122 const SMDS_MeshNode * n4,
124 virtual SMDS_MeshFace* AddFace(const SMDS_MeshNode * n1,
125 const SMDS_MeshNode * n2,
126 const SMDS_MeshNode * n3,
127 const SMDS_MeshNode * n4);
130 // 2d order triangle of 6 nodes
131 virtual SMDS_MeshFace* AddFaceWithID(smIdType n1, smIdType n2, smIdType n3,
132 smIdType n12,smIdType n23,smIdType n31, smIdType ID);
133 virtual SMDS_MeshFace* AddFaceWithID(const SMDS_MeshNode * n1,
134 const SMDS_MeshNode * n2,
135 const SMDS_MeshNode * n3,
136 const SMDS_MeshNode * n12,
137 const SMDS_MeshNode * n23,
138 const SMDS_MeshNode * n31,
140 virtual SMDS_MeshFace* AddFace(const SMDS_MeshNode * n1,
141 const SMDS_MeshNode * n2,
142 const SMDS_MeshNode * n3,
143 const SMDS_MeshNode * n12,
144 const SMDS_MeshNode * n23,
145 const SMDS_MeshNode * n31);
147 // 2d order triangle of 7 nodes
148 virtual SMDS_MeshFace* AddFaceWithID(smIdType n1, smIdType n2, smIdType n3,
149 smIdType n12,smIdType n23,smIdType n31, smIdType nCenter, smIdType ID);
150 virtual SMDS_MeshFace* AddFaceWithID(const SMDS_MeshNode * n1,
151 const SMDS_MeshNode * n2,
152 const SMDS_MeshNode * n3,
153 const SMDS_MeshNode * n12,
154 const SMDS_MeshNode * n23,
155 const SMDS_MeshNode * n31,
156 const SMDS_MeshNode * nCenter,
158 virtual SMDS_MeshFace* AddFace(const SMDS_MeshNode * n1,
159 const SMDS_MeshNode * n2,
160 const SMDS_MeshNode * n3,
161 const SMDS_MeshNode * n12,
162 const SMDS_MeshNode * n23,
163 const SMDS_MeshNode * n31,
164 const SMDS_MeshNode * nCenter);
166 // 2d order quadrangle
167 virtual SMDS_MeshFace* AddFaceWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4,
168 smIdType n12,smIdType n23,smIdType n34,smIdType n41, smIdType ID);
169 virtual SMDS_MeshFace* AddFaceWithID(const SMDS_MeshNode * n1,
170 const SMDS_MeshNode * n2,
171 const SMDS_MeshNode * n3,
172 const SMDS_MeshNode * n4,
173 const SMDS_MeshNode * n12,
174 const SMDS_MeshNode * n23,
175 const SMDS_MeshNode * n34,
176 const SMDS_MeshNode * n41,
178 virtual SMDS_MeshFace* AddFace(const SMDS_MeshNode * n1,
179 const SMDS_MeshNode * n2,
180 const SMDS_MeshNode * n3,
181 const SMDS_MeshNode * n4,
182 const SMDS_MeshNode * n12,
183 const SMDS_MeshNode * n23,
184 const SMDS_MeshNode * n34,
185 const SMDS_MeshNode * n41);
187 virtual SMDS_MeshFace* AddFaceWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4,
188 smIdType n12,smIdType n23,smIdType n34,smIdType n41, smIdType nCenter, smIdType ID);
189 virtual SMDS_MeshFace* AddFaceWithID(const SMDS_MeshNode * n1,
190 const SMDS_MeshNode * n2,
191 const SMDS_MeshNode * n3,
192 const SMDS_MeshNode * n4,
193 const SMDS_MeshNode * n12,
194 const SMDS_MeshNode * n23,
195 const SMDS_MeshNode * n34,
196 const SMDS_MeshNode * n41,
197 const SMDS_MeshNode * nCenter,
199 virtual SMDS_MeshFace* AddFace(const SMDS_MeshNode * n1,
200 const SMDS_MeshNode * n2,
201 const SMDS_MeshNode * n3,
202 const SMDS_MeshNode * n4,
203 const SMDS_MeshNode * n12,
204 const SMDS_MeshNode * n23,
205 const SMDS_MeshNode * n34,
206 const SMDS_MeshNode * n41,
207 const SMDS_MeshNode * nCenter);
209 virtual SMDS_MeshVolume* AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4, smIdType ID);
210 virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
211 const SMDS_MeshNode * n2,
212 const SMDS_MeshNode * n3,
213 const SMDS_MeshNode * n4,
215 virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
216 const SMDS_MeshNode * n2,
217 const SMDS_MeshNode * n3,
218 const SMDS_MeshNode * n4);
220 virtual SMDS_MeshVolume* AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4,
221 smIdType n5, smIdType ID);
222 virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
223 const SMDS_MeshNode * n2,
224 const SMDS_MeshNode * n3,
225 const SMDS_MeshNode * n4,
226 const SMDS_MeshNode * n5,
228 virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
229 const SMDS_MeshNode * n2,
230 const SMDS_MeshNode * n3,
231 const SMDS_MeshNode * n4,
232 const SMDS_MeshNode * n5);
234 virtual SMDS_MeshVolume* AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4,
235 smIdType n5, smIdType n6, smIdType ID);
236 virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
237 const SMDS_MeshNode * n2,
238 const SMDS_MeshNode * n3,
239 const SMDS_MeshNode * n4,
240 const SMDS_MeshNode * n5,
241 const SMDS_MeshNode * n6,
243 virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
244 const SMDS_MeshNode * n2,
245 const SMDS_MeshNode * n3,
246 const SMDS_MeshNode * n4,
247 const SMDS_MeshNode * n5,
248 const SMDS_MeshNode * n6);
250 virtual SMDS_MeshVolume* AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4,
251 smIdType n5, smIdType n6, smIdType n7, smIdType n8, smIdType ID);
252 virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
253 const SMDS_MeshNode * n2,
254 const SMDS_MeshNode * n3,
255 const SMDS_MeshNode * n4,
256 const SMDS_MeshNode * n5,
257 const SMDS_MeshNode * n6,
258 const SMDS_MeshNode * n7,
259 const SMDS_MeshNode * n8,
261 virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
262 const SMDS_MeshNode * n2,
263 const SMDS_MeshNode * n3,
264 const SMDS_MeshNode * n4,
265 const SMDS_MeshNode * n5,
266 const SMDS_MeshNode * n6,
267 const SMDS_MeshNode * n7,
268 const SMDS_MeshNode * n8);
272 virtual SMDS_MeshVolume* AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4, smIdType n5, smIdType n6,
273 smIdType n7, smIdType n8, smIdType n9, smIdType n10, smIdType n11, smIdType n12,
275 virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
276 const SMDS_MeshNode * n2,
277 const SMDS_MeshNode * n3,
278 const SMDS_MeshNode * n4,
279 const SMDS_MeshNode * n5,
280 const SMDS_MeshNode * n6,
281 const SMDS_MeshNode * n7,
282 const SMDS_MeshNode * n8,
283 const SMDS_MeshNode * n9,
284 const SMDS_MeshNode * n10,
285 const SMDS_MeshNode * n11,
286 const SMDS_MeshNode * n12,
288 virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
289 const SMDS_MeshNode * n2,
290 const SMDS_MeshNode * n3,
291 const SMDS_MeshNode * n4,
292 const SMDS_MeshNode * n5,
293 const SMDS_MeshNode * n6,
294 const SMDS_MeshNode * n7,
295 const SMDS_MeshNode * n8,
296 const SMDS_MeshNode * n9,
297 const SMDS_MeshNode * n10,
298 const SMDS_MeshNode * n11,
299 const SMDS_MeshNode * n12);
301 // 2d order tetrahedron of 10 nodes
302 virtual SMDS_MeshVolume* AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4,
303 smIdType n12,smIdType n23,smIdType n31,
304 smIdType n14,smIdType n24,smIdType n34, smIdType ID);
305 virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
306 const SMDS_MeshNode * n2,
307 const SMDS_MeshNode * n3,
308 const SMDS_MeshNode * n4,
309 const SMDS_MeshNode * n12,
310 const SMDS_MeshNode * n23,
311 const SMDS_MeshNode * n31,
312 const SMDS_MeshNode * n14,
313 const SMDS_MeshNode * n24,
314 const SMDS_MeshNode * n34,
316 virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
317 const SMDS_MeshNode * n2,
318 const SMDS_MeshNode * n3,
319 const SMDS_MeshNode * n4,
320 const SMDS_MeshNode * n12,
321 const SMDS_MeshNode * n23,
322 const SMDS_MeshNode * n31,
323 const SMDS_MeshNode * n14,
324 const SMDS_MeshNode * n24,
325 const SMDS_MeshNode * n34);
327 // 2d order pyramid of 13 nodes
328 virtual SMDS_MeshVolume* AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4, smIdType n5,
329 smIdType n12,smIdType n23,smIdType n34,smIdType n41,
330 smIdType n15,smIdType n25,smIdType n35,smIdType n45,
332 virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
333 const SMDS_MeshNode * n2,
334 const SMDS_MeshNode * n3,
335 const SMDS_MeshNode * n4,
336 const SMDS_MeshNode * n5,
337 const SMDS_MeshNode * n12,
338 const SMDS_MeshNode * n23,
339 const SMDS_MeshNode * n34,
340 const SMDS_MeshNode * n41,
341 const SMDS_MeshNode * n15,
342 const SMDS_MeshNode * n25,
343 const SMDS_MeshNode * n35,
344 const SMDS_MeshNode * n45,
346 virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
347 const SMDS_MeshNode * n2,
348 const SMDS_MeshNode * n3,
349 const SMDS_MeshNode * n4,
350 const SMDS_MeshNode * n5,
351 const SMDS_MeshNode * n12,
352 const SMDS_MeshNode * n23,
353 const SMDS_MeshNode * n34,
354 const SMDS_MeshNode * n41,
355 const SMDS_MeshNode * n15,
356 const SMDS_MeshNode * n25,
357 const SMDS_MeshNode * n35,
358 const SMDS_MeshNode * n45);
360 // 2d order Pentahedron with 15 nodes
361 virtual SMDS_MeshVolume* AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3,
362 smIdType n4, smIdType n5, smIdType n6,
363 smIdType n12,smIdType n23,smIdType n31,
364 smIdType n45,smIdType n56,smIdType n64,
365 smIdType n14,smIdType n25,smIdType n36,
367 virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
368 const SMDS_MeshNode * n2,
369 const SMDS_MeshNode * n3,
370 const SMDS_MeshNode * n4,
371 const SMDS_MeshNode * n5,
372 const SMDS_MeshNode * n6,
373 const SMDS_MeshNode * n12,
374 const SMDS_MeshNode * n23,
375 const SMDS_MeshNode * n31,
376 const SMDS_MeshNode * n45,
377 const SMDS_MeshNode * n56,
378 const SMDS_MeshNode * n64,
379 const SMDS_MeshNode * n14,
380 const SMDS_MeshNode * n25,
381 const SMDS_MeshNode * n36,
383 virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
384 const SMDS_MeshNode * n2,
385 const SMDS_MeshNode * n3,
386 const SMDS_MeshNode * n4,
387 const SMDS_MeshNode * n5,
388 const SMDS_MeshNode * n6,
389 const SMDS_MeshNode * n12,
390 const SMDS_MeshNode * n23,
391 const SMDS_MeshNode * n31,
392 const SMDS_MeshNode * n45,
393 const SMDS_MeshNode * n56,
394 const SMDS_MeshNode * n64,
395 const SMDS_MeshNode * n14,
396 const SMDS_MeshNode * n25,
397 const SMDS_MeshNode * n36);
399 // 2d order Pentahedron with 18 nodes
400 virtual SMDS_MeshVolume* AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3,
401 smIdType n4, smIdType n5, smIdType n6,
402 smIdType n12,smIdType n23,smIdType n31,
403 smIdType n45,smIdType n56,smIdType n64,
404 smIdType n14,smIdType n25,smIdType n36,
405 smIdType n1245, smIdType n2356, smIdType n1346,
407 virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
408 const SMDS_MeshNode * n2,
409 const SMDS_MeshNode * n3,
410 const SMDS_MeshNode * n4,
411 const SMDS_MeshNode * n5,
412 const SMDS_MeshNode * n6,
413 const SMDS_MeshNode * n12,
414 const SMDS_MeshNode * n23,
415 const SMDS_MeshNode * n31,
416 const SMDS_MeshNode * n45,
417 const SMDS_MeshNode * n56,
418 const SMDS_MeshNode * n64,
419 const SMDS_MeshNode * n14,
420 const SMDS_MeshNode * n25,
421 const SMDS_MeshNode * n36,
422 const SMDS_MeshNode * n1245,
423 const SMDS_MeshNode * n2356,
424 const SMDS_MeshNode * n1346,
426 virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
427 const SMDS_MeshNode * n2,
428 const SMDS_MeshNode * n3,
429 const SMDS_MeshNode * n4,
430 const SMDS_MeshNode * n5,
431 const SMDS_MeshNode * n6,
432 const SMDS_MeshNode * n12,
433 const SMDS_MeshNode * n23,
434 const SMDS_MeshNode * n31,
435 const SMDS_MeshNode * n45,
436 const SMDS_MeshNode * n56,
437 const SMDS_MeshNode * n64,
438 const SMDS_MeshNode * n14,
439 const SMDS_MeshNode * n25,
440 const SMDS_MeshNode * n36,
441 const SMDS_MeshNode * n1245,
442 const SMDS_MeshNode * n2356,
443 const SMDS_MeshNode * n1346);
446 // 2d oreder Hexahedrons with 20 nodes
447 virtual SMDS_MeshVolume* AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4,
448 smIdType n5, smIdType n6, smIdType n7, smIdType n8,
449 smIdType n12,smIdType n23,smIdType n34,smIdType n41,
450 smIdType n56,smIdType n67,smIdType n78,smIdType n85,
451 smIdType n15,smIdType n26,smIdType n37,smIdType n48,
453 virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
454 const SMDS_MeshNode * n2,
455 const SMDS_MeshNode * n3,
456 const SMDS_MeshNode * n4,
457 const SMDS_MeshNode * n5,
458 const SMDS_MeshNode * n6,
459 const SMDS_MeshNode * n7,
460 const SMDS_MeshNode * n8,
461 const SMDS_MeshNode * n12,
462 const SMDS_MeshNode * n23,
463 const SMDS_MeshNode * n34,
464 const SMDS_MeshNode * n41,
465 const SMDS_MeshNode * n56,
466 const SMDS_MeshNode * n67,
467 const SMDS_MeshNode * n78,
468 const SMDS_MeshNode * n85,
469 const SMDS_MeshNode * n15,
470 const SMDS_MeshNode * n26,
471 const SMDS_MeshNode * n37,
472 const SMDS_MeshNode * n48,
474 virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
475 const SMDS_MeshNode * n2,
476 const SMDS_MeshNode * n3,
477 const SMDS_MeshNode * n4,
478 const SMDS_MeshNode * n5,
479 const SMDS_MeshNode * n6,
480 const SMDS_MeshNode * n7,
481 const SMDS_MeshNode * n8,
482 const SMDS_MeshNode * n12,
483 const SMDS_MeshNode * n23,
484 const SMDS_MeshNode * n34,
485 const SMDS_MeshNode * n41,
486 const SMDS_MeshNode * n56,
487 const SMDS_MeshNode * n67,
488 const SMDS_MeshNode * n78,
489 const SMDS_MeshNode * n85,
490 const SMDS_MeshNode * n15,
491 const SMDS_MeshNode * n26,
492 const SMDS_MeshNode * n37,
493 const SMDS_MeshNode * n48);
495 // 2d oreder Hexahedrons with 27 nodes
496 virtual SMDS_MeshVolume* AddVolumeWithID(smIdType n1, smIdType n2, smIdType n3, smIdType n4,
497 smIdType n5, smIdType n6, smIdType n7, smIdType n8,
498 smIdType n12,smIdType n23,smIdType n34,smIdType n41,
499 smIdType n56,smIdType n67,smIdType n78,smIdType n85,
500 smIdType n15,smIdType n26,smIdType n37,smIdType n48,
501 smIdType n1234,smIdType n1256,smIdType n2367,smIdType n3478,
502 smIdType n1458,smIdType n5678,smIdType nCenter,
504 virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
505 const SMDS_MeshNode * n2,
506 const SMDS_MeshNode * n3,
507 const SMDS_MeshNode * n4,
508 const SMDS_MeshNode * n5,
509 const SMDS_MeshNode * n6,
510 const SMDS_MeshNode * n7,
511 const SMDS_MeshNode * n8,
512 const SMDS_MeshNode * n12,
513 const SMDS_MeshNode * n23,
514 const SMDS_MeshNode * n34,
515 const SMDS_MeshNode * n41,
516 const SMDS_MeshNode * n56,
517 const SMDS_MeshNode * n67,
518 const SMDS_MeshNode * n78,
519 const SMDS_MeshNode * n85,
520 const SMDS_MeshNode * n15,
521 const SMDS_MeshNode * n26,
522 const SMDS_MeshNode * n37,
523 const SMDS_MeshNode * n48,
524 const SMDS_MeshNode * n1234,
525 const SMDS_MeshNode * n1256,
526 const SMDS_MeshNode * n2367,
527 const SMDS_MeshNode * n3478,
528 const SMDS_MeshNode * n1458,
529 const SMDS_MeshNode * n5678,
530 const SMDS_MeshNode * nCenter,
532 virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
533 const SMDS_MeshNode * n2,
534 const SMDS_MeshNode * n3,
535 const SMDS_MeshNode * n4,
536 const SMDS_MeshNode * n5,
537 const SMDS_MeshNode * n6,
538 const SMDS_MeshNode * n7,
539 const SMDS_MeshNode * n8,
540 const SMDS_MeshNode * n12,
541 const SMDS_MeshNode * n23,
542 const SMDS_MeshNode * n34,
543 const SMDS_MeshNode * n41,
544 const SMDS_MeshNode * n56,
545 const SMDS_MeshNode * n67,
546 const SMDS_MeshNode * n78,
547 const SMDS_MeshNode * n85,
548 const SMDS_MeshNode * n15,
549 const SMDS_MeshNode * n26,
550 const SMDS_MeshNode * n37,
551 const SMDS_MeshNode * n48,
552 const SMDS_MeshNode * n1234,
553 const SMDS_MeshNode * n1256,
554 const SMDS_MeshNode * n2367,
555 const SMDS_MeshNode * n3478,
556 const SMDS_MeshNode * n1458,
557 const SMDS_MeshNode * n5678,
558 const SMDS_MeshNode * nCenter);
560 virtual SMDS_MeshFace* AddPolygonalFaceWithID (const std::vector<smIdType> & nodes_ids,
563 virtual SMDS_MeshFace* AddPolygonalFaceWithID (const std::vector<const SMDS_MeshNode*> & nodes,
566 virtual SMDS_MeshFace* AddPolygonalFace (const std::vector<const SMDS_MeshNode*> & nodes);
568 virtual SMDS_MeshFace* AddQuadPolygonalFaceWithID(const std::vector<smIdType> & nodes_ids,
571 virtual SMDS_MeshFace* AddQuadPolygonalFaceWithID(const std::vector<const SMDS_MeshNode*> & nodes,
574 virtual SMDS_MeshFace* AddQuadPolygonalFace(const std::vector<const SMDS_MeshNode*> & nodes);
576 virtual SMDS_MeshVolume* AddPolyhedralVolumeWithID
577 (const std::vector<smIdType> & nodes_ids,
578 const std::vector<int> & quantities,
581 virtual SMDS_MeshVolume* AddPolyhedralVolumeWithID
582 (const std::vector<const SMDS_MeshNode*> & nodes,
583 const std::vector<int> & quantities,
586 virtual SMDS_MeshVolume* AddPolyhedralVolume
587 (const std::vector<const SMDS_MeshNode*> & nodes,
588 const std::vector<int> & quantities);
590 virtual SMDS_MeshVolume* AddVolumeFromVtkIds(const std::vector<vtkIdType>& vtkNodeIds);
592 virtual SMDS_MeshFace* AddFaceFromVtkIds(const std::vector<vtkIdType>& vtkNodeIds);
594 virtual void MoveNode(const SMDS_MeshNode *n, double x, double y, double z);
596 virtual void RemoveElement(const SMDS_MeshElement * elem,
597 std::vector<const SMDS_MeshElement *>& removedElems,
598 std::vector<const SMDS_MeshElement *>& removedNodes,
599 const bool removenodes = false);
600 virtual void RemoveElement(const SMDS_MeshElement * elem, bool removenodes = false);
601 virtual void RemoveNode(const SMDS_MeshNode * node);
603 /*! Remove only the given element and only if it is free.
604 * Method does not work for meshes with descendants.
605 * Implemented for fast cleaning of meshes.
607 virtual void RemoveFreeElement(const SMDS_MeshElement * elem);
609 virtual void Clear();
611 virtual bool RemoveFromParent();
612 virtual bool RemoveSubMesh(const SMDS_Mesh * aMesh);
614 bool ChangeElementNodes(const SMDS_MeshElement * elem,
615 const SMDS_MeshNode * nodes[],
617 bool ChangePolyhedronNodes(const SMDS_MeshElement * elem,
618 const std::vector<const SMDS_MeshNode*>& nodes,
619 const std::vector<int>& quantities);
621 void SetAllNodesNotMarked();
622 void SetAllCellsNotMarked();
624 //virtual void Renumber (const bool isNodes, const int startID = 1, const int deltaID = 1);
625 // Renumber all nodes or elements.
627 virtual void CompactMesh();
628 virtual bool IsCompacted();
629 virtual bool HasNumerationHoles();
631 template<class ELEMTYPE>
632 static const ELEMTYPE* DownCast( const SMDS_MeshElement* e )
634 return (( e && !e->IsNull() && ELEMTYPE::Type() == e->GetType() ) ?
635 static_cast<const ELEMTYPE*>(e) : 0 );
638 const SMDS_MeshNode *FindNode(smIdType idnode) const;
639 const SMDS_MeshNode *FindNodeVtk(vtkIdType idnode) const;
640 const SMDS_MeshElement *FindElementVtk(vtkIdType IDelem) const;
641 virtual const SMDS_MeshElement * FindElement(smIdType IDelem) const;
642 static const SMDS_Mesh0DElement* Find0DElement(const SMDS_MeshNode * n);
643 static const SMDS_BallElement* FindBall(const SMDS_MeshNode * n);
644 static const SMDS_MeshEdge* FindEdge(const SMDS_MeshNode * n1,
645 const SMDS_MeshNode * n2);
646 static const SMDS_MeshEdge* FindEdge(const SMDS_MeshNode * n1,
647 const SMDS_MeshNode * n2,
648 const SMDS_MeshNode * n3);
649 static const SMDS_MeshFace* FindFace(const SMDS_MeshNode *n1,
650 const SMDS_MeshNode *n2,
651 const SMDS_MeshNode *n3);
652 static const SMDS_MeshFace* FindFace(const SMDS_MeshNode *n1,
653 const SMDS_MeshNode *n2,
654 const SMDS_MeshNode *n3,
655 const SMDS_MeshNode *n4);
656 static const SMDS_MeshFace* FindFace(const SMDS_MeshNode *n1,
657 const SMDS_MeshNode *n2,
658 const SMDS_MeshNode *n3,
659 const SMDS_MeshNode *n4,
660 const SMDS_MeshNode *n5,
661 const SMDS_MeshNode *n6);
662 static const SMDS_MeshFace* FindFace(const SMDS_MeshNode *n1,
663 const SMDS_MeshNode *n2,
664 const SMDS_MeshNode *n3,
665 const SMDS_MeshNode *n4,
666 const SMDS_MeshNode *n5,
667 const SMDS_MeshNode *n6,
668 const SMDS_MeshNode *n7,
669 const SMDS_MeshNode *n8);
671 static const SMDS_MeshFace* FindFace (const std::vector<const SMDS_MeshNode *>& nodes);
672 static const SMDS_MeshElement* FindElement(const std::vector<const SMDS_MeshNode *>& nodes,
673 const SMDSAbs_ElementType type=SMDSAbs_All,
674 const bool noMedium=true);
675 static int GetElementsByNodes(const std::vector<const SMDS_MeshNode *>& nodes,
676 std::vector<const SMDS_MeshElement *>& foundElems,
677 const SMDSAbs_ElementType type=SMDSAbs_All);
679 virtual bool Contains( const SMDS_MeshElement* elem ) const;
682 * \brief Raise an exception if free memory (ram+swap) too low
683 * \param doNotRaise - if true, suppress exception, just return free memory size
684 * \retval int - amount of available memory in MB or negative number in failure case
686 static int CheckMemory(const bool doNotRaise=false);
688 virtual smIdType MaxNodeID() const;
689 virtual smIdType MinNodeID() const;
690 virtual smIdType MaxElementID() const;
691 virtual smIdType MinElementID() const;
693 const SMDS_MeshInfo& GetMeshInfo() const { return myInfo; }
695 virtual smIdType NbNodes() const;
696 virtual smIdType NbElements() const;
697 virtual smIdType Nb0DElements() const;
698 virtual smIdType NbBalls() const;
699 virtual smIdType NbEdges() const;
700 virtual smIdType NbFaces() const;
701 virtual smIdType NbVolumes() const;
702 virtual smIdType NbSubMesh() const;
704 virtual ~SMDS_Mesh();
707 smIdType FromVtkToSmds(vtkIdType vtkid) const;
709 void dumpGrid(std::string ficdump="dumpGrid");
710 static int chunkSize;
712 //! low level modification: add, change or remove node or element
713 inline void setMyModified() { this->myModified = true; }
716 vtkMTimeType GetMTime() const;
719 SMDS_Mesh(SMDS_Mesh * parent);
721 void addChildrenWithNodes(std::set<const SMDS_MeshElement*>& setOfChildren,
722 const SMDS_MeshElement * element,
723 std::set<const SMDS_MeshElement*>& nodes);
725 inline void adjustBoundingBox(double x, double y, double z)
727 if (x > xmax) xmax = x;
728 else if (x < xmin) xmin = x;
729 if (y > ymax) ymax = y;
730 else if (y < ymin) ymin = y;
731 if (z > zmax) zmax = z;
732 else if (z < zmin) zmin = z;
735 void updateInverseElements( const SMDS_MeshElement * element,
736 const SMDS_MeshNode* const* nodes,
738 std::set<const SMDS_MeshNode*>& oldNodes );
740 void setNbShapes( size_t nbShapes );
744 //! actual nodes coordinates, cells definition and reverse connectivity are stored in a vtkUnstructuredGrid
745 SMDS_UnstructuredGrid* myGrid;
747 //! Small objects like SMDS_MeshNode are allocated by chunks to limit memory costs of new
748 SMDS_NodeFactory* myNodeFactory;
749 SMDS_ElementFactory* myCellFactory;
751 SMDS_Mesh * myParent;
752 std::list<SMDS_Mesh *> myChildren;
753 SMDS_MeshInfo myInfo;
755 //! any add, remove or change of node or cell
757 //! use a counter to keep track of modifications
758 unsigned long myModifTime, myCompactTime;
760 friend class SMDS_ElementHolder;
761 std::set< SMDS_ElementHolder* > myElemHolders;