src/GEOMImpl/GEOMImpl_GlueDriver.cxx

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