0023881: BRepAlgoAPI_Section HasAncestorFaceOn1 returned False on the boundary(6...
[occt.git] / src / BRepAlgoAPI / BRepAlgoAPI_Section.cxx
1 // Created on: 1994-02-18
2 // Created by: Remi LEQUETTE
3 // Copyright (c) 1994-1999 Matra Datavision
4 // Copyright (c) 1999-2012 OPEN CASCADE SAS
5 //
6 // The content of this file is subject to the Open CASCADE Technology Public
7 // License Version 6.5 (the "License"). You may not use the content of this file
8 // except in compliance with the License. Please obtain a copy of the License
9 // at http://www.opencascade.org and read it completely before using this file.
10 //
11 // The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
12 // main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
13 //
14 // The Original Code and all software distributed under the License is
15 // distributed on an "AS IS" basis, without warranty of any kind, and the
16 // Initial Developer hereby disclaims all such warranties, including without
17 // limitation, any warranties of merchantability, fitness for a particular
18 // purpose or non-infringement. Please see the License for the specific terms
19 // and conditions governing the rights and limitations under the License.
20
21 // modified by Michael KLOKOV  Wed Mar  6 15:01:25 2002
22 // modified by  Eugeny MALTCHIKOV Wed Jul 04 11:13:01 2012 
23
24 #include <BRepAlgoAPI_Section.ixx>
25
26
27 #include <BRepBuilderAPI_MakeFace.hxx>
28 #include <BRepBuilderAPI_MakeShell.hxx>
29 #include <Geom_Plane.hxx>
30 #include <Geom2d_TrimmedCurve.hxx>
31
32 #include <TopExp.hxx>
33 #include <TopTools_IndexedMapOfShape.hxx>
34 #include <TopTools_MapOfShape.hxx>
35 #include <TopoDS.hxx>
36 #include <TopoDS_Face.hxx>
37 #include <BRep_Tool.hxx>
38
39 #include <BOPAlgo_PaveFiller.hxx>
40 #include <BOPAlgo_BOP.hxx>
41 #include <BOPDS_DS.hxx>
42
43 static TopoDS_Shape MakeShape(const Handle(Geom_Surface)& S)
44 {
45   GeomAbs_Shape c = S->Continuity();
46   if (c >= GeomAbs_C2) return BRepBuilderAPI_MakeFace(S, Precision::Confusion());
47   else return BRepBuilderAPI_MakeShell(S);
48 }
49
50 static Standard_Boolean HasAncestorFaces(const BOPAlgo_PPaveFiller& theDSFiller, 
51                                          const TopoDS_Shape&      E,
52                                          TopoDS_Shape&            F1,
53                                          TopoDS_Shape&            F2);
54
55 //=======================================================================
56 //function : Constructor
57 //purpose  : 
58 //=======================================================================
59 BRepAlgoAPI_Section::BRepAlgoAPI_Section(const TopoDS_Shape& Sh1,
60                                          const TopoDS_Shape& Sh2,
61                                          const Standard_Boolean PerformNow)
62 : BRepAlgoAPI_BooleanOperation(Sh1, Sh2, BOPAlgo_SECTION)
63 {
64   InitParameters();
65   myparameterschanged = Standard_True;
66
67   if(myS1.IsNull() || myS2.IsNull()) {
68     //     StdFail_NotDone::Raise("BRepAlgoAPI_Section : Shape NULL");
69     myshapeisnull = Standard_True;
70   }
71   if (PerformNow) {
72     Build();
73   }
74 }
75
76 BRepAlgoAPI_Section::BRepAlgoAPI_Section(const TopoDS_Shape&      aS1, 
77                                          const TopoDS_Shape&      aS2,
78                                          const BOPAlgo_PaveFiller& aDSF,
79                                          const Standard_Boolean   PerformNow)
80 : BRepAlgoAPI_BooleanOperation(aS1, aS2, aDSF, BOPAlgo_SECTION)
81 {
82   InitParameters();
83   myparameterschanged = Standard_True;
84
85   if(myS1.IsNull() || myS2.IsNull()) {
86     //     StdFail_NotDone::Raise("BRepAlgoAPI_Section : Shape NULL");
87     myshapeisnull = Standard_True;
88   }
89
90   if (PerformNow) {
91     Build();
92   }
93 }
94
95
96 //=======================================================================
97 //function : Constructor
98 //purpose  : 
99 //=======================================================================
100 BRepAlgoAPI_Section::BRepAlgoAPI_Section(const TopoDS_Shape&    Sh,
101                                          const gp_Pln&          Pl,
102                                          const Standard_Boolean PerformNow)
103 : BRepAlgoAPI_BooleanOperation(Sh, MakeShape(new Geom_Plane(Pl)), BOPAlgo_SECTION)
104 {
105   InitParameters();
106   myparameterschanged = Standard_True;
107
108   if(Sh.IsNull() || myS2.IsNull()) {
109     //     StdFail_NotDone::Raise("BRepAlgoAPI_Section : Shape NULL");
110     myshapeisnull = Standard_True;
111   }
112   if (PerformNow) {
113     Build();
114   }
115 }
116
117 //=======================================================================
118 //function : Constructor
119 //purpose  : 
120 //=======================================================================
121 BRepAlgoAPI_Section::BRepAlgoAPI_Section(const TopoDS_Shape&         Sh,
122                                                const Handle(Geom_Surface)& Sf,
123                                                const Standard_Boolean      PerformNow)
124 : BRepAlgoAPI_BooleanOperation(Sh, MakeShape(Sf), BOPAlgo_SECTION)
125 {
126   InitParameters();
127   myparameterschanged = Standard_True;
128
129   if(Sh.IsNull() || myS2.IsNull()) {
130 //     StdFail_NotDone::Raise("BRepAlgoAPI_Section : Shape NULL");
131     myshapeisnull = Standard_True;
132   }
133   if (PerformNow) {
134     Build();
135   }
136 }
137
138 //=======================================================================
139 //function : Constructor
140 //purpose  : 
141 //=======================================================================
142 BRepAlgoAPI_Section::BRepAlgoAPI_Section(const Handle(Geom_Surface)& Sf,
143                                                const TopoDS_Shape&         Sh,
144                                                const Standard_Boolean      PerformNow)
145 : BRepAlgoAPI_BooleanOperation(MakeShape(Sf), Sh, BOPAlgo_SECTION)
146 {
147   InitParameters();
148   myparameterschanged = Standard_True;
149
150   if(myS1.IsNull() || Sh.IsNull()) {
151 //     StdFail_NotDone::Raise("BRepAlgoAPI_Section : Shape NULL");
152     myshapeisnull = Standard_True;
153   }
154   if (PerformNow) {
155     Build();
156   }
157 }
158
159 //=======================================================================
160 //function : Constructor
161 //purpose  : 
162 //=======================================================================
163 BRepAlgoAPI_Section::BRepAlgoAPI_Section(const Handle(Geom_Surface)& Sf1,
164                                                const Handle(Geom_Surface)& Sf2,
165                                                const Standard_Boolean      PerformNow)
166 : BRepAlgoAPI_BooleanOperation(MakeShape(Sf1), MakeShape(Sf2), BOPAlgo_SECTION)
167 {
168   InitParameters();
169   myparameterschanged = Standard_True;
170
171   if(myS1.IsNull() || myS2.IsNull()) {
172 //     StdFail_NotDone::Raise("BRepAlgoAPI_Section : Shape NULL");
173     myshapeisnull = Standard_True;
174   }
175
176   if (PerformNow) {
177     Build();
178   }
179 }
180
181 //=======================================================================
182 //function : Init1
183 //purpose  : 
184 //=======================================================================
185 void BRepAlgoAPI_Section::Init1(const TopoDS_Shape& S1) 
186 {
187   if(!S1.IsNull()) {
188     if (!S1.IsEqual(myS1)) {
189       myS1 = S1;
190
191       if(!myS2.IsNull()) {
192         myshapeisnull = Standard_False;
193       }
194       myparameterschanged = Standard_True;
195     }
196   }
197   else {
198     if(!myS1.IsNull()) {
199       myS1 = S1;
200       myshapeisnull = Standard_True;
201       myparameterschanged = Standard_True;
202     }
203   }
204
205   if(myparameterschanged)
206     NotDone();
207 }
208
209 //=======================================================================
210 //function : Init1
211 //purpose  : 
212 //=======================================================================
213 void BRepAlgoAPI_Section::Init1(const gp_Pln& Pl) 
214 {
215   Init1(MakeShape(new Geom_Plane(Pl)));
216 }
217
218 //=======================================================================
219 //function : Init1
220 //purpose  : 
221 //=======================================================================
222 void BRepAlgoAPI_Section::Init1(const Handle(Geom_Surface)& Sf) 
223 {
224   Init1(MakeShape(Sf));
225 }
226
227 //=======================================================================
228 //function : Init2
229 //purpose  : 
230 //=======================================================================
231 void BRepAlgoAPI_Section::Init2(const TopoDS_Shape& S2) 
232 {
233   if(!S2.IsNull()) {
234     if (!S2.IsEqual(myS2)) {
235       myS2 = S2;
236
237       if(!myS1.IsNull()) {
238         myshapeisnull = Standard_False;
239       }
240       myparameterschanged = Standard_True;
241     }
242   }
243   else {
244     if(!myS2.IsNull()) {
245       myS2 = S2;
246       myshapeisnull = Standard_True;
247       myparameterschanged = Standard_True;
248     }
249   }
250
251   if(myparameterschanged) {
252     NotDone();
253   }
254 }
255
256 //=======================================================================
257 //function : Init2
258 //purpose  : 
259 //=======================================================================
260 void BRepAlgoAPI_Section::Init2(const gp_Pln& Pl) 
261 {
262   Init2(MakeShape(new Geom_Plane(Pl)));
263 }
264
265 //=======================================================================
266 //function : Init2
267 //purpose  : 
268 //=======================================================================
269 void BRepAlgoAPI_Section::Init2(const Handle(Geom_Surface)& Sf) 
270 {
271   Init2(MakeShape(Sf));
272 }
273
274 //=======================================================================
275 //function : Approximation
276 //purpose  : 
277 //=======================================================================
278 void BRepAlgoAPI_Section::Approximation(const Standard_Boolean B) 
279 {
280   if(myApprox != B) {
281     myApprox = B;
282     myparameterschanged = Standard_True;
283   }
284 }
285
286 //=======================================================================
287 //function : ComputePCurveOn1
288 //purpose  : 
289 //=======================================================================
290 void BRepAlgoAPI_Section::ComputePCurveOn1(const Standard_Boolean B) 
291 {
292   if(myComputePCurve1 != B) {
293     myComputePCurve1 = B;
294     myparameterschanged = Standard_True;
295   }
296 }
297
298 //=======================================================================
299 //function : ComputePCurveOn2
300 //purpose  : 
301 //=======================================================================
302 void BRepAlgoAPI_Section::ComputePCurveOn2(const Standard_Boolean B) 
303 {
304   if(myComputePCurve2 != B) {
305     myComputePCurve2 = B;
306     myparameterschanged = Standard_True;
307   }
308 }
309
310 //=======================================================================
311 //function : Build
312 //purpose  : 
313 //=======================================================================
314 void BRepAlgoAPI_Section::Build() 
315 {
316   if(myshapeisnull) {
317     myErrorStatus = 2;
318     NotDone();
319     return;
320   }
321
322   if(myparameterschanged) {
323     Standard_Boolean bIsNewFiller = PrepareFiller();
324     //
325     if (myErrorStatus!=1) {
326       NotDone();
327       // there were errors during the preparation 
328       return;
329     }
330     //
331     if (bIsNewFiller) {
332       BOPAlgo_SectionAttribute theSecAttr(myApprox, myComputePCurve1, myComputePCurve2);
333       myDSFiller->SetSectionAttribute(theSecAttr);
334       myDSFiller->Perform();
335     }
336     //
337     BRepAlgoAPI_BooleanOperation::Build();
338     //
339     myparameterschanged = Standard_False;
340   }
341 }
342
343 //=======================================================================
344 //function : HasAncestorFaceOn1
345 //purpose  : 
346 //=======================================================================
347 Standard_Boolean BRepAlgoAPI_Section::HasAncestorFaceOn1(const TopoDS_Shape& E, TopoDS_Shape& F) const
348 {
349   Standard_Boolean aResult = Standard_False;
350   if(E.IsNull()) {
351     return aResult;
352   }
353
354   if(E.ShapeType() != TopAbs_EDGE) {
355     return aResult;
356   }
357   TopoDS_Shape F1, F2;
358   aResult = HasAncestorFaces(myDSFiller, E, F1, F2);
359
360   if(F1.IsNull()) {
361     return Standard_False;
362   }
363   F = F1;
364   return aResult;
365 }
366
367 //=======================================================================
368 //function : HasAncestorFaceOn2
369 //purpose  : 
370 //=======================================================================
371 Standard_Boolean BRepAlgoAPI_Section::HasAncestorFaceOn2(const TopoDS_Shape& E,TopoDS_Shape& F) const
372 {
373   Standard_Boolean aResult = Standard_False;
374   if(E.IsNull()) {
375     return aResult;
376   }
377
378   if(E.ShapeType() != TopAbs_EDGE) {
379     return aResult;
380   }
381   TopoDS_Shape F1, F2;
382   aResult = HasAncestorFaces(myDSFiller, E, F1, F2);
383
384   if(F2.IsNull()) {
385     return Standard_False;
386   }
387   F = F2;
388   return aResult;
389 }
390
391 //=======================================================================
392 //function : InitParameters
393 //purpose  : 
394 //=======================================================================
395 void BRepAlgoAPI_Section::InitParameters() 
396 {
397   myparameterschanged = Standard_False;
398   myshapeisnull = Standard_False;
399   myApprox = Standard_False;
400   myComputePCurve1 = Standard_False;
401   myComputePCurve2 = Standard_False;
402 }
403
404 // ------------------------------------------------------------------------
405 // static function : HasAncestorFaces
406 // purpose  : 
407 // ------------------------------------------------------------------------
408 static Standard_Boolean HasAncestorFaces(const BOPAlgo_PPaveFiller& theDSFiller, 
409                                          const TopoDS_Shape&      E,
410                                          TopoDS_Shape&            F1,
411                                          TopoDS_Shape&            F2) {
412
413   Standard_Integer aNb, i, j, nE, nF1, nF2, aNbCurves;
414   //
415   const BOPDS_PDS& pDS = theDSFiller->PDS();
416   BOPDS_VectorOfInterfFF& aFFs=pDS->InterfFF();
417   //
418   aNb=aFFs.Extent();
419   //section edges
420   for (i = 0; i < aNb; i++) {
421     BOPDS_InterfFF& aFFi=aFFs(i);
422     aFFi.Indices(nF1, nF2);
423     //
424     const BOPDS_VectorOfCurve& aSeqOfCurve=aFFi.Curves();
425     aNbCurves=aSeqOfCurve.Extent();
426     for (j=0; j<aNbCurves; j++) {
427       const BOPDS_Curve& aCurve=aSeqOfCurve(j);
428
429       const BOPDS_ListOfPaveBlock& aSectEdges = aCurve.PaveBlocks();
430
431       BOPDS_ListIteratorOfListOfPaveBlock anIt;
432       anIt.Initialize(aSectEdges);
433
434       for(; anIt.More(); anIt.Next()) {
435         const Handle(BOPDS_PaveBlock)& aPB = anIt.Value();
436         nE = aPB->Edge();
437         if(nE < 0)  continue;
438         //
439         if(E.IsSame(pDS->Shape(nE))) {
440           F1 = pDS->Shape(nF1);
441           F2 = pDS->Shape(nF2);
442           return Standard_True;
443         }
444       }
445     }
446   }
447   return Standard_False;
448 }