src/GEOMAlgo/GEOMAlgo_Gluer2.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.cxx
  24 // Author:      Peter KURNEV
  25
  26 #include <GEOMAlgo_Gluer2.hxx>
  27
  28 #include <TopAbs_ShapeEnum.hxx>
  29
  30 #include <TopoDS_Compound.hxx>
  31 #include <TopoDS_Iterator.hxx>
  32 #include <TopoDS_Shape.hxx>
  33 #include <TopoDS_Vertex.hxx>
  34
  35 #include <BRep_Builder.hxx>
  36 #include <TopExp.hxx>
  37 #include <BRepLib.hxx>
  38
  39 #include <TopTools_MapOfShape.hxx>
  40 #include <TopTools_MapIteratorOfMapOfShape.hxx>
  41 #include <TopTools_DataMapOfShapeShape.hxx>
  42 #include <TopTools_ListOfShape.hxx>
  43 #include <TopTools_ListIteratorOfListOfShape.hxx>
  44 #include <TopTools_DataMapIteratorOfDataMapOfShapeListOfShape.hxx>
  45 #include <TopTools_IndexedMapOfShape.hxx>
  46
  47 #include <BOPTools_AlgoTools.hxx>
  48
  49 #include <GEOMAlgo_GlueDetector.hxx>
  50 #include <GEOMAlgo_AlgoTools.hxx>
  51
  52 //=======================================================================
  53 //function : GEOMAlgo_Gluer2
  54 //purpose  :
  55 //=======================================================================
  56 GEOMAlgo_Gluer2::GEOMAlgo_Gluer2()
  57 :
  58   GEOMAlgo_GluerAlgo(),
  59   GEOMAlgo_BuilderShape()
  60 {
  61   myTolerance=0.0001;
  62 }
  63 //=======================================================================
  64 //function : ~GEOMAlgo_Gluer2
  65 //purpose  :
  66 //=======================================================================
  67 GEOMAlgo_Gluer2::~GEOMAlgo_Gluer2()
  68 {
  69 }
  70 //=======================================================================
  71 //function : Clear
  72 //purpose  :
  73 //=======================================================================
  74 void GEOMAlgo_Gluer2::Clear()
  75 {
  76   myErrorStatus=0;
  77   myWarningStatus=0;
  78   //
  79   GEOMAlgo_GluerAlgo::Clear();
  80   //
  81   myImagesDetected.Clear();
  82   myOriginsDetected.Clear();
  83   myShapesToGlue.Clear();
  84   myImagesToWork.Clear();
  85   myOriginsToWork.Clear();
  86   myKeepNonSolids=Standard_False;
  87   myDetector.Clear();
  88 }
  89 //=======================================================================
  90 //function : StickedShapes
  91 //purpose  :
  92 //=======================================================================
  93 const TopTools_IndexedDataMapOfShapeListOfShape&
  94   GEOMAlgo_Gluer2::StickedShapes()
  95 {
  96   return myDetector.StickedShapes();
  97 }
  98 //=======================================================================
  99 //function : SetShapesToGlue
 100 //purpose  :
 101 //=======================================================================
 102 void GEOMAlgo_Gluer2::SetShapesToGlue
 103   (const TopTools_DataMapOfShapeListOfShape& aM)
 104 {
 105   myShapesToGlue=aM;
 106 }
 107 //=======================================================================
 108 //function : ShapesToGlue
 109 //purpose  :
 110 //=======================================================================
 111 const TopTools_DataMapOfShapeListOfShape& GEOMAlgo_Gluer2::ShapesToGlue()const
 112 {
 113   return myShapesToGlue;
 114 }
 115 //=======================================================================
 116 //function : SetKeepNonSolids
 117 //purpose  :
 118 //=======================================================================
 119 void GEOMAlgo_Gluer2::SetKeepNonSolids(const Standard_Boolean aFlag)
 120 {
 121   myKeepNonSolids=aFlag;
 122 }
 123 //=======================================================================
 124 //function : KeepNonSolids
 125 //purpose  :
 126 //=======================================================================
 127 Standard_Boolean GEOMAlgo_Gluer2::KeepNonSolids()const
 128 {
 129   return myKeepNonSolids;
 130 }
 131 //=======================================================================
 132 //function : ShapesDetected
 133 //purpose  :
 134 //=======================================================================
 135 const TopTools_DataMapOfShapeListOfShape& GEOMAlgo_Gluer2::ShapesDetected()const
 136 {
 137   return myImagesDetected;
 138 }
 139 //=======================================================================
 140 //function : ImagesToWork
 141 //purpose  :
 142 //=======================================================================
 143 const TopTools_DataMapOfShapeListOfShape& GEOMAlgo_Gluer2::ImagesToWork()const
 144 {
 145   return myImagesToWork;
 146 }
 147 //=======================================================================
 148 //function : Perform
 149 //purpose  :
 150 //=======================================================================
 151 void GEOMAlgo_Gluer2::Perform()
 152 {
 153   myErrorStatus=0;
 154   myWarningStatus=0;
 155   //
 156   CheckData();
 157   if (myErrorStatus) {
 158     return;
 159   }
 160   //
 161   // Initialize the context
 162   GEOMAlgo_GluerAlgo::Perform();
 163   //
 164   PerformShapesToWork();
 165   if (myErrorStatus) {
 166     return;
 167   }
 168   if (myWarningStatus==1) {
 169     // no shapes to glue
 170     myShape=myArgument;
 171     return;
 172   }
 173   //
 174   FillVertices();
 175   if (myErrorStatus) {
 176     return;
 177   }
 178   //
 179   FillEdges();
 180   if (myErrorStatus) {
 181     return;
 182   }
 183   //
 184   FillWires();
 185   if (myErrorStatus) {
 186     return;
 187   }
 188   //
 189   FillFaces();
 190   if (myErrorStatus) {
 191     return;
 192   }
 193   //
 194   FillShells();
 195   if (myErrorStatus) {
 196     return;
 197   }
 198   //
 199   FillSolids();
 200   if (myErrorStatus) {
 201     return;
 202   }
 203   //
 204   FillCompSolids();
 205   if (myErrorStatus) {
 206     return;
 207   }
 208   //
 209   FillCompounds();
 210   if (myErrorStatus) {
 211     return;
 212   }
 213   //
 214   BuildResult();
 215   if (myErrorStatus) {
 216     return;
 217   }
 218   //
 219   PrepareHistory();
 220   if (myErrorStatus) {
 221     return;
 222   }
 223   //
 224   BRepLib::SameParameter(myShape, myTolerance, Standard_True);
 225 }
 226 //=======================================================================
 227 //function : CheckData
 228 //purpose  :
 229 //=======================================================================
 230 void GEOMAlgo_Gluer2::CheckData()
 231 {
 232   Standard_Integer aNbSG, i;
 233   TopAbs_ShapeEnum aType, aTypeX;
 234   TopTools_ListIteratorOfListOfShape aItLS;
 235   TopTools_DataMapIteratorOfDataMapOfShapeListOfShape aItDMSLS;
 236   //
 237   myErrorStatus=0;
 238   myWarningStatus=0;
 239   //
 240   aNbSG=myShapesToGlue.Extent();
 241   aType=TopAbs_SHAPE;
 242   if (aNbSG) {
 243     // Check myShapesToGlue
 244     aItDMSLS.Initialize(myShapesToGlue);
 245     for (; aItDMSLS.More(); aItDMSLS.Next()) {
 246       //const TopoDS_Shape& aSkey=aItDMSLS.Key();
 247       const TopTools_ListOfShape& aLSG=aItDMSLS.Value();
 248       aItLS.Initialize(aLSG);
 249       for (i=0; aItLS.More(); aItLS.Next(), ++i) {
 250         const TopoDS_Shape& aSG=aItLS.Value();
 251         aTypeX=aSG.ShapeType();
 252         if (!i) {
 253           aType=aTypeX;
 254           if (!(aType==TopAbs_VERTEX ||
 255                 aType==TopAbs_EDGE ||
 256                 aType==TopAbs_FACE)) {
 257             myErrorStatus=21;// non-brep shapes
 258             return;
 259           }
 260           continue;
 261         }
 262         if (aTypeX!=aType) {
 263           myErrorStatus=20;// non-homogeneous shapes
 264           return;
 265         }
 266       }
 267     }
 268   }// if (aNbSG) {
 269 }
 270 //=======================================================================
 271 //function : FillEdges
 272 //purpose  :
 273 //=======================================================================
 274 void GEOMAlgo_Gluer2::FillEdges()
 275 {
 276   FillBRepShapes(TopAbs_EDGE);
 277 }
 278 //=======================================================================
 279 //function : FillFaces
 280 //purpose  :
 281 //=======================================================================
 282 void GEOMAlgo_Gluer2::FillFaces()
 283 {
 284   FillBRepShapes(TopAbs_FACE);
 285 }
 286 //=======================================================================
 287 //function : FillWires
 288 //purpose  :
 289 //=======================================================================
 290 void GEOMAlgo_Gluer2::FillWires()
 291 {
 292   FillContainers(TopAbs_WIRE);
 293 }
 294 //=======================================================================
 295 //function : FillShells
 296 //purpose  :
 297 //=======================================================================
 298 void GEOMAlgo_Gluer2::FillShells()
 299 {
 300   FillContainers(TopAbs_SHELL);
 301 }
 302 //=======================================================================
 303 //function : FillSolids
 304 //purpose  :
 305 //=======================================================================
 306 void GEOMAlgo_Gluer2::FillSolids()
 307 {
 308   FillContainers(TopAbs_SOLID);
 309 }
 310 //=======================================================================
 311 //function : FillCompSolids
 312 //purpose  :
 313 //=======================================================================
 314 void GEOMAlgo_Gluer2::FillCompSolids()
 315 {
 316   FillContainers(TopAbs_COMPSOLID);
 317 }
 318 //=======================================================================
 319 //function : FillVertices
 320 //purpose  :
 321 //=======================================================================
 322 void GEOMAlgo_Gluer2::FillVertices()
 323 {
 324   TopAbs_ShapeEnum aType;
 325   TopoDS_Vertex aVnew;
 326   TopTools_ListIteratorOfListOfShape aItLS;
 327   TopTools_DataMapIteratorOfDataMapOfShapeListOfShape aItDMSLS;
 328   //
 329   myErrorStatus=0;
 330   myWarningStatus=0;
 331   //
 332   aItDMSLS.Initialize(myImagesToWork);
 333   for (; aItDMSLS.More(); aItDMSLS.Next()) {
 334     const TopoDS_Shape& aSkey=aItDMSLS.Key();
 335     aType=aSkey.ShapeType();
 336     if (aType!=TopAbs_VERTEX) {
 337       continue;
 338     }
 339     //
 340     const TopTools_ListOfShape& aLSD=aItDMSLS.Value();
 341     //
 342     GEOMAlgo_Gluer2::MakeVertex(aLSD, aVnew);
 343     //
 344     myImages.Bind(aVnew, aLSD);
 345     //
 346     aItLS.Initialize(aLSD);
 347     for (; aItLS.More(); aItLS.Next()) {
 348       const TopoDS_Shape& aV=aItLS.Value();
 349       myOrigins.Bind(aV, aVnew);
 350     }
 351   }
 352 }
 353 //=======================================================================
 354 //function : FillBRepShapes
 355 //purpose  :
 356 //=======================================================================
 357 void GEOMAlgo_Gluer2::FillBRepShapes(const TopAbs_ShapeEnum theType)
 358 {
 359   Standard_Boolean bHasImage, bIsToWork;
 360   Standard_Integer i, aNbE;
 361   TopoDS_Iterator aItS;
 362   TopoDS_Shape aEnew;
 363   TopTools_IndexedMapOfShape aME;
 364   TopTools_MapOfShape aMFence;
 365   TopTools_ListIteratorOfListOfShape aItLS;
 366   //
 367   myErrorStatus=0;
 368   myWarningStatus=0;
 369   //
 370   TopExp::MapShapes(myArgument, theType, aME);
 371   //
 372   aNbE=aME.Extent();
 373   for (i=1; i<=aNbE; ++i) {
 374     const TopoDS_Shape& aE=aME(i);
 375     //
 376     if (!aMFence.Add(aE)) {
 377       continue;
 378     }
 379     //
 380     bIsToWork=myOriginsToWork.IsBound(aE);
 381     bHasImage=HasImage(aE);
 382     if (!bHasImage && !bIsToWork) {
 383       continue;
 384     }
 385     //
 386     MakeBRepShapes(aE, aEnew);
 387     //
 388     //myImages / myOrigins
 389     if (bIsToWork) {
 390       const TopoDS_Shape& aSkey=myOriginsToWork.Find(aE);
 391       const TopTools_ListOfShape& aLSD=myImagesToWork.Find(aSkey);
 392       //
 393       myImages.Bind(aEnew, aLSD);
 394       //
 395       aItLS.Initialize(aLSD);
 396       for (; aItLS.More(); aItLS.Next()) {
 397         const TopoDS_Shape& aEx=aItLS.Value();
 398         myOrigins.Bind(aEx, aEnew);
 399         //
 400         aMFence.Add(aEx);
 401       }
 402     }
 403     else {
 404       TopTools_ListOfShape aLSD;
 405       //
 406       aLSD.Append(aE);
 407       myImages.Bind(aEnew, aLSD);
 408       myOrigins.Bind(aE, aEnew);
 409     }
 410   }//for (i=1; i<=aNbF; ++i) {
 411 }
 412 //=======================================================================
 413 //function : FillContainers
 414 //purpose  :
 415 //=======================================================================
 416 void GEOMAlgo_Gluer2::FillContainers(const TopAbs_ShapeEnum aType)
 417 {
 418   Standard_Boolean bHasImage, bToReverse;
 419   Standard_Integer i, aNbW;
 420   TopoDS_Shape aWnew, aEnew;
 421   TopoDS_Iterator aItS;
 422   BRep_Builder aBB;
 423   TopTools_IndexedMapOfShape aMW;
 424   TopTools_MapOfShape aMFence;
 425   //
 426   myErrorStatus=0;
 427   myWarningStatus=0;
 428   //
 429   TopExp::MapShapes(myArgument, aType, aMW);
 430   //
 431   aNbW=aMW.Extent();
 432   for (i=1; i<=aNbW; ++i) {
 433     const TopoDS_Shape& aW=aMW(i);
 434     //
 435     if (!aMFence.Add(aW)) {
 436       continue;
 437     }
 438     //
 439     bHasImage=HasImage(aW);
 440     if (!bHasImage) {
 441       continue;
 442     }
 443     //
 444     GEOMAlgo_AlgoTools::MakeContainer(aType, aWnew);
 445     aWnew.Orientation(aW.Orientation());
 446     //
 447     aItS.Initialize(aW);
 448     for (; aItS.More(); aItS.Next()) {
 449       const TopoDS_Shape& aE=aItS.Value();
 450       if (myOrigins.IsBound(aE)) {
 451         aEnew=myOrigins.Find(aE);
 452         //
 453         bToReverse=BOPTools_AlgoTools::IsSplitToReverse(aEnew, aE, myContext);
 454         if (bToReverse) {
 455           aEnew.Reverse();
 456         }
 457         //
 458         aBB.Add(aWnew, aEnew);
 459       }
 460       else {
 461         aBB.Add(aWnew, aE);
 462       }
 463     }
 464     //
 465     //myImages / myOrigins
 466     TopTools_ListOfShape aLSD;
 467     //
 468     aLSD.Append(aW);
 469     myImages.Bind(aWnew, aLSD);
 470     myOrigins.Bind(aW, aWnew);
 471     //
 472   }//for (i=1; i<=aNbE; ++i) {
 473 }
 474 //=======================================================================
 475 //function : FillCompounds
 476 //purpose  :
 477 //=======================================================================
 478 void GEOMAlgo_Gluer2::FillCompounds()
 479 {
 480   TopAbs_ShapeEnum aType;
 481   TopoDS_Iterator aItC;
 482   //
 483   myErrorStatus=0;
 484   myWarningStatus=0;
 485   //
 486   aItC.Initialize(myArgument);
 487   for (; aItC.More(); aItC.Next()) {
 488     const TopoDS_Shape& aCx=aItC.Value();
 489     aType=aCx.ShapeType();
 490     if (aType==TopAbs_COMPOUND) {
 491       FillCompound(aCx);
 492     }
 493   }
 494 }
 495 //=======================================================================
 496 //function : FillCompound
 497 //purpose  :
 498 //=======================================================================
 499 void GEOMAlgo_Gluer2::FillCompound(const TopoDS_Shape& aC)
 500 {
 501   Standard_Boolean bHasImage;
 502   TopAbs_ShapeEnum aType;
 503   TopoDS_Shape aCnew, aCXnew;
 504   TopoDS_Iterator aItC;
 505   BRep_Builder aBB;
 506   //
 507   bHasImage=HasImage(aC);
 508   if (!bHasImage) {
 509     return;
 510   }
 511   //
 512   GEOMAlgo_AlgoTools::MakeContainer(TopAbs_COMPOUND, aCnew);
 513   //
 514   aItC.Initialize(aC);
 515   for (; aItC.More(); aItC.Next()) {
 516     const TopoDS_Shape& aCX=aItC.Value();
 517     aType=aCX.ShapeType();
 518     //
 519     if (aType==TopAbs_COMPOUND) {
 520       FillCompound(aCX);
 521     }
 522     //
 523     if (myOrigins.IsBound(aCX)) {
 524       aCXnew=myOrigins.Find(aCX);
 525       aCXnew.Orientation(aCX.Orientation());
 526       aBB.Add(aCnew, aCXnew);
 527     }
 528     else {
 529       aBB.Add(aCnew, aCX);
 530     }
 531   }
 532   //
 533   //myImages / myOrigins
 534   TopTools_ListOfShape aLSD;
 535   //
 536   aLSD.Append(aC);
 537   myImages.Bind(aCnew, aLSD);
 538   myOrigins.Bind(aC, aCnew);
 539 }
 540 //=======================================================================
 541 //function : HasImage
 542 //purpose  :
 543 //=======================================================================
 544 Standard_Boolean GEOMAlgo_Gluer2::HasImage(const TopoDS_Shape& aC)
 545 {
 546   Standard_Boolean bRet;
 547   TopAbs_ShapeEnum aType;
 548   TopoDS_Iterator aItC;
 549   //
 550   bRet=Standard_False;
 551   aItC.Initialize(aC);
 552   for (; aItC.More(); aItC.Next()) {
 553     const TopoDS_Shape& aCx=aItC.Value();
 554     aType=aCx.ShapeType();
 555     //
 556     if (aType==TopAbs_COMPOUND) {
 557       bRet=HasImage(aCx);
 558       if (bRet) {
 559         return bRet;
 560       }
 561     }
 562     else {
 563       bRet=myOrigins.IsBound(aCx);
 564       if (bRet) {
 565         return bRet;
 566       }
 567     }
 568   }
 569   //
 570   bRet=myOrigins.IsBound(aC);
 571   //
 572   return bRet;
 573 }
 574 //=======================================================================
 575 //function : BuildResult
 576 //purpose  :
 577 //=======================================================================
 578 void GEOMAlgo_Gluer2::BuildResult()
 579 {
 580   Standard_Boolean bHasImage;
 581   TopoDS_Shape aCnew, aCXnew;
 582   TopoDS_Iterator aItC;
 583   BRep_Builder aBB;
 584   //
 585   myErrorStatus=0;
 586   myWarningStatus=0;
 587   //
 588   bHasImage=Standard_False;
 589   aItC.Initialize(myArgument);
 590   for (; aItC.More(); aItC.Next()) {
 591     const TopoDS_Shape& aCx=aItC.Value();
 592     bHasImage=HasImage(aCx);
 593     if (bHasImage) {
 594       break;
 595     }
 596   }
 597   //
 598   if (!bHasImage) {
 599     myShape=myArgument;
 600     return;
 601   }
 602   //
 603   GEOMAlgo_AlgoTools::MakeContainer(TopAbs_COMPOUND, aCnew);
 604   //
 605   aItC.Initialize(myArgument);
 606   for (; aItC.More(); aItC.Next()) {
 607     const TopoDS_Shape& aCX=aItC.Value();
 608     if (myOrigins.IsBound(aCX)) {
 609       aCXnew=myOrigins.Find(aCX);
 610       aCXnew.Orientation(aCX.Orientation());
 611       aBB.Add(aCnew, aCXnew);
 612     }
 613     else {
 614       aBB.Add(aCnew, aCX);
 615     }
 616   }
 617   //
 618   if (!myKeepNonSolids) {
 619     Standard_Integer i, aNb;
 620     TopoDS_Shape aCnew1;
 621     TopTools_IndexedMapOfShape aM;
 622     //
 623     GEOMAlgo_AlgoTools::MakeContainer(TopAbs_COMPOUND, aCnew1);
 624     //
 625     TopExp::MapShapes(aCnew, TopAbs_SOLID, aM);
 626
 627     aNb=aM.Extent();
 628     for (i=1; i<=aNb; ++i) {
 629       const TopoDS_Shape& aS=aM(i);
 630       aBB.Add(aCnew1, aS);
 631     }
 632     aCnew=aCnew1;
 633   }
 634   //
 635   myShape=aCnew;
 636 }
 637 //--------------------------------------------------------
 638 //
 639 // ErrorStatus
 640 // 11   - GEOMAlgo_GlueDetector failed
 641 // 13   - PerformImagesToWork failed
 642 // 14   - PerformImagesToWork failed
 643 //
 644 // WarningStatus
 645 // 1   - no shapes to glue
 646 // 2  - sticked shapes are detected.
 647 //      The value of myTolerance is so large that
 648 //      subshapes of a shape becomes intefere
 649 //      (e.g. vertices of an edge).
 650 //      In the case
 651 //      the result is can not be obtained