projects / occt.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
0023648: Add tool for checking shapes on validity for boolean operations.
[occt.git] / src / BRepAlgoAPI / BRepAlgoAPI_Section.cxx
CommitLineData
b311480e 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
7fd59977 21// modified by Michael KLOKOV Wed Mar 6 15:01:25 2002
4e57c75e 22// modified by Eugeny MALTCHIKOV Wed Jul 04 11:13:01 2012
7fd59977 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>
4e57c75e 31
7fd59977 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
4e57c75e 39#include <BOPAlgo_PaveFiller.hxx>
40#include <BOPAlgo_BOP.hxx>
41#include <BOPDS_DS.hxx>
7fd59977 42
43static TopoDS_Shape MakeShape(const Handle(Geom_Surface)& S)
44{
45 GeomAbs_Shape c = S->Continuity();
1c72dff6 46 if (c >= GeomAbs_C2) return BRepBuilderAPI_MakeFace(S, Precision::Confusion());
7fd59977 47 else return BRepBuilderAPI_MakeShell(S);
48}
49
4e57c75e 50static Standard_Boolean HasAncestorFaces(const BOPAlgo_PPaveFiller& theDSFiller,
51 const TopoDS_Shape& E,
52 TopoDS_Shape& F1,
53 TopoDS_Shape& F2);
7fd59977 54
55//=======================================================================
56//function : Constructor
57//purpose :
58//=======================================================================
59BRepAlgoAPI_Section::BRepAlgoAPI_Section(const TopoDS_Shape& Sh1,
4e57c75e 60 const TopoDS_Shape& Sh2,
61 const Standard_Boolean PerformNow)
62: BRepAlgoAPI_BooleanOperation(Sh1, Sh2, BOPAlgo_SECTION)
7fd59977 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 }
4e57c75e 71 if (PerformNow) {
7fd59977 72 Build();
4e57c75e 73 }
7fd59977 74}
75
76BRepAlgoAPI_Section::BRepAlgoAPI_Section(const TopoDS_Shape& aS1,
4e57c75e 77 const TopoDS_Shape& aS2,
78 const BOPAlgo_PaveFiller& aDSF,
79 const Standard_Boolean PerformNow)
80: BRepAlgoAPI_BooleanOperation(aS1, aS2, aDSF, BOPAlgo_SECTION)
7fd59977 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
4e57c75e 90 if (PerformNow) {
7fd59977 91 Build();
92 }
93}
94
95
96//=======================================================================
97//function : Constructor
98//purpose :
99//=======================================================================
100BRepAlgoAPI_Section::BRepAlgoAPI_Section(const TopoDS_Shape& Sh,
4e57c75e 101 const gp_Pln& Pl,
102 const Standard_Boolean PerformNow)
103: BRepAlgoAPI_BooleanOperation(Sh, MakeShape(new Geom_Plane(Pl)), BOPAlgo_SECTION)
7fd59977 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 }
4e57c75e 112 if (PerformNow) {
7fd59977 113 Build();
4e57c75e 114 }
7fd59977 115}
116
117//=======================================================================
118//function : Constructor
119//purpose :
120//=======================================================================
121BRepAlgoAPI_Section::BRepAlgoAPI_Section(const TopoDS_Shape& Sh,
4e57c75e 122 const Handle(Geom_Surface)& Sf,
123 const Standard_Boolean PerformNow)
124: BRepAlgoAPI_BooleanOperation(Sh, MakeShape(Sf), BOPAlgo_SECTION)
7fd59977 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 }
4e57c75e 133 if (PerformNow) {
7fd59977 134 Build();
4e57c75e 135 }
7fd59977 136}
137
138//=======================================================================
139//function : Constructor
140//purpose :
141//=======================================================================
142BRepAlgoAPI_Section::BRepAlgoAPI_Section(const Handle(Geom_Surface)& Sf,
4e57c75e 143 const TopoDS_Shape& Sh,
144 const Standard_Boolean PerformNow)
145: BRepAlgoAPI_BooleanOperation(MakeShape(Sf), Sh, BOPAlgo_SECTION)
7fd59977 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 }
4e57c75e 154 if (PerformNow) {
7fd59977 155 Build();
4e57c75e 156 }
7fd59977 157}
158
159//=======================================================================
160//function : Constructor
161//purpose :
162//=======================================================================
163BRepAlgoAPI_Section::BRepAlgoAPI_Section(const Handle(Geom_Surface)& Sf1,
4e57c75e 164 const Handle(Geom_Surface)& Sf2,
165 const Standard_Boolean PerformNow)
166: BRepAlgoAPI_BooleanOperation(MakeShape(Sf1), MakeShape(Sf2), BOPAlgo_SECTION)
7fd59977 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
4e57c75e 176 if (PerformNow) {
7fd59977 177 Build();
4e57c75e 178 }
7fd59977 179}
180
181//=======================================================================
182//function : Init1
183//purpose :
184//=======================================================================
185void 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()) {
4e57c75e 192 myshapeisnull = Standard_False;
7fd59977 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//=======================================================================
213void BRepAlgoAPI_Section::Init1(const gp_Pln& Pl)
214{
215 Init1(MakeShape(new Geom_Plane(Pl)));
216}
217
218//=======================================================================
219//function : Init1
220//purpose :
221//=======================================================================
222void BRepAlgoAPI_Section::Init1(const Handle(Geom_Surface)& Sf)
223{
224 Init1(MakeShape(Sf));
225}
226
227//=======================================================================
228//function : Init2
229//purpose :
230//=======================================================================
231void 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()) {
4e57c75e 238 myshapeisnull = Standard_False;
7fd59977 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//=======================================================================
260void BRepAlgoAPI_Section::Init2(const gp_Pln& Pl)
261{
262 Init2(MakeShape(new Geom_Plane(Pl)));
263}
264
265//=======================================================================
266//function : Init2
267//purpose :
268//=======================================================================
269void BRepAlgoAPI_Section::Init2(const Handle(Geom_Surface)& Sf)
270{
271 Init2(MakeShape(Sf));
272}
273
274//=======================================================================
275//function : Approximation
276//purpose :
277//=======================================================================
278void 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//=======================================================================
290void 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//=======================================================================
302void 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//=======================================================================
314void BRepAlgoAPI_Section::Build()
315{
316 if(myshapeisnull) {
317 myErrorStatus = 2;
318 NotDone();
319 return;
320 }
321
322 if(myparameterschanged) {
7fd59977 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) {
4e57c75e 332 BOPAlgo_SectionAttribute theSecAttr(myApprox, myComputePCurve1, myComputePCurve2);
333 myDSFiller->SetSectionAttribute(theSecAttr);
334 myDSFiller->Perform();
7fd59977 335 }
4e57c75e 336 //
337 BRepAlgoAPI_BooleanOperation::Build();
338 //
7fd59977 339 myparameterschanged = Standard_False;
340 }
341}
342
343//=======================================================================
344//function : HasAncestorFaceOn1
345//purpose :
346//=======================================================================
347Standard_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;
4e57c75e 358 aResult = HasAncestorFaces(myDSFiller, E, F1, F2);
7fd59977 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//=======================================================================
371Standard_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;
4e57c75e 382 aResult = HasAncestorFaces(myDSFiller, E, F1, F2);
7fd59977 383
384 if(F2.IsNull()) {
385 return Standard_False;
386 }
387 F = F2;
388 return aResult;
389}
390
391//=======================================================================
392//function : PCurveOn1
393//purpose :
394//=======================================================================
395Handle(Geom2d_Curve) BRepAlgoAPI_Section::PCurveOn1(const TopoDS_Shape& E) const
396{
397 Handle(Geom2d_Curve) aResult;
398
399 if(myComputePCurve1) {
400 TopoDS_Shape aShape;
401
402 if(HasAncestorFaceOn1(E, aShape)) {
403 const TopoDS_Edge& anEdge = TopoDS::Edge(E);
404 const TopoDS_Face& aFace = TopoDS::Face(aShape);
405 Standard_Real f, l;
406 aResult = BRep_Tool::CurveOnSurface(anEdge, aFace, f, l);
407
408 if(!aResult->IsKind(STANDARD_TYPE(Geom2d_TrimmedCurve))) {
4e57c75e 409 aResult = new Geom2d_TrimmedCurve(aResult, f, l);
7fd59977 410 }
411 }
412 }
413 return aResult;
414}
415
416//=======================================================================
417//function : PCurveOn2
418//purpose :
419//=======================================================================
420Handle(Geom2d_Curve) BRepAlgoAPI_Section::PCurveOn2(const TopoDS_Shape& E) const
421{
422 Handle(Geom2d_Curve) aResult;
423
424 if(myComputePCurve2) {
425 TopoDS_Shape aShape;
426
427 if(HasAncestorFaceOn2(E, aShape)) {
428 const TopoDS_Edge& anEdge = TopoDS::Edge(E);
429 const TopoDS_Face& aFace = TopoDS::Face(aShape);
430 Standard_Real f, l;
431 aResult = BRep_Tool::CurveOnSurface(anEdge, aFace, f, l);
432
433 if(!aResult->IsKind(STANDARD_TYPE(Geom2d_TrimmedCurve))) {
4e57c75e 434 aResult = new Geom2d_TrimmedCurve(aResult, f, l);
7fd59977 435 }
436 }
437 }
438 return aResult;
439}
440
441//=======================================================================
442//function : InitParameters
443//purpose :
444//=======================================================================
445void BRepAlgoAPI_Section::InitParameters()
446{
447 myparameterschanged = Standard_False;
448 myshapeisnull = Standard_False;
449 myApprox = Standard_False;
450 myComputePCurve1 = Standard_False;
451 myComputePCurve2 = Standard_False;
452}
453
454// ------------------------------------------------------------------------
455// static function : HasAncestorFaces
456// purpose :
457// ------------------------------------------------------------------------
4e57c75e 458static Standard_Boolean HasAncestorFaces(const BOPAlgo_PPaveFiller& theDSFiller,
459 const TopoDS_Shape& E,
460 TopoDS_Shape& F1,
461 TopoDS_Shape& F2) {
462
463 Standard_Integer aNb, i, j, nE, nF1, nF2, aNbCurves;;
464 //
465 const BOPDS_PDS& pDS = theDSFiller->PDS();
466 BOPDS_VectorOfInterfFF& aFFs=pDS->InterfFF();
467 //
468 aNb=aFFs.Extent();
469
470 for (i = 0; i < aNb; i++) {
471 BOPDS_InterfFF& aFFi=aFFs(i);
472 aFFi.Indices(nF1, nF2);
473 //
474 //section edges
475 const BOPDS_VectorOfCurve& aSeqOfCurve=aFFi.Curves();
476 aNbCurves=aSeqOfCurve.Extent();
477 for (j=0; j<aNbCurves; j++) {
478 const BOPDS_Curve& aCurve=aSeqOfCurve(j);
479
480 const BOPDS_ListOfPaveBlock& aSectEdges = aCurve.PaveBlocks();
481
482 BOPDS_ListIteratorOfListOfPaveBlock anIt;
483 anIt.Initialize(aSectEdges);
7fd59977 484
485 for(; anIt.More(); anIt.Next()) {
4e57c75e 486 const Handle(BOPDS_PaveBlock)& aPB = anIt.Value();
487 nE = aPB->Edge();
488 if(nE < 0) continue;
489 //
490 if(E.IsSame(pDS->Shape(nE))) {
491 F1 = pDS->Shape(nF1);
492 F2 = pDS->Shape(nF2);
493 return Standard_True;
494 }
7fd59977 495 }
496 }
497
4e57c75e 498 //existing pave blocks
499 BOPCol_ListOfInteger aLSE;
500 BOPCol_ListIteratorOfListOfInteger aItLSE;
501 //
502 pDS->SharedEdges(nF1, nF2, aLSE, theDSFiller->Allocator());
503 aItLSE.Initialize(aLSE);
504 for (; aItLSE.More(); aItLSE.Next()) {
505 nE = aItLSE.Value();
506 //
507 if(E.IsSame(pDS->Shape(nE))) {
508 F1 = pDS->Shape(nF1);
509 F2 = pDS->Shape(nF2);
510 return Standard_True;
7fd59977 511 }
512 }
513 }
514 return Standard_False;
515}
Open CASCADE Technology repository