1 // Copyright (C) 2007-2010 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
23 // SMESH SMDS : implementaion of Salome mesh data structure
26 #pragma warning(disable:4786)
29 #include "SMDS_MeshElement.hxx"
30 #include "SMDS_MeshNode.hxx"
31 #include "SMDS_MeshEdge.hxx"
32 #include "SMDS_MeshFace.hxx"
33 #include "SMDS_MeshVolume.hxx"
34 #include "utilities.h"
38 SMDS_MeshElement::SMDS_MeshElement(int ID):myID(ID), myMeshId(-1), myShapeId(0), myIdInShape(-1)
42 SMDS_MeshElement::SMDS_MeshElement(int id, ShortType meshId, ShortType shapeId):
43 myID(id), myMeshId(meshId), myShapeId(shapeId), myIdInShape(-1)
47 void SMDS_MeshElement::Print(ostream & OS) const
49 OS << "dump of mesh element" << endl;
52 ostream & operator <<(ostream & OS, const SMDS_MeshElement * ME)
58 ///////////////////////////////////////////////////////////////////////////////
59 /// Create an iterator which iterate on nodes owned by the element.
60 /// This method call elementsIterator().
61 ///////////////////////////////////////////////////////////////////////////////
62 SMDS_ElemIteratorPtr SMDS_MeshElement::nodesIterator() const
64 return elementsIterator(SMDSAbs_Node);
67 ///////////////////////////////////////////////////////////////////////////////
68 /// Create an iterator which iterate on edges linked with or owned by the element.
69 /// This method call elementsIterator().
70 ///////////////////////////////////////////////////////////////////////////////
71 SMDS_ElemIteratorPtr SMDS_MeshElement::edgesIterator() const
73 return elementsIterator(SMDSAbs_Edge);
76 ///////////////////////////////////////////////////////////////////////////////
77 /// Create an iterator which iterate on faces linked with or owned by the element.
78 /// This method call elementsIterator().
79 ///////////////////////////////////////////////////////////////////////////////
80 SMDS_ElemIteratorPtr SMDS_MeshElement::facesIterator() const
82 return elementsIterator(SMDSAbs_Face);
85 ///////////////////////////////////////////////////////////////////////////////
86 ///Return The number of nodes owned by the current element
87 ///////////////////////////////////////////////////////////////////////////////
88 int SMDS_MeshElement::NbNodes() const
91 SMDS_ElemIteratorPtr it=nodesIterator();
100 ///////////////////////////////////////////////////////////////////////////////
101 ///Return the number of edges owned by or linked with the current element
102 ///////////////////////////////////////////////////////////////////////////////
103 int SMDS_MeshElement::NbEdges() const
106 SMDS_ElemIteratorPtr it=edgesIterator();
115 ///////////////////////////////////////////////////////////////////////////////
116 ///Return the number of faces owned by or linked with the current element
117 ///////////////////////////////////////////////////////////////////////////////
118 int SMDS_MeshElement::NbFaces() const
121 SMDS_ElemIteratorPtr it=facesIterator();
130 ///////////////////////////////////////////////////////////////////////////////
131 ///Create an iterator which iterate on elements linked with the current element.
132 ///@param type The of elements on which you want to iterate
133 ///@return A smart pointer to iterator, you are not to take care of freeing memory
134 ///////////////////////////////////////////////////////////////////////////////
135 class SMDS_MeshElement_MyIterator:public SMDS_ElemIterator
137 const SMDS_MeshElement * myElement;
140 SMDS_MeshElement_MyIterator(const SMDS_MeshElement * element):
141 myElement(element),myMore(true) {}
148 const SMDS_MeshElement* next()
155 SMDS_ElemIteratorPtr SMDS_MeshElement::
156 elementsIterator(SMDSAbs_ElementType type) const
158 /** @todo Check that iterator in the child classes return elements
159 in the same order for each different implementation (i.e: SMDS_VolumeOfNodes
160 and SMDS_VolumeOfFaces */
163 return SMDS_ElemIteratorPtr(new SMDS_MeshElement_MyIterator(this));
166 MESSAGE("Iterator not implemented");
167 return SMDS_ElemIteratorPtr((SMDS_ElemIterator*)NULL);
171 //! virtual, redefined in vtkEdge, vtkFace and vtkVolume classes
172 SMDS_ElemIteratorPtr SMDS_MeshElement::nodesIteratorToUNV() const
174 MESSAGE("Iterator not implemented");
175 return SMDS_ElemIteratorPtr((SMDS_ElemIterator*) NULL);
178 //! virtual, redefined in vtkEdge, vtkFace and vtkVolume classes
179 SMDS_ElemIteratorPtr SMDS_MeshElement::interlacedNodesElemIterator() const
181 MESSAGE("Iterator not implemented");
182 return SMDS_ElemIteratorPtr((SMDS_ElemIterator*) NULL);
185 bool operator<(const SMDS_MeshElement& e1, const SMDS_MeshElement& e2)
187 if(e1.GetType()!=e2.GetType()) return false;
191 return static_cast<const SMDS_MeshNode &>(e1) <
192 static_cast<const SMDS_MeshNode &>(e2);
195 return static_cast<const SMDS_MeshEdge &>(e1) <
196 static_cast<const SMDS_MeshEdge &>(e2);
199 return static_cast<const SMDS_MeshFace &>(e1) <
200 static_cast<const SMDS_MeshFace &>(e2);
203 return static_cast<const SMDS_MeshVolume &>(e1) <
204 static_cast<const SMDS_MeshVolume &>(e2);
206 default : MESSAGE("Internal Error");
211 bool SMDS_MeshElement::IsValidIndex(const int ind) const
213 return ( ind>-1 && ind<NbNodes() );
216 const SMDS_MeshNode* SMDS_MeshElement::GetNode(const int ind) const
219 SMDS_ElemIteratorPtr it = nodesIterator();
220 for ( int i = 0; i < ind; ++i )
223 return static_cast<const SMDS_MeshNode*> (it->next());
228 bool SMDS_MeshElement::IsQuadratic() const
233 bool SMDS_MeshElement::IsMediumNode(const SMDS_MeshNode* node) const
238 //================================================================================
240 * \brief Return number of nodes excluding medium ones
242 //================================================================================
244 int SMDS_MeshElement::NbCornerNodes() const
246 return IsQuadratic() ? NbNodes() - NbEdges() : NbNodes();
249 //================================================================================
251 * \brief Check if a node belongs to the element
252 * \param node - the node to check
253 * \retval int - node index within the element, -1 if not found
255 //================================================================================
257 int SMDS_MeshElement::GetNodeIndex( const SMDS_MeshNode* node ) const
259 SMDS_ElemIteratorPtr nIt = nodesIterator();
260 for ( int i = 0; nIt->more(); ++i )
261 if ( nIt->next() == node )