-// SMESH SMDS : implementaion of Salome mesh data structure
+// Copyright (C) 2007-2011 CEA/DEN, EDF R&D, OPEN CASCADE
+//
+// Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
+// CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
//
-// Copyright (C) 2003 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
-// CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
-//
-// This library is free software; you can redistribute it and/or
-// modify it under the terms of the GNU Lesser General Public
-// License as published by the Free Software Foundation; either
-// version 2.1 of the License.
-//
-// This library is distributed in the hope that it will be useful,
-// but WITHOUT ANY WARRANTY; without even the implied warranty of
-// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
-// Lesser General Public License for more details.
-//
-// You should have received a copy of the GNU Lesser General Public
-// License along with this library; if not, write to the Free Software
-// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
-//
-// See http://www.opencascade.org/SALOME/ or email : webmaster.salome@opencascade.org
+// This library is free software; you can redistribute it and/or
+// modify it under the terms of the GNU Lesser General Public
+// License as published by the Free Software Foundation; either
+// version 2.1 of the License.
//
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
+// Lesser General Public License for more details.
//
+// You should have received a copy of the GNU Lesser General Public
+// License along with this library; if not, write to the Free Software
+// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
//
+// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
+//
+
+// SMESH SMDS : implementaion of Salome mesh data structure
// File : SMDS_VolumeTool.hxx
// Module : SMESH
// Created : Tue Jul 13 11:27:17 2004
// Author : Edward AGAPOV (eap)
-
-
+//
#ifndef SMDS_VolumeTool_HeaderFile
#define SMDS_VolumeTool_HeaderFile
+#include "SMESH_SMDS.hxx"
+
class SMDS_MeshElement;
class SMDS_MeshNode;
+class SMDS_VtkVolume;
+class SMDS_MeshVolume;
#include <vector>
#include <set>
//
// =========================================================================
-class SMDS_VolumeTool
+class SMDS_EXPORT SMDS_VolumeTool
{
public:
+ enum VolumeType { UNKNOWN = -1, TETRA = 0, PYRAM, PENTA, HEXA, QUAD_TETRA,
+ QUAD_PYRAM, QUAD_PENTA, QUAD_HEXA, POLYHEDA };
+
SMDS_VolumeTool ();
~SMDS_VolumeTool ();
SMDS_VolumeTool (const SMDS_MeshElement* theVolume);
// Set volume.
// Return false if theVolume is not of type SMDSAbs_Volume
+ const SMDS_MeshVolume* Element() const;
+ // return element
+
+ int ID() const;
+ // return element ID
+
// -----------------------
// general info
// -----------------------
+ VolumeType GetVolumeType() const;
+
bool IsForward() const { return myVolForward; }
// Check volume orientation. can be changed by Inverse().
- // See node order of forward volumes at file bottom
+ // See node order of forward volumes at the file bottom
void Inverse();
// Change nodes order as if the volume changes its orientation:
bool GetBaryCenter (double & X, double & Y, double & Z) const;
+ bool IsOut(double X, double Y, double Z, double tol) const;
+ // Classify a point
// -----------------------
// info on node connection
// -----------------------
bool IsLinked (const SMDS_MeshNode* theNode1,
- const SMDS_MeshNode* theNode2) const;
+ const SMDS_MeshNode* theNode2,
+ const bool theIgnoreMediumNodes=false) const;
// Return true if theNode1 is linked with theNode2.
+ // If theIgnoreMediumNodes then corner nodes of quadratic cell are considered linked as well
bool IsLinked (const int theNode1Index,
- const int theNode2Index) const;
+ const int theNode2Index,
+ bool theIgnoreMediumNodes=false) const;
// Return true if the node with theNode1Index is linked
// with the node with theNode2Index
+ // If theIgnoreMediumNodes then corner nodes of quadratic cell are considered linked as well
int GetNodeIndex(const SMDS_MeshNode* theNode) const;
// Return an index of theNode
+ int GetAllExistingEdges(std::vector<const SMDS_MeshElement*> & edges) const;
+ // Fill vector with boundary edges existing in the mesh
+
+ double MinLinearSize2() const;
+ // Return minimal square distance between connected corner nodes
+
// -------------
// info on faces
// -------------
- void SetForwardOrientation ();
- // Node order in faces will be as for forward orientation
-
void SetExternalNormal ();
// Node order in faces will be so that faces normals are external.
- // It overrides SetForwardOrientation()
int NbFaces() const { return myNbFaces; }
// Return number of faces of the volume. In the following
// methods 0 <= faceIndex < NbFaces()
- int NbFaceNodes( int faceIndex );
+ int NbFaceNodes( int faceIndex ) const;
// Return number of nodes in the array of face nodes
- const int* GetFaceNodesIndices( int faceIndex );
+ const int* GetFaceNodesIndices( int faceIndex ) const;
// Return the array of face nodes indices
// To comfort link iteration, the array
// length == NbFaceNodes( faceIndex ) + 1 and
// the last node index == the first one.
+ // NOTE: for the quadratic volume, node indoces are in the order the nodes encounter
+ // in face boundary and not the order they are in the mesh face
- const SMDS_MeshNode** GetFaceNodes( int faceIndex );
+ const SMDS_MeshNode** GetFaceNodes( int faceIndex ) const;
// Return the array of face nodes.
// To comfort link iteration, the array
// length == NbFaceNodes( faceIndex ) + 1 and
// the last node == the first one.
+ // NOTE: for the quadratic volume, nodes are in the order they encounter in face boundary
+ // and not the order they are in the mesh face
// WARNING: do not modify the array, some methods
// work basing on its contents
bool GetFaceNodes (int faceIndex,
- std::set<const SMDS_MeshNode*>& theFaceNodes );
+ std::set<const SMDS_MeshNode*>& theFaceNodes ) const;
// Return a set of face nodes.
- bool IsFaceExternal( int faceIndex );
+ bool IsFaceExternal( int faceIndex ) const;
// Check normal orientation of a face.
- // SetForwardOrientation() and SetForwardOrientation() are taken
- // into account.
+ // SetExternalNormal() is taken into account.
- bool IsFreeFace( int faceIndex );
+ bool IsFreeFace( int faceIndex, const SMDS_MeshElement** otherVol=0 ) const;
// Check that all volumes built on the face nodes lays on one side
+ // otherVol returns another volume sharing the given facet
- bool GetFaceNormal (int faceIndex, double & X, double & Y, double & Z);
+ bool GetFaceNormal (int faceIndex, double & X, double & Y, double & Z) const;
// Return a normal to a face
- double GetFaceArea( int faceIndex );
+ bool GetFaceBaryCenter (int faceIndex, double & X, double & Y, double & Z) const;
+ // Return barycenter of a face
+
+ double GetFaceArea( int faceIndex ) const;
// Return face area
int GetOppFaceIndex( int faceIndex ) const;
// Return index of the opposite face if it exists, else -1.
- int GetFaceIndex( const std::set<const SMDS_MeshNode*>& theFaceNodes );
+ int GetFaceIndex( const std::set<const SMDS_MeshNode*>& theFaceNodes ) const;
// Return index of a face formed by theFaceNodes.
// Return -1 if a face not found
- int GetFaceIndex( const std::set<int>& theFaceNodesIndices );
+ //int GetFaceIndex( const std::set<int>& theFaceNodesIndices );
// Return index of a face formed by theFaceNodesIndices
// Return -1 if a face not found
+ int GetAllExistingFaces(std::vector<const SMDS_MeshElement*> & faces) const;
+ // Fill vector with boundary faces existing in the mesh
+
+ // ------------------------
+ // static methods for faces
+ // ------------------------
+
+ static VolumeType GetType(int nbNodes);
+ // return VolumeType by nb of nodes in a volume
- private:
+ static int NbFaces( VolumeType type );
+ // return nb of faces by volume type
- bool setFace( int faceIndex );
+ static const int* GetFaceNodesIndices(VolumeType type,
+ int faceIndex,
+ bool external);
+ // Return the array of face nodes indices
+ // To comfort link iteration, the array
+ // length == NbFaceNodes( faceIndex ) + 1 and
+ // the last node index == the first one.
+
+ static int NbFaceNodes(VolumeType type,
+ int faceIndex );
+ // Return number of nodes in the array of face nodes
+
+ static int NbCornerNodes(VolumeType type);
+ // Useful to know nb of corner nodes of a quadratic volume
+
+private:
+
+ bool setFace( int faceIndex ) const;
const SMDS_MeshElement* myVolume;
+ const SMDS_VtkVolume* myPolyedre;
+
bool myVolForward;
int myNbFaces;
int myVolumeNbNodes;
- const SMDS_MeshNode* myVolumeNodes[ 8 ];
+ const SMDS_MeshNode** myVolumeNodes;
+ std::vector< int > myPolyIndices;
+
+ mutable bool myExternalFaces;
- bool myForwardFaces;
- bool myExternalFaces;
- int* myFaceNodeIndices;
- int* myFaceNbNodes;
- const SMDS_MeshNode* myFaceNodes[ 5 ];
- int myCurFace;
+ mutable int myCurFace;
+ mutable int myFaceNbNodes;
+ mutable int* myFaceNodeIndices;
+ mutable const SMDS_MeshNode** myFaceNodes;
};
#endif
// /|\
// / | \
// / | \
-// N0 +---|---+ N2 TETRAHEDRON
+// N0 +---|---+ N1 TETRAHEDRON
// \ | /
// \ | /
// \ | /
// \|/
// +
-// N1
+// N2
// + N4
// /|\
// |/ \|
// N0 +---------+ N2
-// N7+----------+N6
+// N5+----------+N6
// /| /|
// / | / |
// / | / |
-// N4+----------+N5 |
+// N4+----------+N7 |
// | | | | HEXAHEDRON
// | | | |
// | | | |
-// | N3+------|---+N2
+// | N1+------|---+N2
// | / | /
// | / | /
// |/ |/
-// N0+----------+N1
+// N0+----------+N3
//
*/