+++ /dev/null
-// File: StdMeshers_Helper.cxx
-// Created: 15.02.06 15:22:41
-// Author: Sergey KUUL
-// Copyright: Open CASCADE 2006
-
-
-#include "StdMeshers_Helper.hxx"
-
-#include "SMDS_FacePosition.hxx"
-#include "SMDS_EdgePosition.hxx"
-#include "SMESH_subMesh.hxx"
-#include "SMESH_MeshEditor.hxx"
-
-#include <BRepAdaptor_Surface.hxx>
-#include <BRepTools.hxx>
-#include <BRep_Tool.hxx>
-#include <Geom2d_Curve.hxx>
-#include <Geom_Curve.hxx>
-#include <Geom_Surface.hxx>
-#include <TopExp_Explorer.hxx>
-#include <TopTools_MapOfShape.hxx>
-#include <gp_Pnt2d.hxx>
-
-
-//=======================================================================
-//function : CheckShape
-//purpose :
-//=======================================================================
-
-bool StdMeshers_Helper::IsQuadraticSubMesh(const TopoDS_Shape& aSh)
-{
- SMESHDS_Mesh* meshDS = GetMesh()->GetMeshDS();
- myShapeID = meshDS->ShapeToIndex(aSh);
- // we can create quadratic elements only if all elements
- // created on subshapes of given shape are quadratic
- // also we have to fill myNLinkNodeMap
- myCreateQuadratic = true;
- mySeamShapeIds.clear();
- TopAbs_ShapeEnum subType( aSh.ShapeType()==TopAbs_FACE ? TopAbs_EDGE : TopAbs_FACE );
- SMDSAbs_ElementType elemType( subType==TopAbs_FACE ? SMDSAbs_Face : SMDSAbs_Edge );
-
- TopExp_Explorer exp( aSh, subType );
- for (; exp.More() && myCreateQuadratic; exp.Next()) {
- if ( SMESHDS_SubMesh * subMesh = meshDS->MeshElements( exp.Current() )) {
- if ( SMDS_ElemIteratorPtr it = subMesh->GetElements() ) {
- while(it->more()) {
- const SMDS_MeshElement* e = it->next();
- if ( e->GetType() != elemType || !e->IsQuadratic() ) {
- myCreateQuadratic = false;
- break;
- }
- else {
- // fill NLinkNodeMap
- switch ( e->NbNodes() ) {
- case 3:
- AddNLinkNode(e->GetNode(0),e->GetNode(1),e->GetNode(2)); break;
- case 6:
- AddNLinkNode(e->GetNode(0),e->GetNode(1),e->GetNode(3));
- AddNLinkNode(e->GetNode(1),e->GetNode(2),e->GetNode(4));
- AddNLinkNode(e->GetNode(2),e->GetNode(0),e->GetNode(5)); break;
- case 8:
- AddNLinkNode(e->GetNode(0),e->GetNode(1),e->GetNode(4));
- AddNLinkNode(e->GetNode(1),e->GetNode(2),e->GetNode(5));
- AddNLinkNode(e->GetNode(2),e->GetNode(3),e->GetNode(6));
- AddNLinkNode(e->GetNode(3),e->GetNode(0),e->GetNode(7));
- break;
- default:
- myCreateQuadratic = false;
- break;
- }
- }
- }
- }
- }
- }
-
- if(!myCreateQuadratic) {
- myNLinkNodeMap.clear();
- }
- else {
- // treatment of periodic faces
- if ( aSh.ShapeType() == TopAbs_FACE ) {
- const TopoDS_Face& face = TopoDS::Face( aSh );
- BRepAdaptor_Surface surface( face );
- if ( surface.IsUPeriodic() || surface.IsVPeriodic() ) {
- // look for a seam edge
- for ( exp.Init( face, TopAbs_EDGE ); exp.More(); exp.Next()) {
- const TopoDS_Edge& edge = TopoDS::Edge( exp.Current() );
- if ( BRep_Tool::IsClosed( edge, face )) {
- // initialize myPar1, myPar2 and myParIndex
- if ( mySeamShapeIds.empty() ) {
- gp_Pnt2d uv1, uv2;
- BRep_Tool::UVPoints( edge, face, uv1, uv2 );
- if ( Abs( uv1.Coord(1) - uv2.Coord(1) ) < Abs( uv1.Coord(2) - uv2.Coord(2) ))
- {
- myParIndex = 1; // U periodic
- myPar1 = surface.FirstUParameter();
- myPar2 = surface.LastUParameter();
- }
- else {
- myParIndex = 2; // V periodic
- myPar1 = surface.FirstVParameter();
- myPar2 = surface.LastVParameter();
- }
- }
- // store shapes indices
- mySeamShapeIds.insert( meshDS->ShapeToIndex( exp.Current() ));
- for ( TopExp_Explorer v( exp.Current(), TopAbs_VERTEX ); v.More(); v.Next() )
- mySeamShapeIds.insert( meshDS->ShapeToIndex( v.Current() ));
- }
- }
- }
- }
- }
- return myCreateQuadratic;
-}
-
-
-//=======================================================================
-//function : IsMedium
-//purpose :
-//=======================================================================
-
-bool StdMeshers_Helper::IsMedium(const SMDS_MeshNode* node,
- const SMDSAbs_ElementType typeToCheck)
-{
- return SMESH_MeshEditor::IsMedium( node, typeToCheck );
-}
-
-//=======================================================================
-//function : AddNLinkNode
-//purpose :
-//=======================================================================
-/*!
- * Auxilary function for filling myNLinkNodeMap
- */
-void StdMeshers_Helper::AddNLinkNode(const SMDS_MeshNode* n1,
- const SMDS_MeshNode* n2,
- const SMDS_MeshNode* n12)
-{
- NLink link( n1, n2 );
- if ( n1 > n2 ) link = NLink( n2, n1 );
- // add new record to map
- myNLinkNodeMap.insert( make_pair(link,n12));
-}
-
-//=======================================================================
-/*!
- * \brief Select UV on either of 2 pcurves of a seam edge, closest to the given UV
- * \param uv1 - UV on the seam
- * \param uv2 - UV within a face
- * \retval gp_Pnt2d - selected UV
- */
-//=======================================================================
-
-gp_Pnt2d StdMeshers_Helper::GetUVOnSeam( const gp_Pnt2d& uv1, const gp_Pnt2d& uv2 ) const
-{
- double p1 = uv1.Coord( myParIndex );
- double p2 = uv2.Coord( myParIndex );
- double p3 = ( Abs( p1 - myPar1 ) < Abs( p1 - myPar2 )) ? myPar2 : myPar1;
- if ( Abs( p2 - p1 ) > Abs( p2 - p3 ))
- p1 = p3;
- gp_Pnt2d result = uv1;
- result.SetCoord( myParIndex, p1 );
- return result;
-}
-
-//=======================================================================
-/*!
- * \brief Return node UV on face
- * \param F - the face
- * \param n - the node
- * \param n2 - a medium node will be placed between n and n2
- * \retval gp_XY - resulting UV
- *
- * Auxilary function called form GetMediumNode()
- */
-//=======================================================================
-
-gp_XY StdMeshers_Helper::GetNodeUV(const TopoDS_Face& F,
- const SMDS_MeshNode* n,
- const SMDS_MeshNode* n2)
-{
- gp_Pnt2d uv;
- const SMDS_PositionPtr Pos = n->GetPosition();
- if(Pos->GetTypeOfPosition()==SMDS_TOP_FACE) {
- // node has position on face
- const SMDS_FacePosition* fpos =
- static_cast<const SMDS_FacePosition*>(n->GetPosition().get());
- uv = gp_Pnt2d(fpos->GetUParameter(),fpos->GetVParameter());
- }
- else if(Pos->GetTypeOfPosition()==SMDS_TOP_EDGE) {
- // node has position on edge => it is needed to find
- // corresponding edge from face, get pcurve for this
- // edge and recieve value from this pcurve
- const SMDS_EdgePosition* epos =
- static_cast<const SMDS_EdgePosition*>(n->GetPosition().get());
- SMESHDS_Mesh* meshDS = GetMesh()->GetMeshDS();
- int edgeID = Pos->GetShapeId();
- TopoDS_Edge E = TopoDS::Edge(meshDS->IndexToShape(edgeID));
- double f, l;
- TopLoc_Location loc;
- Handle(Geom2d_Curve) C2d = BRep_Tool::CurveOnSurface(E, F, f, l);
- uv = C2d->Value( epos->GetUParameter() );
- // for a node on a seam edge select one of UVs on 2 pcurves
- if ( n2 && mySeamShapeIds.find( edgeID ) != mySeamShapeIds.end() )
- uv = GetUVOnSeam( uv, GetNodeUV( F, n2, 0 ));
- }
- else if(Pos->GetTypeOfPosition()==SMDS_TOP_VERTEX) {
- SMESHDS_Mesh * meshDS = GetMesh()->GetMeshDS();
- int vertexID = n->GetPosition()->GetShapeId();
- const TopoDS_Vertex& V = TopoDS::Vertex(meshDS->IndexToShape(vertexID));
- uv = BRep_Tool::Parameters( V, F );
- if ( n2 && mySeamShapeIds.find( vertexID ) != mySeamShapeIds.end() )
- uv = GetUVOnSeam( uv, GetNodeUV( F, n2, 0 ));
- }
- return uv.XY();
-}
-
-
-//=======================================================================
-//function : GetMediumNode
-//purpose :
-//=======================================================================
-/*!
- * Special function for search or creation medium node
- */
-const SMDS_MeshNode* StdMeshers_Helper::GetMediumNode(const SMDS_MeshNode* n1,
- const SMDS_MeshNode* n2,
- const bool force3d)
-{
-//cout<<"n1: "<<n1;
-//cout<<"n2: "<<n2;
- NLink link(( n1 < n2 ? n1 : n2 ), ( n1 < n2 ? n2 : n1 ));
- ItNLinkNode itLN = myNLinkNodeMap.find( link );
- if ( itLN != myNLinkNodeMap.end() ) {
- return (*itLN).second;
- }
- else {
- // create medium node
- SMDS_MeshNode* n12;
- SMESHDS_Mesh* meshDS = GetMesh()->GetMeshDS();
- if(!force3d) {
- // we try to create medium node using UV parameters of
- // nodes, else - medium between corresponding 3d points
- const SMDS_PositionPtr Pos1 = n1->GetPosition();
- const SMDS_PositionPtr Pos2 = n2->GetPosition();
- int faceID = -1;
- if( Pos1->GetTypeOfPosition()==SMDS_TOP_FACE ) {
- faceID = Pos1->GetShapeId();
- }
- else if( Pos2->GetTypeOfPosition()==SMDS_TOP_FACE ) {
- faceID = Pos2->GetShapeId();
- }
- if(faceID>-1) {
- TopoDS_Face F = TopoDS::Face(meshDS->IndexToShape(faceID));
- gp_XY p1 = GetNodeUV(F,n1,n2);
- gp_XY p2 = GetNodeUV(F,n2,n1);
- double u = (p1.X()+p2.X())/2.;
- double v = (p1.Y()+p2.Y())/2.;
- Handle(Geom_Surface) S = BRep_Tool::Surface(F);
- gp_Pnt P = S->Value(u, v);
- n12 = meshDS->AddNode(P.X(), P.Y(), P.Z());
- meshDS->SetNodeOnFace(n12, faceID, u, v);
- myNLinkNodeMap.insert(NLinkNodeMap::value_type(link,n12));
- return n12;
- }
- }
- // 3d variant
- double x = ( n1->X() + n2->X() )/2.;
- double y = ( n1->Y() + n2->Y() )/2.;
- double z = ( n1->Z() + n2->Z() )/2.;
- n12 = meshDS->AddNode(x,y,z);
- meshDS->SetNodeInVolume(n12, myShapeID);
- myNLinkNodeMap.insert(NLinkNodeMap::value_type(link,n12));
- return n12;
- }
-}
-
-
-//=======================================================================
-//function : AddFace
-//purpose :
-//=======================================================================
-/*!
- * Special function for creation quadratic triangle
- */
-SMDS_MeshFace* StdMeshers_Helper::AddFace(const SMDS_MeshNode* n1,
- const SMDS_MeshNode* n2,
- const SMDS_MeshNode* n3)
-{
- SMESHDS_Mesh * meshDS = GetMesh()->GetMeshDS();
- if(!myCreateQuadratic) {
- return meshDS->AddFace(n1, n2, n3);
- }
-
- const SMDS_MeshNode* n12 = GetMediumNode(n1,n2,false);
- const SMDS_MeshNode* n23 = GetMediumNode(n2,n3,false);
- const SMDS_MeshNode* n31 = GetMediumNode(n3,n1,false);
-
- return meshDS->AddFace(n1, n2, n3, n12, n23, n31);
-}
-
-
-//=======================================================================
-//function : AddFace
-//purpose :
-//=======================================================================
-/*!
- * Special function for creation quadratic quadrangle
- */
-SMDS_MeshFace* StdMeshers_Helper::AddFace(const SMDS_MeshNode* n1,
- const SMDS_MeshNode* n2,
- const SMDS_MeshNode* n3,
- const SMDS_MeshNode* n4)
-{
- SMESHDS_Mesh * meshDS = GetMesh()->GetMeshDS();
- if(!myCreateQuadratic) {
- return meshDS->AddFace(n1, n2, n3, n4);
- }
-
- const SMDS_MeshNode* n12 = GetMediumNode(n1,n2,false);
- const SMDS_MeshNode* n23 = GetMediumNode(n2,n3,false);
- const SMDS_MeshNode* n34 = GetMediumNode(n3,n4,false);
- const SMDS_MeshNode* n41 = GetMediumNode(n4,n1,false);
-
- return meshDS->AddFace(n1, n2, n3, n4, n12, n23, n34, n41);
-}
-
-
-//=======================================================================
-//function : AddVolume
-//purpose :
-//=======================================================================
-/*!
- * Special function for creation quadratic volume
- */
-SMDS_MeshVolume* StdMeshers_Helper::AddVolume(const SMDS_MeshNode* n1,
- const SMDS_MeshNode* n2,
- const SMDS_MeshNode* n3,
- const SMDS_MeshNode* n4,
- const SMDS_MeshNode* n5,
- const SMDS_MeshNode* n6)
-{
- SMESHDS_Mesh * meshDS = GetMesh()->GetMeshDS();
- if(!myCreateQuadratic) {
- return meshDS->AddVolume(n1, n2, n3, n4, n5, n6);
- }
-
- const SMDS_MeshNode* n12 = GetMediumNode(n1,n2,true);
- const SMDS_MeshNode* n23 = GetMediumNode(n2,n3,true);
- const SMDS_MeshNode* n31 = GetMediumNode(n3,n1,true);
-
- const SMDS_MeshNode* n45 = GetMediumNode(n4,n5,true);
- const SMDS_MeshNode* n56 = GetMediumNode(n5,n6,true);
- const SMDS_MeshNode* n64 = GetMediumNode(n6,n4,true);
-
- const SMDS_MeshNode* n14 = GetMediumNode(n1,n4,true);
- const SMDS_MeshNode* n25 = GetMediumNode(n2,n5,true);
- const SMDS_MeshNode* n36 = GetMediumNode(n3,n6,true);
-
- return meshDS->AddVolume(n1, n2, n3, n4, n5, n6,
- n12, n23, n31, n45, n56, n64, n14, n25, n36);
-}
-
-
-//=======================================================================
-//function : AddVolume
-//purpose :
-//=======================================================================
-/*!
- * Special function for creation quadratic volume
- */
-SMDS_MeshVolume* StdMeshers_Helper::AddVolume(const SMDS_MeshNode* n1,
- const SMDS_MeshNode* n2,
- const SMDS_MeshNode* n3,
- const SMDS_MeshNode* n4)
-{
- SMESHDS_Mesh * meshDS = GetMesh()->GetMeshDS();
- if(!myCreateQuadratic) {
- return meshDS->AddVolume(n1, n2, n3, n4);
- }
-
- const SMDS_MeshNode* n12 = GetMediumNode(n1,n2,true);
- const SMDS_MeshNode* n23 = GetMediumNode(n2,n3,true);
- const SMDS_MeshNode* n31 = GetMediumNode(n3,n1,true);
-
- const SMDS_MeshNode* n14 = GetMediumNode(n1,n4,true);
- const SMDS_MeshNode* n24 = GetMediumNode(n2,n4,true);
- const SMDS_MeshNode* n34 = GetMediumNode(n3,n4,true);
-
- return meshDS->AddVolume(n1, n2, n3, n4, n12, n23, n31, n14, n24, n34);
-}
-
-
-//=======================================================================
-//function : AddVolume
-//purpose :
-//=======================================================================
-/*!
- * Special function for creation quadratic volume
- */
-SMDS_MeshVolume* StdMeshers_Helper::AddVolume(const SMDS_MeshNode* n1,
- const SMDS_MeshNode* n2,
- const SMDS_MeshNode* n3,
- const SMDS_MeshNode* n4,
- const SMDS_MeshNode* n5,
- const SMDS_MeshNode* n6,
- const SMDS_MeshNode* n7,
- const SMDS_MeshNode* n8)
-{
- SMESHDS_Mesh * meshDS = GetMesh()->GetMeshDS();
- if(!myCreateQuadratic) {
- return meshDS->AddVolume(n1, n2, n3, n4, n5, n6, n7, n8);
- }
-
- const SMDS_MeshNode* n12 = GetMediumNode(n1,n2,true);
- const SMDS_MeshNode* n23 = GetMediumNode(n2,n3,true);
- const SMDS_MeshNode* n34 = GetMediumNode(n3,n4,true);
- const SMDS_MeshNode* n41 = GetMediumNode(n4,n1,true);
-
- const SMDS_MeshNode* n56 = GetMediumNode(n5,n6,true);
- const SMDS_MeshNode* n67 = GetMediumNode(n6,n7,true);
- const SMDS_MeshNode* n78 = GetMediumNode(n7,n8,true);
- const SMDS_MeshNode* n85 = GetMediumNode(n8,n5,true);
-
- const SMDS_MeshNode* n15 = GetMediumNode(n1,n5,true);
- const SMDS_MeshNode* n26 = GetMediumNode(n2,n6,true);
- const SMDS_MeshNode* n37 = GetMediumNode(n3,n7,true);
- const SMDS_MeshNode* n48 = GetMediumNode(n4,n8,true);
-
- return meshDS->AddVolume(n1, n2, n3, n4, n5, n6, n7, n8,
- n12, n23, n34, n41, n56, n67,
- n78, n85, n15, n26, n37, n48);
-}
-
-
+++ /dev/null
-// File: StdMeshers_Helper.hxx
-// Created: 15.02.06 14:48:09
-// Author: Sergey KUUL
-// Copyright: Open CASCADE 2006
-
-
-#ifndef StdMeshers_Helper_HeaderFile
-#define StdMeshers_Helper_HeaderFile
-
-#include <SMESH_Mesh.hxx>
-#include <TopoDS_Shape.hxx>
-#include <SMDS_MeshNode.hxx>
-#include <TopoDS_Face.hxx>
-#include <gp_Pnt2d.hxx>
-
-#include <map>
-
-typedef pair<const SMDS_MeshNode*, const SMDS_MeshNode*> NLink;
-typedef map<NLink, const SMDS_MeshNode*> NLinkNodeMap;
-typedef map<NLink, const SMDS_MeshNode*>::iterator ItNLinkNode;
-
-/*!
- * \brief It helps meshers to add elements
- *
- * It allow meshers not to care about creation of medium nodes
- * when filling a quadratic mesh. Helper does it itself.
- * It defines degree of elements to create when IsQuadraticSubMesh()
- * is called.
- */
-
-class StdMeshers_Helper
-{
- public:
- // ---------- PUBLIC METHODS ----------
-
- /// Empty constructor
- StdMeshers_Helper(SMESH_Mesh& theMesh)
- { myMesh=(void *)&theMesh; myCreateQuadratic = false; }
-
- SMESH_Mesh* GetMesh() const
- { return (SMESH_Mesh*)myMesh; }
-
- /// Copy constructor
- //Standard_EXPORT StdMeshers_Helper (const StdMeshers_Helper& theOther);
-
- /// Destructor
- //Standard_EXPORT virtual ~StdMeshers_Helper ();
-
- /**
- * Check submesh for given shape
- * Check if all elements on this shape
- * are quadratic, if yes => set true to myCreateQuadratic
- * (default value is false). Also fill myNLinkNodeMap
- * Returns myCreateQuadratic
- */
- bool IsQuadraticSubMesh(const TopoDS_Shape& theShape);
-
- /*!
- * \brief Returns true if given node is medium
- * \param n - node to check
- * \param typeToCheck - type of elements containing the node to ask about node status
- * \retval bool - check result
- */
- static bool IsMedium(const SMDS_MeshNode* node,
- const SMDSAbs_ElementType typeToCheck = SMDSAbs_All);
-
- /**
- * Auxilary function for filling myNLinkNodeMap
- */
- void AddNLinkNode(const SMDS_MeshNode* n1,
- const SMDS_MeshNode* n2,
- const SMDS_MeshNode* n12);
-
- /**
- * Auxilary function for filling myNLinkNodeMap
- */
- void AddNLinkNodeMap(const NLinkNodeMap& aMap)
- { myNLinkNodeMap.insert(aMap.begin(), aMap.end()); }
-
- /**
- * Returns myNLinkNodeMap
- */
- const NLinkNodeMap& GetNLinkNodeMap() { return myNLinkNodeMap; }
-
- /*!
- * \brief Return node UV on face
- * \param F - the face
- * \param n - the node
- * \param n2 - a medium node will be placed between n and n2
- * \retval gp_XY - resulting UV
- *
- * Auxilary function called form GetMediumNode()
- */
- gp_XY GetNodeUV(const TopoDS_Face& F,
- const SMDS_MeshNode* n,
- const SMDS_MeshNode* n2=0);
-
- /**
- * Special function for search or creation medium node
- */
- const SMDS_MeshNode* GetMediumNode(const SMDS_MeshNode* n1,
- const SMDS_MeshNode* n2,
- const bool force3d);
- /**
- * Special function for creation quadratic triangle
- */
- SMDS_MeshFace* AddFace(const SMDS_MeshNode* n1,
- const SMDS_MeshNode* n2,
- const SMDS_MeshNode* n3);
-
- /**
- * Special function for creation quadratic quadrangle
- */
- SMDS_MeshFace* AddFace(const SMDS_MeshNode* n1,
- const SMDS_MeshNode* n2,
- const SMDS_MeshNode* n3,
- const SMDS_MeshNode* n4);
-
- /**
- * Special function for creation quadratic tetraahedron
- */
- SMDS_MeshVolume* AddVolume(const SMDS_MeshNode* n1,
- const SMDS_MeshNode* n2,
- const SMDS_MeshNode* n3,
- const SMDS_MeshNode* n4);
-
- /**
- * Special function for creation quadratic pentahedron
- */
- SMDS_MeshVolume* AddVolume(const SMDS_MeshNode* n1,
- const SMDS_MeshNode* n2,
- const SMDS_MeshNode* n3,
- const SMDS_MeshNode* n4,
- const SMDS_MeshNode* n5,
- const SMDS_MeshNode* n6);
-
- /**
- * Special function for creation quadratic hexahedron
- */
- SMDS_MeshVolume* AddVolume(const SMDS_MeshNode* n1,
- const SMDS_MeshNode* n2,
- const SMDS_MeshNode* n3,
- const SMDS_MeshNode* n4,
- const SMDS_MeshNode* n5,
- const SMDS_MeshNode* n6,
- const SMDS_MeshNode* n7,
- const SMDS_MeshNode* n8);
-
-
- protected:
-
- /*!
- * \brief Select UV on either of 2 pcurves of a seam edge, closest to the given UV
- * \param uv1 - UV on the seam
- * \param uv2 - UV within a face
- * \retval gp_Pnt2d - selected UV
- */
- gp_Pnt2d GetUVOnSeam( const gp_Pnt2d& uv1, const gp_Pnt2d& uv2 ) const;
-
- private:
-
- void* myMesh;
-
- int myShapeID;
-
- // Key for creation quadratic faces
- bool myCreateQuadratic;
-
- // special map for using during creation quadratic faces
- NLinkNodeMap myNLinkNodeMap;
-
- std::set< int > mySeamShapeIds;
- double myPar1, myPar2; // bounds of a closed periodic surface
- int myParIndex; // bounds' index (1-U, 2-V)
-
-};
-
-
-#endif
