Salome HOME
Issue #2565: CEA 2018-1 Intersection
[modules/shaper.git] / src / GeomAlgoAPI / GeomAlgoAPI_Intersection.cpp
1 // Copyright (C) 2014-2017  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
18 // email : webmaster.salome@opencascade.com<mailto:webmaster.salome@opencascade.com>
19 //
20
21 #include "GeomAlgoAPI_Intersection.h"
22
23 #include <GeomAlgoAPI_DFLoader.h>
24
25 #include <BOPAlgo_PaveFiller.hxx>
26 #include <BOPAlgo_Section.hxx>
27
28 //==================================================================================================
29 GeomAlgoAPI_Intersection::GeomAlgoAPI_Intersection(const ListOfShape& theObjects)
30 {
31   build(theObjects);
32 }
33
34 //==================================================================================================
35 void GeomAlgoAPI_Intersection::build(const ListOfShape& theObjects)
36 {
37   if (theObjects.empty()) {
38     return;
39   }
40
41   // Creating partition operation.
42   BOPAlgo_Section* anOperation = new BOPAlgo_Section;
43   this->setImpl(anOperation);
44   this->setBuilderType(OCCT_BOPAlgo_Builder);
45
46   // Getting objects.
47   TopTools_ListOfShape anObjects;
48   for (ListOfShape::const_iterator
49     anObjectsIt = theObjects.begin(); anObjectsIt != theObjects.end(); anObjectsIt++) {
50     const TopoDS_Shape& aShape = (*anObjectsIt)->impl<TopoDS_Shape>();
51     if(!aShape.IsNull()) {
52       anObjects.Append(aShape);
53     }
54   }
55
56   BOPAlgo_PaveFiller aDSFiller;
57   aDSFiller.SetArguments(anObjects);
58
59   aDSFiller.SetRunParallel(false);
60   aDSFiller.SetNonDestructive(false);
61   aDSFiller.SetGlue(BOPAlgo_GlueOff);
62
63   // optimization for the issue #2399
64   BOPAlgo_SectionAttribute theSecAttr(Standard_True,
65                                       Standard_True,
66                                       Standard_True);
67   aDSFiller.SetSectionAttribute(theSecAttr);
68
69   aDSFiller.Perform();
70   if (aDSFiller.HasErrors()) {
71     return;
72   }
73
74   anOperation->SetArguments(anObjects);
75   anOperation->SetRunParallel(false);
76   anOperation->SetCheckInverted(true);
77
78   anOperation->PerformWithFiller(aDSFiller);
79   if(anOperation->HasErrors()) {
80     return;
81   }
82   TopoDS_Shape aResult = anOperation->Shape();
83
84   if(aResult.ShapeType() == TopAbs_COMPOUND) {
85     aResult = GeomAlgoAPI_DFLoader::refineResult(aResult);
86   }
87
88   std::shared_ptr<GeomAPI_Shape> aShape(new GeomAPI_Shape());
89   aShape->setImpl(new TopoDS_Shape(aResult));
90   this->setShape(aShape);
91   this->setDone(true);
92 }