StdMeshers_FaceSide.hxx

Go to the documentation of this file.

// Copyright (C) 2007-2026  CEA, EDF, OPEN CASCADE
//
// Copyright (C) 2003-2007  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, or (at your option) any later version.
//
// 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
//
 
// File      : StdMeshers_FaceSide.hxx
// Created   : Wed Jan 31 18:41:25 2007
// Author    : Edward AGAPOV (eap)
// Module    : SMESH
//
#ifndef StdMeshers_FaceSide_HeaderFile
#define StdMeshers_FaceSide_HeaderFile
 
#include "SMESH_StdMeshers.hxx"
 
#include "SMESH_ProxyMesh.hxx"
 
#include <Geom2d_Curve.hxx>
#include <GeomAdaptor_Curve.hxx>
#include <TopoDS_Edge.hxx>
#include <TopoDS_Face.hxx>
#include <TopoDS_Vertex.hxx>
#include <gp_Pnt2d.hxx>
 
#include <vector>
#include <list>
#include <boost/shared_ptr.hpp>
 
class Adaptor2d_Curve2d;
class Adaptor3d_Curve;
class BRepAdaptor_CompCurve;
class SMDS_MeshNode;
class SMESH_Mesh;
class SMESH_MesherHelper;
class StdMeshers_FaceSide;
struct SMESH_ComputeError;
 
typedef boost::shared_ptr< SMESH_ComputeError >  TError;
typedef boost::shared_ptr< StdMeshers_FaceSide > StdMeshers_FaceSidePtr;
typedef std::vector< StdMeshers_FaceSidePtr >    TSideVector;
 
//================================================================================
//================================================================================
 
class STDMESHERS_EXPORT StdMeshers_FaceSide
{
public:
 
  enum { ALL_EDGES = -1, LAST_EDGE = -1 }; 
 
  StdMeshers_FaceSide(const TopoDS_Face&   theFace,
                      const TopoDS_Edge&   theEdge,
                      SMESH_Mesh*          theMesh,
                      const bool           theIsForward,
                      const bool           theIgnoreMediumNodes,
                      SMESH_MesherHelper*  theFaceHelper = NULL,
                      SMESH_ProxyMesh::Ptr theProxyMesh = SMESH_ProxyMesh::Ptr());
  StdMeshers_FaceSide(const TopoDS_Face&            theFace,
                      const std::list<TopoDS_Edge>& theEdges,
                      SMESH_Mesh*                   theMesh,
                      const bool                    theIsForward,
                      const bool                    theIgnoreMediumNodes,
                      SMESH_MesherHelper*           theFaceHelper = NULL,
                      SMESH_ProxyMesh::Ptr          theProxyMesh = SMESH_ProxyMesh::Ptr());
  StdMeshers_FaceSide(const StdMeshers_FaceSide*  theSide,
                      const SMDS_MeshNode*        theNode,
                      const gp_Pnt2d*             thePnt2d1,
                      const gp_Pnt2d*             thePnt2d2 = NULL,
                      const Handle(Geom2d_Curve)& theC2d = NULL,
                      const double                theUFirst = 0.,
                      const double                theULast = 1.);
  StdMeshers_FaceSide(UVPtStructVec&     theSideNodes,
                      const TopoDS_Face& theFace = TopoDS_Face(),
                      const TopoDS_Edge& theEdge = TopoDS_Edge(),
                      SMESH_Mesh*        theMesh = 0);
 
  ~StdMeshers_FaceSide();
 
  // static "constructors"
  static StdMeshers_FaceSidePtr New(const TopoDS_Face&   Face,
                                    const TopoDS_Edge&   Edge,
                                    SMESH_Mesh*          Mesh,
                                    const bool           IsForward,
                                    const bool           IgnoreMediumNodes,
                                    SMESH_MesherHelper*  FaceHelper = NULL,
                                    SMESH_ProxyMesh::Ptr ProxyMesh = SMESH_ProxyMesh::Ptr())
  { return StdMeshers_FaceSidePtr
      ( new StdMeshers_FaceSide( Face,Edge,Mesh,IsForward,IgnoreMediumNodes,FaceHelper,ProxyMesh ));
  }
  static StdMeshers_FaceSidePtr New (const TopoDS_Face&            Face,
                                     const std::list<TopoDS_Edge>& Edges,
                                     SMESH_Mesh*                   Mesh,
                                     const bool                    IsForward,
                                     const bool                    IgnoreMediumNodes,
                                     SMESH_MesherHelper*           FaceHelper = NULL,
                                     SMESH_ProxyMesh::Ptr          ProxyMesh = SMESH_ProxyMesh::Ptr())
  { return StdMeshers_FaceSidePtr
      ( new StdMeshers_FaceSide( Face,Edges,Mesh,IsForward,IgnoreMediumNodes,FaceHelper,ProxyMesh ));
  }
  static StdMeshers_FaceSidePtr New (const StdMeshers_FaceSide*  Side,
                                     const SMDS_MeshNode*        Node,
                                     const gp_Pnt2d*             Pnt2d1,
                                     const gp_Pnt2d*             Pnt2d2 = NULL,
                                     const Handle(Geom2d_Curve)& C2d = NULL,
                                     const double                UFirst = 0.,
                                     const double                ULast = 1.)
  { return StdMeshers_FaceSidePtr
      ( new StdMeshers_FaceSide( Side,Node,Pnt2d1,Pnt2d2,C2d,UFirst,ULast ));
  }
  static StdMeshers_FaceSidePtr New (UVPtStructVec&     theSideNodes,
                                     const TopoDS_Face& theFace = TopoDS_Face(),
                                     const TopoDS_Edge& theEdge = TopoDS_Edge(),
                                     SMESH_Mesh*        theMesh = 0)
  {
    return StdMeshers_FaceSidePtr( new StdMeshers_FaceSide( theSideNodes, theFace, theEdge, theMesh ));
  }
 
  static TSideVector GetFaceWires(const TopoDS_Face&   theFace,
                                  SMESH_Mesh &         theMesh,
                                  const bool           theIgnoreMediumNodes,
                                  TError &             theError,
                                  SMESH_MesherHelper*  theFaceHelper = NULL,
                                  SMESH_ProxyMesh::Ptr theProxyMesh = SMESH_ProxyMesh::Ptr(),
                                  const bool           theCheckVertexNodes = true);
  void Reverse();
  void SetIgnoreMediumNodes(bool toIgnore);
 
  smIdType NbPoints(const bool update = false) const;
  smIdType NbSegments(const bool update = false) const;
  SMESH_Mesh* GetMesh() const { return myProxyMesh->GetMesh(); }
  bool MissVertexNode() const { return myMissingVertexNodes; }
 
  const UVPtStructVec& GetUVPtStruct( bool isXConst = 0, double constValue = 0 ) const;
  const UVPtStructVec& SimulateUVPtStruct(int    nbSeg,
                                          bool   isXConst   = 0,
                                          double constValue = 0) const;
  std::vector<const SMDS_MeshNode*> GetOrderedNodes( int iE = ALL_EDGES ) const;
 
  bool GetEdgeNodes(const size_t                       i,
                    std::vector<const SMDS_MeshNode*>& nodes,
                    bool                               inlude1stVertex=true,
                    bool                               inludeLastVertex=true) const;
 
  const SMDS_MeshNode* VertexNode(std::size_t i, bool* isMoved = 0) const;
 
  /*
   * \brief Return edge and parameter on edge by normalized parameter
   */
  inline double Parameter(double U, TopoDS_Edge & edge) const;
  /*
   * \brief Return edge ID and parameter on edge by normalized parameter
   */
  inline double Parameter(double U, int & edgeID) const;
  gp_Pnt2d Value2d(double U) const;
  gp_Pnt   Value3d(double U) const;
  Adaptor2d_Curve2d* GetCurve2d() const;
  BRepAdaptor_CompCurve* GetCurve3d() const;
  int NbEdges() const { return (int)myEdge.size(); }
  const TopoDS_Edge& Edge(int i) const { return myEdge[i]; }
  const std::vector<TopoDS_Edge>& Edges() const { return myEdge; }
  const TopoDS_Face& Face() const { return myFace; }
  TopoDS_Vertex FirstVertex(int i=0) const;
  TopoDS_Vertex LastVertex(int i = LAST_EDGE) const;
  bool IsClosed() const;
  double Length() const { return myLength; }
  inline int EdgeIndex( double U ) const;
 
  void dump(const char* msg=0) const;
  
  inline int EdgeID(int i) const;
  inline Handle(Geom2d_Curve) Curve2d(int i) const;
  inline double FirstParameter(int i) const;
  inline double LastParameter(int i) const;
  inline double FirstU(int i) const;
  inline double LastU(int i) const;
  inline double EdgeLength(int i) const;
  inline bool IsReversed(int i) const;
  SMESH_MesherHelper* FaceHelper() const;
 
protected:
 
  void reverseProxySubmesh( const TopoDS_Edge& E );
 
  // DON't FORGET to update Reverse() when adding one more vector!
  TopoDS_Face                       myFace;
  std::vector<uvPtStruct>           myPoints, myFalsePoints;
  std::vector<TopoDS_Edge>          myEdge;
  std::vector<int>                  myEdgeID;
  std::vector<Handle(Geom2d_Curve)> myC2d;
  std::vector<GeomAdaptor_Curve>    myC3dAdaptor;
  std::vector<double>               myFirst, myLast;
  std::vector<double>               myNormPar;
  std::vector<double>               myEdgeLength;
  std::vector<int>                  myIsUniform;
  double                            myLength;
  smIdType                          myNbPonits, myNbSegments;
  SMESH_ProxyMesh::Ptr              myProxyMesh;
  bool                              myMissingVertexNodes, myIgnoreMediumNodes;
  gp_Pnt2d                          myDefaultPnt2d;
  SMESH_MesherHelper*               myHelper;
};
 
 
//================================================================================
//================================================================================
 
inline int StdMeshers_FaceSide::EdgeIndex( double U ) const
{
  int i = (int)myNormPar.size() - 1;
  while ( i > 0 && U < myNormPar[ i-1 ] ) --i;
  return i;
}
 
//================================================================================
//================================================================================
 
inline double StdMeshers_FaceSide::Parameter(double U, TopoDS_Edge & edge) const
{
  int i = EdgeIndex( U );
  edge = myEdge[ i ];
  double prevU = i ? myNormPar[ i-1 ] : 0;
  double r = ( U - prevU )/ ( myNormPar[ i ] - prevU );
  return myFirst[i] * ( 1 - r ) + myLast[i] * r;
}
 
//================================================================================
//================================================================================
 
inline double StdMeshers_FaceSide::Parameter(double U, int & edgeID) const
{
  int i = EdgeIndex( U );
  edgeID = myEdgeID[ i ];
  double prevU = i ? myNormPar[ i-1 ] : 0;
  double r = ( U - prevU )/ ( myNormPar[ i ] - prevU );
  return myFirst[i] * ( 1 - r ) + myLast[i] * r;
}
 
//================================================================================
//================================================================================
 
inline double StdMeshers_FaceSide::FirstParameter(int i) const
{
  return i==0 ? 0. : i<(int)myNormPar.size() ? myNormPar[i-1] : 1.;
}
 
//================================================================================
//================================================================================
 
inline double StdMeshers_FaceSide::LastParameter(int i) const
{
  return i < (int)myNormPar.size() ? myNormPar[i] : 1;
}
 
//================================================================================
//================================================================================
 
inline double StdMeshers_FaceSide::FirstU(int i) const
{
  return myFirst[ i % myFirst.size() ];
}
 
//================================================================================
//================================================================================
 
inline double StdMeshers_FaceSide::LastU(int i) const
{
  return myLast[ i % myLast.size() ];
}
 
//================================================================================
//================================================================================
 
inline int StdMeshers_FaceSide::EdgeID(int i) const
{
  return myEdgeID[ i % myEdgeID.size() ];
}
 
//================================================================================
//================================================================================
 
inline Handle(Geom2d_Curve) StdMeshers_FaceSide::Curve2d(int i) const
{
  return myC2d[ i % myC2d.size() ];
}
 
//================================================================================
 //================================================================================
 
inline double StdMeshers_FaceSide::EdgeLength(int i) const
{
  return myEdgeLength[ i % myEdgeLength.size() ];
}
 
//================================================================================
 //================================================================================
 
inline bool StdMeshers_FaceSide::IsReversed(int i) const
{
  return myFirst[i] > myLast[i];
}
 
#endif