src/GEOMImpl/GEOMImpl_GlueDriver.cxx

   1 // Copyright (C) 2007-2012  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.
  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 #include <Standard_Stream.hxx>
  24
  25 #include <GEOMImpl_GlueDriver.hxx>
  26 #include <GEOMImpl_IGlue.hxx>
  27 #include <GEOMImpl_Types.hxx>
  28
  29 #include <GEOM_Object.hxx>
  30 #include <GEOM_Function.hxx>
  31
  32 #include "GEOMAlgo_Gluer2.hxx"
  33 #include "GEOMAlgo_ListIteratorOfListOfCoupleOfShapes.hxx"
  34 #include "GEOMAlgo_CoupleOfShapes.hxx"
  35 #include "GEOMAlgo_ListOfCoupleOfShapes.hxx"
  36
  37 #include <Basics_OCCTVersion.hxx>
  38
  39 #include "utilities.h"
  40
  41 #include <TDataStd_IntegerArray.hxx>
  42
  43 #include <TopExp.hxx>
  44 #include <TopoDS_Shape.hxx>
  45 #include <TopTools_DataMapIteratorOfDataMapOfShapeListOfShape.hxx>
  46 #include <TopTools_IndexedMapOfShape.hxx>
  47 #include <TopTools_ListOfShape.hxx>
  48 #include <TopTools_ListIteratorOfListOfShape.hxx>
  49
  50 #include <ShapeFix_Shape.hxx>
  51
  52 #include <Standard_NullObject.hxx>
  53 #include <Standard_Failure.hxx>
  54
  55 #define MSG_BAD_TOLERANCE "Tolerance is too big"
  56 #define MSG_BAD_ARG_SHAPE "Argument shape is not a compound of hexahedral solids"
  57
  58 //=======================================================================
  59 //function : GEOMImpl_GlueDriver
  60 //purpose  :
  61 //=======================================================================
  62 GEOMImpl_GlueDriver::GEOMImpl_GlueDriver()
  63 {
  64 }
  65
  66 //=======================================================================
  67 //function : GetID
  68 //purpose  :
  69 //=======================================================================
  70 const Standard_GUID& GEOMImpl_GlueDriver::GetID()
  71 {
  72   static Standard_GUID aGlueDriver("FF1BBB63-5D14-4df2-980B-3A668264EA16");
  73   return aGlueDriver;
  74 }
  75
  76 //=======================================================================
  77 //function : GlueFacesWithWarnings
  78 //purpose  :
  79 //=======================================================================
  80 /*
  81 TopoDS_Shape GEOMImpl_GlueDriver::GlueFacesWithWarnings (const TopoDS_Shape& theShape,
  82                                                          const Standard_Real theTolerance,
  83                                                          const Standard_Boolean doKeepNonSolids,
  84                                                          TCollection_AsciiString& theWarning) const
  85 {
  86   Standard_Integer iErr, iWrn;
  87   TopoDS_Shape aRes;
  88   GEOMAlgo_Gluer aGluer;
  89
  90   aGluer.SetShape(theShape);
  91   aGluer.SetTolerance(theTolerance);
  92   aGluer.SetCheckGeometry(Standard_True);
  93   aGluer.SetKeepNonSolids(doKeepNonSolids);
  94
  95   aGluer.Perform();
  96
  97   iErr = aGluer.ErrorStatus();
  98   if (iErr) {
  99     switch (iErr) {
 100     case 2:
 101       Standard_Failure::Raise("No vertices found in source shape");
 102       break;
 103     case 3:
 104     case 4:
 105       Standard_Failure::Raise(MSG_BAD_TOLERANCE " or " MSG_BAD_ARG_SHAPE);
 106       break;
 107     case 5:
 108       Standard_Failure::Raise("Source shape is Null");
 109       break;
 110     case 6:
 111       Standard_Failure::Raise("Result shape is Null");
 112       break;
 113     case 100:
 114       Standard_Failure::Raise(MSG_BAD_TOLERANCE);
 115       break;
 116     case 101:
 117     case 102:
 118       Standard_Failure::Raise(MSG_BAD_ARG_SHAPE);
 119       break;
 120     case 200:
 121       Standard_Failure::Raise("Error occured during check of geometric coincidence");
 122       break;
 123     default:
 124       {
 125         // description of all errors see in GEOMAlgo_Gluer.cxx
 126         TCollection_AsciiString aMsg ("Error in GEOMAlgo_Gluer with code ");
 127         aMsg += TCollection_AsciiString(iErr);
 128         Standard_Failure::Raise(aMsg.ToCString());
 129         break;
 130       }
 131     }
 132     return aRes;
 133   }
 134
 135   iWrn = aGluer.WarningStatus();
 136   if (iWrn) {
 137     switch (iWrn) {
 138     case 1:
 139       {
 140         Standard_Integer nbAlone = aGluer.AloneShapes();
 141         theWarning = TCollection_AsciiString(nbAlone);
 142         theWarning += " solid(s) can not be glued by faces";
 143       }
 144       break;
 145     default:
 146       // description of all warnings see in GEOMAlgo_Gluer.cxx
 147       theWarning = "Warning in GEOMAlgo_Gluer with code ";
 148       theWarning += TCollection_AsciiString(iWrn);
 149       break;
 150     }
 151   }
 152
 153   aRes = aGluer.Result();
 154
 155   // SKL 18.01.2010 - patch for 20662
 156   Handle(ShapeFix_Shape) aSfs = new ShapeFix_Shape(aRes);
 157   aSfs->SetPrecision(Precision::Confusion());
 158   aSfs->Perform();
 159   aRes = aSfs->Shape();
 160
 161   // Fill history to be used by GetInPlace functionality
 162   TopTools_IndexedMapOfShape aResIndices;
 163   TopExp::MapShapes(aRes, aResIndices);
 164
 165   Handle(GEOM_Function) aFunction = GEOM_Function::GetFunction(Label());
 166
 167   // history for all argument shapes
 168   TDF_LabelSequence aLabelSeq;
 169   aFunction->GetDependency(aLabelSeq);
 170   Standard_Integer nbArg = aLabelSeq.Length();
 171
 172   for (Standard_Integer iarg = 1; iarg <= nbArg; iarg++) {
 173
 174     TDF_Label anArgumentRefLabel = aLabelSeq.Value(iarg);
 175
 176     Handle(GEOM_Object) anArgumentObject = GEOM_Object::GetReferencedObject(anArgumentRefLabel);
 177     TopoDS_Shape anArgumentShape = anArgumentObject->GetValue();
 178
 179     TopTools_IndexedMapOfShape anArgumentIndices;
 180     TopExp::MapShapes(anArgumentShape, anArgumentIndices);
 181     Standard_Integer nbArgumentEntities = anArgumentIndices.Extent();
 182
 183     // Find corresponding label in history
 184     TDF_Label anArgumentHistoryLabel =
 185       aFunction->GetArgumentHistoryEntry(anArgumentRefLabel, Standard_True);
 186
 187     for (Standard_Integer ie = 1; ie <= nbArgumentEntities; ie++) {
 188       TopoDS_Shape anEntity = anArgumentIndices.FindKey(ie);
 189       const TopTools_ListOfShape& aModified = aGluer.Modified(anEntity);
 190       Standard_Integer nbModified = aModified.Extent();
 191
 192       if (nbModified > 0) {
 193         TDF_Label aWhatHistoryLabel = anArgumentHistoryLabel.FindChild(ie, Standard_True);
 194         Handle(TDataStd_IntegerArray) anAttr =
 195           TDataStd_IntegerArray::Set(aWhatHistoryLabel, 1, nbModified);
 196
 197         TopTools_ListIteratorOfListOfShape itM (aModified);
 198         for (int im = 1; itM.More(); itM.Next(), ++im) {
 199           int id = aResIndices.FindIndex(itM.Value());
 200           anAttr->SetValue(im, id);
 201         }
 202       }
 203     }
 204   }
 205
 206   return aRes;
 207 }
 208
 209 //=======================================================================
 210 //function : GlueFaces
 211 //purpose  :
 212 //=======================================================================
 213 TopoDS_Shape GEOMImpl_GlueDriver::GlueFaces (const TopoDS_Shape& theShape,
 214                                              const Standard_Real theTolerance,
 215                                              const Standard_Boolean doKeepNonSolids)
 216 {
 217   TopoDS_Shape aRes;
 218   GEOMAlgo_Gluer aGluer;
 219
 220   aGluer.SetShape(theShape);
 221   aGluer.SetTolerance(theTolerance);
 222   aGluer.SetCheckGeometry(Standard_True);
 223   aGluer.SetKeepNonSolids(doKeepNonSolids);
 224
 225   aGluer.Perform();
 226
 227   Standard_Integer iErr = aGluer.ErrorStatus();
 228   if (iErr) {
 229     switch (iErr) {
 230     case 2:
 231       Standard_Failure::Raise("No vertices found in source shape");
 232       break;
 233     case 5:
 234       Standard_Failure::Raise("Source shape is Null");
 235       break;
 236     case 6:
 237       Standard_Failure::Raise("Result shape is Null");
 238       break;
 239     case 200:
 240       Standard_Failure::Raise("Error occured during check of geometric coincidence");
 241       break;
 242     default:
 243       {
 244         // description of all errors see in GEOMAlgo_Gluer.cxx
 245         TCollection_AsciiString aMsg ("Error in GEOMAlgo_Gluer with code ");
 246         aMsg += TCollection_AsciiString(iErr);
 247         Standard_Failure::Raise(aMsg.ToCString());
 248         break;
 249       }
 250     }
 251     return aRes;
 252   }
 253
 254   Standard_Integer iWrn = aGluer.WarningStatus();
 255   if (iWrn) {
 256     switch (iWrn) {
 257     case 1:
 258       MESSAGE("Some shapes can not be glued by faces");
 259       break;
 260     default:
 261       // description of all warnings see in GEOMAlgo_Gluer.cxx
 262       MESSAGE("Warning in GEOMAlgo_Gluer with code " << iWrn);
 263       break;
 264     }
 265   }
 266
 267   aRes = aGluer.Result();
 268
 269   return aRes;
 270 }
 271
 272 //=======================================================================
 273 //function : GlueFacesByList
 274 //purpose  :
 275 //=======================================================================
 276 TopoDS_Shape GEOMImpl_GlueDriver::GlueFacesByList (const TopoDS_Shape& theShape,
 277                                                    const Standard_Real theTolerance,
 278                                                    const Standard_Boolean doKeepNonSolids,
 279                                                    const TopTools_MapOfShape& aFaces)
 280 {
 281   TopoDS_Shape aRes;
 282
 283   GEOMAlgo_Gluer1 aGluer;
 284   GEOMAlgo_ListIteratorOfListOfCoupleOfShapes aItCS;
 285   GEOMAlgo_CoupleOfShapes aCS;
 286   GEOMAlgo_ListOfCoupleOfShapes aLCS;
 287
 288   aGluer.SetShape(theShape);
 289   aGluer.SetTolerance(theTolerance);
 290   aGluer.SetKeepNonSolids(doKeepNonSolids);
 291   aGluer.Perform();
 292   Standard_Integer iErr = aGluer.ErrorStatus();
 293   if (iErr) return aRes;
 294
 295   TopTools_ListOfShape listShape;
 296   const GEOMAlgo_ListOfCoupleOfShapes& aLCSG = aGluer.GluedFaces();
 297   // Access to faces
 298   aItCS.Initialize(aLCSG);
 299   for (; aItCS.More(); aItCS.Next()) {
 300     const GEOMAlgo_CoupleOfShapes& aCSG = aItCS.Value();
 301     const TopoDS_Shape& aF1 = aCSG.Shape1();
 302     const TopoDS_Shape& aF2 = aCSG.Shape2();
 303     if (aFaces.Contains(aF1) || aFaces.Contains(aF2))
 304       continue;
 305     aCS.SetShapes(aF1,aF2);
 306     aLCS.Append(aCS);
 307   }
 308
 309   //cout<<"aLCS.Extent() = "<<aLCS.Extent()<<endl;
 310   if (aLCS.Extent()>0) {
 311     aGluer.SetFacesToUnglue(aLCS);
 312     aGluer.UnglueFaces();
 313     iErr = aGluer.ErrorStatus();
 314     if (iErr) return aRes;
 315   }
 316
 317   aRes = aGluer.Result();
 318
 319   return aRes;
 320 }
 321 */
 322
 323 //=======================================================================
 324 //function : GlueFaces
 325 //purpose  :
 326 //=======================================================================
 327 TopoDS_Shape GEOMImpl_GlueDriver::GlueFaces (const TopoDS_Shape& theShape,
 328                                              const Standard_Real theTolerance,
 329                                              const Standard_Boolean doKeepNonSolids)
 330 {
 331   TopoDS_Shape aRes;
 332
 333   GEOMAlgo_Gluer2 aGA;
 334
 335   // 1. Initialization
 336   aGA.SetArgument(theShape);
 337   aGA.SetTolerance(theTolerance);
 338   aGA.SetKeepNonSolids(doKeepNonSolids);
 339
 340   // 2. Detect interferred shapes
 341   aGA.Detect();
 342   Standard_Integer iErr = aGA.ErrorStatus();
 343   if (iErr) {
 344     switch (iErr) {
 345     case 11:
 346       Standard_Failure::Raise("GEOMAlgo_GlueDetector failed");
 347       break;
 348     case 13:
 349     case 14:
 350       Standard_Failure::Raise("PerformImagesToWork failed");
 351       break;
 352     default:
 353       {
 354         // description of all errors see in GEOMAlgo_Gluer2.cxx
 355         TCollection_AsciiString aMsg ("Error in GEOMAlgo_Gluer2 with code ");
 356         aMsg += TCollection_AsciiString(iErr);
 357         Standard_Failure::Raise(aMsg.ToCString());
 358         break;
 359       }
 360     }
 361     return aRes;
 362   }
 363
 364   // 3. Set shapes to glue. If the operator is absent, the whole gluing will be done
 365   //aGA.SetShapesToGlue(aMSG);
 366
 367   // 4. Gluing
 368   aGA.Perform();
 369   iErr = aGA.ErrorStatus();
 370   if (iErr) {
 371     switch (iErr) {
 372     case 11:
 373       Standard_Failure::Raise("GEOMAlgo_GlueDetector failed");
 374       break;
 375     case 13:
 376     case 14:
 377       Standard_Failure::Raise("PerformImagesToWork failed");
 378       break;
 379     default:
 380       {
 381         // description of all errors see in GEOMAlgo_Gluer2.cxx
 382         TCollection_AsciiString aMsg ("Error in GEOMAlgo_Gluer2 with code ");
 383         aMsg += TCollection_AsciiString(iErr);
 384         Standard_Failure::Raise(aMsg.ToCString());
 385         break;
 386       }
 387     }
 388     return aRes;
 389   }
 390
 391   Standard_Integer iWrn = aGA.WarningStatus();
 392   if (iWrn) {
 393     switch (iWrn) {
 394     case 1:
 395       MESSAGE("No shapes to glue");
 396       break;
 397     default:
 398       // description of all warnings see in GEOMAlgo_Gluer2.cxx
 399       MESSAGE("Warning in GEOMAlgo_Gluer2 with code " << iWrn);
 400       break;
 401     }
 402   }
 403
 404   // 5. Result
 405   aRes = aGA.Shape();
 406
 407   return aRes;
 408 }
 409
 410 //=======================================================================
 411 //function : GlueWithWarnings
 412 //purpose  :
 413 //=======================================================================
 414 TopoDS_Shape GEOMImpl_GlueDriver::GlueWithWarnings (const TopoDS_Shape& theShape,
 415                                                     const Standard_Real theTolerance,
 416                                                     const TopAbs_ShapeEnum theShapeType,
 417                                                     const Standard_Boolean doKeepNonSolids,
 418                                                     TCollection_AsciiString& theWarning) const
 419 {
 420   TopoDS_Shape aRes;
 421
 422   GEOMAlgo_Gluer2 aGA;
 423
 424   // 1. Initialization
 425   aGA.SetArgument(theShape);
 426   aGA.SetTolerance(theTolerance);
 427   aGA.SetKeepNonSolids(doKeepNonSolids);
 428
 429   // 2. Detect interferred shapes
 430   aGA.Detect();
 431
 432   //modified by NIZNHY-PKV Tue Mar 13 14:07:12 2012f
 433 #if OCC_VERSION_LARGE > 0x06050200
 434   Standard_Integer iWrnDetect = aGA.WarningStatus();
 435   if (iWrnDetect == 2) {
 436     /*
 437     TopTools_ListIteratorOfListOfShape aItLS;
 438
 439     // Sticked shapes are detected
 440     const TopTools_IndexedDataMapOfShapeListOfShape& aIDMSS = pGluer2->StickedShapes();
 441
 442     Standard_Integer i, aNb = aIDMSS.Extent();
 443     for (i = 1; i <= aNb; ++i) {
 444       // ancestor aSa (edge, wire face,..)
 445       const TopoDS_Shape& aSa = aIDMSS.FindKey(i);
 446
 447       // successors aSs (vertex, edge, ...)
 448       // of the ancestor that are sticked
 449       // for given value of the tolerance
 450       const TopTools_ListOfShape& aLSS = aIDMSS.FindFromIndex(i);
 451       aItLS.Initialize(aLSS);
 452       for (; aItLS.More(); aItLS.Next()) {
 453         const TopoDS_Shape& aSs = aItLS.Value();
 454       }
 455     }
 456     */
 457     Standard_Failure::Raise("GLUE_ERROR_STICKED_SHAPES");
 458     //Standard_Failure::Raise("Sticked shapes are detected. The tolerance value is too big.");
 459     //theWarning = "Sticked shapes are detected. The tolerance value is too big.";
 460   }
 461 #endif
 462   //modified by NIZNHY-PKV Tue Mar 13 14:07:14 2012t
 463
 464   Standard_Integer iErr = aGA.ErrorStatus();
 465   if (iErr) {
 466     switch (iErr) {
 467     case 11:
 468       Standard_Failure::Raise("GEOMAlgo_GlueDetector failed");
 469       break;
 470     case 13:
 471     case 14:
 472       Standard_Failure::Raise("PerformImagesToWork failed");
 473       break;
 474     default:
 475       {
 476         // description of all errors see in GEOMAlgo_Gluer2.cxx
 477         TCollection_AsciiString aMsg ("Error in GEOMAlgo_Gluer2 with code ");
 478         aMsg += TCollection_AsciiString(iErr);
 479         Standard_Failure::Raise(aMsg.ToCString());
 480         break;
 481       }
 482     }
 483     return aRes;
 484   }
 485
 486   if (theShapeType != TopAbs_FACE) {
 487     // 3. Fill shapes to glue aMSG
 488     TopTools_DataMapOfShapeListOfShape aMSG;
 489     const TopTools_DataMapOfShapeListOfShape& aMSD = aGA.ShapesDetected();
 490     TopTools_DataMapIteratorOfDataMapOfShapeListOfShape aItMSD;
 491     aItMSD.Initialize(aMSD);
 492     for (; aItMSD.More(); aItMSD.Next()) {
 493       const TopoDS_Shape& aSx = aItMSD.Key();
 494       const TopTools_ListOfShape& aLSD = aItMSD.Value();
 495       if (aSx.ShapeType() == theShapeType) {
 496         aMSG.Bind(aSx, aLSD);
 497       }
 498     }
 499
 500     // 4. Set shapes to glue. If the operator is absent, the whole gluing will be done
 501     aGA.SetShapesToGlue(aMSG);
 502   }
 503
 504   // 5. Gluing
 505   aGA.Perform();
 506   iErr = aGA.ErrorStatus();
 507   if (iErr) {
 508     switch (iErr) {
 509     case 11:
 510       Standard_Failure::Raise("GEOMAlgo_GlueDetector failed");
 511       break;
 512     case 13:
 513     case 14:
 514       Standard_Failure::Raise("PerformImagesToWork failed");
 515       break;
 516     default:
 517       {
 518         // description of all errors see in GEOMAlgo_Gluer2.cxx
 519         TCollection_AsciiString aMsg ("Error in GEOMAlgo_Gluer2 with code ");
 520         aMsg += TCollection_AsciiString(iErr);
 521         Standard_Failure::Raise(aMsg.ToCString());
 522         break;
 523       }
 524     }
 525     return aRes;
 526   }
 527
 528   Standard_Integer iWrn = aGA.WarningStatus();
 529   if (iWrn) {
 530     switch (iWrn) {
 531     case 1:
 532       theWarning = "No shapes to glue";
 533       break;
 534     default:
 535       // description of all warnings see in GEOMAlgo_Gluer2.cxx
 536       theWarning = "Warning in GEOMAlgo_Gluer2 with code ";
 537       theWarning += TCollection_AsciiString(iWrn);
 538       break;
 539     }
 540   }
 541
 542   // 6. Result
 543   aRes = aGA.Shape();
 544
 545   // 7. Fill history to be used by GetInPlace functionality
 546   TopTools_IndexedMapOfShape aResIndices;
 547   TopExp::MapShapes(aRes, aResIndices);
 548
 549   Handle(GEOM_Function) aFunction = GEOM_Function::GetFunction(Label());
 550
 551   // history for all argument shapes
 552   TDF_LabelSequence aLabelSeq;
 553   aFunction->GetDependency(aLabelSeq);
 554   Standard_Integer nbArg = aLabelSeq.Length();
 555
 556   for (Standard_Integer iarg = 1; iarg <= nbArg; iarg++) {
 557
 558     TDF_Label anArgumentRefLabel = aLabelSeq.Value(iarg);
 559
 560     Handle(GEOM_Object) anArgumentObject = GEOM_Object::GetReferencedObject(anArgumentRefLabel);
 561     TopoDS_Shape anArgumentShape = anArgumentObject->GetValue();
 562
 563     TopTools_IndexedMapOfShape anArgumentIndices;
 564     TopExp::MapShapes(anArgumentShape, anArgumentIndices);
 565     Standard_Integer nbArgumentEntities = anArgumentIndices.Extent();
 566
 567     // Find corresponding label in history
 568     TDF_Label anArgumentHistoryLabel =
 569       aFunction->GetArgumentHistoryEntry(anArgumentRefLabel, Standard_True);
 570
 571     for (Standard_Integer ie = 1; ie <= nbArgumentEntities; ie++) {
 572       TopoDS_Shape anEntity = anArgumentIndices.FindKey(ie);
 573       const TopTools_ListOfShape& aModified = aGA.Modified(anEntity);
 574       Standard_Integer nbModified = aModified.Extent();
 575
 576       if (nbModified > 0) {
 577         TDF_Label aWhatHistoryLabel = anArgumentHistoryLabel.FindChild(ie, Standard_True);
 578         Handle(TDataStd_IntegerArray) anAttr =
 579           TDataStd_IntegerArray::Set(aWhatHistoryLabel, 1, nbModified);
 580
 581         TopTools_ListIteratorOfListOfShape itM (aModified);
 582         for (int im = 1; itM.More(); itM.Next(), ++im) {
 583           int id = aResIndices.FindIndex(itM.Value());
 584           anAttr->SetValue(im, id);
 585         }
 586       }
 587     }
 588   }
 589
 590   return aRes;
 591 }
 592
 593 //=======================================================================
 594 //function : GlueByList
 595 //purpose  :
 596 //=======================================================================
 597 TopoDS_Shape GEOMImpl_GlueDriver::GlueByList (const TopoDS_Shape& theShape,
 598                                               const Standard_Real theTolerance,
 599                                               const Standard_Boolean doKeepNonSolids,
 600                                               const TopTools_MapOfShape& theShapesList,
 601                                               const Standard_Boolean doGlueAllEdges)
 602 {
 603   TopoDS_Shape aRes;
 604
 605   GEOMAlgo_Gluer2 aGA;
 606
 607   // 1. Initialization
 608   aGA.SetArgument(theShape);
 609   aGA.SetTolerance(theTolerance);
 610   aGA.SetKeepNonSolids(doKeepNonSolids);
 611
 612   // 2. Detect interferred shapes
 613   aGA.Detect();
 614   Standard_Integer iErr = aGA.ErrorStatus();
 615   if (iErr) {
 616     switch (iErr) {
 617     case 11:
 618       Standard_Failure::Raise("GEOMAlgo_GlueDetector failed");
 619       break;
 620     case 13:
 621     case 14:
 622       Standard_Failure::Raise("PerformImagesToWork failed");
 623       break;
 624     default:
 625       {
 626         // description of all errors see in GEOMAlgo_Gluer2.cxx
 627         TCollection_AsciiString aMsg ("Error in GEOMAlgo_Gluer2 with code ");
 628         aMsg += TCollection_AsciiString(iErr);
 629         Standard_Failure::Raise(aMsg.ToCString());
 630         break;
 631       }
 632     }
 633     return aRes;
 634   }
 635
 636   // 3. Fill shapes to glue aMSG
 637   TopTools_DataMapOfShapeListOfShape aMSG;
 638   const TopTools_DataMapOfShapeListOfShape& aMSD = aGA.ShapesDetected();
 639   TopTools_DataMapIteratorOfDataMapOfShapeListOfShape aItMSD;
 640   aItMSD.Initialize(aMSD);
 641   for (; aItMSD.More(); aItMSD.Next()) {
 642     const TopoDS_Shape& aSx = aItMSD.Key();
 643     const TopTools_ListOfShape& aLSD = aItMSD.Value();
 644     TopTools_ListIteratorOfListOfShape anItLSD (aLSD);
 645     bool isToGlue = false;
 646     if (doGlueAllEdges && aSx.ShapeType() == TopAbs_EDGE) {
 647       isToGlue = true;
 648     }
 649     else {
 650       for (; anItLSD.More() && !isToGlue; anItLSD.Next()) {
 651         if (theShapesList.Contains(anItLSD.Value())) {
 652           isToGlue = true;
 653         }
 654       }
 655     }
 656     if (isToGlue) {
 657       aMSG.Bind(aSx, aLSD);
 658     }
 659   }
 660
 661   // 4. Set shapes to glue. If the operator is absent, the whole gluing will be done
 662   aGA.SetShapesToGlue(aMSG);
 663
 664   // 5. Gluing
 665   aGA.Perform();
 666   iErr = aGA.ErrorStatus();
 667   if (iErr) {
 668     switch (iErr) {
 669     case 11:
 670       Standard_Failure::Raise("GEOMAlgo_GlueDetector failed");
 671       break;
 672     case 13:
 673     case 14:
 674       Standard_Failure::Raise("PerformImagesToWork failed");
 675       break;
 676     default:
 677       {
 678         // description of all errors see in GEOMAlgo_Gluer2.cxx
 679         TCollection_AsciiString aMsg ("Error in GEOMAlgo_Gluer2 with code ");
 680         aMsg += TCollection_AsciiString(iErr);
 681         Standard_Failure::Raise(aMsg.ToCString());
 682         break;
 683       }
 684     }
 685     return aRes;
 686   }
 687
 688   Standard_Integer iWrn = aGA.WarningStatus();
 689   if (iWrn) {
 690     switch (iWrn) {
 691     case 1:
 692       MESSAGE("No shapes to glue");
 693       break;
 694     default:
 695       // description of all warnings see in GEOMAlgo_Gluer2.cxx
 696       MESSAGE("Warning in GEOMAlgo_Gluer2 with code " << iWrn);
 697       break;
 698     }
 699   }
 700
 701   // 6. Result
 702   aRes = aGA.Shape();
 703
 704   return aRes;
 705 }
 706
 707 //=======================================================================
 708 //function : Execute
 709 //purpose  :
 710 //=======================================================================
 711 Standard_Integer GEOMImpl_GlueDriver::Execute(TFunction_Logbook& log) const
 712 {
 713   if (Label().IsNull()) return 0;
 714   Handle(GEOM_Function) aFunction = GEOM_Function::GetFunction(Label());
 715
 716   GEOMImpl_IGlue aCI (aFunction);
 717   Standard_Integer aType = aFunction->GetType();
 718
 719   TopoDS_Shape aShape;
 720   TCollection_AsciiString aWrn;
 721
 722   Handle(GEOM_Function) aRefBase = aCI.GetBase();
 723   TopoDS_Shape aShapeBase = aRefBase->GetValue();
 724   if (aShapeBase.IsNull()) {
 725     Standard_NullObject::Raise("Shape for gluing is null");
 726   }
 727
 728   Standard_Real tol3d = aCI.GetTolerance();
 729
 730   Standard_Boolean aKeepNonSolids = aCI.GetKeepNonSolids();
 731
 732   if (aType == GLUE_FACES) {
 733     //aShape = GlueFacesWithWarnings(aShapeBase, tol3d, aKeepNonSolids, aWrn);
 734     aShape = GlueWithWarnings(aShapeBase, tol3d, TopAbs_FACE, aKeepNonSolids, aWrn);
 735   }
 736   else if (aType == GLUE_EDGES) {
 737     aShape = GlueWithWarnings(aShapeBase, tol3d, TopAbs_EDGE, aKeepNonSolids, aWrn);
 738   }
 739   else if (aType == GLUE_FACES_BY_LIST || aType == GLUE_EDGES_BY_LIST) {
 740     Handle(TColStd_HSequenceOfTransient) SF = aCI.GetFaces();
 741     TopTools_MapOfShape aFaces;
 742     int i = 1;
 743     for (; i <= SF->Length(); i++) {
 744       Handle(Standard_Transient) anItem = SF->Value(i);
 745       if (anItem.IsNull())
 746         continue;
 747       Handle(GEOM_Function) aRefSh = Handle(GEOM_Function)::DownCast(anItem);
 748       if (aRefSh.IsNull())
 749         continue;
 750       TopoDS_Shape aFace = aRefSh->GetValue();
 751       if (aFace.IsNull())
 752         continue;
 753       if (!aFaces.Contains(aFace))
 754         aFaces.Add(aFace);
 755     }
 756
 757     Standard_Boolean aGlueAllEdges = Standard_False;
 758     if (aType == GLUE_FACES_BY_LIST)
 759       aGlueAllEdges = aCI.GetGlueAllEdges();
 760
 761     //aShape = GlueFacesByList(aShapeBase, tol3d, aKeepNonSolids, aFaces);
 762     aShape = GlueByList(aShapeBase, tol3d, aKeepNonSolids, aFaces, aGlueAllEdges);
 763   }
 764
 765   if (aShape.IsNull()) return 0;
 766
 767   aFunction->SetValue(aShape);
 768
 769   log.SetTouched(Label());
 770
 771   if (!aWrn.IsEmpty()) {
 772     Standard_Failure::Raise(aWrn.ToCString());
 773   }
 774
 775   return 1;
 776 }
 777
 778 //=======================================================================
 779 //function :  GEOMImpl_GlueDriver_Type_
 780 //purpose  :
 781 //=======================================================================
 782 Standard_EXPORT Handle_Standard_Type& GEOMImpl_GlueDriver_Type_()
 783 {
 784
 785   static Handle_Standard_Type aType1 = STANDARD_TYPE(TFunction_Driver);
 786   if (aType1.IsNull()) aType1 = STANDARD_TYPE(TFunction_Driver);
 787   static Handle_Standard_Type aType2 = STANDARD_TYPE(MMgt_TShared);
 788   if (aType2.IsNull()) aType2 = STANDARD_TYPE(MMgt_TShared);
 789   static Handle_Standard_Type aType3 = STANDARD_TYPE(Standard_Transient);
 790   if (aType3.IsNull()) aType3 = STANDARD_TYPE(Standard_Transient);
 791
 792
 793   static Handle_Standard_Transient _Ancestors[]= {aType1,aType2,aType3,NULL};
 794   static Handle_Standard_Type _aType = new Standard_Type("GEOMImpl_GlueDriver",
 795                                                          sizeof(GEOMImpl_GlueDriver),
 796                                                          1,
 797                                                          (Standard_Address)_Ancestors,
 798                                                          (Standard_Address)NULL);
 799
 800   return _aType;
 801 }
 802
 803 //=======================================================================
 804 //function : DownCast
 805 //purpose  :
 806 //=======================================================================
 807 const Handle(GEOMImpl_GlueDriver) Handle(GEOMImpl_GlueDriver)::DownCast(const Handle(Standard_Transient)& AnObject)
 808 {
 809   Handle(GEOMImpl_GlueDriver) _anOtherObject;
 810
 811   if (!AnObject.IsNull()) {
 812      if (AnObject->IsKind(STANDARD_TYPE(GEOMImpl_GlueDriver))) {
 813        _anOtherObject = Handle(GEOMImpl_GlueDriver)((Handle(GEOMImpl_GlueDriver)&)AnObject);
 814      }
 815   }
 816
 817   return _anOtherObject;
 818 }