Salome HOME
IMP 23373: [CEA 1170] Optimization of a 3D mesh using MG-Tetra
[modules/smesh.git] / src / SMDS / SMDS_MeshElementIDFactory.cxx
1 // Copyright (C) 2007-2016  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, or (at your option) any later version.
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
23 //  SMESH SMDS : implementaion of Salome mesh data structure
24 //  File   : SMDS_MeshElementIDFactory.cxx
25 //  Author : Jean-Michel BOULCOURT
26 //  Module : SMESH
27 //
28 #ifdef _MSC_VER
29 #pragma warning(disable:4786)
30 #endif
31
32 #include "SMDS_MeshElementIDFactory.hxx"
33 #include "SMDS_MeshElement.hxx"
34 #include "SMDS_Mesh.hxx"
35
36 #include "utilities.h"
37
38 #include "SMDS_UnstructuredGrid.hxx"
39 #include <vtkCellType.h>
40
41 #include <climits>
42
43 using namespace std;
44
45 //=======================================================================
46 //function : SMDS_MeshElementIDFactory
47 //purpose  : 
48 //=======================================================================
49 SMDS_MeshElementIDFactory::SMDS_MeshElementIDFactory():
50   SMDS_MeshNodeIDFactory()
51 {
52 }
53
54 int SMDS_MeshElementIDFactory::SetInVtkGrid(SMDS_MeshElement * elem)
55 {
56   // --- retrieve nodes ID
57
58   SMDS_MeshCell *cell = dynamic_cast<SMDS_MeshCell*>(elem);
59   assert(cell);
60   vector<vtkIdType> nodeIds( elem->NbNodes() );
61   SMDS_ElemIteratorPtr it = elem->nodesIterator();
62   for( int i = 0; it->more(); ++i )
63   {
64     int nodeId = (static_cast<const SMDS_MeshNode*>(it->next()))->getVtkId();
65     nodeIds[i] = nodeId;
66   }
67
68   // --- insert cell in vtkUnstructuredGrid
69
70   int typ = VTK_VERTEX;
71   int cellId = myMesh->getGrid()->InsertNextLinkedCell(typ, nodeIds.size(), &nodeIds[0]);
72   cell->setVtkId(cellId);
73   return cellId;
74 }
75
76 //=======================================================================
77 //function : BindID
78 //purpose  :
79 //=======================================================================
80
81 bool SMDS_MeshElementIDFactory::BindID(int ID, SMDS_MeshElement * elem)
82 {
83   MESSAGE("SMDS_MeshElementIDFactory::BindID " << ID);
84   SetInVtkGrid(elem);
85   return myMesh->registerElement(ID, elem);
86 }
87
88 //=======================================================================
89 //function : MeshElement
90 //purpose  : 
91 //=======================================================================
92 SMDS_MeshElement* SMDS_MeshElementIDFactory::MeshElement(int ID)
93 {
94   if ( ID<1 || ID >= (int) myMesh->myCells.size() )
95     return NULL;
96   const SMDS_MeshElement* elem = GetMesh()->FindElement(ID);
97   return (SMDS_MeshElement*)(elem);
98 }
99
100 //=======================================================================
101 //function : GetFreeID
102 //purpose  : 
103 //=======================================================================
104
105 int SMDS_MeshElementIDFactory::GetFreeID()
106 {
107   int ID;
108   do {
109     ID = SMDS_MeshIDFactory::GetFreeID();
110   } while ( MeshElement( ID ));
111   return ID;
112 }
113
114 //=======================================================================
115 //function : ReleaseID
116 //purpose  : 
117 //=======================================================================
118 void SMDS_MeshElementIDFactory::ReleaseID(int ID, int vtkId)
119 {
120   if (ID < 1) // TODO check case ID == O
121   {
122     MESSAGE("~~~~~~~~~~~~~~ SMDS_MeshElementIDFactory::ReleaseID ID = " << ID);
123     return;
124   }
125   if (vtkId >= 0)
126   {
127     assert(vtkId < (int)myMesh->myCellIdVtkToSmds.size());
128     myMesh->myCellIdVtkToSmds[vtkId] = -1;
129     myMesh->setMyModified();
130   }
131   SMDS_MeshIDFactory::ReleaseID(ID);
132   if (ID == myMax)
133     myMax = 0;
134   if (ID == myMin)
135     myMax = 0;
136 }
137
138 //=======================================================================
139 //function : updateMinMax
140 //purpose  :
141 //=======================================================================
142
143 void SMDS_MeshElementIDFactory::updateMinMax() const
144 {
145   myMin = INT_MAX;
146   myMax = 0;
147   for (size_t i = 0; i < myMesh->myCells.size(); i++)
148   {
149     if (myMesh->myCells[i])
150     {
151       int id = myMesh->myCells[i]->GetID();
152       if (id > myMax)
153         myMax = id;
154       if (id < myMin)
155         myMin = id;
156     }
157   }
158   if (myMin == INT_MAX)
159     myMin = 0;
160 }
161
162 //=======================================================================
163 //function : elementsIterator
164 //purpose  : Return an iterator on elements of the factory
165 //=======================================================================
166
167 SMDS_ElemIteratorPtr SMDS_MeshElementIDFactory::elementsIterator() const
168 {
169   return myMesh->elementsIterator(SMDSAbs_All);
170 }
171
172 void SMDS_MeshElementIDFactory::Clear()
173 {
174   myMesh->myCellIdVtkToSmds.clear();
175   myMin = myMax = 0;
176   SMDS_MeshIDFactory::Clear();
177 }