src/GeomAlgoAPI/GeomAlgoAPI_Intersection.cpp

   1 // Copyright (C) 2014-2022  CEA/DEN, EDF R&D
   2 //
   3 // This library is free software; you can redistribute it and/or
   4 // modify it under the terms of the GNU Lesser General Public
   5 // License as published by the Free Software Foundation; either
   6 // version 2.1 of the License, or (at your option) any later version.
   7 //
   8 // This library is distributed in the hope that it will be useful,
   9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11 // Lesser General Public License for more details.
  12 //
  13 // You should have received a copy of the GNU Lesser General Public
  14 // License along with this library; if not, write to the Free Software
  15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  16 //
  17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  18 //
  19
  20 #include "GeomAlgoAPI_Intersection.h"
  21
  22 #include <GeomAlgoAPI_DFLoader.h>
  23
  24 #include <BOPAlgo_PaveFiller.hxx>
  25 #include <BOPAlgo_Section.hxx>
  26
  27
  28 //==================================================================================================
  29 GeomAlgoAPI_Intersection::GeomAlgoAPI_Intersection(const ListOfShape& theObjects, const double theFuzzy)
  30   : myFiller(0)
  31 {
  32   build(theObjects, theFuzzy);
  33 }
  34
  35 //==================================================================================================
  36 GeomAlgoAPI_Intersection::~GeomAlgoAPI_Intersection() {
  37   if (myFiller)
  38     delete (BOPAlgo_PaveFiller*)myFiller;
  39 }
  40
  41 //==================================================================================================
  42 void GeomAlgoAPI_Intersection::build(const ListOfShape& theObjects, const double theFuzzy)
  43 {
  44   if (theObjects.empty()) {
  45     return;
  46   }
  47
  48   // Creating partition operation.
  49   BOPAlgo_Section* anOperation = new BOPAlgo_Section;
  50   this->setImpl(anOperation);
  51   this->setBuilderType(OCCT_BOPAlgo_Builder);
  52
  53   // Getting objects.
  54   TopTools_ListOfShape anObjects;
  55   for (ListOfShape::const_iterator
  56     anObjectsIt = theObjects.begin(); anObjectsIt != theObjects.end(); anObjectsIt++) {
  57     const TopoDS_Shape& aShape = (*anObjectsIt)->impl<TopoDS_Shape>();
  58     if(!aShape.IsNull()) {
  59       anObjects.Append(aShape);
  60     }
  61   }
  62
  63   BOPAlgo_PaveFiller* aDSFiller = new BOPAlgo_PaveFiller;
  64   myFiller = aDSFiller;
  65   aDSFiller->SetArguments(anObjects);
  66
  67   aDSFiller->SetRunParallel(false);
  68   aDSFiller->SetNonDestructive(false);
  69   aDSFiller->SetGlue(BOPAlgo_GlueOff);
  70   if (theFuzzy > 0) aDSFiller->SetFuzzyValue(theFuzzy);
  71
  72   // optimization for the issue #2399
  73   BOPAlgo_SectionAttribute theSecAttr(Standard_True,
  74                                       Standard_True,
  75                                       Standard_True);
  76   aDSFiller->SetSectionAttribute(theSecAttr);
  77
  78   aDSFiller->Perform();
  79   if (aDSFiller->HasErrors()) {
  80     return;
  81   }
  82
  83   anOperation->SetArguments(anObjects);
  84   anOperation->SetRunParallel(false);
  85   anOperation->SetCheckInverted(true);
  86   if (theFuzzy > 0) anOperation->SetFuzzyValue(theFuzzy);
  87
  88   anOperation->PerformWithFiller(*aDSFiller); // it references a filler fields, so keep the filler
  89   myFiller = 0;
  90   if(anOperation->HasErrors()) {
  91     return;
  92   }
  93   TopoDS_Shape aResult = anOperation->Shape();
  94
  95   if(aResult.ShapeType() == TopAbs_COMPOUND) {
  96     aResult = GeomAlgoAPI_DFLoader::refineResult(aResult);
  97   }
  98
  99   std::shared_ptr<GeomAPI_Shape> aShape(new GeomAPI_Shape());
 100   aShape->setImpl(new TopoDS_Shape(aResult));
 101   this->setShape(aShape);
 102   this->setDone(true);
 103 }