src/StdMeshers/StdMeshers_FaceSide.hxx

   1 //  SMESH SMESH : implementaion of SMESH idl descriptions
   2 //
   3 //  Copyright (C) 2003  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
   4 //  CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
   5 //
   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.
  10 //
  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.
  15 //
  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
  19 //
  20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  21 //
  22 //
  23 //
  24 // File      : StdMeshers_FaceSide.hxx
  25 // Created   : Wed Jan 31 18:41:25 2007
  26 // Author    : Edward AGAPOV (eap)
  27 // Module    : SMESH
  28
  29 #ifndef StdMeshers_FaceSide_HeaderFile
  30 #define StdMeshers_FaceSide_HeaderFile
  31
  32 #include <gp_Pnt2d.hxx>
  33 #include <TopoDS_Edge.hxx>
  34 #include <TopoDS_Vertex.hxx>
  35 #include <Geom2d_Curve.hxx>
  36 #include <TopExp.hxx>
  37
  38 #include <vector>
  39 #include <list>
  40 #include <boost/shared_ptr.hpp>
  41
  42 class SMDS_MeshNode;
  43 class SMESH_Mesh;
  44 class Adaptor2d_Curve2d;
  45 class Adaptor3d_Curve;
  46 class BRepAdaptor_CompCurve;
  47 class TopoDS_Face;
  48
  49 typedef struct uvPtStruct
  50 {
  51   double param;
  52   //int    curvIndex;
  53   double normParam;
  54   double u; // original 2d parameter
  55   double v;
  56   double x; // 2d parameter, normalized [0,1]
  57   double y;
  58   const SMDS_MeshNode * node;
  59 } UVPtStruct;
  60
  61 class StdMeshers_FaceSide;
  62 typedef boost::shared_ptr< StdMeshers_FaceSide > StdMeshers_FaceSidePtr;
  63 typedef boost::shared_ptr< uvPtStruct > UVPtStructPtr;
  64
  65 //================================================================================
  66 /*!
  67  * \brief Represents a side of a quasi quadrilateral face.
  68  * It can be composed of several edges. Gives access to geometry and 1D mesh of a side.
  69  */
  70 //================================================================================
  71
  72 class StdMeshers_FaceSide
  73 {
  74 public:
  75   /*!
  76    * \brief Wrap one edge
  77    */
  78   StdMeshers_FaceSide(const TopoDS_Face& theFace,
  79                       const TopoDS_Edge& theEdge,
  80                       SMESH_Mesh*        theMesh,
  81                       const bool         theIsForward);
  82   /*!
  83    * \brief Wrap several edges. Edges must be properly ordered and oriented.
  84    */
  85   StdMeshers_FaceSide(const TopoDS_Face& theFace,
  86                       list<TopoDS_Edge>& theEdges,
  87                       SMESH_Mesh*        theMesh,
  88                       const bool         theIsForward);
  89   /*!
  90    * \brief Change orientation of side geometry
  91    */
  92   void Reverse();
  93   /*!
  94    * \brief Return nb nodes on edges including ones on vertices
  95    */
  96   int NbPoints() const { return myNbPonits; }
  97   /*!
  98    * \brief Return nb nodes on edges including ones on vertices
  99    */
 100   int NbSegments() const { return myNbSegments; }
 101   /*!
 102    * \brief Return mesh
 103    */
 104   SMESH_Mesh* GetMesh() const { return myMesh; }
 105   /*!
 106    * \brief Return true if there vertices without nodes
 107    */
 108   bool MissVertexNode() const { return myMissingVertexNodes; }
 109   /*!
 110    * \brief Return detailed data on nodes
 111     * \param isXConst - true if normalized parameter X is constant
 112     * \param constValue - constant parameter value
 113     *
 114     * Missing nodes are allowed only on internal vertices
 115    */
 116   const vector<UVPtStruct>& GetUVPtStruct(bool isXConst, double constValue) const;
 117   /*!
 118    * \brief Simulates detailed data on nodes
 119     * \param isXConst - true if normalized parameter X is constant
 120     * \param constValue - constant parameter value
 121    */
 122   const vector<UVPtStruct>& SimulateUVPtStruct(int    nbSeg,
 123                                                bool   isXConst,
 124                                                double constValue) const;
 125   /*!
 126    * \brief Return edge and parameter on edge by normalized parameter
 127    */
 128   inline double Parameter(double U, TopoDS_Edge & edge) const;
 129   /*!
 130    * \brief Return UV by normalized parameter
 131    */
 132   gp_Pnt2d Value2d(double U) const;
 133   /*!
 134    * \brief Creates a Adaptor2d_Curve2d to be used in SMESH_Block
 135    */
 136   Adaptor2d_Curve2d* GetCurve2d() const;
 137   /*!
 138    * \brief Creates a fully functional Adaptor_Curve
 139    */
 140   BRepAdaptor_CompCurve* GetCurve3d() const;
 141   /*!
 142    * \brief Return nb of wrapped edges
 143    */
 144   int NbEdges() const { return myEdge.size(); }
 145   /*!
 146    * \brief Return i-th wrapped edge (count starts from zero)
 147    */
 148   const TopoDS_Edge& Edge(int i) const { return myEdge[i]; }
 149   /*!
 150    * \brief Return 1st vertex of the i-the edge (count starts from zero)
 151    */
 152   inline TopoDS_Vertex FirstVertex(int i=0) const;
 153   /*!
 154    * \brief Return last vertex of the i-the edge (count starts from zero)
 155    */
 156   inline TopoDS_Vertex LastVertex(int i=-1) const;
 157   /*!
 158    * \brief Return first normalized parameter of the i-the edge (count starts from zero)
 159    */
 160   inline double FirstParameter(int i) const;
 161   /*!
 162    * \brief Return ast normalized parameter of the i-the edge (count starts from zero)
 163    */
 164   inline double LastParameter(int i) const;
 165   /*!
 166    * \brief Return side length
 167    */
 168   double Length() const { return myLength; }
 169   /*!
 170    * \brief Return edge index corresponding to normalized parameter
 171    */
 172   inline int EdgeIndex( double U ) const;
 173
 174   //virtual gp_Pnt Value(double U) const;
 175
 176   void dump(const char* msg=0) const;
 177
 178
 179 protected:
 180   vector<uvPtStruct>           myPoints, myFalsePoints;
 181   vector<TopoDS_Edge>          myEdge;
 182   vector<Handle(Geom2d_Curve)> myC2d;
 183   vector<double>               myFirst;
 184   vector<double>               myLast;
 185   vector<double>               myNormPar;
 186   double                       myLength;
 187   int                          myNbPonits, myNbSegments;
 188   SMESH_Mesh*                  myMesh;
 189   bool                         myMissingVertexNodes;
 190 };
 191
 192
 193 //================================================================================
 194 /*!
 195  * \brief Return edge index corresponding to normalized parameter
 196   * \param U - the parameter
 197   * \retval int - index
 198  */
 199 //================================================================================
 200
 201 inline int StdMeshers_FaceSide::EdgeIndex( double U ) const
 202 {
 203   int i = myNormPar.size() - 1;
 204   while ( i > 0 && U < myNormPar[ i-1 ] ) --i;
 205   return i;
 206 }
 207
 208 //================================================================================
 209 /*!
 210  * \brief Return edge and parameter on edge by normalized parameter
 211   * \param U - the parameter
 212   * \retval double - pameter on a curve
 213  */
 214 //================================================================================
 215
 216 inline double StdMeshers_FaceSide::Parameter(double U, TopoDS_Edge & edge) const
 217 {
 218   int i = EdgeIndex( U );
 219   edge = myEdge[ i ];
 220   double prevU = i ? myNormPar[ i-1 ] : 0;
 221   double r = ( U - prevU )/ ( myNormPar[ i ] - prevU );
 222   return myFirst[i] * ( 1 - r ) + myLast[i] * r;
 223 }
 224
 225 //================================================================================
 226 /*!
 227  * \brief Return 1st vertex of the i-the edge
 228  */
 229 //================================================================================
 230
 231 inline TopoDS_Vertex StdMeshers_FaceSide::FirstVertex(int i) const
 232 {
 233   return i < myEdge.size() ? TopExp::FirstVertex( myEdge[i], 1 ) : TopoDS_Vertex();
 234 }
 235
 236 //================================================================================
 237 /*!
 238  * \brief Return last vertex of the i-the edge
 239  */
 240 //================================================================================
 241
 242 inline TopoDS_Vertex StdMeshers_FaceSide::LastVertex(int i) const
 243 {
 244   return i<0 ? TopExp::LastVertex( myEdge.back(), 1) : i<myEdge.size() ? TopExp::LastVertex( myEdge[i], 1 ) : TopoDS_Vertex();
 245 }
 246
 247 //================================================================================
 248 /*!
 249  * \brief Return first normalized parameter of the i-the edge
 250  */
 251 //================================================================================
 252
 253 inline double StdMeshers_FaceSide::FirstParameter(int i) const
 254 {
 255   return i==0 ? 0. : i<myNormPar.size() ? myNormPar[i-1] : 1.;
 256 }
 257
 258 //================================================================================
 259 /*!
 260  * \brief Return ast normalized parameter of the i-the edge
 261  */
 262 //================================================================================
 263
 264 inline double StdMeshers_FaceSide::LastParameter(int i) const
 265 {
 266   return i<myNormPar.size() ? myNormPar[i] : 1;
 267 }
 268
 269 #endif