Salome HOME
PAL12470: Cannot use boolean operations on COMPOUNDs and COMPSOLIDs. Workaround OCCT...
[modules/geom.git] / src / GEOMImpl / GEOMImpl_BooleanDriver.cxx
1 // Copyright (C) 2005  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
2 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
3 //
4 // This library is free software; you can redistribute it and/or
5 // modify it under the terms of the GNU Lesser General Public
6 // License as published by the Free Software Foundation; either
7 // version 2.1 of the License.
8 //
9 // This library is distributed in the hope that it will be useful
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
12 // Lesser General Public License for more details.
13 //
14 // You should have received a copy of the GNU Lesser General Public
15 // License along with this library; if not, write to the Free Software
16 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
17 //
18 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
19 //
20
21 #include <GEOMImpl_BooleanDriver.hxx>
22 #include <GEOMImpl_IBoolean.hxx>
23 #include <GEOMImpl_Types.hxx>
24 #include <GEOMImpl_GlueDriver.hxx>
25 #include <GEOM_Function.hxx>
26
27 #include <BRep_Builder.hxx>
28 #include <BRepAlgo.hxx>
29 #include <BRepAlgoAPI_Common.hxx>
30 #include <BRepAlgoAPI_Cut.hxx>
31 #include <BRepAlgoAPI_Fuse.hxx>
32 #include <BRepAlgoAPI_Section.hxx>
33 #include <TopoDS_Shape.hxx>
34 #include <TopoDS_Compound.hxx>
35 #include <TopoDS_Iterator.hxx>
36 #include <TopTools_MapOfShape.hxx>
37 #include <TopTools_ListOfShape.hxx>
38 #include <TopTools_ListIteratorOfListOfShape.hxx>
39 #include <Precision.hxx>
40
41 #include <Standard_ConstructionError.hxx>
42 #include <StdFail_NotDone.hxx>
43
44 //=======================================================================
45 //function : GetID
46 //purpose  :
47 //=======================================================================
48 const Standard_GUID& GEOMImpl_BooleanDriver::GetID()
49 {
50   static Standard_GUID aBooleanDriver("FF1BBB21-5D14-4df2-980B-3A668264EA16");
51   return aBooleanDriver;
52 }
53
54
55 //=======================================================================
56 //function : GEOMImpl_BooleanDriver
57 //purpose  :
58 //=======================================================================
59 GEOMImpl_BooleanDriver::GEOMImpl_BooleanDriver()
60 {
61 }
62
63 void AddSimpleShapes(TopoDS_Shape theShape, TopTools_ListOfShape& theList)
64 {
65   if (theShape.ShapeType() != TopAbs_COMPOUND &&
66       theShape.ShapeType() != TopAbs_COMPSOLID) {
67     theList.Append(theShape);
68     return;
69   }
70
71   TopTools_MapOfShape mapShape;
72   TopoDS_Iterator It (theShape, Standard_True, Standard_True);
73
74   for (; It.More(); It.Next()) {
75     TopoDS_Shape aShape_i = It.Value();
76     if (mapShape.Add(aShape_i)) {
77       if (aShape_i.ShapeType() == TopAbs_COMPOUND ||
78           aShape_i.ShapeType() == TopAbs_COMPSOLID) {
79         AddSimpleShapes(aShape_i, theList);
80       } else {
81         theList.Append(aShape_i);
82       }
83     }
84   }
85 }
86
87 //=======================================================================
88 //function : Execute
89 //purpose  :
90 //=======================================================================
91 Standard_Integer GEOMImpl_BooleanDriver::Execute(TFunction_Logbook& log) const
92 {
93   if (Label().IsNull()) return 0;
94   Handle(GEOM_Function) aFunction = GEOM_Function::GetFunction(Label());
95
96   GEOMImpl_IBoolean aCI (aFunction);
97   Standard_Integer aType = aFunction->GetType();
98
99   TopoDS_Shape aShape;
100
101   Handle(GEOM_Function) aRefShape1 = aCI.GetShape1();
102   Handle(GEOM_Function) aRefShape2 = aCI.GetShape2();
103   TopoDS_Shape aShape1 = aRefShape1->GetValue();
104   TopoDS_Shape aShape2 = aRefShape2->GetValue();
105   if (!aShape1.IsNull() && !aShape2.IsNull()) {
106
107     // perform COMMON operation
108     if (aType == BOOLEAN_COMMON) {
109       BRep_Builder B;
110       TopoDS_Compound C;
111       B.MakeCompound(C);
112
113       TopTools_ListOfShape listShape1, listShape2;
114       AddSimpleShapes(aShape1, listShape1);
115       AddSimpleShapes(aShape2, listShape2);
116
117       Standard_Boolean isCompound =
118         (listShape1.Extent() > 1 || listShape2.Extent() > 1);
119
120       TopTools_ListIteratorOfListOfShape itSub1 (listShape1);
121       for (; itSub1.More(); itSub1.Next()) {
122         TopoDS_Shape aValue1 = itSub1.Value();
123         TopTools_ListIteratorOfListOfShape itSub2 (listShape2);
124         for (; itSub2.More(); itSub2.Next()) {
125           TopoDS_Shape aValue2 = itSub2.Value();
126           BRepAlgoAPI_Common BO (aValue1, aValue2);
127           if (!BO.IsDone()) {
128             StdFail_NotDone::Raise("Common operation can not be performed on the given shapes");
129           }
130           if (isCompound)
131             B.Add(C, BO.Shape());
132           else
133             aShape = BO.Shape();
134         }
135       }
136
137       if (isCompound)
138         aShape = GEOMImpl_GlueDriver::GlueFaces(C, Precision::Confusion());
139     }
140
141     // perform CUT operation
142     else if (aType == BOOLEAN_CUT) {
143       BRep_Builder B;
144       TopoDS_Compound C;
145       B.MakeCompound(C);
146
147       TopTools_ListOfShape listShapes, listTools;
148       AddSimpleShapes(aShape1, listShapes);
149       AddSimpleShapes(aShape2, listTools);
150
151       Standard_Boolean isCompound = (listShapes.Extent() > 1);
152
153       TopTools_ListIteratorOfListOfShape itSub1 (listShapes);
154       for (; itSub1.More(); itSub1.Next()) {
155         TopoDS_Shape aCut = itSub1.Value();
156         // tools
157         TopTools_ListIteratorOfListOfShape itSub2 (listTools);
158         for (; itSub2.More(); itSub2.Next()) {
159           TopoDS_Shape aTool = itSub2.Value();
160           BRepAlgoAPI_Cut BO (aCut, aTool);
161           if (!BO.IsDone()) {
162             StdFail_NotDone::Raise("Cut operation can not be performed on the given shapes");
163           }
164           aCut = BO.Shape();
165         }
166         if (isCompound)
167           B.Add(C, aCut);
168         else
169           aShape = aCut;
170       }
171
172       if (isCompound)
173         aShape = GEOMImpl_GlueDriver::GlueFaces(C, Precision::Confusion());
174     }
175
176     // perform FUSE operation
177     else if (aType == BOOLEAN_FUSE) {
178       // check shapes to provide special processing for compounds and compsolids
179       if (aShape1.ShapeType() == TopAbs_COMPOUND ||
180           aShape1.ShapeType() == TopAbs_COMPSOLID ||
181           aShape2.ShapeType() == TopAbs_COMPOUND ||
182           aShape2.ShapeType() == TopAbs_COMPSOLID) {
183         // at least one of arguments is complex
184         StdFail_NotDone::Raise("Fuse operation aborted: one of arguments is a complex shape");
185       } else {
186         // simple arguments
187         BRepAlgoAPI_Fuse BO (aShape1, aShape2);
188         if (!BO.IsDone()) {
189           StdFail_NotDone::Raise("Fuse operation can not be performed on the given shapes");
190         }
191         aShape = BO.Shape();
192       }
193     }
194
195     // perform SECTION operation
196     else if (aType == BOOLEAN_SECTION) {
197       BRep_Builder B;
198       TopoDS_Compound C;
199       B.MakeCompound(C);
200
201       TopTools_ListOfShape listShape1, listShape2;
202       AddSimpleShapes(aShape1, listShape1);
203       AddSimpleShapes(aShape2, listShape2);
204
205       Standard_Boolean isCompound =
206         (listShape1.Extent() > 1 || listShape2.Extent() > 1);
207
208       TopTools_ListIteratorOfListOfShape itSub1 (listShape1);
209       for (; itSub1.More(); itSub1.Next()) {
210         TopoDS_Shape aValue1 = itSub1.Value();
211         TopTools_ListIteratorOfListOfShape itSub2 (listShape2);
212         for (; itSub2.More(); itSub2.Next()) {
213           TopoDS_Shape aValue2 = itSub2.Value();
214           BRepAlgoAPI_Section BO (aValue1, aValue2);
215           if (!BO.IsDone()) {
216             StdFail_NotDone::Raise("Section operation can not be performed on the given shapes");
217           }
218           if (isCompound)
219             B.Add(C, BO.Shape());
220           else
221             aShape = BO.Shape();
222         }
223       }
224
225       if (isCompound)
226         aShape = C;
227     }
228
229     // UNKNOWN operation
230     else {
231     }
232   }
233
234   if (aShape.IsNull()) return 0;
235   if (!BRepAlgo::IsValid(aShape)) {
236     Standard_ConstructionError::Raise("Boolean operation aborted : non valid shape result");
237   }
238
239   aFunction->SetValue(aShape);
240
241   log.SetTouched(Label());
242
243   return 1;
244 }
245
246
247 //=======================================================================
248 //function :  GEOMImpl_BooleanDriver_Type_
249 //purpose  :
250 //=======================================================================
251 Standard_EXPORT Handle_Standard_Type& GEOMImpl_BooleanDriver_Type_()
252 {
253
254   static Handle_Standard_Type aType1 = STANDARD_TYPE(TFunction_Driver);
255   if ( aType1.IsNull()) aType1 = STANDARD_TYPE(TFunction_Driver);
256   static Handle_Standard_Type aType2 = STANDARD_TYPE(MMgt_TShared);
257   if ( aType2.IsNull()) aType2 = STANDARD_TYPE(MMgt_TShared);
258   static Handle_Standard_Type aType3 = STANDARD_TYPE(Standard_Transient);
259   if ( aType3.IsNull()) aType3 = STANDARD_TYPE(Standard_Transient);
260
261
262   static Handle_Standard_Transient _Ancestors[]= {aType1,aType2,aType3,NULL};
263   static Handle_Standard_Type _aType = new Standard_Type("GEOMImpl_BooleanDriver",
264                                                          sizeof(GEOMImpl_BooleanDriver),
265                                                          1,
266                                                          (Standard_Address)_Ancestors,
267                                                          (Standard_Address)NULL);
268
269   return _aType;
270 }
271
272 //=======================================================================
273 //function : DownCast
274 //purpose  :
275 //=======================================================================
276 const Handle(GEOMImpl_BooleanDriver) Handle(GEOMImpl_BooleanDriver)::DownCast(const Handle(Standard_Transient)& AnObject)
277 {
278   Handle(GEOMImpl_BooleanDriver) _anOtherObject;
279
280   if (!AnObject.IsNull()) {
281      if (AnObject->IsKind(STANDARD_TYPE(GEOMImpl_BooleanDriver))) {
282        _anOtherObject = Handle(GEOMImpl_BooleanDriver)((Handle(GEOMImpl_BooleanDriver)&)AnObject);
283      }
284   }
285
286   return _anOtherObject ;
287 }