Salome HOME
a little optimize FindFace (std::vector<const SMDS_MeshNode *> nodes)
[modules/smesh.git] / src / SMDS / SMDS_Mesh.cxx
1 //  Copyright (C) 2007-2008  CEA/DEN, EDF R&D, OPEN CASCADE
2 //
3 //  Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 //  CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
5 //
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.
10 //
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.
15 //
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
19 //
20 //  See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
21 //
22 //  SMESH SMDS : implementaion of Salome mesh data structure
23
24 #ifdef _MSC_VER
25 #pragma warning(disable:4786)
26 #endif
27
28 #include "utilities.h"
29 #include "SMDS_Mesh.hxx"
30 #include "SMDS_VolumeOfNodes.hxx"
31 #include "SMDS_VolumeOfFaces.hxx"
32 #include "SMDS_FaceOfNodes.hxx"
33 #include "SMDS_FaceOfEdges.hxx"
34 #include "SMDS_PolyhedralVolumeOfNodes.hxx"
35 #include "SMDS_PolygonalFaceOfNodes.hxx"
36 #include "SMDS_QuadraticEdge.hxx"
37 #include "SMDS_QuadraticFaceOfNodes.hxx"
38 #include "SMDS_QuadraticVolumeOfNodes.hxx"
39
40 #include <algorithm>
41 #include <map>
42 using namespace std;
43
44 #ifndef WIN32
45 #include <sys/sysinfo.h>
46 #endif
47
48 // number of added entitis to check memory after
49 #define CHECKMEMORY_INTERVAL 1000
50
51 //================================================================================
52 /*!
53  * \brief Raise an exception if free memory (ram+swap) too low
54  * \param doNotRaise - if true, suppres exception, just return free memory size
55  * \retval int - amount of available memory in MB or negative number in failure case
56  */
57 //================================================================================
58
59 int SMDS_Mesh::CheckMemory(const bool doNotRaise) throw (std::bad_alloc)
60 {
61 #ifndef WIN32
62   struct sysinfo si;
63   int err = sysinfo( &si );
64   if ( err )
65     return -1;
66
67   static int limit = -1;
68   if ( limit < 0 ) {
69     int status = system("SMDS_MemoryLimit"); // it returns lower limit of free RAM
70     if (status >= 0 ) {
71       limit = WEXITSTATUS(status);
72     }
73     if ( limit < 20 )
74       limit = 20;
75     else
76       limit = int( limit * 1.5 );
77 #ifdef _DEBUG_
78     MESSAGE ( "SMDS_Mesh::CheckMemory() memory limit = " << limit << " MB" );
79 #endif
80   }
81
82   const unsigned long Mbyte = 1024 * 1024;
83   // compute separately to avoid overflow
84   int freeMb =
85     ( si.freeram  * si.mem_unit ) / Mbyte +
86     ( si.freeswap * si.mem_unit ) / Mbyte;
87
88   if ( freeMb > limit )
89     return freeMb - limit;
90
91   if ( doNotRaise )
92     return 0;
93 #ifdef _DEBUG_
94   MESSAGE ("SMDS_Mesh::CheckMemory() throws as free memory too low: " << freeMb <<" MB" );
95 #endif
96   throw std::bad_alloc();
97 #else
98   return -1;
99 #endif
100 }
101
102 ///////////////////////////////////////////////////////////////////////////////
103 /// Create a new mesh object
104 ///////////////////////////////////////////////////////////////////////////////
105 SMDS_Mesh::SMDS_Mesh()
106         :myParent(NULL),
107         myNodeIDFactory(new SMDS_MeshElementIDFactory()),
108         myElementIDFactory(new SMDS_MeshElementIDFactory()),
109         myHasConstructionEdges(false), myHasConstructionFaces(false),
110         myHasInverseElements(true)
111 {
112 }
113
114 ///////////////////////////////////////////////////////////////////////////////
115 /// Create a new child mesh
116 /// Note that the tree structure of SMDS_Mesh seems to be unused in this version
117 /// (2003-09-08) of SMESH
118 ///////////////////////////////////////////////////////////////////////////////
119 SMDS_Mesh::SMDS_Mesh(SMDS_Mesh * parent)
120         :myParent(parent), myNodeIDFactory(parent->myNodeIDFactory),
121         myElementIDFactory(parent->myElementIDFactory),
122         myHasConstructionEdges(false), myHasConstructionFaces(false),
123         myHasInverseElements(true)
124 {
125 }
126
127 ///////////////////////////////////////////////////////////////////////////////
128 ///Create a submesh and add it to the current mesh
129 ///////////////////////////////////////////////////////////////////////////////
130
131 SMDS_Mesh *SMDS_Mesh::AddSubMesh()
132 {
133         SMDS_Mesh *submesh = new SMDS_Mesh(this);
134         myChildren.insert(myChildren.end(), submesh);
135         return submesh;
136 }
137
138 ///////////////////////////////////////////////////////////////////////////////
139 ///create a MeshNode and add it to the current Mesh
140 ///An ID is automatically assigned to the node.
141 ///@return : The created node
142 ///////////////////////////////////////////////////////////////////////////////
143
144 SMDS_MeshNode * SMDS_Mesh::AddNode(double x, double y, double z)
145 {
146   return SMDS_Mesh::AddNodeWithID(x,y,z,myNodeIDFactory->GetFreeID());
147 }
148
149 ///////////////////////////////////////////////////////////////////////////////
150 ///create a MeshNode and add it to the current Mesh
151 ///@param ID : The ID of the MeshNode to create
152 ///@return : The created node or NULL if a node with this ID already exists
153 ///////////////////////////////////////////////////////////////////////////////
154 SMDS_MeshNode * SMDS_Mesh::AddNodeWithID(double x, double y, double z, int ID)
155 {
156   // find the MeshNode corresponding to ID
157   const SMDS_MeshElement *node = myNodeIDFactory->MeshElement(ID);
158   if(!node){
159     if ( myNodes.Extent() % CHECKMEMORY_INTERVAL == 0 ) CheckMemory();
160     SMDS_MeshNode * node=new SMDS_MeshNode(x, y, z);
161     myNodes.Add(node);
162     myNodeIDFactory->BindID(ID,node);
163     myInfo.myNbNodes++;
164     return node;
165   }else
166     return NULL;
167 }
168
169 ///////////////////////////////////////////////////////////////////////////////
170 /// create a Mesh0DElement and add it to the current Mesh
171 /// @return : The created Mesh0DElement
172 ///////////////////////////////////////////////////////////////////////////////
173 SMDS_Mesh0DElement* SMDS_Mesh::Add0DElementWithID(int idnode, int ID)
174 {
175   SMDS_MeshNode * node = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode);
176   if (!node) return NULL;
177   return SMDS_Mesh::Add0DElementWithID(node, ID);
178 }
179
180 ///////////////////////////////////////////////////////////////////////////////
181 /// create a Mesh0DElement and add it to the current Mesh
182 /// @return : The created Mesh0DElement
183 ///////////////////////////////////////////////////////////////////////////////
184 SMDS_Mesh0DElement* SMDS_Mesh::Add0DElement(const SMDS_MeshNode * node)
185 {
186   return SMDS_Mesh::Add0DElementWithID(node, myElementIDFactory->GetFreeID());
187 }
188
189 ///////////////////////////////////////////////////////////////////////////////
190 /// Create a new Mesh0DElement and at it to the mesh
191 /// @param idnode ID of the node
192 /// @param ID ID of the 0D element to create
193 /// @return The created 0D element or NULL if an element with this
194 ///         ID already exists or if input node is not found.
195 ///////////////////////////////////////////////////////////////////////////////
196 SMDS_Mesh0DElement* SMDS_Mesh::Add0DElementWithID(const SMDS_MeshNode * n, int ID)
197 {
198   if (!n) return 0;
199
200   if (my0DElements.Extent() % CHECKMEMORY_INTERVAL == 0) CheckMemory();
201
202   SMDS_Mesh0DElement * el0d = new SMDS_Mesh0DElement(n);
203   if (myElementIDFactory->BindID(ID, el0d)) {
204     SMDS_MeshNode *node = const_cast<SMDS_MeshNode*>(n);
205     node->AddInverseElement(el0d);
206     my0DElements.Add(el0d);
207     myInfo.myNb0DElements++;
208     return el0d;
209   }
210
211   delete el0d;
212   return NULL;
213 }
214
215 ///////////////////////////////////////////////////////////////////////////////
216 /// create a MeshEdge and add it to the current Mesh
217 /// @return : The created MeshEdge
218 ///////////////////////////////////////////////////////////////////////////////
219
220 SMDS_MeshEdge* SMDS_Mesh::AddEdgeWithID(int idnode1, int idnode2, int ID)
221 {
222   SMDS_MeshNode * node1 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode1);
223   SMDS_MeshNode * node2 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode2);
224   if(!node1 || !node2) return NULL;
225   return SMDS_Mesh::AddEdgeWithID(node1, node2, ID);
226 }
227
228 ///////////////////////////////////////////////////////////////////////////////
229 /// create a MeshEdge and add it to the current Mesh
230 /// @return : The created MeshEdge
231 ///////////////////////////////////////////////////////////////////////////////
232
233 SMDS_MeshEdge* SMDS_Mesh::AddEdge(const SMDS_MeshNode * node1,
234                                   const SMDS_MeshNode * node2)
235 {
236   return SMDS_Mesh::AddEdgeWithID(node1, node2, myElementIDFactory->GetFreeID());
237 }
238
239 ///////////////////////////////////////////////////////////////////////////////
240 /// Create a new edge and at it to the mesh
241 /// @param idnode1 ID of the first node
242 /// @param idnode2 ID of the second node
243 /// @param ID ID of the edge to create
244 /// @return The created edge or NULL if an element with this ID already exists or
245 /// if input nodes are not found.
246 ///////////////////////////////////////////////////////////////////////////////
247
248 SMDS_MeshEdge* SMDS_Mesh::AddEdgeWithID(const SMDS_MeshNode * n1,
249                                         const SMDS_MeshNode * n2,
250                                         int ID)
251 {
252   if ( !n1 || !n2 ) return 0;
253
254   if ( myEdges.Extent() % CHECKMEMORY_INTERVAL == 0 ) CheckMemory();
255
256   SMDS_MeshEdge * edge=new SMDS_MeshEdge(n1,n2);
257   if(myElementIDFactory->BindID(ID, edge)) {
258     SMDS_MeshNode *node1,*node2;
259     node1=const_cast<SMDS_MeshNode*>(n1);
260     node2=const_cast<SMDS_MeshNode*>(n2);
261     node1->AddInverseElement(edge);
262     node2->AddInverseElement(edge);
263     myEdges.Add(edge);
264     myInfo.myNbEdges++;
265     return edge;
266   }
267   else {
268     delete edge;
269     return NULL;
270   }
271 }
272
273 ///////////////////////////////////////////////////////////////////////////////
274 /// Add a triangle defined by its nodes. An ID is automatically affected to the
275 /// Created face
276 ///////////////////////////////////////////////////////////////////////////////
277
278 SMDS_MeshFace* SMDS_Mesh::AddFace(const SMDS_MeshNode * n1,
279                                   const SMDS_MeshNode * n2,
280                                   const SMDS_MeshNode * n3)
281 {
282   return SMDS_Mesh::AddFaceWithID(n1,n2,n3, myElementIDFactory->GetFreeID());
283 }
284
285 ///////////////////////////////////////////////////////////////////////////////
286 /// Add a triangle defined by its nodes IDs
287 ///////////////////////////////////////////////////////////////////////////////
288
289 SMDS_MeshFace* SMDS_Mesh::AddFaceWithID(int idnode1, int idnode2, int idnode3, int ID)
290 {
291   SMDS_MeshNode * node1 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode1);
292   SMDS_MeshNode * node2 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode2);
293   SMDS_MeshNode * node3 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode3);
294   if(!node1 || !node2 || !node3) return NULL;
295   return SMDS_Mesh::AddFaceWithID(node1, node2, node3, ID);
296 }
297
298 ///////////////////////////////////////////////////////////////////////////////
299 /// Add a triangle defined by its nodes
300 ///////////////////////////////////////////////////////////////////////////////
301
302 SMDS_MeshFace* SMDS_Mesh::AddFaceWithID(const SMDS_MeshNode * n1,
303                                         const SMDS_MeshNode * n2,
304                                         const SMDS_MeshNode * n3,
305                                         int ID)
306 {
307   SMDS_MeshFace * face=createTriangle(n1, n2, n3);
308
309   if (face && !registerElement(ID, face)) {
310     RemoveElement(face, false);
311     face = NULL;
312   }
313   return face;
314 }
315
316 ///////////////////////////////////////////////////////////////////////////////
317 /// Add a quadrangle defined by its nodes. An ID is automatically affected to the
318 /// created face
319 ///////////////////////////////////////////////////////////////////////////////
320
321 SMDS_MeshFace* SMDS_Mesh::AddFace(const SMDS_MeshNode * n1,
322                                   const SMDS_MeshNode * n2,
323                                   const SMDS_MeshNode * n3,
324                                   const SMDS_MeshNode * n4)
325 {
326   return SMDS_Mesh::AddFaceWithID(n1,n2,n3, n4, myElementIDFactory->GetFreeID());
327 }
328
329 ///////////////////////////////////////////////////////////////////////////////
330 /// Add a quadrangle defined by its nodes IDs
331 ///////////////////////////////////////////////////////////////////////////////
332
333 SMDS_MeshFace* SMDS_Mesh::AddFaceWithID(int idnode1,
334                                         int idnode2,
335                                         int idnode3,
336                                         int idnode4,
337                                         int ID)
338 {
339   SMDS_MeshNode *node1, *node2, *node3, *node4;
340   node1 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode1);
341   node2 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode2);
342   node3 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode3);
343   node4 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode4);
344   if(!node1 || !node2 || !node3 || !node4) return NULL;
345   return SMDS_Mesh::AddFaceWithID(node1, node2, node3, node4, ID);
346 }
347
348 ///////////////////////////////////////////////////////////////////////////////
349 /// Add a quadrangle defined by its nodes
350 ///////////////////////////////////////////////////////////////////////////////
351
352 SMDS_MeshFace* SMDS_Mesh::AddFaceWithID(const SMDS_MeshNode * n1,
353                                         const SMDS_MeshNode * n2,
354                                         const SMDS_MeshNode * n3,
355                                         const SMDS_MeshNode * n4,
356                                         int ID)
357 {
358   SMDS_MeshFace * face=createQuadrangle(n1, n2, n3, n4);
359
360   if (face && !registerElement(ID, face)) {
361     RemoveElement(face, false);
362     face = NULL;
363   }
364   return face;
365 }
366
367 ///////////////////////////////////////////////////////////////////////////////
368 /// Add a triangle defined by its edges. An ID is automatically assigned to the
369 /// Created face
370 ///////////////////////////////////////////////////////////////////////////////
371
372 SMDS_MeshFace* SMDS_Mesh::AddFace(const SMDS_MeshEdge * e1,
373                                   const SMDS_MeshEdge * e2,
374                                   const SMDS_MeshEdge * e3)
375 {
376   if (!hasConstructionEdges())
377     return NULL;
378   return AddFaceWithID(e1,e2,e3, myElementIDFactory->GetFreeID());
379 }
380
381 ///////////////////////////////////////////////////////////////////////////////
382 /// Add a triangle defined by its edges
383 ///////////////////////////////////////////////////////////////////////////////
384
385 SMDS_MeshFace* SMDS_Mesh::AddFaceWithID(const SMDS_MeshEdge * e1,
386                                         const SMDS_MeshEdge * e2,
387                                         const SMDS_MeshEdge * e3,
388                                         int ID)
389 {
390   if (!hasConstructionEdges())
391     return NULL;
392   if ( !e1 || !e2 || !e3 ) return 0;
393
394   if ( myFaces.Extent() % CHECKMEMORY_INTERVAL == 0 ) CheckMemory();
395
396   SMDS_MeshFace * face = new SMDS_FaceOfEdges(e1,e2,e3);
397   myFaces.Add(face);
398   myInfo.myNbTriangles++;
399
400   if (!registerElement(ID, face)) {
401     RemoveElement(face, false);
402     face = NULL;
403   }
404   return face;
405 }
406
407 ///////////////////////////////////////////////////////////////////////////////
408 /// Add a quadrangle defined by its edges. An ID is automatically assigned to the
409 /// Created face
410 ///////////////////////////////////////////////////////////////////////////////
411
412 SMDS_MeshFace* SMDS_Mesh::AddFace(const SMDS_MeshEdge * e1,
413                                   const SMDS_MeshEdge * e2,
414                                   const SMDS_MeshEdge * e3,
415                                   const SMDS_MeshEdge * e4)
416 {
417   if (!hasConstructionEdges())
418     return NULL;
419   return AddFaceWithID(e1,e2,e3,e4, myElementIDFactory->GetFreeID());
420 }
421
422 ///////////////////////////////////////////////////////////////////////////////
423 /// Add a quadrangle defined by its edges
424 ///////////////////////////////////////////////////////////////////////////////
425
426 SMDS_MeshFace* SMDS_Mesh::AddFaceWithID(const SMDS_MeshEdge * e1,
427                                         const SMDS_MeshEdge * e2,
428                                         const SMDS_MeshEdge * e3,
429                                         const SMDS_MeshEdge * e4,
430                                         int ID)
431 {
432   if (!hasConstructionEdges())
433     return NULL;
434   if ( !e1 || !e2 || !e3 || !e4 ) return 0;
435   if ( myFaces.Extent() % CHECKMEMORY_INTERVAL == 0 ) CheckMemory();
436   SMDS_MeshFace * face = new SMDS_FaceOfEdges(e1,e2,e3,e4);
437   myFaces.Add(face);
438   myInfo.myNbQuadrangles++;
439
440   if (!registerElement(ID, face))
441   {
442     RemoveElement(face, false);
443     face = NULL;
444   }
445   return face;
446 }
447
448 ///////////////////////////////////////////////////////////////////////////////
449 ///Create a new tetrahedron and add it to the mesh.
450 ///@return The created tetrahedron
451 ///////////////////////////////////////////////////////////////////////////////
452
453 SMDS_MeshVolume* SMDS_Mesh::AddVolume(const SMDS_MeshNode * n1,
454                                       const SMDS_MeshNode * n2,
455                                       const SMDS_MeshNode * n3,
456                                       const SMDS_MeshNode * n4)
457 {
458   int ID = myElementIDFactory->GetFreeID();
459   SMDS_MeshVolume * v = SMDS_Mesh::AddVolumeWithID(n1, n2, n3, n4, ID);
460   if(v==NULL) myElementIDFactory->ReleaseID(ID);
461   return v;
462 }
463
464 ///////////////////////////////////////////////////////////////////////////////
465 ///Create a new tetrahedron and add it to the mesh.
466 ///@param ID The ID of the new volume
467 ///@return The created tetrahedron or NULL if an element with this ID already exists
468 ///or if input nodes are not found.
469 ///////////////////////////////////////////////////////////////////////////////
470
471 SMDS_MeshVolume * SMDS_Mesh::AddVolumeWithID(int idnode1,
472                                              int idnode2,
473                                              int idnode3,
474                                              int idnode4,
475                                              int ID)
476 {
477   SMDS_MeshNode *node1, *node2, *node3, *node4;
478   node1 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode1);
479   node2 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode2);
480   node3 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode3);
481   node4 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode4);
482   if(!node1 || !node2 || !node3 || !node4) return NULL;
483   return SMDS_Mesh::AddVolumeWithID(node1, node2, node3, node4, ID);
484 }
485
486 ///////////////////////////////////////////////////////////////////////////////
487 ///Create a new tetrahedron and add it to the mesh.
488 ///@param ID The ID of the new volume
489 ///@return The created tetrahedron
490 ///////////////////////////////////////////////////////////////////////////////
491
492 SMDS_MeshVolume* SMDS_Mesh::AddVolumeWithID(const SMDS_MeshNode * n1,
493                                             const SMDS_MeshNode * n2,
494                                             const SMDS_MeshNode * n3,
495                                             const SMDS_MeshNode * n4,
496                                             int ID)
497 {
498   SMDS_MeshVolume* volume = 0;
499   if ( !n1 || !n2 || !n3 || !n4) return volume;
500   if ( myVolumes.Extent() % CHECKMEMORY_INTERVAL == 0 ) CheckMemory();
501   if(hasConstructionFaces()) {
502     SMDS_MeshFace * f1=FindFaceOrCreate(n1,n2,n3);
503     SMDS_MeshFace * f2=FindFaceOrCreate(n1,n2,n4);
504     SMDS_MeshFace * f3=FindFaceOrCreate(n1,n3,n4);
505     SMDS_MeshFace * f4=FindFaceOrCreate(n2,n3,n4);
506     volume=new SMDS_VolumeOfFaces(f1,f2,f3,f4);
507     myVolumes.Add(volume);
508     myInfo.myNbTetras++;
509   }
510   else if(hasConstructionEdges()) {
511     MESSAGE("Error : Not implemented");
512     return NULL;
513   }
514   else {
515     volume=new SMDS_VolumeOfNodes(n1,n2,n3,n4);
516     myVolumes.Add(volume);
517     myInfo.myNbTetras++;
518   }
519
520   if (!registerElement(ID, volume)) {
521     RemoveElement(volume, false);
522     volume = NULL;
523   }
524   return volume;
525 }
526
527 ///////////////////////////////////////////////////////////////////////////////
528 ///Create a new pyramid and add it to the mesh.
529 ///Nodes 1,2,3 and 4 define the base of the pyramid
530 ///@return The created pyramid
531 ///////////////////////////////////////////////////////////////////////////////
532
533 SMDS_MeshVolume* SMDS_Mesh::AddVolume(const SMDS_MeshNode * n1,
534                                       const SMDS_MeshNode * n2,
535                                       const SMDS_MeshNode * n3,
536                                       const SMDS_MeshNode * n4,
537                                       const SMDS_MeshNode * n5)
538 {
539   int ID = myElementIDFactory->GetFreeID();
540   SMDS_MeshVolume * v = SMDS_Mesh::AddVolumeWithID(n1, n2, n3, n4, n5, ID);
541   if(v==NULL) myElementIDFactory->ReleaseID(ID);
542   return v;
543 }
544
545 ///////////////////////////////////////////////////////////////////////////////
546 ///Create a new pyramid and add it to the mesh.
547 ///Nodes 1,2,3 and 4 define the base of the pyramid
548 ///@param ID The ID of the new volume
549 ///@return The created pyramid or NULL if an element with this ID already exists
550 ///or if input nodes are not found.
551 ///////////////////////////////////////////////////////////////////////////////
552
553 SMDS_MeshVolume * SMDS_Mesh::AddVolumeWithID(int idnode1,
554                                              int idnode2,
555                                              int idnode3,
556                                              int idnode4,
557                                              int idnode5,
558                                              int ID)
559 {
560   SMDS_MeshNode *node1, *node2, *node3, *node4, *node5;
561   node1 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode1);
562   node2 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode2);
563   node3 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode3);
564   node4 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode4);
565   node5 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode5);
566   if(!node1 || !node2 || !node3 || !node4 || !node5) return NULL;
567   return SMDS_Mesh::AddVolumeWithID(node1, node2, node3, node4, node5, ID);
568 }
569
570 ///////////////////////////////////////////////////////////////////////////////
571 ///Create a new pyramid and add it to the mesh.
572 ///Nodes 1,2,3 and 4 define the base of the pyramid
573 ///@param ID The ID of the new volume
574 ///@return The created pyramid
575 ///////////////////////////////////////////////////////////////////////////////
576
577 SMDS_MeshVolume* SMDS_Mesh::AddVolumeWithID(const SMDS_MeshNode * n1,
578                                             const SMDS_MeshNode * n2,
579                                             const SMDS_MeshNode * n3,
580                                             const SMDS_MeshNode * n4,
581                                             const SMDS_MeshNode * n5,
582                                             int ID)
583 {
584   SMDS_MeshVolume* volume = 0;
585   if ( !n1 || !n2 || !n3 || !n4 || !n5) return volume;
586   if ( myVolumes.Extent() % CHECKMEMORY_INTERVAL == 0 ) CheckMemory();
587   if(hasConstructionFaces()) {
588     SMDS_MeshFace * f1=FindFaceOrCreate(n1,n2,n3,n4);
589     SMDS_MeshFace * f2=FindFaceOrCreate(n1,n2,n5);
590     SMDS_MeshFace * f3=FindFaceOrCreate(n2,n3,n5);
591     SMDS_MeshFace * f4=FindFaceOrCreate(n3,n4,n5);
592     volume=new SMDS_VolumeOfFaces(f1,f2,f3,f4);
593     myVolumes.Add(volume);
594     myInfo.myNbPyramids++;
595   }
596   else if(hasConstructionEdges()) {
597     MESSAGE("Error : Not implemented");
598     return NULL;
599   }
600   else {
601     volume=new SMDS_VolumeOfNodes(n1,n2,n3,n4,n5);
602     myVolumes.Add(volume);
603     myInfo.myNbPyramids++;
604   }
605
606   if (!registerElement(ID, volume)) {
607     RemoveElement(volume, false);
608     volume = NULL;
609   }
610   return volume;
611 }
612
613 ///////////////////////////////////////////////////////////////////////////////
614 ///Create a new prism and add it to the mesh.
615 ///Nodes 1,2,3 is a triangle and 1,2,5,4 a quadrangle.
616 ///@return The created prism
617 ///////////////////////////////////////////////////////////////////////////////
618
619 SMDS_MeshVolume* SMDS_Mesh::AddVolume(const SMDS_MeshNode * n1,
620                                       const SMDS_MeshNode * n2,
621                                       const SMDS_MeshNode * n3,
622                                       const SMDS_MeshNode * n4,
623                                       const SMDS_MeshNode * n5,
624                                       const SMDS_MeshNode * n6)
625 {
626   int ID = myElementIDFactory->GetFreeID();
627   SMDS_MeshVolume * v = SMDS_Mesh::AddVolumeWithID(n1, n2, n3, n4, n5, n6, ID);
628   if(v==NULL) myElementIDFactory->ReleaseID(ID);
629   return v;
630 }
631
632 ///////////////////////////////////////////////////////////////////////////////
633 ///Create a new prism and add it to the mesh.
634 ///Nodes 1,2,3 is a triangle and 1,2,5,4 a quadrangle.
635 ///@param ID The ID of the new volume
636 ///@return The created prism or NULL if an element with this ID already exists
637 ///or if input nodes are not found.
638 ///////////////////////////////////////////////////////////////////////////////
639
640 SMDS_MeshVolume * SMDS_Mesh::AddVolumeWithID(int idnode1,
641                                              int idnode2,
642                                              int idnode3,
643                                              int idnode4,
644                                              int idnode5,
645                                              int idnode6,
646                                              int ID)
647 {
648   SMDS_MeshNode *node1, *node2, *node3, *node4, *node5, *node6;
649   node1 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode1);
650   node2 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode2);
651   node3 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode3);
652   node4 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode4);
653   node5 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode5);
654   node6 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode6);
655   if(!node1 || !node2 || !node3 || !node4 || !node5 || !node6) return NULL;
656   return SMDS_Mesh::AddVolumeWithID(node1, node2, node3, node4, node5, node6, ID);
657 }
658
659 ///////////////////////////////////////////////////////////////////////////////
660 ///Create a new prism and add it to the mesh.
661 ///Nodes 1,2,3 is a triangle and 1,2,5,4 a quadrangle.
662 ///@param ID The ID of the new volume
663 ///@return The created prism
664 ///////////////////////////////////////////////////////////////////////////////
665
666 SMDS_MeshVolume* SMDS_Mesh::AddVolumeWithID(const SMDS_MeshNode * n1,
667                                             const SMDS_MeshNode * n2,
668                                             const SMDS_MeshNode * n3,
669                                             const SMDS_MeshNode * n4,
670                                             const SMDS_MeshNode * n5,
671                                             const SMDS_MeshNode * n6,
672                                             int ID)
673 {
674   SMDS_MeshVolume* volume = 0;
675   if ( !n1 || !n2 || !n3 || !n4 || !n5 || !n6) return volume;
676   if ( myVolumes.Extent() % CHECKMEMORY_INTERVAL == 0 ) CheckMemory();
677   if(hasConstructionFaces()) {
678     SMDS_MeshFace * f1=FindFaceOrCreate(n1,n2,n3);
679     SMDS_MeshFace * f2=FindFaceOrCreate(n4,n5,n6);
680     SMDS_MeshFace * f3=FindFaceOrCreate(n1,n4,n5,n2);
681     SMDS_MeshFace * f4=FindFaceOrCreate(n2,n5,n6,n3);
682     SMDS_MeshFace * f5=FindFaceOrCreate(n3,n6,n4,n1);
683     volume=new SMDS_VolumeOfFaces(f1,f2,f3,f4,f5);
684     myVolumes.Add(volume);
685     myInfo.myNbPrisms++;
686   }
687   else if(hasConstructionEdges()) {
688     MESSAGE("Error : Not implemented");
689     return NULL;
690   }
691   else {
692     volume=new SMDS_VolumeOfNodes(n1,n2,n3,n4,n5,n6);
693     myVolumes.Add(volume);
694     myInfo.myNbPrisms++;
695   }
696
697   if (!registerElement(ID, volume)) {
698     RemoveElement(volume, false);
699     volume = NULL;
700   }
701   return volume;
702 }
703
704 ///////////////////////////////////////////////////////////////////////////////
705 ///Create a new hexahedron and add it to the mesh.
706 ///Nodes 1,2,3,4 and 5,6,7,8 are quadrangle and 5,1 and 7,3 are an edges.
707 ///@return The created hexahedron
708 ///////////////////////////////////////////////////////////////////////////////
709
710 SMDS_MeshVolume* SMDS_Mesh::AddVolume(const SMDS_MeshNode * n1,
711                                       const SMDS_MeshNode * n2,
712                                       const SMDS_MeshNode * n3,
713                                       const SMDS_MeshNode * n4,
714                                       const SMDS_MeshNode * n5,
715                                       const SMDS_MeshNode * n6,
716                                       const SMDS_MeshNode * n7,
717                                       const SMDS_MeshNode * n8)
718 {
719   int ID = myElementIDFactory->GetFreeID();
720   SMDS_MeshVolume * v = SMDS_Mesh::AddVolumeWithID(n1, n2, n3, n4, n5, n6, n7, n8, ID);
721   if(v==NULL) myElementIDFactory->ReleaseID(ID);
722   return v;
723 }
724
725 ///////////////////////////////////////////////////////////////////////////////
726 ///Create a new hexahedron and add it to the mesh.
727 ///Nodes 1,2,3,4 and 5,6,7,8 are quadrangle and 5,1 and 7,3 are an edges.
728 ///@param ID The ID of the new volume
729 ///@return The created hexahedron or NULL if an element with this ID already
730 ///exists or if input nodes are not found.
731 ///////////////////////////////////////////////////////////////////////////////
732
733 SMDS_MeshVolume * SMDS_Mesh::AddVolumeWithID(int idnode1,
734                                              int idnode2,
735                                              int idnode3,
736                                              int idnode4,
737                                              int idnode5,
738                                              int idnode6,
739                                              int idnode7,
740                                              int idnode8,
741                                              int ID)
742 {
743   SMDS_MeshNode *node1, *node2, *node3, *node4, *node5, *node6, *node7, *node8;
744   node1 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode1);
745   node2 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode2);
746   node3 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode3);
747   node4 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode4);
748   node5 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode5);
749   node6 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode6);
750   node7 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode7);
751   node8 = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(idnode8);
752   if(!node1 || !node2 || !node3 || !node4 || !node5 || !node6 || !node7 || !node8)
753     return NULL;
754   return SMDS_Mesh::AddVolumeWithID(node1, node2, node3, node4, node5, node6,
755                                     node7, node8, ID);
756 }
757
758 ///////////////////////////////////////////////////////////////////////////////
759 ///Create a new hexahedron and add it to the mesh.
760 ///Nodes 1,2,3,4 and 5,6,7,8 are quadrangle and 5,1 and 7,3 are an edges.
761 ///@param ID The ID of the new volume
762 ///@return The created prism or NULL if an element with this ID already exists
763 ///or if input nodes are not found.
764 ///////////////////////////////////////////////////////////////////////////////
765
766 SMDS_MeshVolume* SMDS_Mesh::AddVolumeWithID(const SMDS_MeshNode * n1,
767                                             const SMDS_MeshNode * n2,
768                                             const SMDS_MeshNode * n3,
769                                             const SMDS_MeshNode * n4,
770                                             const SMDS_MeshNode * n5,
771                                             const SMDS_MeshNode * n6,
772                                             const SMDS_MeshNode * n7,
773                                             const SMDS_MeshNode * n8,
774                                             int ID)
775 {
776   SMDS_MeshVolume* volume = 0;
777   if ( !n1 || !n2 || !n3 || !n4 || !n5 || !n6 || !n7 || !n8) return volume;
778   if ( myVolumes.Extent() % CHECKMEMORY_INTERVAL == 0 ) CheckMemory();
779   if(hasConstructionFaces()) {
780     SMDS_MeshFace * f1=FindFaceOrCreate(n1,n2,n3,n4);
781     SMDS_MeshFace * f2=FindFaceOrCreate(n5,n6,n7,n8);
782     SMDS_MeshFace * f3=FindFaceOrCreate(n1,n4,n8,n5);
783     SMDS_MeshFace * f4=FindFaceOrCreate(n1,n2,n6,n5);
784     SMDS_MeshFace * f5=FindFaceOrCreate(n2,n3,n7,n6);
785     SMDS_MeshFace * f6=FindFaceOrCreate(n3,n4,n8,n7);
786     volume=new SMDS_VolumeOfFaces(f1,f2,f3,f4,f5,f6);
787     myVolumes.Add(volume);
788     myInfo.myNbHexas++;
789   }
790   else if(hasConstructionEdges()) {
791     MESSAGE("Error : Not implemented");
792     return NULL;
793   }
794   else {
795 //    volume=new SMDS_HexahedronOfNodes(n1,n2,n3,n4,n5,n6,n7,n8);
796     volume=new SMDS_VolumeOfNodes(n1,n2,n3,n4,n5,n6,n7,n8);
797     myVolumes.Add(volume);
798     myInfo.myNbHexas++;
799   }
800
801   if (!registerElement(ID, volume)) {
802     RemoveElement(volume, false);
803     volume = NULL;
804   }
805   return volume;
806 }
807
808 ///////////////////////////////////////////////////////////////////////////////
809 ///Create a new tetrahedron defined by its faces and add it to the mesh.
810 ///@return The created tetrahedron
811 ///////////////////////////////////////////////////////////////////////////////
812
813 SMDS_MeshVolume* SMDS_Mesh::AddVolume(const SMDS_MeshFace * f1,
814                                       const SMDS_MeshFace * f2,
815                                       const SMDS_MeshFace * f3,
816                                       const SMDS_MeshFace * f4)
817 {
818   if (!hasConstructionFaces())
819     return NULL;
820   return AddVolumeWithID(f1,f2,f3,f4, myElementIDFactory->GetFreeID());
821 }
822
823 ///////////////////////////////////////////////////////////////////////////////
824 ///Create a new tetrahedron defined by its faces and add it to the mesh.
825 ///@param ID The ID of the new volume
826 ///@return The created tetrahedron
827 ///////////////////////////////////////////////////////////////////////////////
828
829 SMDS_MeshVolume* SMDS_Mesh::AddVolumeWithID(const SMDS_MeshFace * f1,
830                                             const SMDS_MeshFace * f2,
831                                             const SMDS_MeshFace * f3,
832                                             const SMDS_MeshFace * f4,
833                                             int ID)
834 {
835   if (!hasConstructionFaces())
836     return NULL;
837   if ( !f1 || !f2 || !f3 || !f4) return 0;
838   if ( myVolumes.Extent() % CHECKMEMORY_INTERVAL == 0 ) CheckMemory();
839   SMDS_MeshVolume * volume = new SMDS_VolumeOfFaces(f1,f2,f3,f4);
840   myVolumes.Add(volume);
841   myInfo.myNbTetras++;
842
843   if (!registerElement(ID, volume)) {
844     RemoveElement(volume, false);
845     volume = NULL;
846   }
847   return volume;
848 }
849
850 ///////////////////////////////////////////////////////////////////////////////
851 ///Create a new pyramid defined by its faces and add it to the mesh.
852 ///@return The created pyramid
853 ///////////////////////////////////////////////////////////////////////////////
854
855 SMDS_MeshVolume* SMDS_Mesh::AddVolume(const SMDS_MeshFace * f1,
856                                       const SMDS_MeshFace * f2,
857                                       const SMDS_MeshFace * f3,
858                                       const SMDS_MeshFace * f4,
859                                       const SMDS_MeshFace * f5)
860 {
861   if (!hasConstructionFaces())
862     return NULL;
863   return AddVolumeWithID(f1,f2,f3,f4,f5, myElementIDFactory->GetFreeID());
864 }
865
866 ///////////////////////////////////////////////////////////////////////////////
867 ///Create a new pyramid defined by its faces and add it to the mesh.
868 ///@param ID The ID of the new volume
869 ///@return The created pyramid
870 ///////////////////////////////////////////////////////////////////////////////
871
872 SMDS_MeshVolume* SMDS_Mesh::AddVolumeWithID(const SMDS_MeshFace * f1,
873                                             const SMDS_MeshFace * f2,
874                                             const SMDS_MeshFace * f3,
875                                             const SMDS_MeshFace * f4,
876                                             const SMDS_MeshFace * f5,
877                                             int ID)
878 {
879   if (!hasConstructionFaces())
880     return NULL;
881   if ( !f1 || !f2 || !f3 || !f4 || !f5) return 0;
882   if ( myVolumes.Extent() % CHECKMEMORY_INTERVAL == 0 ) CheckMemory();
883   SMDS_MeshVolume * volume = new SMDS_VolumeOfFaces(f1,f2,f3,f4,f5);
884   myVolumes.Add(volume);
885   myInfo.myNbPyramids++;
886
887   if (!registerElement(ID, volume)) {
888     RemoveElement(volume, false);
889     volume = NULL;
890   }
891   return volume;
892 }
893
894 ///////////////////////////////////////////////////////////////////////////////
895 ///Create a new prism defined by its faces and add it to the mesh.
896 ///@return The created prism
897 ///////////////////////////////////////////////////////////////////////////////
898
899 SMDS_MeshVolume* SMDS_Mesh::AddVolume(const SMDS_MeshFace * f1,
900                                       const SMDS_MeshFace * f2,
901                                       const SMDS_MeshFace * f3,
902                                       const SMDS_MeshFace * f4,
903                                       const SMDS_MeshFace * f5,
904                                       const SMDS_MeshFace * f6)
905 {
906   if (!hasConstructionFaces())
907     return NULL;
908   return AddVolumeWithID(f1,f2,f3,f4,f5,f6, myElementIDFactory->GetFreeID());
909 }
910
911 ///////////////////////////////////////////////////////////////////////////////
912 ///Create a new prism defined by its faces and add it to the mesh.
913 ///@param ID The ID of the new volume
914 ///@return The created prism
915 ///////////////////////////////////////////////////////////////////////////////
916
917 SMDS_MeshVolume* SMDS_Mesh::AddVolumeWithID(const SMDS_MeshFace * f1,
918                                             const SMDS_MeshFace * f2,
919                                             const SMDS_MeshFace * f3,
920                                             const SMDS_MeshFace * f4,
921                                             const SMDS_MeshFace * f5,
922                                             const SMDS_MeshFace * f6,
923                                             int ID)
924 {
925   if (!hasConstructionFaces())
926     return NULL;
927   if ( !f1 || !f2 || !f3 || !f4 || !f5 || !f6) return 0;
928   if ( myVolumes.Extent() % CHECKMEMORY_INTERVAL == 0 ) CheckMemory();
929   SMDS_MeshVolume * volume = new SMDS_VolumeOfFaces(f1,f2,f3,f4,f5,f6);
930   myVolumes.Add(volume);
931   myInfo.myNbPrisms++;
932
933   if (!registerElement(ID, volume)) {
934     RemoveElement(volume, false);
935     volume = NULL;
936   }
937   return volume;
938 }
939
940 ///////////////////////////////////////////////////////////////////////////////
941 /// Add a polygon defined by its nodes IDs
942 ///////////////////////////////////////////////////////////////////////////////
943
944 SMDS_MeshFace* SMDS_Mesh::AddPolygonalFaceWithID (std::vector<int> nodes_ids,
945                                                   const int        ID)
946 {
947   int nbNodes = nodes_ids.size();
948   std::vector<const SMDS_MeshNode*> nodes (nbNodes);
949   for (int i = 0; i < nbNodes; i++) {
950     nodes[i] = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(nodes_ids[i]);
951     if (!nodes[i]) return NULL;
952   }
953   return SMDS_Mesh::AddPolygonalFaceWithID(nodes, ID);
954 }
955
956 ///////////////////////////////////////////////////////////////////////////////
957 /// Add a polygon defined by its nodes
958 ///////////////////////////////////////////////////////////////////////////////
959
960 SMDS_MeshFace* SMDS_Mesh::AddPolygonalFaceWithID
961                           (std::vector<const SMDS_MeshNode*> nodes,
962                            const int                         ID)
963 {
964   SMDS_MeshFace * face;
965
966   if ( myFaces.Extent() % CHECKMEMORY_INTERVAL == 0 ) CheckMemory();
967   if (hasConstructionEdges())
968   {
969     MESSAGE("Error : Not implemented");
970     return NULL;
971   }
972   else
973   {
974     for ( int i = 0; i < nodes.size(); ++i )
975       if ( !nodes[ i ] ) return 0;
976     face = new SMDS_PolygonalFaceOfNodes(nodes);
977     myFaces.Add(face);
978     myInfo.myNbPolygons++;
979   }
980
981   if (!registerElement(ID, face)) {
982     RemoveElement(face, false);
983     face = NULL;
984   }
985   return face;
986 }
987
988 ///////////////////////////////////////////////////////////////////////////////
989 /// Add a polygon defined by its nodes.
990 /// An ID is automatically affected to the created face.
991 ///////////////////////////////////////////////////////////////////////////////
992
993 SMDS_MeshFace* SMDS_Mesh::AddPolygonalFace (std::vector<const SMDS_MeshNode*> nodes)
994 {
995   return SMDS_Mesh::AddPolygonalFaceWithID(nodes, myElementIDFactory->GetFreeID());
996 }
997
998 ///////////////////////////////////////////////////////////////////////////////
999 /// Create a new polyhedral volume and add it to the mesh.
1000 /// @param ID The ID of the new volume
1001 /// @return The created volume or NULL if an element with this ID already exists
1002 /// or if input nodes are not found.
1003 ///////////////////////////////////////////////////////////////////////////////
1004
1005 SMDS_MeshVolume * SMDS_Mesh::AddPolyhedralVolumeWithID
1006                              (std::vector<int> nodes_ids,
1007                               std::vector<int> quantities,
1008                               const int        ID)
1009 {
1010   int nbNodes = nodes_ids.size();
1011   std::vector<const SMDS_MeshNode*> nodes (nbNodes);
1012   for (int i = 0; i < nbNodes; i++) {
1013     nodes[i] = (SMDS_MeshNode *)myNodeIDFactory->MeshElement(nodes_ids[i]);
1014     if (!nodes[i]) return NULL;
1015   }
1016   return SMDS_Mesh::AddPolyhedralVolumeWithID(nodes, quantities, ID);
1017 }
1018
1019 ///////////////////////////////////////////////////////////////////////////////
1020 /// Create a new polyhedral volume and add it to the mesh.
1021 /// @param ID The ID of the new volume
1022 /// @return The created  volume
1023 ///////////////////////////////////////////////////////////////////////////////
1024
1025 SMDS_MeshVolume* SMDS_Mesh::AddPolyhedralVolumeWithID
1026                             (std::vector<const SMDS_MeshNode*> nodes,
1027                              std::vector<int>                  quantities,
1028                              const int                         ID)
1029 {
1030   SMDS_MeshVolume* volume;
1031   if ( myVolumes.Extent() % CHECKMEMORY_INTERVAL == 0 ) CheckMemory();
1032   if (hasConstructionFaces()) {
1033     MESSAGE("Error : Not implemented");
1034     return NULL;
1035   } else if (hasConstructionEdges()) {
1036     MESSAGE("Error : Not implemented");
1037     return NULL;
1038   } else {
1039     for ( int i = 0; i < nodes.size(); ++i )
1040       if ( !nodes[ i ] ) return 0;
1041     volume = new SMDS_PolyhedralVolumeOfNodes(nodes, quantities);
1042     myVolumes.Add(volume);
1043     myInfo.myNbPolyhedrons++;
1044   }
1045
1046   if (!registerElement(ID, volume)) {
1047     RemoveElement(volume, false);
1048     volume = NULL;
1049   }
1050   return volume;
1051 }
1052
1053 ///////////////////////////////////////////////////////////////////////////////
1054 /// Create a new polyhedral volume and add it to the mesh.
1055 /// @return The created  volume
1056 ///////////////////////////////////////////////////////////////////////////////
1057
1058 SMDS_MeshVolume* SMDS_Mesh::AddPolyhedralVolume
1059                             (std::vector<const SMDS_MeshNode*> nodes,
1060                              std::vector<int>                  quantities)
1061 {
1062   int ID = myElementIDFactory->GetFreeID();
1063   SMDS_MeshVolume * v = SMDS_Mesh::AddPolyhedralVolumeWithID(nodes, quantities, ID);
1064   if (v == NULL) myElementIDFactory->ReleaseID(ID);
1065   return v;
1066 }
1067
1068 ///////////////////////////////////////////////////////////////////////////////
1069 /// Registers element with the given ID, maintains inverse connections
1070 ///////////////////////////////////////////////////////////////////////////////
1071 bool SMDS_Mesh::registerElement(int ID, SMDS_MeshElement * element)
1072 {
1073   if (myElementIDFactory->BindID(ID, element)) {
1074     SMDS_ElemIteratorPtr it = element->nodesIterator();
1075     while (it->more()) {
1076       SMDS_MeshNode *node = static_cast<SMDS_MeshNode*>
1077         (const_cast<SMDS_MeshElement*>(it->next()));
1078       node->AddInverseElement(element);
1079     }
1080     return true;
1081   }
1082   return false;
1083 }
1084
1085 ///////////////////////////////////////////////////////////////////////////////
1086 /// Return the node whose ID is 'ID'.
1087 ///////////////////////////////////////////////////////////////////////////////
1088 const SMDS_MeshNode * SMDS_Mesh::FindNode(int ID) const
1089 {
1090   return (const SMDS_MeshNode *)myNodeIDFactory->MeshElement(ID);
1091 }
1092
1093 ///////////////////////////////////////////////////////////////////////////////
1094 ///Create a triangle and add it to the current mesh. This methode do not bind a
1095 ///ID to the create triangle.
1096 ///////////////////////////////////////////////////////////////////////////////
1097 SMDS_MeshFace * SMDS_Mesh::createTriangle(const SMDS_MeshNode * node1,
1098                                           const SMDS_MeshNode * node2,
1099                                           const SMDS_MeshNode * node3)
1100 {
1101   if ( !node1 || !node2 || !node3) return 0;
1102   if ( myFaces.Extent() % CHECKMEMORY_INTERVAL == 0 ) CheckMemory();
1103   if(hasConstructionEdges())
1104   {
1105     SMDS_MeshEdge *edge1, *edge2, *edge3;
1106     edge1=FindEdgeOrCreate(node1,node2);
1107     edge2=FindEdgeOrCreate(node2,node3);
1108     edge3=FindEdgeOrCreate(node3,node1);
1109
1110     SMDS_MeshFace * face = new SMDS_FaceOfEdges(edge1,edge2,edge3);
1111     myFaces.Add(face);
1112     myInfo.myNbTriangles++;
1113     return face;
1114   }
1115   else
1116   {
1117     SMDS_MeshFace * face = new SMDS_FaceOfNodes(node1,node2,node3);
1118     myFaces.Add(face);
1119     myInfo.myNbTriangles++;
1120     return face;
1121   }
1122 }
1123
1124 ///////////////////////////////////////////////////////////////////////////////
1125 ///Create a quadrangle and add it to the current mesh. This methode do not bind
1126 ///a ID to the create triangle.
1127 ///////////////////////////////////////////////////////////////////////////////
1128 SMDS_MeshFace * SMDS_Mesh::createQuadrangle(const SMDS_MeshNode * node1,
1129                                             const SMDS_MeshNode * node2,
1130                                             const SMDS_MeshNode * node3,
1131                                             const SMDS_MeshNode * node4)
1132 {
1133   if ( !node1 || !node2 || !node3 || !node4 ) return 0;
1134   if ( myFaces.Extent() % CHECKMEMORY_INTERVAL == 0 ) CheckMemory();
1135   if(hasConstructionEdges())
1136   {
1137     SMDS_MeshEdge *edge1, *edge2, *edge3, *edge4;
1138     edge1=FindEdgeOrCreate(node1,node2);
1139     edge2=FindEdgeOrCreate(node2,node3);
1140     edge3=FindEdgeOrCreate(node3,node4);
1141     edge4=FindEdgeOrCreate(node4,node1);
1142
1143     SMDS_MeshFace * face = new SMDS_FaceOfEdges(edge1,edge2,edge3,edge4);
1144     myFaces.Add(face);
1145     myInfo.myNbQuadrangles++;
1146     return face;
1147   }
1148   else
1149   {
1150     SMDS_MeshFace * face = new SMDS_FaceOfNodes(node1,node2,node3,node4);
1151     myFaces.Add(face);
1152     myInfo.myNbQuadrangles++;
1153     return face;
1154   }
1155 }
1156
1157 ///////////////////////////////////////////////////////////////////////////////
1158 /// Remove a node and all the elements which own this node
1159 ///////////////////////////////////////////////////////////////////////////////
1160
1161 void SMDS_Mesh::RemoveNode(const SMDS_MeshNode * node)
1162 {
1163         RemoveElement(node, true);
1164 }
1165
1166 ///////////////////////////////////////////////////////////////////////////////
1167 /// Remove an edge and all the elements which own this edge
1168 ///////////////////////////////////////////////////////////////////////////////
1169
1170 void SMDS_Mesh::Remove0DElement(const SMDS_Mesh0DElement * elem0d)
1171 {
1172   RemoveElement(elem0d,true);
1173 }
1174
1175 ///////////////////////////////////////////////////////////////////////////////
1176 /// Remove an edge and all the elements which own this edge
1177 ///////////////////////////////////////////////////////////////////////////////
1178
1179 void SMDS_Mesh::RemoveEdge(const SMDS_MeshEdge * edge)
1180 {
1181         RemoveElement(edge,true);
1182 }
1183
1184 ///////////////////////////////////////////////////////////////////////////////
1185 /// Remove an face and all the elements which own this face
1186 ///////////////////////////////////////////////////////////////////////////////
1187
1188 void SMDS_Mesh::RemoveFace(const SMDS_MeshFace * face)
1189 {
1190         RemoveElement(face, true);
1191 }
1192
1193 ///////////////////////////////////////////////////////////////////////////////
1194 /// Remove a volume
1195 ///////////////////////////////////////////////////////////////////////////////
1196
1197 void SMDS_Mesh::RemoveVolume(const SMDS_MeshVolume * volume)
1198 {
1199         RemoveElement(volume, true);
1200 }
1201
1202 //=======================================================================
1203 //function : RemoveFromParent
1204 //purpose  :
1205 //=======================================================================
1206
1207 bool SMDS_Mesh::RemoveFromParent()
1208 {
1209         if (myParent==NULL) return false;
1210         else return (myParent->RemoveSubMesh(this));
1211 }
1212
1213 //=======================================================================
1214 //function : RemoveSubMesh
1215 //purpose  :
1216 //=======================================================================
1217
1218 bool SMDS_Mesh::RemoveSubMesh(const SMDS_Mesh * aMesh)
1219 {
1220         bool found = false;
1221
1222         list<SMDS_Mesh *>::iterator itmsh=myChildren.begin();
1223         for (; itmsh!=myChildren.end() && !found; itmsh++)
1224         {
1225                 SMDS_Mesh * submesh = *itmsh;
1226                 if (submesh == aMesh)
1227                 {
1228                         found = true;
1229                         myChildren.erase(itmsh);
1230                 }
1231         }
1232
1233         return found;
1234 }
1235
1236 //=======================================================================
1237 //function : ChangeElementNodes
1238 //purpose  :
1239 //=======================================================================
1240
1241 bool SMDS_Mesh::ChangeElementNodes(const SMDS_MeshElement * element,
1242                                    const SMDS_MeshNode    * nodes[],
1243                                    const int                nbnodes)
1244 {
1245   // keep current nodes of elem
1246   set<const SMDS_MeshElement*> oldNodes;
1247   SMDS_ElemIteratorPtr itn = element->nodesIterator();
1248   while(itn->more())
1249     oldNodes.insert(  itn->next() );
1250
1251   if ( !element->IsPoly() )
1252     myInfo.remove( element ); // element may change type
1253
1254   // change nodes
1255   bool Ok = false;
1256   SMDS_MeshElement* elem = const_cast<SMDS_MeshElement*>(element);
1257   switch ( elem->GetType() )
1258   {
1259   case SMDSAbs_0DElement: {
1260     if ( SMDS_Mesh0DElement* elem0d = dynamic_cast<SMDS_Mesh0DElement*>( elem ))
1261       Ok = elem0d->ChangeNode( nodes[0] );
1262     break;
1263   }
1264   case SMDSAbs_Edge: {
1265     if ( nbnodes == 2 ) {
1266       if ( SMDS_MeshEdge* edge = dynamic_cast<SMDS_MeshEdge*>( elem ))
1267         Ok = edge->ChangeNodes( nodes[0], nodes[1] );
1268     }
1269     else if ( nbnodes == 3 ) {
1270       if ( SMDS_QuadraticEdge* edge = dynamic_cast<SMDS_QuadraticEdge*>( elem ))
1271         Ok = edge->ChangeNodes( nodes[0], nodes[1], nodes[2] );
1272     }
1273     break;
1274   }
1275   case SMDSAbs_Face: {
1276     if ( SMDS_FaceOfNodes* face = dynamic_cast<SMDS_FaceOfNodes*>( elem ))
1277       Ok = face->ChangeNodes( nodes, nbnodes );
1278     else
1279       if ( SMDS_QuadraticFaceOfNodes* QF = dynamic_cast<SMDS_QuadraticFaceOfNodes*>( elem ))
1280         Ok = QF->ChangeNodes( nodes, nbnodes );
1281       else
1282         if (SMDS_PolygonalFaceOfNodes* face = dynamic_cast<SMDS_PolygonalFaceOfNodes*>(elem))
1283           Ok = face->ChangeNodes(nodes, nbnodes);
1284     break;
1285   }
1286   case SMDSAbs_Volume: {
1287     if ( SMDS_VolumeOfNodes* vol = dynamic_cast<SMDS_VolumeOfNodes*>( elem ))
1288       Ok = vol->ChangeNodes( nodes, nbnodes );
1289     else
1290       if ( SMDS_QuadraticVolumeOfNodes* QV = dynamic_cast<SMDS_QuadraticVolumeOfNodes*>( elem ))
1291         Ok = QV->ChangeNodes( nodes, nbnodes );
1292     break;
1293   }
1294   default:
1295     MESSAGE ( "WRONG ELEM TYPE");
1296   }
1297
1298   if ( Ok ) { // update InverseElements
1299
1300     set<const SMDS_MeshElement*>::iterator it;
1301
1302     // AddInverseElement to new nodes
1303     for ( int i = 0; i < nbnodes; i++ ) {
1304       it = oldNodes.find( nodes[i] );
1305       if ( it == oldNodes.end() )
1306         // new node
1307         const_cast<SMDS_MeshNode*>( nodes[i] )->AddInverseElement( elem );
1308       else
1309         // remove from oldNodes a node that remains in elem
1310         oldNodes.erase( it );
1311     }
1312     // RemoveInverseElement from the nodes removed from elem
1313     for ( it = oldNodes.begin(); it != oldNodes.end(); it++ )
1314     {
1315       SMDS_MeshNode * n = static_cast<SMDS_MeshNode *>
1316         (const_cast<SMDS_MeshElement *>( *it ));
1317       n->RemoveInverseElement( elem );
1318     }
1319   }
1320
1321   if ( !element->IsPoly() )
1322     myInfo.add( element ); // element may change type
1323
1324   return Ok;
1325 }
1326
1327 //=======================================================================
1328 //function : ChangePolyhedronNodes
1329 //purpose  : to change nodes of polyhedral volume
1330 //=======================================================================
1331 bool SMDS_Mesh::ChangePolyhedronNodes (const SMDS_MeshElement *            elem,
1332                                        const vector<const SMDS_MeshNode*>& nodes,
1333                                        const vector<int>                 & quantities)
1334 {
1335   if (elem->GetType() != SMDSAbs_Volume) {
1336     MESSAGE("WRONG ELEM TYPE");
1337     return false;
1338   }
1339
1340   const SMDS_PolyhedralVolumeOfNodes* vol = dynamic_cast<const SMDS_PolyhedralVolumeOfNodes*>(elem);
1341   if (!vol) {
1342     return false;
1343   }
1344
1345   // keep current nodes of elem
1346   set<const SMDS_MeshElement*> oldNodes;
1347   SMDS_ElemIteratorPtr itn = elem->nodesIterator();
1348   while (itn->more()) {
1349     oldNodes.insert(itn->next());
1350   }
1351
1352   // change nodes
1353   bool Ok = const_cast<SMDS_PolyhedralVolumeOfNodes*>(vol)->ChangeNodes(nodes, quantities);
1354   if (!Ok) {
1355     return false;
1356   }
1357
1358   // update InverseElements
1359
1360   // AddInverseElement to new nodes
1361   int nbnodes = nodes.size();
1362   set<const SMDS_MeshElement*>::iterator it;
1363   for (int i = 0; i < nbnodes; i++) {
1364     it = oldNodes.find(nodes[i]);
1365     if (it == oldNodes.end()) {
1366       // new node
1367       const_cast<SMDS_MeshNode*>(nodes[i])->AddInverseElement(elem);
1368     } else {
1369       // remove from oldNodes a node that remains in elem
1370       oldNodes.erase(it);
1371     }
1372   }
1373
1374   // RemoveInverseElement from the nodes removed from elem
1375   for (it = oldNodes.begin(); it != oldNodes.end(); it++) {
1376     SMDS_MeshNode * n = static_cast<SMDS_MeshNode *>
1377       (const_cast<SMDS_MeshElement *>( *it ));
1378     n->RemoveInverseElement(elem);
1379   }
1380
1381   return Ok;
1382 }
1383
1384
1385 //=======================================================================
1386 //function : Find0DElement
1387 //purpose  :
1388 //=======================================================================
1389 const SMDS_Mesh0DElement* SMDS_Mesh::Find0DElement(int idnode) const
1390 {
1391   const SMDS_MeshNode * node = FindNode(idnode);
1392   if(node == NULL) return NULL;
1393   return Find0DElement(node);
1394 }
1395
1396 const SMDS_Mesh0DElement* SMDS_Mesh::Find0DElement(const SMDS_MeshNode * node)
1397 {
1398   if (!node) return 0;
1399   const SMDS_Mesh0DElement* toReturn = NULL;
1400   SMDS_ElemIteratorPtr it1 = node->GetInverseElementIterator(SMDSAbs_0DElement);
1401   while (it1->more() && (toReturn == NULL)) {
1402     const SMDS_MeshElement* e = it1->next();
1403     if (e->NbNodes() == 1) {
1404       toReturn = static_cast<const SMDS_Mesh0DElement*>(e);
1405     }
1406   }
1407   return toReturn;
1408 }
1409
1410 //=======================================================================
1411 //function : Find0DElementOrCreate
1412 //purpose  :
1413 //=======================================================================
1414 SMDS_Mesh0DElement* SMDS_Mesh::Find0DElementOrCreate(const SMDS_MeshNode * node)
1415 {
1416   if (!node) return 0;
1417   SMDS_Mesh0DElement * toReturn = NULL;
1418   toReturn = const_cast<SMDS_Mesh0DElement*>(Find0DElement(node));
1419   if (toReturn == NULL) {
1420     if (my0DElements.Extent() % CHECKMEMORY_INTERVAL == 0) CheckMemory();
1421     toReturn = new SMDS_Mesh0DElement(node);
1422     my0DElements.Add(toReturn);
1423     myInfo.myNb0DElements++;
1424   }
1425   return toReturn;
1426 }
1427
1428
1429 //=======================================================================
1430 //function : FindEdge
1431 //purpose  :
1432 //=======================================================================
1433
1434 const SMDS_MeshEdge* SMDS_Mesh::FindEdge(int idnode1, int idnode2) const
1435 {
1436   const SMDS_MeshNode * node1=FindNode(idnode1);
1437   const SMDS_MeshNode * node2=FindNode(idnode2);
1438   if((node1==NULL)||(node2==NULL)) return NULL;
1439   return FindEdge(node1,node2);
1440 }
1441
1442 //#include "Profiler.h"
1443 const SMDS_MeshEdge* SMDS_Mesh::FindEdge(const SMDS_MeshNode * node1,
1444                                          const SMDS_MeshNode * node2)
1445 {
1446   if ( !node1 ) return 0;
1447   const SMDS_MeshEdge * toReturn=NULL;
1448   //PROFILER_Init();
1449   //PROFILER_Set();
1450   SMDS_ElemIteratorPtr it1=node1->GetInverseElementIterator(SMDSAbs_Edge);
1451   //PROFILER_Get(0);
1452   //PROFILER_Set();
1453   while(it1->more()) {
1454     const SMDS_MeshElement * e = it1->next();
1455     if ( e->NbNodes() == 2 && e->GetNodeIndex( node2 ) >= 0 ) {
1456       toReturn = static_cast<const SMDS_MeshEdge*>( e );
1457       break;
1458     }
1459   }
1460   //PROFILER_Get(1);
1461   return toReturn;
1462 }
1463
1464
1465 //=======================================================================
1466 //function : FindEdgeOrCreate
1467 //purpose  :
1468 //=======================================================================
1469
1470 SMDS_MeshEdge* SMDS_Mesh::FindEdgeOrCreate(const SMDS_MeshNode * node1,
1471                                            const SMDS_MeshNode * node2)
1472 {
1473   if ( !node1 || !node2) return 0;
1474   SMDS_MeshEdge * toReturn=NULL;
1475   toReturn=const_cast<SMDS_MeshEdge*>(FindEdge(node1,node2));
1476   if(toReturn==NULL) {
1477     if ( myEdges.Extent() % CHECKMEMORY_INTERVAL == 0 ) CheckMemory();
1478     toReturn=new SMDS_MeshEdge(node1,node2);
1479     myEdges.Add(toReturn);
1480     myInfo.myNbEdges++;
1481   }
1482   return toReturn;
1483 }
1484
1485
1486 //=======================================================================
1487 //function : FindEdge
1488 //purpose  :
1489 //=======================================================================
1490
1491 const SMDS_MeshEdge* SMDS_Mesh::FindEdge(int idnode1, int idnode2,
1492                                          int idnode3) const
1493 {
1494   const SMDS_MeshNode * node1=FindNode(idnode1);
1495   const SMDS_MeshNode * node2=FindNode(idnode2);
1496   const SMDS_MeshNode * node3=FindNode(idnode3);
1497   return FindEdge(node1,node2,node3);
1498 }
1499
1500 const SMDS_MeshEdge* SMDS_Mesh::FindEdge(const SMDS_MeshNode * node1,
1501                                          const SMDS_MeshNode * node2,
1502                                          const SMDS_MeshNode * node3)
1503 {
1504   if ( !node1 ) return 0;
1505   SMDS_ElemIteratorPtr it1 = node1->GetInverseElementIterator(SMDSAbs_Edge);
1506   while(it1->more()) {
1507     const SMDS_MeshElement * e = it1->next();
1508     if ( e->NbNodes() == 3 ) {
1509       SMDS_ElemIteratorPtr it2 = e->nodesIterator();
1510       while(it2->more()) {
1511         const SMDS_MeshElement* n = it2->next();
1512         if( n!=node1 &&
1513             n!=node2 &&
1514             n!=node3 )
1515         {
1516           e = 0;
1517           break;
1518         }
1519       }
1520       if ( e )
1521         return static_cast<const SMDS_MeshEdge *> (e);
1522     }
1523   }
1524   return 0;
1525 }
1526
1527
1528 //=======================================================================
1529 //function : FindFace
1530 //purpose  :
1531 //=======================================================================
1532
1533 const SMDS_MeshFace* SMDS_Mesh::FindFace(int idnode1, int idnode2,
1534         int idnode3) const
1535 {
1536   const SMDS_MeshNode * node1=FindNode(idnode1);
1537   const SMDS_MeshNode * node2=FindNode(idnode2);
1538   const SMDS_MeshNode * node3=FindNode(idnode3);
1539   return FindFace(node1, node2, node3);
1540 }
1541
1542 const SMDS_MeshFace* SMDS_Mesh::FindFace(const SMDS_MeshNode *node1,
1543                                          const SMDS_MeshNode *node2,
1544                                          const SMDS_MeshNode *node3)
1545 {
1546   if ( !node1 ) return 0;
1547   SMDS_ElemIteratorPtr it1 = node1->GetInverseElementIterator(SMDSAbs_Face);
1548   while(it1->more()) {
1549     const SMDS_MeshElement * e = it1->next();
1550     if ( e->NbNodes() == 3 ) {
1551       SMDS_ElemIteratorPtr it2 = e->nodesIterator();
1552       while(it2->more()) {
1553         const SMDS_MeshElement* n = it2->next();
1554         if( n!=node1 &&
1555             n!=node2 &&
1556             n!=node3 )
1557         {
1558           e = 0;
1559           break;
1560         }
1561       }
1562       if ( e )
1563         return static_cast<const SMDS_MeshFace *> (e);
1564     }
1565   }
1566   return 0;
1567 }
1568
1569 SMDS_MeshFace* SMDS_Mesh::FindFaceOrCreate(const SMDS_MeshNode *node1,
1570                                            const SMDS_MeshNode *node2,
1571                                            const SMDS_MeshNode *node3)
1572 {
1573   SMDS_MeshFace * toReturn=NULL;
1574   toReturn = const_cast<SMDS_MeshFace*>(FindFace(node1,node2,node3));
1575   if(toReturn==NULL) {
1576     toReturn = createTriangle(node1,node2,node3);
1577   }
1578   return toReturn;
1579 }
1580
1581
1582 //=======================================================================
1583 //function : FindFace
1584 //purpose  :
1585 //=======================================================================
1586
1587 const SMDS_MeshFace* SMDS_Mesh::FindFace(int idnode1, int idnode2,
1588                                          int idnode3, int idnode4) const
1589 {
1590   const SMDS_MeshNode * node1=FindNode(idnode1);
1591   const SMDS_MeshNode * node2=FindNode(idnode2);
1592   const SMDS_MeshNode * node3=FindNode(idnode3);
1593   const SMDS_MeshNode * node4=FindNode(idnode4);
1594   return FindFace(node1, node2, node3, node4);
1595 }
1596
1597 const SMDS_MeshFace* SMDS_Mesh::FindFace(const SMDS_MeshNode *node1,
1598                                          const SMDS_MeshNode *node2,
1599                                          const SMDS_MeshNode *node3,
1600                                          const SMDS_MeshNode *node4)
1601 {
1602   if ( !node1 ) return 0;
1603   SMDS_ElemIteratorPtr it1 = node1->GetInverseElementIterator(SMDSAbs_Face);
1604   while(it1->more()) {
1605     const SMDS_MeshElement * e = it1->next();
1606     if ( e->NbNodes() == 4 ) {
1607       SMDS_ElemIteratorPtr it2 = e->nodesIterator();
1608       while(it2->more()) {
1609         const SMDS_MeshElement* n = it2->next();
1610         if( n!=node1 &&
1611             n!=node2 &&
1612             n!=node3 &&
1613             n!=node4 )
1614         {
1615           e = 0;
1616           break;
1617         }
1618       }
1619       if ( e )
1620         return static_cast<const SMDS_MeshFace *> (e);
1621     }
1622   }
1623   return 0;
1624 }
1625
1626 SMDS_MeshFace* SMDS_Mesh::FindFaceOrCreate(const SMDS_MeshNode *node1,
1627                                            const SMDS_MeshNode *node2,
1628                                            const SMDS_MeshNode *node3,
1629                                            const SMDS_MeshNode *node4)
1630 {
1631   SMDS_MeshFace * toReturn=NULL;
1632   toReturn=const_cast<SMDS_MeshFace*>(FindFace(node1,node2,node3,node4));
1633   if(toReturn==NULL) {
1634     toReturn=createQuadrangle(node1,node2,node3,node4);
1635   }
1636   return toReturn;
1637 }
1638
1639
1640 //=======================================================================
1641 //function : FindFace
1642 //purpose  :quadratic triangle
1643 //=======================================================================
1644
1645 const SMDS_MeshFace* SMDS_Mesh::FindFace(int idnode1, int idnode2,
1646                                          int idnode3, int idnode4,
1647                                          int idnode5, int idnode6) const
1648 {
1649   const SMDS_MeshNode * node1 = FindNode(idnode1);
1650   const SMDS_MeshNode * node2 = FindNode(idnode2);
1651   const SMDS_MeshNode * node3 = FindNode(idnode3);
1652   const SMDS_MeshNode * node4 = FindNode(idnode4);
1653   const SMDS_MeshNode * node5 = FindNode(idnode5);
1654   const SMDS_MeshNode * node6 = FindNode(idnode6);
1655   return FindFace(node1, node2, node3, node4, node5, node6);
1656 }
1657
1658 const SMDS_MeshFace* SMDS_Mesh::FindFace(const SMDS_MeshNode *node1,
1659                                          const SMDS_MeshNode *node2,
1660                                          const SMDS_MeshNode *node3,
1661                                          const SMDS_MeshNode *node4,
1662                                          const SMDS_MeshNode *node5,
1663                                          const SMDS_MeshNode *node6)
1664 {
1665   if ( !node1 ) return 0;
1666   SMDS_ElemIteratorPtr it1 = node1->GetInverseElementIterator(SMDSAbs_Face);
1667   while(it1->more()) {
1668     const SMDS_MeshElement * e = it1->next();
1669     if ( e->NbNodes() == 6 ) {
1670       SMDS_ElemIteratorPtr it2 = e->nodesIterator();
1671       while(it2->more()) {
1672         const SMDS_MeshElement* n = it2->next();
1673         if( n!=node1 &&
1674             n!=node2 &&
1675             n!=node3 &&
1676             n!=node4 &&
1677             n!=node5 &&
1678             n!=node6 )
1679         {
1680           e = 0;
1681           break;
1682         }
1683       }
1684       if ( e )
1685         return static_cast<const SMDS_MeshFace *> (e);
1686     }
1687   }
1688   return 0;
1689 }
1690
1691
1692 //=======================================================================
1693 //function : FindFace
1694 //purpose  : quadratic quadrangle
1695 //=======================================================================
1696
1697 const SMDS_MeshFace* SMDS_Mesh::FindFace(int idnode1, int idnode2,
1698                                          int idnode3, int idnode4,
1699                                          int idnode5, int idnode6,
1700                                          int idnode7, int idnode8) const
1701 {
1702   const SMDS_MeshNode * node1 = FindNode(idnode1);
1703   const SMDS_MeshNode * node2 = FindNode(idnode2);
1704   const SMDS_MeshNode * node3 = FindNode(idnode3);
1705   const SMDS_MeshNode * node4 = FindNode(idnode4);
1706   const SMDS_MeshNode * node5 = FindNode(idnode5);
1707   const SMDS_MeshNode * node6 = FindNode(idnode6);
1708   const SMDS_MeshNode * node7 = FindNode(idnode7);
1709   const SMDS_MeshNode * node8 = FindNode(idnode8);
1710   return FindFace(node1, node2, node3, node4, node5, node6, node7, node8);
1711 }
1712
1713 const SMDS_MeshFace* SMDS_Mesh::FindFace(const SMDS_MeshNode *node1,
1714                                          const SMDS_MeshNode *node2,
1715                                          const SMDS_MeshNode *node3,
1716                                          const SMDS_MeshNode *node4,
1717                                          const SMDS_MeshNode *node5,
1718                                          const SMDS_MeshNode *node6,
1719                                          const SMDS_MeshNode *node7,
1720                                          const SMDS_MeshNode *node8)
1721 {
1722   if ( !node1 ) return 0;
1723   SMDS_ElemIteratorPtr it1 = node1->GetInverseElementIterator(SMDSAbs_Face);
1724   while(it1->more()) {
1725     const SMDS_MeshElement * e = it1->next();
1726     if ( e->NbNodes() == 8 ) {
1727       SMDS_ElemIteratorPtr it2 = e->nodesIterator();
1728       while(it2->more()) {
1729         const SMDS_MeshElement* n = it2->next();
1730         if( n!=node1 &&
1731             n!=node2 &&
1732             n!=node3 &&
1733             n!=node4 &&
1734             n!=node5 &&
1735             n!=node6 &&
1736             n!=node7 &&
1737             n!=node8 )
1738         {
1739           e = 0;
1740           break;
1741         }
1742       }
1743       if ( e )
1744         return static_cast<const SMDS_MeshFace *> (e);
1745     }
1746   }
1747   return 0;
1748 }
1749
1750
1751 //=======================================================================
1752 //function : FindElement
1753 //purpose  :
1754 //=======================================================================
1755
1756 const SMDS_MeshElement* SMDS_Mesh::FindElement(int IDelem) const
1757 {
1758   return myElementIDFactory->MeshElement(IDelem);
1759 }
1760
1761 //=======================================================================
1762 //function : FindFace
1763 //purpose  : find polygon
1764 //=======================================================================
1765
1766 const SMDS_MeshFace* SMDS_Mesh::FindFace (std::vector<int> nodes_ids) const
1767 {
1768   int nbnodes = nodes_ids.size();
1769   std::vector<const SMDS_MeshNode *> poly_nodes (nbnodes);
1770   for (int inode = 0; inode < nbnodes; inode++) {
1771     const SMDS_MeshNode * node = FindNode(nodes_ids[inode]);
1772     if (node == NULL) return NULL;
1773     poly_nodes[inode] = node;
1774   }
1775   return FindFace(poly_nodes);
1776 }
1777
1778 const SMDS_MeshFace* SMDS_Mesh::FindFace (std::vector<const SMDS_MeshNode *> nodes)
1779 {
1780   if ( nodes.size() > 2 && nodes[0] ) {
1781     SMDS_ElemIteratorPtr itF = nodes[0]->GetInverseElementIterator(SMDSAbs_Face);
1782     while (itF->more()) {
1783       const SMDS_MeshElement* f = itF->next();
1784       if ( f->NbNodes() == nodes.size() ) {
1785         for ( int i = 1; i < nodes.size(); ++ i )
1786           if ( f->GetNodeIndex( nodes[ i ]) < 0 )
1787             return NULL;
1788         return static_cast<const SMDS_MeshFace *> (f);
1789       }
1790     }
1791   }
1792   return NULL;
1793 }
1794
1795 //=======================================================================
1796 //function : DumpNodes
1797 //purpose  :
1798 //=======================================================================
1799
1800 void SMDS_Mesh::DumpNodes() const
1801 {
1802         MESSAGE("dump nodes of mesh : ");
1803         SMDS_NodeIteratorPtr itnode=nodesIterator();
1804         while(itnode->more()) MESSAGE(itnode->next());
1805 }
1806
1807 //=======================================================================
1808 //function : Dump0DElements
1809 //purpose  :
1810 //=======================================================================
1811 void SMDS_Mesh::Dump0DElements() const
1812 {
1813   MESSAGE("dump 0D elements of mesh : ");
1814   SMDS_0DElementIteratorPtr it0d = elements0dIterator();
1815   while(it0d->more()) MESSAGE(it0d->next());
1816 }
1817
1818 //=======================================================================
1819 //function : DumpEdges
1820 //purpose  :
1821 //=======================================================================
1822
1823 void SMDS_Mesh::DumpEdges() const
1824 {
1825         MESSAGE("dump edges of mesh : ");
1826         SMDS_EdgeIteratorPtr itedge=edgesIterator();
1827         while(itedge->more()) MESSAGE(itedge->next());
1828 }
1829
1830 //=======================================================================
1831 //function : DumpFaces
1832 //purpose  :
1833 //=======================================================================
1834
1835 void SMDS_Mesh::DumpFaces() const
1836 {
1837         MESSAGE("dump faces of mesh : ");
1838         SMDS_FaceIteratorPtr itface=facesIterator();
1839         while(itface->more()) MESSAGE(itface->next());
1840 }
1841
1842 //=======================================================================
1843 //function : DumpVolumes
1844 //purpose  :
1845 //=======================================================================
1846
1847 void SMDS_Mesh::DumpVolumes() const
1848 {
1849         MESSAGE("dump volumes of mesh : ");
1850         SMDS_VolumeIteratorPtr itvol=volumesIterator();
1851         while(itvol->more()) MESSAGE(itvol->next());
1852 }
1853
1854 //=======================================================================
1855 //function : DebugStats
1856 //purpose  :
1857 //=======================================================================
1858
1859 void SMDS_Mesh::DebugStats() const
1860 {
1861   MESSAGE("Debug stats of mesh : ");
1862
1863   MESSAGE("===== NODES ====="<<NbNodes());
1864   MESSAGE("===== 0DELEMS ====="<<Nb0DElements());
1865   MESSAGE("===== EDGES ====="<<NbEdges());
1866   MESSAGE("===== FACES ====="<<NbFaces());
1867   MESSAGE("===== VOLUMES ====="<<NbVolumes());
1868
1869   MESSAGE("End Debug stats of mesh ");
1870
1871   //#ifdef DEB
1872
1873   SMDS_NodeIteratorPtr itnode=nodesIterator();
1874   int sizeofnodes = 0;
1875   int sizeoffaces = 0;
1876
1877   while(itnode->more())
1878   {
1879     const SMDS_MeshNode *node = itnode->next();
1880
1881     sizeofnodes += sizeof(*node);
1882
1883     SMDS_ElemIteratorPtr it = node->GetInverseElementIterator();
1884     while(it->more())
1885     {
1886       const SMDS_MeshElement *me = it->next();
1887       sizeofnodes += sizeof(me);
1888     }
1889   }
1890
1891   SMDS_FaceIteratorPtr itface=facesIterator();
1892   while(itface->more())
1893   {
1894     const SMDS_MeshElement *face = itface->next();
1895     sizeoffaces += sizeof(*face);
1896   }
1897
1898   MESSAGE("total size of node elements = " << sizeofnodes);;
1899   MESSAGE("total size of face elements = " << sizeoffaces);;
1900
1901   //#endif
1902 }
1903
1904 ///////////////////////////////////////////////////////////////////////////////
1905 /// Return the number of nodes
1906 ///////////////////////////////////////////////////////////////////////////////
1907 int SMDS_Mesh::NbNodes() const
1908 {
1909         return myNodes.Size();
1910 }
1911
1912 ///////////////////////////////////////////////////////////////////////////////
1913 /// Return the number of 0D elements
1914 ///////////////////////////////////////////////////////////////////////////////
1915 int SMDS_Mesh::Nb0DElements() const
1916 {
1917   return my0DElements.Size();
1918 }
1919
1920 ///////////////////////////////////////////////////////////////////////////////
1921 /// Return the number of edges (including construction edges)
1922 ///////////////////////////////////////////////////////////////////////////////
1923 int SMDS_Mesh::NbEdges() const
1924 {
1925         return myEdges.Size();
1926 }
1927
1928 ///////////////////////////////////////////////////////////////////////////////
1929 /// Return the number of faces (including construction faces)
1930 ///////////////////////////////////////////////////////////////////////////////
1931 int SMDS_Mesh::NbFaces() const
1932 {
1933         return myFaces.Size();
1934 }
1935
1936 ///////////////////////////////////////////////////////////////////////////////
1937 /// Return the number of volumes
1938 ///////////////////////////////////////////////////////////////////////////////
1939 int SMDS_Mesh::NbVolumes() const
1940 {
1941         return myVolumes.Size();
1942 }
1943
1944 ///////////////////////////////////////////////////////////////////////////////
1945 /// Return the number of child mesh of this mesh.
1946 /// Note that the tree structure of SMDS_Mesh seems to be unused in this version
1947 /// (2003-09-08) of SMESH
1948 ///////////////////////////////////////////////////////////////////////////////
1949 int SMDS_Mesh::NbSubMesh() const
1950 {
1951         return myChildren.size();
1952 }
1953
1954 ///////////////////////////////////////////////////////////////////////////////
1955 /// Destroy the mesh and all its elements
1956 /// All pointer on elements owned by this mesh become illegals.
1957 ///////////////////////////////////////////////////////////////////////////////
1958 SMDS_Mesh::~SMDS_Mesh()
1959 {
1960   list<SMDS_Mesh*>::iterator itc=myChildren.begin();
1961   while(itc!=myChildren.end())
1962   {
1963     delete *itc;
1964     itc++;
1965   }
1966
1967   if(myParent==NULL)
1968   {
1969     delete myNodeIDFactory;
1970     delete myElementIDFactory;
1971   }
1972   else
1973   {
1974     SMDS_ElemIteratorPtr eIt = elementsIterator();
1975     while ( eIt->more() )
1976       myElementIDFactory->ReleaseID(eIt->next()->GetID());
1977     SMDS_NodeIteratorPtr itn = nodesIterator();
1978     while (itn->more())
1979       myNodeIDFactory->ReleaseID(itn->next()->GetID());
1980   }
1981
1982   SetOfNodes::Iterator itn(myNodes);
1983   for (; itn.More(); itn.Next())
1984     delete itn.Value();
1985
1986   SetOf0DElements::Iterator it0d (my0DElements);
1987   for (; it0d.More(); it0d.Next())
1988   {
1989     SMDS_MeshElement* elem = it0d.Value();
1990     delete elem;
1991   }
1992
1993   SetOfEdges::Iterator ite(myEdges);
1994   for (; ite.More(); ite.Next())
1995   {
1996     SMDS_MeshElement* elem = ite.Value();
1997     delete elem;
1998   }
1999
2000   SetOfFaces::Iterator itf(myFaces);
2001   for (; itf.More(); itf.Next())
2002   {
2003     SMDS_MeshElement* elem = itf.Value();
2004     delete elem;
2005   }
2006
2007   SetOfVolumes::Iterator itv(myVolumes);
2008   for (; itv.More(); itv.Next())
2009   {
2010     SMDS_MeshElement* elem = itv.Value();
2011     delete elem;
2012   }
2013 }
2014
2015 //================================================================================
2016 /*!
2017  * \brief Clear all data
2018  */
2019 //================================================================================
2020
2021 void SMDS_Mesh::Clear()
2022 {
2023   if (myParent!=NULL) {
2024     SMDS_ElemIteratorPtr eIt = elementsIterator();
2025     while ( eIt->more() )
2026       myElementIDFactory->ReleaseID(eIt->next()->GetID());
2027     SMDS_NodeIteratorPtr itn = nodesIterator();
2028     while (itn->more())
2029       myNodeIDFactory->ReleaseID(itn->next()->GetID());
2030   }
2031   else {
2032     myNodeIDFactory->Clear();
2033     myElementIDFactory->Clear();
2034   }
2035
2036   SMDS_VolumeIteratorPtr itv = volumesIterator();
2037   while (itv->more())
2038     delete itv->next();
2039   myVolumes.Clear();
2040
2041   SMDS_FaceIteratorPtr itf = facesIterator();
2042   while (itf->more())
2043     delete itf->next();
2044   myFaces.Clear();
2045
2046   SMDS_EdgeIteratorPtr ite = edgesIterator();
2047   while (ite->more())
2048     delete ite->next();
2049   myEdges.Clear();
2050
2051   SMDS_0DElementIteratorPtr it0d = elements0dIterator();
2052   while (it0d->more())
2053     delete it0d->next();
2054   my0DElements.Clear();
2055
2056   SMDS_NodeIteratorPtr itn = nodesIterator();
2057   while (itn->more())
2058     delete itn->next();
2059   myNodes.Clear();
2060
2061   list<SMDS_Mesh*>::iterator itc=myChildren.begin();
2062   while(itc!=myChildren.end())
2063     (*itc)->Clear();
2064
2065   myInfo.Clear();
2066 }
2067
2068 ///////////////////////////////////////////////////////////////////////////////
2069 /// Return true if this mesh create faces with edges.
2070 /// A false returned value mean that faces are created with nodes. A concequence
2071 /// is, iteration on edges (SMDS_Element::edgesIterator) will be unavailable.
2072 ///////////////////////////////////////////////////////////////////////////////
2073 bool SMDS_Mesh::hasConstructionEdges()
2074 {
2075         return myHasConstructionEdges;
2076 }
2077
2078 ///////////////////////////////////////////////////////////////////////////////
2079 /// Return true if this mesh create volumes with faces
2080 /// A false returned value mean that volumes are created with nodes or edges.
2081 /// (see hasConstructionEdges)
2082 /// A concequence is, iteration on faces (SMDS_Element::facesIterator) will be
2083 /// unavailable.
2084 ///////////////////////////////////////////////////////////////////////////////
2085 bool SMDS_Mesh::hasConstructionFaces()
2086 {
2087         return myHasConstructionFaces;
2088 }
2089
2090 ///////////////////////////////////////////////////////////////////////////////
2091 /// Return true if nodes are linked to the finit elements, they are belonging to.
2092 /// Currently, It always return true.
2093 ///////////////////////////////////////////////////////////////////////////////
2094 bool SMDS_Mesh::hasInverseElements()
2095 {
2096         return myHasInverseElements;
2097 }
2098
2099 ///////////////////////////////////////////////////////////////////////////////
2100 /// Make this mesh creating construction edges (see hasConstructionEdges)
2101 /// @param b true to have construction edges, else false.
2102 ///////////////////////////////////////////////////////////////////////////////
2103 void SMDS_Mesh::setConstructionEdges(bool b)
2104 {
2105         myHasConstructionEdges=b;
2106 }
2107
2108 ///////////////////////////////////////////////////////////////////////////////
2109 /// Make this mesh creating construction faces (see hasConstructionFaces)
2110 /// @param b true to have construction faces, else false.
2111 ///////////////////////////////////////////////////////////////////////////////
2112 void SMDS_Mesh::setConstructionFaces(bool b)
2113 {
2114          myHasConstructionFaces=b;
2115 }
2116
2117 ///////////////////////////////////////////////////////////////////////////////
2118 /// Make this mesh creating link from nodes to elements (see hasInverseElements)
2119 /// @param b true to link nodes to elements, else false.
2120 ///////////////////////////////////////////////////////////////////////////////
2121 void SMDS_Mesh::setInverseElements(bool b)
2122 {
2123   if(!b) MESSAGE("Error : inverseElement=false not implemented");
2124   myHasInverseElements=b;
2125 }
2126
2127 namespace {
2128
2129   ///////////////////////////////////////////////////////////////////////////////
2130   ///Iterator on NCollection_Map
2131   ///////////////////////////////////////////////////////////////////////////////
2132   template <class MAP, typename ELEM=const SMDS_MeshElement*, class FATHER=SMDS_ElemIterator>
2133   struct MYNCollection_Map_Iterator: public FATHER
2134   {
2135     typename MAP::Iterator myIterator;
2136
2137     MYNCollection_Map_Iterator(const MAP& map):myIterator(map){}
2138
2139     bool more()
2140     {
2141       while(myIterator.More())
2142       {
2143         if(myIterator.Value()->GetID()!=-1)
2144           return true;
2145         myIterator.Next();
2146       }
2147       return false;
2148     }
2149
2150     ELEM next()
2151     {
2152       ELEM current = (ELEM) myIterator.Value();
2153       myIterator.Next();
2154       return current;
2155     }
2156   };
2157   //================================================================================
2158   /*!
2159    * \brief Iterator on elements in id increasing order
2160    */
2161   //================================================================================
2162
2163   template <typename ELEM=const SMDS_MeshElement*>
2164   class IdSortedIterator : public SMDS_Iterator<ELEM>
2165   {
2166     const SMDS_MeshElementIDFactory& myIDFact;
2167     int                              myID, myMaxID;
2168     ELEM                             myElem;
2169
2170   public:
2171     IdSortedIterator(const SMDS_MeshElementIDFactory& fact)
2172       :myIDFact( fact ), myID(1), myMaxID( myIDFact.GetMaxID() ), myElem(0)
2173     {
2174       next();
2175     }
2176     bool more()
2177     {
2178       return myElem;
2179     }
2180     ELEM next()
2181     {
2182       ELEM current = myElem;
2183       for ( myElem = 0; myID <= myMaxID && !myElem; ++myID )
2184         myElem = (ELEM) myIDFact.MeshElement( myID );
2185       return current;
2186     }
2187   };
2188 }
2189
2190 ///////////////////////////////////////////////////////////////////////////////
2191 /// Return an iterator on nodes of the current mesh factory
2192 ///////////////////////////////////////////////////////////////////////////////
2193
2194 SMDS_NodeIteratorPtr SMDS_Mesh::nodesIterator(bool idInceasingOrder) const
2195 {
2196   typedef MYNCollection_Map_Iterator
2197     < SetOfNodes, const SMDS_MeshNode*, SMDS_NodeIterator > TIterator;
2198   typedef IdSortedIterator< const SMDS_MeshNode* >          TSortedIterator;
2199   return ( idInceasingOrder ?
2200            SMDS_NodeIteratorPtr( new TSortedIterator( *myNodeIDFactory )) :
2201            SMDS_NodeIteratorPtr( new TIterator(myNodes)));
2202 }
2203
2204 ///////////////////////////////////////////////////////////////////////////////
2205 ///Return an iterator on 0D elements of the current mesh.
2206 ///////////////////////////////////////////////////////////////////////////////
2207
2208 SMDS_0DElementIteratorPtr SMDS_Mesh::elements0dIterator() const
2209 {
2210   typedef MYNCollection_Map_Iterator
2211     < SetOf0DElements, const SMDS_Mesh0DElement*, SMDS_0DElementIterator > TIterator;
2212   return SMDS_0DElementIteratorPtr(new TIterator(my0DElements));
2213 }
2214
2215 ///////////////////////////////////////////////////////////////////////////////
2216 ///Return an iterator on edges of the current mesh.
2217 ///////////////////////////////////////////////////////////////////////////////
2218
2219 SMDS_EdgeIteratorPtr SMDS_Mesh::edgesIterator() const
2220 {
2221   typedef MYNCollection_Map_Iterator
2222     < SetOfEdges, const SMDS_MeshEdge*, SMDS_EdgeIterator > TIterator;
2223   return SMDS_EdgeIteratorPtr(new TIterator(myEdges));
2224 }
2225
2226 ///////////////////////////////////////////////////////////////////////////////
2227 ///Return an iterator on faces of the current mesh.
2228 ///////////////////////////////////////////////////////////////////////////////
2229
2230 SMDS_FaceIteratorPtr SMDS_Mesh::facesIterator() const
2231 {
2232   typedef MYNCollection_Map_Iterator
2233     < SetOfFaces, const SMDS_MeshFace*, SMDS_FaceIterator > TIterator;
2234   return SMDS_FaceIteratorPtr(new TIterator(myFaces));
2235 }
2236
2237 ///////////////////////////////////////////////////////////////////////////////
2238 ///Return an iterator on volumes of the current mesh.
2239 ///////////////////////////////////////////////////////////////////////////////
2240
2241 SMDS_VolumeIteratorPtr SMDS_Mesh::volumesIterator() const
2242 {
2243   typedef MYNCollection_Map_Iterator
2244     < SetOfVolumes, const SMDS_MeshVolume*, SMDS_VolumeIterator > TIterator;
2245   return SMDS_VolumeIteratorPtr(new TIterator(myVolumes));
2246 }
2247
2248 ///////////////////////////////////////////////////////////////////////////////
2249 /// Return an iterator on elements of the current mesh factory
2250 ///////////////////////////////////////////////////////////////////////////////
2251 SMDS_ElemIteratorPtr SMDS_Mesh::elementsIterator(SMDSAbs_ElementType type) const
2252 {
2253   switch (type) {
2254   case SMDSAbs_All:
2255     break;
2256   case SMDSAbs_Volume:
2257     return SMDS_ElemIteratorPtr (new MYNCollection_Map_Iterator< SetOfVolumes >(myVolumes));
2258   case SMDSAbs_Face:
2259     return SMDS_ElemIteratorPtr (new MYNCollection_Map_Iterator< SetOfFaces >(myFaces));
2260   case SMDSAbs_Edge:
2261     return SMDS_ElemIteratorPtr (new MYNCollection_Map_Iterator< SetOfEdges >(myEdges));
2262   case SMDSAbs_0DElement:
2263     return SMDS_ElemIteratorPtr (new MYNCollection_Map_Iterator< SetOf0DElements >(my0DElements));
2264   case SMDSAbs_Node:
2265     return myNodeIDFactory->elementsIterator();
2266   default:;
2267   }
2268   return myElementIDFactory->elementsIterator();
2269 }
2270
2271 ///////////////////////////////////////////////////////////////////////////////
2272 /// Do intersection of sets (more than 2)
2273 ///////////////////////////////////////////////////////////////////////////////
2274 static set<const SMDS_MeshElement*> * intersectionOfSets(
2275         set<const SMDS_MeshElement*> vs[], int numberOfSets)
2276 {
2277         set<const SMDS_MeshElement*>* rsetA=new set<const SMDS_MeshElement*>(vs[0]);
2278         set<const SMDS_MeshElement*>* rsetB;
2279
2280         for(int i=0; i<numberOfSets-1; i++)
2281         {
2282                 rsetB=new set<const SMDS_MeshElement*>();
2283                 set_intersection(
2284                         rsetA->begin(), rsetA->end(),
2285                         vs[i+1].begin(), vs[i+1].end(),
2286                         inserter(*rsetB, rsetB->begin()));
2287                 delete rsetA;
2288                 rsetA=rsetB;
2289         }
2290         return rsetA;
2291 }
2292
2293 ///////////////////////////////////////////////////////////////////////////////
2294 /// Return the list of finit elements owning the given element
2295 ///////////////////////////////////////////////////////////////////////////////
2296 static set<const SMDS_MeshElement*> * getFinitElements(const SMDS_MeshElement * element)
2297 {
2298         int numberOfSets=element->NbNodes();
2299         set<const SMDS_MeshElement*> *initSet = new set<const SMDS_MeshElement*>[numberOfSets];
2300
2301         SMDS_ElemIteratorPtr itNodes=element->nodesIterator();
2302
2303         int i=0;
2304         while(itNodes->more())
2305         {
2306                 const SMDS_MeshNode * n=static_cast<const SMDS_MeshNode*>(itNodes->next());
2307                 SMDS_ElemIteratorPtr itFe = n->GetInverseElementIterator();
2308
2309                 //initSet[i]=set<const SMDS_MeshElement*>();
2310                 while(itFe->more())
2311                   initSet[i].insert(itFe->next());
2312
2313                 i++;
2314         }
2315         set<const SMDS_MeshElement*> *retSet=intersectionOfSets(initSet, numberOfSets);
2316         delete [] initSet;
2317         return retSet;
2318 }
2319
2320 ///////////////////////////////////////////////////////////////////////////////
2321 /// Return the list of nodes used only by the given elements
2322 ///////////////////////////////////////////////////////////////////////////////
2323 static set<const SMDS_MeshElement*> * getExclusiveNodes(
2324         set<const SMDS_MeshElement*>& elements)
2325 {
2326         set<const SMDS_MeshElement*> * toReturn=new set<const SMDS_MeshElement*>();
2327         set<const SMDS_MeshElement*>::iterator itElements=elements.begin();
2328
2329         while(itElements!=elements.end())
2330         {
2331                 SMDS_ElemIteratorPtr itNodes = (*itElements)->nodesIterator();
2332                 itElements++;
2333
2334                 while(itNodes->more())
2335                 {
2336                         const SMDS_MeshNode * n=static_cast<const SMDS_MeshNode*>(itNodes->next());
2337                         SMDS_ElemIteratorPtr itFe = n->GetInverseElementIterator();
2338                         set<const SMDS_MeshElement*> s;
2339                         while(itFe->more())
2340                           s.insert(itFe->next());
2341                         if(s==elements) toReturn->insert(n);
2342                 }
2343         }
2344         return toReturn;
2345 }
2346
2347 ///////////////////////////////////////////////////////////////////////////////
2348 ///Find the children of an element that are made of given nodes
2349 ///@param setOfChildren The set in which matching children will be inserted
2350 ///@param element The element were to search matching children
2351 ///@param nodes The nodes that the children must have to be selected
2352 ///////////////////////////////////////////////////////////////////////////////
2353 void SMDS_Mesh::addChildrenWithNodes(set<const SMDS_MeshElement*>& setOfChildren,
2354                                      const SMDS_MeshElement *      element,
2355                                      set<const SMDS_MeshElement*>& nodes)
2356 {
2357   switch(element->GetType())
2358     {
2359     case SMDSAbs_Node:
2360       MESSAGE("Internal Error: This should not happend");
2361       break;
2362     case SMDSAbs_0DElement:
2363       {
2364       }
2365       break;
2366     case SMDSAbs_Edge:
2367         {
2368                 SMDS_ElemIteratorPtr itn=element->nodesIterator();
2369                 while(itn->more())
2370                 {
2371                         const SMDS_MeshElement * e=itn->next();
2372                         if(nodes.find(e)!=nodes.end())
2373                         {
2374                           setOfChildren.insert(element);
2375                           break;
2376                         }
2377                 }
2378         } break;
2379     case SMDSAbs_Face:
2380         {
2381                 SMDS_ElemIteratorPtr itn=element->nodesIterator();
2382                 while(itn->more())
2383                 {
2384                         const SMDS_MeshElement * e=itn->next();
2385                         if(nodes.find(e)!=nodes.end())
2386                         {
2387                           setOfChildren.insert(element);
2388                           break;
2389                         }
2390                 }
2391                 if(hasConstructionEdges())
2392                 {
2393                         SMDS_ElemIteratorPtr ite=element->edgesIterator();
2394                         while(ite->more())
2395                                 addChildrenWithNodes(setOfChildren, ite->next(), nodes);
2396                 }
2397         } break;
2398     case SMDSAbs_Volume:
2399         {
2400                 if(hasConstructionFaces())
2401                 {
2402                         SMDS_ElemIteratorPtr ite=element->facesIterator();
2403                         while(ite->more())
2404                                 addChildrenWithNodes(setOfChildren, ite->next(), nodes);
2405                 }
2406                 else if(hasConstructionEdges())
2407                 {
2408                         SMDS_ElemIteratorPtr ite=element->edgesIterator();
2409                         while(ite->more())
2410                                 addChildrenWithNodes(setOfChildren, ite->next(), nodes);
2411                 }
2412         }
2413     }
2414 }
2415
2416 ///////////////////////////////////////////////////////////////////////////////
2417 ///@param elem The element to delete
2418 ///@param removenodes if true remaining nodes will be removed
2419 ///////////////////////////////////////////////////////////////////////////////
2420 void SMDS_Mesh::RemoveElement(const SMDS_MeshElement * elem,
2421                               const bool removenodes)
2422 {
2423   list<const SMDS_MeshElement *> removedElems;
2424   list<const SMDS_MeshElement *> removedNodes;
2425   RemoveElement( elem, removedElems, removedNodes, removenodes );
2426 }
2427
2428 ///////////////////////////////////////////////////////////////////////////////
2429 ///@param elem The element to delete
2430 ///@param removedElems contains all removed elements
2431 ///@param removedNodes contains all removed nodes
2432 ///@param removenodes if true remaining nodes will be removed
2433 ///////////////////////////////////////////////////////////////////////////////
2434 void SMDS_Mesh::RemoveElement(const SMDS_MeshElement *        elem,
2435                               list<const SMDS_MeshElement *>& removedElems,
2436                               list<const SMDS_MeshElement *>& removedNodes,
2437                               bool                            removenodes)
2438 {
2439   // get finite elements built on elem
2440   set<const SMDS_MeshElement*> * s1;
2441   if (elem->GetType() == SMDSAbs_0DElement ||
2442       elem->GetType() == SMDSAbs_Edge && !hasConstructionEdges() ||
2443       elem->GetType() == SMDSAbs_Face && !hasConstructionFaces() ||
2444       elem->GetType() == SMDSAbs_Volume)
2445   {
2446     s1 = new set<const SMDS_MeshElement*>();
2447     s1->insert(elem);
2448   }
2449   else
2450     s1 = getFinitElements(elem);
2451
2452   // get exclusive nodes (which would become free afterwards)
2453   set<const SMDS_MeshElement*> * s2;
2454   if (elem->GetType() == SMDSAbs_Node) // a node is removed
2455   {
2456     // do not remove nodes except elem
2457     s2 = new set<const SMDS_MeshElement*>();
2458     s2->insert(elem);
2459     removenodes = true;
2460   }
2461   else
2462     s2 = getExclusiveNodes(*s1);
2463
2464   // form the set of finite and construction elements to remove
2465   set<const SMDS_MeshElement*> s3;
2466   set<const SMDS_MeshElement*>::iterator it=s1->begin();
2467   while(it!=s1->end())
2468   {
2469     addChildrenWithNodes(s3, *it ,*s2);
2470     s3.insert(*it);
2471     it++;
2472   }
2473   if(elem->GetType()!=SMDSAbs_Node) s3.insert(elem);
2474
2475   // remove finite and construction elements
2476   it=s3.begin();
2477   while(it!=s3.end())
2478   {
2479     // Remove element from <InverseElements> of its nodes
2480     SMDS_ElemIteratorPtr itn=(*it)->nodesIterator();
2481     while(itn->more())
2482     {
2483       SMDS_MeshNode * n = static_cast<SMDS_MeshNode *>
2484         (const_cast<SMDS_MeshElement *>(itn->next()));
2485       n->RemoveInverseElement( (*it) );
2486     }
2487
2488     switch((*it)->GetType())
2489     {
2490     case SMDSAbs_Node:
2491       MESSAGE("Internal Error: This should not happen");
2492       break;
2493     case SMDSAbs_0DElement:
2494       my0DElements.Remove(static_cast<SMDS_Mesh0DElement*>
2495                           (const_cast<SMDS_MeshElement*>(*it)));
2496       //myInfo.Remove0DElement(*it);
2497       myInfo.remove(*it);
2498       break;
2499     case SMDSAbs_Edge:
2500       myEdges.Remove(static_cast<SMDS_MeshEdge*>
2501                     (const_cast<SMDS_MeshElement*>(*it)));
2502       myInfo.RemoveEdge(*it);
2503       break;
2504     case SMDSAbs_Face:
2505       myFaces.Remove(static_cast<SMDS_MeshFace*>
2506                     (const_cast<SMDS_MeshElement*>(*it)));
2507       myInfo.RemoveFace(*it);
2508       break;
2509     case SMDSAbs_Volume:
2510       myVolumes.Remove(static_cast<SMDS_MeshVolume*>
2511                       (const_cast<SMDS_MeshElement*>(*it)));
2512       myInfo.RemoveVolume(*it);
2513       break;
2514     }
2515     //MESSAGE( "SMDS: RM elem " << (*it)->GetID() );
2516     removedElems.push_back( (*it) );
2517     myElementIDFactory->ReleaseID((*it)->GetID());
2518     delete (*it);
2519     it++;
2520   }
2521
2522   // remove exclusive (free) nodes
2523   if(removenodes)
2524   {
2525     it=s2->begin();
2526     while(it!=s2->end())
2527     {
2528       //MESSAGE( "SMDS: RM node " << (*it)->GetID() );
2529       myNodes.Remove(static_cast<SMDS_MeshNode*>
2530                     (const_cast<SMDS_MeshElement*>(*it)));
2531       myInfo.myNbNodes--;
2532       myNodeIDFactory->ReleaseID((*it)->GetID());
2533       removedNodes.push_back( (*it) );
2534       delete *it;
2535       it++;
2536     }
2537   }
2538
2539   delete s2;
2540   delete s1;
2541 }
2542
2543
2544 ///////////////////////////////////////////////////////////////////////////////
2545 ///@param elem The element to delete
2546 ///////////////////////////////////////////////////////////////////////////////
2547 void SMDS_Mesh::RemoveFreeElement(const SMDS_MeshElement * elem)
2548 {
2549   SMDSAbs_ElementType aType = elem->GetType();
2550   if (aType == SMDSAbs_Node) {
2551     // only free node can be removed by this method
2552     const SMDS_MeshNode* n = static_cast<const SMDS_MeshNode*>(elem);
2553     SMDS_ElemIteratorPtr itFe = n->GetInverseElementIterator();
2554     if (!itFe->more()) { // free node
2555       myNodes.Remove(const_cast<SMDS_MeshNode*>(n));
2556       myInfo.myNbNodes--;
2557       myNodeIDFactory->ReleaseID(elem->GetID());
2558       delete elem;
2559     }
2560   } else {
2561     if (hasConstructionEdges() || hasConstructionFaces())
2562       // this methods is only for meshes without descendants
2563       return;
2564
2565     // Remove element from <InverseElements> of its nodes
2566     SMDS_ElemIteratorPtr itn = elem->nodesIterator();
2567     while (itn->more()) {
2568       SMDS_MeshNode * n = static_cast<SMDS_MeshNode *>
2569         (const_cast<SMDS_MeshElement *>(itn->next()));
2570       n->RemoveInverseElement(elem);
2571     }
2572
2573     // in meshes without descendants elements are always free
2574     switch (aType) {
2575     case SMDSAbs_0DElement:
2576       my0DElements.Remove(static_cast<SMDS_Mesh0DElement*>
2577                           (const_cast<SMDS_MeshElement*>(elem)));
2578       //myInfo.Remove0DElement(elem);
2579       myInfo.remove(elem);
2580       break;
2581     case SMDSAbs_Edge:
2582       myEdges.Remove(static_cast<SMDS_MeshEdge*>
2583                      (const_cast<SMDS_MeshElement*>(elem)));
2584       myInfo.RemoveEdge(elem);
2585       break;
2586     case SMDSAbs_Face:
2587       myFaces.Remove(static_cast<SMDS_MeshFace*>
2588                      (const_cast<SMDS_MeshElement*>(elem)));
2589       myInfo.RemoveFace(elem);
2590       break;
2591     case SMDSAbs_Volume:
2592       myVolumes.Remove(static_cast<SMDS_MeshVolume*>
2593                        (const_cast<SMDS_MeshElement*>(elem)));
2594       myInfo.RemoveVolume(elem);
2595       break;
2596     default:
2597       break;
2598     }
2599     myElementIDFactory->ReleaseID(elem->GetID());
2600     delete elem;
2601   }
2602 }
2603
2604 /*!
2605  * Checks if the element is present in mesh.
2606  * Useful to determine dead pointers.
2607  */
2608 bool SMDS_Mesh::Contains (const SMDS_MeshElement* elem) const
2609 {
2610   // we should not imply on validity of *elem, so iterate on containers
2611   // of all types in the hope of finding <elem> somewhere there
2612   SMDS_NodeIteratorPtr itn = nodesIterator();
2613   while (itn->more())
2614     if (elem == itn->next())
2615       return true;
2616   SMDS_0DElementIteratorPtr it0d = elements0dIterator();
2617   while (it0d->more())
2618     if (elem == it0d->next())
2619       return true;
2620   SMDS_EdgeIteratorPtr ite = edgesIterator();
2621   while (ite->more())
2622     if (elem == ite->next())
2623       return true;
2624   SMDS_FaceIteratorPtr itf = facesIterator();
2625   while (itf->more())
2626     if (elem == itf->next())
2627       return true;
2628   SMDS_VolumeIteratorPtr itv = volumesIterator();
2629   while (itv->more())
2630     if (elem == itv->next())
2631       return true;
2632   return false;
2633 }
2634
2635 //=======================================================================
2636 //function : MaxNodeID
2637 //purpose  :
2638 //=======================================================================
2639
2640 int SMDS_Mesh::MaxNodeID() const
2641 {
2642   return myNodeIDFactory->GetMaxID();
2643 }
2644
2645 //=======================================================================
2646 //function : MinNodeID
2647 //purpose  :
2648 //=======================================================================
2649
2650 int SMDS_Mesh::MinNodeID() const
2651 {
2652   return myNodeIDFactory->GetMinID();
2653 }
2654
2655 //=======================================================================
2656 //function : MaxElementID
2657 //purpose  :
2658 //=======================================================================
2659
2660 int SMDS_Mesh::MaxElementID() const
2661 {
2662   return myElementIDFactory->GetMaxID();
2663 }
2664
2665 //=======================================================================
2666 //function : MinElementID
2667 //purpose  :
2668 //=======================================================================
2669
2670 int SMDS_Mesh::MinElementID() const
2671 {
2672   return myElementIDFactory->GetMinID();
2673 }
2674
2675 //=======================================================================
2676 //function : Renumber
2677 //purpose  : Renumber all nodes or elements.
2678 //=======================================================================
2679
2680 void SMDS_Mesh::Renumber (const bool isNodes, const int  startID, const int  deltaID)
2681 {
2682   if ( deltaID == 0 )
2683     return;
2684
2685   SMDS_MeshElementIDFactory * idFactory =
2686     isNodes ? myNodeIDFactory : myElementIDFactory;
2687
2688   // get existing elements in the order of ID increasing
2689   map<int,SMDS_MeshElement*> elemMap;
2690   SMDS_ElemIteratorPtr idElemIt = idFactory->elementsIterator();
2691   while ( idElemIt->more() ) {
2692     SMDS_MeshElement* elem = const_cast<SMDS_MeshElement*>(idElemIt->next());
2693     int id = elem->GetID();
2694     elemMap.insert(map<int,SMDS_MeshElement*>::value_type(id, elem));
2695   }
2696   // release their ids
2697   map<int,SMDS_MeshElement*>::iterator elemIt = elemMap.begin();
2698   idFactory->Clear();
2699 //   for ( ; elemIt != elemMap.end(); elemIt++ )
2700 //   {
2701 //     int id = (*elemIt).first;
2702 //     idFactory->ReleaseID( id );
2703 //   }
2704   // set new IDs
2705   int ID = startID;
2706   elemIt = elemMap.begin();
2707   for ( ; elemIt != elemMap.end(); elemIt++ )
2708   {
2709     idFactory->BindID( ID, (*elemIt).second );
2710     ID += deltaID;
2711   }
2712 }
2713
2714 //=======================================================================
2715 //function : GetElementType
2716 //purpose  : Return type of element or node with id
2717 //=======================================================================
2718
2719 SMDSAbs_ElementType SMDS_Mesh::GetElementType( const int id, const bool iselem ) const
2720 {
2721   SMDS_MeshElement* elem = 0;
2722   if( iselem )
2723     elem = myElementIDFactory->MeshElement( id );
2724   else
2725     elem = myNodeIDFactory->MeshElement( id );
2726
2727   if( !elem )
2728   {
2729     //throw SALOME_Exception(LOCALIZED ("this element isn't exist"));
2730     return SMDSAbs_All;
2731   }
2732   else
2733     return elem->GetType();
2734 }
2735
2736
2737
2738 //********************************************************************
2739 //********************************************************************
2740 //********                                                   *********
2741 //*****       Methods for addition of quadratic elements        ******
2742 //********                                                   *********
2743 //********************************************************************
2744 //********************************************************************
2745
2746 //=======================================================================
2747 //function : AddEdgeWithID
2748 //purpose  :
2749 //=======================================================================
2750 SMDS_MeshEdge* SMDS_Mesh::AddEdgeWithID(int n1, int n2, int n12, int ID)
2751 {
2752   return SMDS_Mesh::AddEdgeWithID
2753     ((SMDS_MeshNode*) myNodeIDFactory->MeshElement(n1),
2754      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n2),
2755      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n12),
2756      ID);
2757 }
2758
2759 //=======================================================================
2760 //function : AddEdge
2761 //purpose  :
2762 //=======================================================================
2763 SMDS_MeshEdge* SMDS_Mesh::AddEdge(const SMDS_MeshNode* n1,
2764                                   const SMDS_MeshNode* n2,
2765                                   const SMDS_MeshNode* n12)
2766 {
2767   return SMDS_Mesh::AddEdgeWithID(n1, n2, n12, myElementIDFactory->GetFreeID());
2768 }
2769
2770 //=======================================================================
2771 //function : AddEdgeWithID
2772 //purpose  :
2773 //=======================================================================
2774 SMDS_MeshEdge* SMDS_Mesh::AddEdgeWithID(const SMDS_MeshNode * n1,
2775                                         const SMDS_MeshNode * n2,
2776                                         const SMDS_MeshNode * n12,
2777                                         int ID)
2778 {
2779   if ( !n1 || !n2 || !n12 ) return 0;
2780   SMDS_QuadraticEdge* edge = new SMDS_QuadraticEdge(n1,n2,n12);
2781   if(myElementIDFactory->BindID(ID, edge)) {
2782     SMDS_MeshNode *node1,*node2, *node12;
2783     node1 = const_cast<SMDS_MeshNode*>(n1);
2784     node2 = const_cast<SMDS_MeshNode*>(n2);
2785     node12 = const_cast<SMDS_MeshNode*>(n12);
2786     node1->AddInverseElement(edge);
2787     node2->AddInverseElement(edge);
2788     node12->AddInverseElement(edge);
2789     myEdges.Add(edge);
2790     myInfo.myNbQuadEdges++;
2791     return edge;
2792   }
2793   else {
2794     delete edge;
2795     return NULL;
2796   }
2797 }
2798
2799
2800 //=======================================================================
2801 //function : AddFace
2802 //purpose  :
2803 //=======================================================================
2804 SMDS_MeshFace* SMDS_Mesh::AddFace(const SMDS_MeshNode * n1,
2805                                   const SMDS_MeshNode * n2,
2806                                   const SMDS_MeshNode * n3,
2807                                   const SMDS_MeshNode * n12,
2808                                   const SMDS_MeshNode * n23,
2809                                   const SMDS_MeshNode * n31)
2810 {
2811   return SMDS_Mesh::AddFaceWithID(n1,n2,n3,n12,n23,n31,
2812                                   myElementIDFactory->GetFreeID());
2813 }
2814
2815 //=======================================================================
2816 //function : AddFaceWithID
2817 //purpose  :
2818 //=======================================================================
2819 SMDS_MeshFace* SMDS_Mesh::AddFaceWithID(int n1, int n2, int n3,
2820                                         int n12,int n23,int n31, int ID)
2821 {
2822   return SMDS_Mesh::AddFaceWithID
2823     ((SMDS_MeshNode *)myNodeIDFactory->MeshElement(n1) ,
2824      (SMDS_MeshNode *)myNodeIDFactory->MeshElement(n2) ,
2825      (SMDS_MeshNode *)myNodeIDFactory->MeshElement(n3) ,
2826      (SMDS_MeshNode *)myNodeIDFactory->MeshElement(n12),
2827      (SMDS_MeshNode *)myNodeIDFactory->MeshElement(n23),
2828      (SMDS_MeshNode *)myNodeIDFactory->MeshElement(n31),
2829      ID);
2830 }
2831
2832 //=======================================================================
2833 //function : AddFaceWithID
2834 //purpose  :
2835 //=======================================================================
2836 SMDS_MeshFace* SMDS_Mesh::AddFaceWithID(const SMDS_MeshNode * n1,
2837                                         const SMDS_MeshNode * n2,
2838                                         const SMDS_MeshNode * n3,
2839                                         const SMDS_MeshNode * n12,
2840                                         const SMDS_MeshNode * n23,
2841                                         const SMDS_MeshNode * n31,
2842                                         int ID)
2843 {
2844   if ( !n1 || !n2 || !n3 || !n12 || !n23 || !n31) return 0;
2845   if(hasConstructionEdges()) {
2846     // creation quadratic edges - not implemented
2847     return 0;
2848   }
2849   SMDS_QuadraticFaceOfNodes* face =
2850     new SMDS_QuadraticFaceOfNodes(n1,n2,n3,n12,n23,n31);
2851   myFaces.Add(face);
2852   myInfo.myNbQuadTriangles++;
2853
2854   if (!registerElement(ID, face)) {
2855     RemoveElement(face, false);
2856     face = NULL;
2857   }
2858   return face;
2859 }
2860
2861
2862 //=======================================================================
2863 //function : AddFace
2864 //purpose  :
2865 //=======================================================================
2866 SMDS_MeshFace* SMDS_Mesh::AddFace(const SMDS_MeshNode * n1,
2867                                   const SMDS_MeshNode * n2,
2868                                   const SMDS_MeshNode * n3,
2869                                   const SMDS_MeshNode * n4,
2870                                   const SMDS_MeshNode * n12,
2871                                   const SMDS_MeshNode * n23,
2872                                   const SMDS_MeshNode * n34,
2873                                   const SMDS_MeshNode * n41)
2874 {
2875   return SMDS_Mesh::AddFaceWithID(n1,n2,n3,n4,n12,n23,n34,n41,
2876                                   myElementIDFactory->GetFreeID());
2877 }
2878
2879 //=======================================================================
2880 //function : AddFaceWithID
2881 //purpose  :
2882 //=======================================================================
2883 SMDS_MeshFace* SMDS_Mesh::AddFaceWithID(int n1, int n2, int n3, int n4,
2884                                         int n12,int n23,int n34,int n41, int ID)
2885 {
2886   return SMDS_Mesh::AddFaceWithID
2887     ((SMDS_MeshNode *)myNodeIDFactory->MeshElement(n1) ,
2888      (SMDS_MeshNode *)myNodeIDFactory->MeshElement(n2) ,
2889      (SMDS_MeshNode *)myNodeIDFactory->MeshElement(n3) ,
2890      (SMDS_MeshNode *)myNodeIDFactory->MeshElement(n4) ,
2891      (SMDS_MeshNode *)myNodeIDFactory->MeshElement(n12),
2892      (SMDS_MeshNode *)myNodeIDFactory->MeshElement(n23),
2893      (SMDS_MeshNode *)myNodeIDFactory->MeshElement(n34),
2894      (SMDS_MeshNode *)myNodeIDFactory->MeshElement(n41),
2895      ID);
2896 }
2897
2898 //=======================================================================
2899 //function : AddFaceWithID
2900 //purpose  :
2901 //=======================================================================
2902 SMDS_MeshFace* SMDS_Mesh::AddFaceWithID(const SMDS_MeshNode * n1,
2903                                         const SMDS_MeshNode * n2,
2904                                         const SMDS_MeshNode * n3,
2905                                         const SMDS_MeshNode * n4,
2906                                         const SMDS_MeshNode * n12,
2907                                         const SMDS_MeshNode * n23,
2908                                         const SMDS_MeshNode * n34,
2909                                         const SMDS_MeshNode * n41,
2910                                         int ID)
2911 {
2912   if ( !n1 || !n2 || !n3 || !n4 || !n12 || !n23 || !n34 || !n41) return 0;
2913   if(hasConstructionEdges()) {
2914     // creation quadratic edges - not implemented
2915   }
2916   SMDS_QuadraticFaceOfNodes* face =
2917     new SMDS_QuadraticFaceOfNodes(n1,n2,n3,n4,n12,n23,n34,n41);
2918   myFaces.Add(face);
2919   myInfo.myNbQuadQuadrangles++;
2920
2921   if (!registerElement(ID, face)) {
2922     RemoveElement(face, false);
2923     face = NULL;
2924   }
2925   return face;
2926 }
2927
2928
2929 //=======================================================================
2930 //function : AddVolume
2931 //purpose  :
2932 //=======================================================================
2933 SMDS_MeshVolume* SMDS_Mesh::AddVolume(const SMDS_MeshNode * n1,
2934                                       const SMDS_MeshNode * n2,
2935                                       const SMDS_MeshNode * n3,
2936                                       const SMDS_MeshNode * n4,
2937                                       const SMDS_MeshNode * n12,
2938                                       const SMDS_MeshNode * n23,
2939                                       const SMDS_MeshNode * n31,
2940                                       const SMDS_MeshNode * n14,
2941                                       const SMDS_MeshNode * n24,
2942                                       const SMDS_MeshNode * n34)
2943 {
2944   int ID = myElementIDFactory->GetFreeID();
2945   SMDS_MeshVolume * v = SMDS_Mesh::AddVolumeWithID(n1, n2, n3, n4, n12, n23,
2946                                                    n31, n14, n24, n34, ID);
2947   if(v==NULL) myElementIDFactory->ReleaseID(ID);
2948   return v;
2949 }
2950
2951 //=======================================================================
2952 //function : AddVolumeWithID
2953 //purpose  :
2954 //=======================================================================
2955 SMDS_MeshVolume* SMDS_Mesh::AddVolumeWithID(int n1, int n2, int n3, int n4,
2956                                             int n12,int n23,int n31,
2957                                             int n14,int n24,int n34, int ID)
2958 {
2959   return SMDS_Mesh::AddVolumeWithID
2960     ((SMDS_MeshNode*) myNodeIDFactory->MeshElement(n1) ,
2961      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n2) ,
2962      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n3) ,
2963      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n4) ,
2964      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n12),
2965      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n23),
2966      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n31),
2967      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n14),
2968      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n24),
2969      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n34),
2970      ID);
2971 }
2972
2973 //=======================================================================
2974 //function : AddVolumeWithID
2975 //purpose  : 2d order tetrahedron of 10 nodes
2976 //=======================================================================
2977 SMDS_MeshVolume* SMDS_Mesh::AddVolumeWithID(const SMDS_MeshNode * n1,
2978                                             const SMDS_MeshNode * n2,
2979                                             const SMDS_MeshNode * n3,
2980                                             const SMDS_MeshNode * n4,
2981                                             const SMDS_MeshNode * n12,
2982                                             const SMDS_MeshNode * n23,
2983                                             const SMDS_MeshNode * n31,
2984                                             const SMDS_MeshNode * n14,
2985                                             const SMDS_MeshNode * n24,
2986                                             const SMDS_MeshNode * n34,
2987                                             int ID)
2988 {
2989   if ( !n1 || !n2 || !n3 || !n4 || !n12 || !n23 || !n31 || !n14 || !n24 || !n34)
2990     return 0;
2991   if(hasConstructionFaces()) {
2992     // creation quadratic faces - not implemented
2993     return 0;
2994   }
2995   SMDS_QuadraticVolumeOfNodes * volume =
2996     new SMDS_QuadraticVolumeOfNodes(n1,n2,n3,n4,n12,n23,n31,n14,n24,n34);
2997   myVolumes.Add(volume);
2998   myInfo.myNbQuadTetras++;
2999
3000   if (!registerElement(ID, volume)) {
3001     RemoveElement(volume, false);
3002     volume = NULL;
3003   }
3004   return volume;
3005 }
3006
3007
3008 //=======================================================================
3009 //function : AddVolume
3010 //purpose  :
3011 //=======================================================================
3012 SMDS_MeshVolume* SMDS_Mesh::AddVolume(const SMDS_MeshNode * n1,
3013                                       const SMDS_MeshNode * n2,
3014                                       const SMDS_MeshNode * n3,
3015                                       const SMDS_MeshNode * n4,
3016                                       const SMDS_MeshNode * n5,
3017                                       const SMDS_MeshNode * n12,
3018                                       const SMDS_MeshNode * n23,
3019                                       const SMDS_MeshNode * n34,
3020                                       const SMDS_MeshNode * n41,
3021                                       const SMDS_MeshNode * n15,
3022                                       const SMDS_MeshNode * n25,
3023                                       const SMDS_MeshNode * n35,
3024                                       const SMDS_MeshNode * n45)
3025 {
3026   int ID = myElementIDFactory->GetFreeID();
3027   SMDS_MeshVolume * v =
3028     SMDS_Mesh::AddVolumeWithID(n1, n2, n3, n4, n5, n12, n23, n34, n41,
3029                                n15, n25, n35, n45, ID);
3030   if(v==NULL) myElementIDFactory->ReleaseID(ID);
3031   return v;
3032 }
3033
3034 //=======================================================================
3035 //function : AddVolumeWithID
3036 //purpose  :
3037 //=======================================================================
3038 SMDS_MeshVolume* SMDS_Mesh::AddVolumeWithID(int n1, int n2, int n3, int n4, int n5,
3039                                             int n12,int n23,int n34,int n41,
3040                                             int n15,int n25,int n35,int n45, int ID)
3041 {
3042   return SMDS_Mesh::AddVolumeWithID
3043     ((SMDS_MeshNode*) myNodeIDFactory->MeshElement(n1) ,
3044      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n2) ,
3045      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n3) ,
3046      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n4) ,
3047      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n5) ,
3048      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n12),
3049      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n23),
3050      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n34),
3051      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n41),
3052      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n15),
3053      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n25),
3054      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n35),
3055      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n45),
3056      ID);
3057 }
3058
3059 //=======================================================================
3060 //function : AddVolumeWithID
3061 //purpose  : 2d order pyramid of 13 nodes
3062 //=======================================================================
3063 SMDS_MeshVolume* SMDS_Mesh::AddVolumeWithID(const SMDS_MeshNode * n1,
3064                                             const SMDS_MeshNode * n2,
3065                                             const SMDS_MeshNode * n3,
3066                                             const SMDS_MeshNode * n4,
3067                                             const SMDS_MeshNode * n5,
3068                                             const SMDS_MeshNode * n12,
3069                                             const SMDS_MeshNode * n23,
3070                                             const SMDS_MeshNode * n34,
3071                                             const SMDS_MeshNode * n41,
3072                                             const SMDS_MeshNode * n15,
3073                                             const SMDS_MeshNode * n25,
3074                                             const SMDS_MeshNode * n35,
3075                                             const SMDS_MeshNode * n45,
3076                                             int ID)
3077 {
3078   if (!n1 || !n2 || !n3 || !n4 || !n5 || !n12 || !n23 ||
3079       !n34 || !n41 || !n15 || !n25 || !n35 || !n45)
3080     return 0;
3081   if(hasConstructionFaces()) {
3082     // creation quadratic faces - not implemented
3083     return 0;
3084   }
3085   SMDS_QuadraticVolumeOfNodes * volume =
3086     new SMDS_QuadraticVolumeOfNodes(n1,n2,n3,n4,n5,n12,n23,
3087                                     n34,n41,n15,n25,n35,n45);
3088   myVolumes.Add(volume);
3089   myInfo.myNbQuadPyramids++;
3090
3091   if (!registerElement(ID, volume)) {
3092     RemoveElement(volume, false);
3093     volume = NULL;
3094   }
3095   return volume;
3096 }
3097
3098
3099 //=======================================================================
3100 //function : AddVolume
3101 //purpose  :
3102 //=======================================================================
3103 SMDS_MeshVolume* SMDS_Mesh::AddVolume(const SMDS_MeshNode * n1,
3104                                       const SMDS_MeshNode * n2,
3105                                       const SMDS_MeshNode * n3,
3106                                       const SMDS_MeshNode * n4,
3107                                       const SMDS_MeshNode * n5,
3108                                       const SMDS_MeshNode * n6,
3109                                       const SMDS_MeshNode * n12,
3110                                       const SMDS_MeshNode * n23,
3111                                       const SMDS_MeshNode * n31,
3112                                       const SMDS_MeshNode * n45,
3113                                       const SMDS_MeshNode * n56,
3114                                       const SMDS_MeshNode * n64,
3115                                       const SMDS_MeshNode * n14,
3116                                       const SMDS_MeshNode * n25,
3117                                       const SMDS_MeshNode * n36)
3118 {
3119   int ID = myElementIDFactory->GetFreeID();
3120   SMDS_MeshVolume * v =
3121     SMDS_Mesh::AddVolumeWithID(n1, n2, n3, n4, n5, n6, n12, n23, n31,
3122                                n45, n56, n64, n14, n25, n36, ID);
3123   if(v==NULL) myElementIDFactory->ReleaseID(ID);
3124   return v;
3125 }
3126
3127 //=======================================================================
3128 //function : AddVolumeWithID
3129 //purpose  :
3130 //=======================================================================
3131 SMDS_MeshVolume* SMDS_Mesh::AddVolumeWithID(int n1, int n2, int n3,
3132                                             int n4, int n5, int n6,
3133                                             int n12,int n23,int n31,
3134                                             int n45,int n56,int n64,
3135                                             int n14,int n25,int n36, int ID)
3136 {
3137   return SMDS_Mesh::AddVolumeWithID
3138     ((SMDS_MeshNode*) myNodeIDFactory->MeshElement(n1) ,
3139      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n2) ,
3140      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n3) ,
3141      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n4) ,
3142      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n5) ,
3143      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n6) ,
3144      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n12),
3145      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n23),
3146      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n31),
3147      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n45),
3148      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n56),
3149      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n64),
3150      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n14),
3151      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n25),
3152      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n36),
3153      ID);
3154 }
3155
3156 //=======================================================================
3157 //function : AddVolumeWithID
3158 //purpose  : 2d order Pentahedron with 15 nodes
3159 //=======================================================================
3160 SMDS_MeshVolume* SMDS_Mesh::AddVolumeWithID(const SMDS_MeshNode * n1,
3161                                             const SMDS_MeshNode * n2,
3162                                             const SMDS_MeshNode * n3,
3163                                             const SMDS_MeshNode * n4,
3164                                             const SMDS_MeshNode * n5,
3165                                             const SMDS_MeshNode * n6,
3166                                             const SMDS_MeshNode * n12,
3167                                             const SMDS_MeshNode * n23,
3168                                             const SMDS_MeshNode * n31,
3169                                             const SMDS_MeshNode * n45,
3170                                             const SMDS_MeshNode * n56,
3171                                             const SMDS_MeshNode * n64,
3172                                             const SMDS_MeshNode * n14,
3173                                             const SMDS_MeshNode * n25,
3174                                             const SMDS_MeshNode * n36,
3175                                             int ID)
3176 {
3177   if (!n1 || !n2 || !n3 || !n4 || !n5 || !n6 || !n12 || !n23 ||
3178       !n31 || !n45 || !n56 || !n64 || !n14 || !n25 || !n36)
3179     return 0;
3180   if(hasConstructionFaces()) {
3181     // creation quadratic faces - not implemented
3182     return 0;
3183   }
3184   SMDS_QuadraticVolumeOfNodes * volume =
3185     new SMDS_QuadraticVolumeOfNodes(n1,n2,n3,n4,n5,n6,n12,n23,n31,
3186                                     n45,n56,n64,n14,n25,n36);
3187   myVolumes.Add(volume);
3188   myInfo.myNbQuadPrisms++;
3189
3190   if (!registerElement(ID, volume)) {
3191     RemoveElement(volume, false);
3192     volume = NULL;
3193   }
3194   return volume;
3195 }
3196
3197
3198 //=======================================================================
3199 //function : AddVolume
3200 //purpose  :
3201 //=======================================================================
3202 SMDS_MeshVolume* SMDS_Mesh::AddVolume(const SMDS_MeshNode * n1,
3203                                       const SMDS_MeshNode * n2,
3204                                       const SMDS_MeshNode * n3,
3205                                       const SMDS_MeshNode * n4,
3206                                       const SMDS_MeshNode * n5,
3207                                       const SMDS_MeshNode * n6,
3208                                       const SMDS_MeshNode * n7,
3209                                       const SMDS_MeshNode * n8,
3210                                       const SMDS_MeshNode * n12,
3211                                       const SMDS_MeshNode * n23,
3212                                       const SMDS_MeshNode * n34,
3213                                       const SMDS_MeshNode * n41,
3214                                       const SMDS_MeshNode * n56,
3215                                       const SMDS_MeshNode * n67,
3216                                       const SMDS_MeshNode * n78,
3217                                       const SMDS_MeshNode * n85,
3218                                       const SMDS_MeshNode * n15,
3219                                       const SMDS_MeshNode * n26,
3220                                       const SMDS_MeshNode * n37,
3221                                       const SMDS_MeshNode * n48)
3222 {
3223   int ID = myElementIDFactory->GetFreeID();
3224   SMDS_MeshVolume * v =
3225     SMDS_Mesh::AddVolumeWithID(n1, n2, n3, n4, n5, n6, n7, n8, n12, n23, n34, n41,
3226                                n56, n67, n78, n85, n15, n26, n37, n48, ID);
3227   if(v==NULL) myElementIDFactory->ReleaseID(ID);
3228   return v;
3229 }
3230
3231 //=======================================================================
3232 //function : AddVolumeWithID
3233 //purpose  :
3234 //=======================================================================
3235 SMDS_MeshVolume* SMDS_Mesh::AddVolumeWithID(int n1, int n2, int n3, int n4,
3236                                             int n5, int n6, int n7, int n8,
3237                                             int n12,int n23,int n34,int n41,
3238                                             int n56,int n67,int n78,int n85,
3239                                             int n15,int n26,int n37,int n48, int ID)
3240 {
3241   return SMDS_Mesh::AddVolumeWithID
3242     ((SMDS_MeshNode*) myNodeIDFactory->MeshElement(n1),
3243      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n2),
3244      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n3),
3245      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n4),
3246      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n5),
3247      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n6),
3248      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n7),
3249      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n8),
3250      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n12),
3251      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n23),
3252      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n34),
3253      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n41),
3254      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n56),
3255      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n67),
3256      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n78),
3257      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n85),
3258      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n15),
3259      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n26),
3260      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n37),
3261      (SMDS_MeshNode*) myNodeIDFactory->MeshElement(n48),
3262      ID);
3263 }
3264
3265 //=======================================================================
3266 //function : AddVolumeWithID
3267 //purpose  : 2d order Hexahedrons with 20 nodes
3268 //=======================================================================
3269 SMDS_MeshVolume* SMDS_Mesh::AddVolumeWithID(const SMDS_MeshNode * n1,
3270                                             const SMDS_MeshNode * n2,
3271                                             const SMDS_MeshNode * n3,
3272                                             const SMDS_MeshNode * n4,
3273                                             const SMDS_MeshNode * n5,
3274                                             const SMDS_MeshNode * n6,
3275                                             const SMDS_MeshNode * n7,
3276                                             const SMDS_MeshNode * n8,
3277                                             const SMDS_MeshNode * n12,
3278                                             const SMDS_MeshNode * n23,
3279                                             const SMDS_MeshNode * n34,
3280                                             const SMDS_MeshNode * n41,
3281                                             const SMDS_MeshNode * n56,
3282                                             const SMDS_MeshNode * n67,
3283                                             const SMDS_MeshNode * n78,
3284                                             const SMDS_MeshNode * n85,
3285                                             const SMDS_MeshNode * n15,
3286                                             const SMDS_MeshNode * n26,
3287                                             const SMDS_MeshNode * n37,
3288                                             const SMDS_MeshNode * n48,
3289                                             int ID)
3290 {
3291   if (!n1 || !n2 || !n3 || !n4 || !n5 || !n6 || !n7 || !n8 || !n12 || !n23 ||
3292       !n34 || !n41 || !n56 || !n67 || !n78 || !n85 || !n15 || !n26 || !n37 || !n48)
3293     return 0;
3294   if(hasConstructionFaces()) {
3295     return 0;
3296     // creation quadratic faces - not implemented
3297   }
3298   SMDS_QuadraticVolumeOfNodes * volume =
3299     new SMDS_QuadraticVolumeOfNodes(n1,n2,n3,n4,n5,n6,n7,n8,n12,n23,n34,n41,
3300                                     n56,n67,n78,n85,n15,n26,n37,n48);
3301   myVolumes.Add(volume);
3302   myInfo.myNbQuadHexas++;
3303
3304   if (!registerElement(ID, volume)) {
3305     RemoveElement(volume, false);
3306     volume = NULL;
3307   }
3308   return volume;
3309 }