src/GeomAlgoImpl/GEOMAlgo_Gluer2_1.cxx

   1 // Copyright (C) 2007-2022  CEA/DEN, EDF R&D, OPEN CASCADE
   2 //
   3 // Copyright (C) 2003-2007  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, or (at your option) any later version.
  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 // File:     GEOMAlgo_Gluer2_1.cxx
  24 // Author:   Peter KURNEV
  25
  26 #include <GEOMAlgo_Gluer2.hxx>
  27
  28 #include <gp_XYZ.hxx>
  29 #include <gp_Pnt.hxx>
  30
  31 #include <TopLoc_Location.hxx>
  32
  33 #include <Geom_Surface.hxx>
  34
  35 #include <TopAbs_ShapeEnum.hxx>
  36
  37 #include <TopoDS_Shape.hxx>
  38 #include <TopoDS_Vertex.hxx>
  39 #include <TopoDS_Edge.hxx>
  40 #include <TopoDS_Wire.hxx>
  41 #include <TopoDS_Face.hxx>
  42 #include <TopoDS_Iterator.hxx>
  43
  44 #include <BRep_Tool.hxx>
  45 #include <BRep_Builder.hxx>
  46
  47 #include <BRepTools.hxx>
  48 #include <TopExp.hxx>
  49
  50 #include <TopTools_MapOfShape.hxx>
  51 #include <TopTools_ListIteratorOfListOfShape.hxx>
  52 #include <TopTools_ListOfShape.hxx>
  53
  54 #include <BOPTools_AlgoTools.hxx>
  55 #include <GEOMAlgo_AlgoTools.hxx>
  56 #include <BOPTools_AlgoTools2D.hxx>
  57
  58
  59 //=======================================================================
  60 //function : MakeBRepShapes
  61 //purpose  :
  62 //=======================================================================
  63 void GEOMAlgo_Gluer2::MakeBRepShapes(const TopoDS_Shape& theS,
  64                                      TopoDS_Shape& theSnew)
  65 {
  66   TopAbs_ShapeEnum aType;
  67   //
  68   aType=theS.ShapeType();
  69   if (aType==TopAbs_EDGE) {
  70     TopoDS_Edge aEE, aEEnew;
  71     //
  72     aEE=*((TopoDS_Edge*)&theS);
  73     MakeEdge(aEE, aEEnew);
  74     if (myErrorStatus) {
  75       return;
  76     }
  77     //
  78     theSnew=aEEnew;
  79   }
  80   else if (aType==TopAbs_FACE) {
  81     TopoDS_Face aFF, aFFnew;
  82     //
  83     aFF=*((TopoDS_Face*)&theS);
  84     MakeFace(aFF, aFFnew);
  85     if (myErrorStatus) {
  86       return;
  87     }
  88     //
  89     theSnew=aFFnew;
  90   }
  91 }
  92
  93 //=======================================================================
  94 //function : MakeFace
  95 //purpose  :
  96 //=======================================================================
  97 void GEOMAlgo_Gluer2::MakeFace(const TopoDS_Face& theF,
  98                                TopoDS_Face& theFnew)
  99 {
 100   Standard_Boolean bIsToReverse, bIsUPeriodic;
 101   Standard_Integer iRet;
 102   Standard_Real aTol, aUMin, aUMax, aVMin, aVMax;
 103   Handle(Geom_Surface) aS;
 104   TopLoc_Location aLoc;
 105   TopoDS_Shape aW, aWr;
 106   TopoDS_Edge aEx;
 107   TopoDS_Face aFF, aFnew;
 108   TopoDS_Iterator aItW, aItE;
 109   BRep_Builder aBB;
 110   TopTools_ListOfShape aLEr;
 111   TopTools_ListIteratorOfListOfShape aItLE;
 112   //
 113   myErrorStatus=0;
 114   //
 115   aFF=theF;
 116   aFF.Orientation(TopAbs_FORWARD);
 117   //
 118   aTol=BRep_Tool::Tolerance(aFF);
 119   aS=BRep_Tool::Surface(aFF, aLoc);
 120   bIsUPeriodic=GEOMAlgo_AlgoTools::IsUPeriodic(aS);
 121   BRepTools::UVBounds(aFF, aUMin, aUMax, aVMin, aVMax);
 122   //
 123   aBB.MakeFace (aFnew, aS, aLoc, aTol);
 124   //
 125   aItW.Initialize(aFF);
 126   for (; aItW.More(); aItW.Next()) {
 127     const TopoDS_Shape& aW=aItW.Value();
 128     //
 129     if (!myOrigins.IsBound(aW)) {
 130       aBB.Add(aFnew, aW);
 131       continue;
 132     }
 133     //
 134     aWr=myOrigins.Find(aW);
 135     //
 136     // clear contents of Wr
 137     aLEr.Clear();
 138     aItE.Initialize(aWr);
 139     for (; aItE.More(); aItE.Next()) {
 140       const TopoDS_Shape& aEr=aItE.Value();
 141       aLEr.Append(aEr);
 142     }
 143     //
 144     aItLE.Initialize(aLEr);
 145     for (; aItLE.More(); aItLE.Next()) {
 146       const TopoDS_Shape& aEr=aItLE.Value();
 147       aBB.Remove(aWr, aEr);
 148     }
 149     //
 150     // refill contents of Wr
 151     aItE.Initialize(aW);
 152     for (; aItE.More(); aItE.Next()) {
 153       const TopoDS_Edge& aE=*((TopoDS_Edge*)(&aItE.Value()));
 154       //
 155       aEx=aE;
 156       if (myOrigins.IsBound(aE)) {
 157         aEx=*((TopoDS_Edge*)(&myOrigins.Find(aE)));
 158       }
 159       //
 160       if (!BRep_Tool::Degenerated(aEx)) {
 161         aEx.Orientation(TopAbs_FORWARD);
 162         TopoDS_Edge aE_forward = aE;
 163         aE_forward.Orientation(TopAbs_FORWARD);
 164         if (bIsUPeriodic) {
 165           GEOMAlgo_AlgoTools::RefinePCurveForEdgeOnFace(aEx, aFF, aUMin, aUMax);
 166         }
 167         //
 168         //modified by NIZNHY-PKV Fri Feb 03 11:18:17 2012f
 169         iRet=GEOMAlgo_AlgoTools::BuildPCurveForEdgeOnFace(aE_forward, aEx, aFF, myContext);
 170         if (iRet) {
 171           continue;
 172         }
 173         //modified by NIZNHY-PKV Fri Feb 03 11:18:20 2012t
 174         //
 175         bIsToReverse=GEOMAlgo_AlgoTools::IsSplitToReverse(aEx, aE_forward, myContext);
 176         //bIsToReverse=BOPTools_AlgoTools::IsSplitToReverse(aEx, aE, myContext);
 177
 178         aEx.Orientation(aE.Orientation());
 179
 180         if (bIsToReverse)
 181           aEx.Reverse();
 182       }
 183       else {
 184         aEx.Orientation(aE.Orientation());
 185       }
 186       aBB.Add(aWr, aEx);
 187     }// for (; aItE.More(); aItE.Next()) {
 188     //
 189     aBB.Add(aFnew, aWr);
 190   }// for (; aItW.More(); aItW.Next()) {
 191   theFnew=aFnew;
 192 }
 193
 194 //=======================================================================
 195 //function : MakeEdge
 196 //purpose  :
 197 //=======================================================================
 198 void GEOMAlgo_Gluer2::MakeEdge(const TopoDS_Edge& aE,
 199                                TopoDS_Edge& aNewEdge)
 200 {
 201   myErrorStatus=0;
 202   //
 203   Standard_Boolean bIsDE;
 204   Standard_Real aT1, aT2;
 205   TopoDS_Vertex aV1, aV2, aVR1, aVR2;
 206   TopoDS_Edge aEx;
 207   //
 208   bIsDE=BRep_Tool::Degenerated(aE);
 209   //
 210   aEx=aE;
 211   aEx.Orientation(TopAbs_FORWARD);
 212   //
 213   TopExp::Vertices(aEx, aV1, aV2);
 214   //
 215   aT1=BRep_Tool::Parameter(aV1, aEx);
 216   aT2=BRep_Tool::Parameter(aV2, aEx);
 217   //
 218   aVR1=aV1;
 219   if (myOrigins.IsBound(aV1)) {
 220     aVR1=*((TopoDS_Vertex*)&myOrigins.Find(aV1));
 221   }
 222   aVR1.Orientation(TopAbs_FORWARD);
 223   //
 224   aVR2=aV2;
 225   if (myOrigins.IsBound(aV2)) {
 226     aVR2=*((TopoDS_Vertex*)&myOrigins.Find(aV2));
 227   }
 228   aVR2.Orientation(TopAbs_REVERSED);
 229   //
 230   if (!bIsDE) {
 231     BOPTools_AlgoTools::MakeSplitEdge(aEx, aVR1, aT1, aVR2, aT2, aNewEdge);
 232   }
 233   else {
 234     Standard_Real aTol;
 235     BRep_Builder aBB;
 236     TopoDS_Edge E;
 237     //
 238     aTol=BRep_Tool::Tolerance(aE);
 239     //
 240     E=aEx;
 241     E.EmptyCopy();
 242     //
 243     aBB.Add  (E, aVR1);
 244     aBB.Add  (E, aVR2);
 245     aBB.Range(E, aT1, aT2);
 246     aBB.Degenerated(E, Standard_True);
 247     aBB.UpdateEdge(E, aTol);
 248     //
 249     aNewEdge=E;
 250   }
 251 }
 252
 253 //=======================================================================
 254 //function : MakeVertex
 255 //purpose  :
 256 //=======================================================================
 257 void GEOMAlgo_Gluer2::MakeVertex(const TopTools_ListOfShape& aLV,
 258                                  TopoDS_Vertex& aNewVertex)
 259 {
 260   Standard_Integer aNbV;
 261   Standard_Real aTolV, aD, aDmax;
 262   gp_XYZ aGC;
 263   gp_Pnt aP3D, aPGC;
 264   TopoDS_Vertex aVx;
 265   BRep_Builder aBB;
 266   TopTools_ListIteratorOfListOfShape aIt;
 267   //
 268   aNbV=aLV.Extent();
 269   if (!aNbV) {
 270     return;
 271   }
 272   //
 273   // center of gravity
 274   aGC.SetCoord(0.,0.,0.);
 275   aIt.Initialize(aLV);
 276   for (; aIt.More(); aIt.Next()) {
 277     aVx=*((TopoDS_Vertex*)(&aIt.Value()));
 278     aP3D=BRep_Tool::Pnt(aVx);
 279     aGC+=aP3D.XYZ();
 280   }
 281   aGC/=(Standard_Real)aNbV;
 282   aPGC.SetXYZ(aGC);
 283   //
 284   // tolerance value
 285   aDmax=-1.;
 286   aIt.Initialize(aLV);
 287   for (; aIt.More(); aIt.Next()) {
 288     aVx=*((TopoDS_Vertex*)(&aIt.Value()));
 289     aP3D=BRep_Tool::Pnt(aVx);
 290     aTolV=BRep_Tool::Tolerance(aVx);
 291     aD=aPGC.Distance(aP3D)+aTolV;
 292     if (aD>aDmax) {
 293       aDmax=aD;
 294     }
 295   }
 296   //
 297   aBB.MakeVertex (aNewVertex, aPGC, aDmax);
 298 }
 299
 300 //=======================================================================
 301 //function : MapBRepShapes
 302 //purpose  :
 303 //=======================================================================
 304 void GEOMAlgo_Gluer2::MapBRepShapes(const TopoDS_Shape& aS,
 305                                     TopTools_MapOfShape& aM)
 306 {
 307   TopAbs_ShapeEnum aType;
 308   TopoDS_Iterator aIt;
 309   //
 310   aType=aS.ShapeType();
 311   if (aType==TopAbs_VERTEX ||
 312       aType==TopAbs_EDGE ||
 313       aType==TopAbs_FACE) {
 314     aM.Add(aS);
 315   }
 316   //
 317   aIt.Initialize(aS);
 318   for(; aIt.More(); aIt.Next()) {
 319     const TopoDS_Shape& aSx=aIt.Value();
 320     aType=aSx.ShapeType();
 321     GEOMAlgo_Gluer2::MapBRepShapes(aSx, aM);
 322   }
 323 }
 324
 325 //
 326 // ErrorStatus
 327 //
 328 // 40   - GEOMAlgo_GlueDetector has failed