src/GEOMAlgo/GEOMAlgo_Gluer2_1.cxx

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