1 // Copyright (C) 2007-2008 CEA/DEN, EDF R&D, OPEN CASCADE
3 // Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
6 // This library is free software; you can redistribute it and/or
7 // modify it under the terms of the GNU Lesser General Public
8 // License as published by the Free Software Foundation; either
9 // version 2.1 of the License.
11 // This library is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 // Lesser General Public License for more details.
16 // You should have received a copy of the GNU Lesser General Public
17 // License along with this library; if not, write to the Free Software
18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
22 // SMESH SMDS : implementaion of Salome mesh data structure
25 #pragma warning(disable:4786)
28 #include "SMDS_VolumeOfNodes.hxx"
29 #include "SMDS_MeshNode.hxx"
30 #include "SMDS_SetIterator.hxx"
31 #include "SMDS_VolumeTool.hxx"
32 #include "SMDS_Mesh.hxx"
33 #include "utilities.h"
41 ///////////////////////////////////////////////////////////////////////////////
42 /// Create an hexahedron. node 1,2,3,4 and 5,6,7,8 are quadrangle and
43 /// 5,1 and 7,3 are an edges.
44 ///////////////////////////////////////////////////////////////////////////////
45 SMDS_VolumeOfNodes::SMDS_VolumeOfNodes(
46 const SMDS_MeshNode * node1,
47 const SMDS_MeshNode * node2,
48 const SMDS_MeshNode * node3,
49 const SMDS_MeshNode * node4,
50 const SMDS_MeshNode * node5,
51 const SMDS_MeshNode * node6,
52 const SMDS_MeshNode * node7,
53 const SMDS_MeshNode * node8)
56 myNodes = new const SMDS_MeshNode* [myNbNodes];
67 SMDS_VolumeOfNodes::SMDS_VolumeOfNodes(
68 const SMDS_MeshNode * node1,
69 const SMDS_MeshNode * node2,
70 const SMDS_MeshNode * node3,
71 const SMDS_MeshNode * node4)
74 myNodes = new const SMDS_MeshNode* [myNbNodes];
81 SMDS_VolumeOfNodes::SMDS_VolumeOfNodes(
82 const SMDS_MeshNode * node1,
83 const SMDS_MeshNode * node2,
84 const SMDS_MeshNode * node3,
85 const SMDS_MeshNode * node4,
86 const SMDS_MeshNode * node5)
89 myNodes = new const SMDS_MeshNode* [myNbNodes];
97 SMDS_VolumeOfNodes::SMDS_VolumeOfNodes(
98 const SMDS_MeshNode * node1,
99 const SMDS_MeshNode * node2,
100 const SMDS_MeshNode * node3,
101 const SMDS_MeshNode * node4,
102 const SMDS_MeshNode * node5,
103 const SMDS_MeshNode * node6)
106 myNodes = new const SMDS_MeshNode* [myNbNodes];
115 bool SMDS_VolumeOfNodes::ChangeNodes(const SMDS_MeshNode* nodes[],
118 if (nbNodes < 4 || nbNodes > 8 || nbNodes == 7)
123 myNodes = new const SMDS_MeshNode* [myNbNodes];
124 for ( int i = 0; i < nbNodes; i++ )
125 myNodes[ i ] = nodes [ i ];
130 SMDS_VolumeOfNodes::~SMDS_VolumeOfNodes()
132 if (myNodes != NULL) {
138 void SMDS_VolumeOfNodes::Print(ostream & OS) const
140 OS << "volume <" << GetID() << "> : ";
142 for (i = 0; i < NbNodes()-1; ++i) OS << myNodes[i] << ",";
143 OS << myNodes[NbNodes()-1]<< ") " << endl;
146 int SMDS_VolumeOfNodes::NbFaces() const
154 default: MESSAGE("invalid number of nodes");
159 int SMDS_VolumeOfNodes::NbNodes() const
164 int SMDS_VolumeOfNodes::NbEdges() const
172 default: MESSAGE("invalid number of nodes");
178 * \brief Iterator on node of volume
180 class SMDS_VolumeOfNodes_MyIterator:public SMDS_NodeArrayElemIterator
183 SMDS_VolumeOfNodes_MyIterator(const SMDS_MeshNode* const* s, int l):
184 SMDS_NodeArrayElemIterator( s, & s[ l ]) {}
188 * \brief Iterator on faces or edges of volume
190 class _MySubIterator : public SMDS_ElemIterator
192 vector< const SMDS_MeshElement* > myElems;
195 _MySubIterator(const SMDS_VolumeOfNodes* vol, SMDSAbs_ElementType type):myIndex(0) {
196 SMDS_VolumeTool vTool(vol);
197 if (type == SMDSAbs_Face)
198 vTool.GetAllExistingFaces( myElems );
200 vTool.GetAllExistingFaces( myElems );
202 /// Return true if and only if there are other object in this iterator
203 virtual bool more() { return myIndex < myElems.size(); }
205 /// Return the current object and step to the next one
206 virtual const SMDS_MeshElement* next() { return myElems[ myIndex++ ]; }
209 SMDS_ElemIteratorPtr SMDS_VolumeOfNodes::elementsIterator(SMDSAbs_ElementType type) const
214 return SMDS_MeshElement::elementsIterator(SMDSAbs_Volume);
216 return SMDS_ElemIteratorPtr(new SMDS_VolumeOfNodes_MyIterator(myNodes,myNbNodes));
218 return SMDS_ElemIteratorPtr(new _MySubIterator(this,SMDSAbs_Face));
220 return SMDS_ElemIteratorPtr(new _MySubIterator(this,SMDSAbs_Edge));
222 MESSAGE("ERROR : Iterator not implemented");
223 return SMDS_ElemIteratorPtr((SMDS_ElemIterator*)NULL);
227 SMDSAbs_ElementType SMDS_VolumeOfNodes::GetType() const
229 return SMDSAbs_Volume;
233 * \brief Return node by its index
234 * \param ind - node index
235 * \retval const SMDS_MeshNode* - the node
237 const SMDS_MeshNode* SMDS_VolumeOfNodes::GetNode(const int ind) const
239 return myNodes[ ind ];
242 SMDSAbs_EntityType SMDS_VolumeOfNodes::GetEntityType() const
244 SMDSAbs_EntityType aType = SMDSEntity_Tetra;
247 case 4: aType = SMDSEntity_Tetra; break;
248 case 5: aType = SMDSEntity_Pyramid; break;
249 case 6: aType = SMDSEntity_Penta; break;
251 default: aType = SMDSEntity_Hexa; break;
257 // *************************************************************************************
260 SMDS_VolumeVtkNodes::SMDS_VolumeVtkNodes()
264 SMDS_VolumeVtkNodes::SMDS_VolumeVtkNodes(std::vector<vtkIdType> nodeIds, SMDS_Mesh* mesh)
269 void SMDS_VolumeVtkNodes::init(std::vector<vtkIdType> nodeIds, SMDS_Mesh* mesh)
271 vtkUnstructuredGrid* grid = mesh->getGrid();
272 myVtkID = grid->InsertNextLinkedCell(GetType(), nodeIds.size(), &nodeIds[0]);
275 bool SMDS_VolumeVtkNodes::ChangeNodes(const SMDS_MeshNode* nodes[],
278 // utilise dans SMDS_Mesh
282 SMDS_VolumeVtkNodes::~SMDS_VolumeVtkNodes()
286 void SMDS_VolumeVtkNodes::Print(ostream & OS) const
288 OS << "volume <" << GetID() << "> : ";
291 int SMDS_VolumeVtkNodes::NbFaces() const
299 default: MESSAGE("invalid number of nodes");
304 int SMDS_VolumeVtkNodes::NbNodes() const
306 vtkUnstructuredGrid* grid =SMDS_Mesh::_meshList[myMeshId]->getGrid();
307 int nbPoints = grid->GetCell(myVtkID)->GetNumberOfPoints();
311 int SMDS_VolumeVtkNodes::NbEdges() const
319 default: MESSAGE("invalid number of nodes");
324 SMDS_ElemIteratorPtr SMDS_VolumeVtkNodes::elementsIterator(SMDSAbs_ElementType type) const
328 // case SMDSAbs_Volume:
329 // return SMDS_MeshElement::elementsIterator(SMDSAbs_Volume);
330 // case SMDSAbs_Node:
331 // return SMDS_ElemIteratorPtr(new SMDS_VolumeVtkNodes_MyIterator(myNodes,myNbNodes));
332 // case SMDSAbs_Face:
333 // return SMDS_ElemIteratorPtr(new _MySubWNIterator(this,SMDSAbs_Face));
334 // case SMDSAbs_Edge:
335 // return SMDS_ElemIteratorPtr(new _MySubWNIterator(this,SMDSAbs_Edge));
337 MESSAGE("ERROR : Iterator not implemented");
338 return SMDS_ElemIteratorPtr((SMDS_ElemIterator*)NULL);
342 SMDSAbs_ElementType SMDS_VolumeVtkNodes::GetType() const
344 return SMDSAbs_Volume;
348 * \brief Return node by its index
349 * \param ind - node index
350 * \retval const SMDS_MeshNode* - the node
352 const SMDS_MeshNode* SMDS_VolumeVtkNodes::GetNode(const int ind) const
357 SMDSAbs_EntityType SMDS_VolumeVtkNodes::GetEntityType() const
359 SMDSAbs_EntityType aType = SMDSEntity_Tetra;
362 case 4: aType = SMDSEntity_Tetra; break;
363 case 5: aType = SMDSEntity_Pyramid; break;
364 case 6: aType = SMDSEntity_Penta; break;
366 default: aType = SMDSEntity_Hexa; break;