0025748: Parallel version of progress indicator
[occt.git] / src / BOPAlgo / BOPAlgo_CheckerSI.cxx
1 // Created by: Peter KURNEV
2 // Copyright (c) 2010-2014 OPEN CASCADE SAS
3 // Copyright (c) 2007-2010 CEA/DEN, EDF R&D, OPEN CASCADE
4 // Copyright (c) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN, CEDRAT,
5 //                         EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
6 //
7 // This file is part of Open CASCADE Technology software library.
8 //
9 // This library is free software; you can redistribute it and/or modify it under
10 // the terms of the GNU Lesser General Public License version 2.1 as published
11 // by the Free Software Foundation, with special exception defined in the file
12 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
13 // distribution for complete text of the license and disclaimer of any warranty.
14 //
15 // Alternatively, this file may be used under the terms of Open CASCADE
16 // commercial license or contractual agreement.
17
18 //
19
20 #include <BOPAlgo_CheckerSI.hxx>
21 #include <BOPAlgo_Alerts.hxx>
22 #include <BOPDS_DS.hxx>
23 #include <BOPDS_Interf.hxx>
24 #include <BOPDS_IteratorSI.hxx>
25 #include <BOPDS_MapOfPair.hxx>
26 #include <BOPDS_Pair.hxx>
27 #include <BOPDS_PIteratorSI.hxx>
28 #include <BOPDS_VectorOfInterfEF.hxx>
29 #include <BOPDS_VectorOfInterfFF.hxx>
30 #include <BOPDS_VectorOfInterfVE.hxx>
31 #include <BOPDS_VectorOfInterfVF.hxx>
32 #include <BOPDS_VectorOfInterfVV.hxx>
33 #include <BRep_Tool.hxx>
34 #include <gp_Torus.hxx>
35 #include <TopExp.hxx>
36 #include <BOPTools_AlgoTools.hxx>
37 #include <BOPTools_Parallel.hxx>
38 #include <BRepBuilderAPI_Copy.hxx>
39 #include <IntTools_Context.hxx>
40 #include <IntTools_Tools.hxx>
41 #include <IntTools_FaceFace.hxx>
42 #include <Standard_ErrorHandler.hxx>
43 #include <Standard_Failure.hxx>
44 #include <TopTools_ListOfShape.hxx>
45 #include <TopTools_MapOfShape.hxx>
46
47 //=======================================================================
48 //class    : BOPAlgo_FaceSelfIntersect
49 //purpose  : 
50 //=======================================================================
51 class BOPAlgo_FaceSelfIntersect : 
52   public IntTools_FaceFace,
53   public BOPAlgo_Algo {
54
55  public:
56   DEFINE_STANDARD_ALLOC
57
58   BOPAlgo_FaceSelfIntersect() : 
59     IntTools_FaceFace(),  
60     BOPAlgo_Algo(),
61     myIF(-1), myTolF(1.e-7) {
62   }
63   //
64   virtual ~BOPAlgo_FaceSelfIntersect() {
65   }
66   //
67   void SetIndex(const Standard_Integer nF) {
68     myIF = nF;
69   }
70   //
71   Standard_Integer IndexOfFace() const {
72     return myIF;
73   }
74   //
75   void SetFace(const TopoDS_Face& aF) {
76     myF = aF;
77   }
78   //
79   const TopoDS_Face& Face()const {
80     return myF;
81   }
82   //
83   void SetTolF(const Standard_Real aTolF) {
84     myTolF = aTolF;
85   }
86   //
87   Standard_Real TolF() const{
88     return myTolF;
89   }
90   //
91   virtual void Perform() {
92     BOPAlgo_Algo::UserBreak();
93     IntTools_FaceFace::Perform(myF, myF);
94   }
95   //
96  protected:
97   Standard_Integer myIF;
98   Standard_Real myTolF;
99   TopoDS_Face myF;
100 };
101 //end of definition of class BOPAlgo_FaceSelfIntersect
102
103 //=======================================================================
104
105 typedef NCollection_Vector<BOPAlgo_FaceSelfIntersect> BOPAlgo_VectorOfFaceSelfIntersect;
106
107 //=======================================================================
108 //function : 
109 //purpose  : 
110 //=======================================================================
111 BOPAlgo_CheckerSI::BOPAlgo_CheckerSI()
112 :
113   BOPAlgo_PaveFiller()
114 {
115   myLevelOfCheck=BOPDS_DS::NbInterfTypes()-1;
116   myNonDestructive=Standard_True;
117   SetAvoidBuildPCurve(Standard_True);
118 }
119 //=======================================================================
120 //function : ~
121 //purpose  : 
122 //=======================================================================
123 BOPAlgo_CheckerSI::~BOPAlgo_CheckerSI()
124 {
125 }
126 //=======================================================================
127 //function : SetLevelOfCheck
128 //purpose  : 
129 //=======================================================================
130 void BOPAlgo_CheckerSI::SetLevelOfCheck(const Standard_Integer theLevel)
131 {
132   Standard_Integer aNbLists;
133   //
134   aNbLists=BOPDS_DS::NbInterfTypes();
135   if (theLevel >= 0 && theLevel < aNbLists) {
136     myLevelOfCheck = theLevel;
137   }
138 }
139 //=======================================================================
140 //function : Init
141 //purpose  : 
142 //=======================================================================
143 void BOPAlgo_CheckerSI::Init()
144 {
145   Clear();
146   //
147   // 1. myDS
148   myDS=new BOPDS_DS(myAllocator);
149   myDS->SetArguments(myArguments);
150   myDS->Init(myFuzzyValue);
151   //
152   // 2 myContext
153   myContext=new IntTools_Context;
154   //
155   // 3.myIterator 
156   BOPDS_PIteratorSI theIterSI=new BOPDS_IteratorSI(myAllocator);
157   theIterSI->SetDS(myDS);
158   theIterSI->Prepare(myContext, myUseOBB, myFuzzyValue);
159   theIterSI->UpdateByLevelOfCheck(myLevelOfCheck);
160   //
161   myIterator=theIterSI;
162 }
163 //=======================================================================
164 //function : Perform
165 //purpose  : 
166 //=======================================================================
167 void BOPAlgo_CheckerSI::Perform()
168 {
169   try {
170     OCC_CATCH_SIGNALS
171     //
172     if (myArguments.Extent() != 1) {
173       AddError (new BOPAlgo_AlertMultipleArguments);
174       return;
175     }
176     //
177     // Perform intersection of sub shapes
178     BOPAlgo_PaveFiller::Perform();
179     //
180     CheckFaceSelfIntersection();
181     
182     // Perform intersection with solids
183     if (!HasErrors())
184       PerformVZ();
185     //
186     if (!HasErrors())
187       PerformEZ();
188     //
189     if (!HasErrors())
190       PerformFZ();
191     //
192     if (!HasErrors())
193       PerformZZ();
194     //
195     // Treat the intersection results
196     PostTreat();
197   }
198   //
199   catch (Standard_Failure const&) {
200     AddError (new BOPAlgo_AlertIntersectionFailed);
201   }
202 }
203 //=======================================================================
204 //function : PostTreat
205 //purpose  : 
206 //=======================================================================
207 void BOPAlgo_CheckerSI::PostTreat()
208 {
209   Standard_Integer i, aNb, n1, n2; 
210   BOPDS_Pair aPK;
211   //
212   BOPDS_MapOfPair& aMPK=
213     *((BOPDS_MapOfPair*)&myDS->Interferences());
214
215   // 0
216   BOPDS_VectorOfInterfVV& aVVs=myDS->InterfVV();
217   aNb=aVVs.Length();
218   for (i=0; i!=aNb; ++i) {
219     const BOPDS_InterfVV& aVV=aVVs(i);
220     aVV.Indices(n1, n2);
221     if (myDS->IsNewShape(n1) || myDS->IsNewShape(n2)) {
222       continue;
223     }
224     aPK.SetIndices(n1, n2);
225     aMPK.Add(aPK);
226   }
227   //
228   // 1
229   BOPDS_VectorOfInterfVE& aVEs=myDS->InterfVE();
230   aNb=aVEs.Length();
231   for (i=0; i!=aNb; ++i) {
232     const BOPDS_InterfVE& aVE=aVEs(i);
233     aVE.Indices(n1, n2);
234     if (myDS->IsNewShape(n1) || myDS->IsNewShape(n2)) {
235       continue;
236     }
237     aPK.SetIndices(n1, n2);
238     aMPK.Add(aPK);
239   }
240   //
241   // 2
242   BOPDS_VectorOfInterfEE& aEEs=myDS->InterfEE();
243   aNb=aEEs.Length();
244   for (i=0; i!=aNb; ++i) {
245     const BOPDS_InterfEE& aEE=aEEs(i);
246     aEE.Indices(n1, n2);
247     if (myDS->IsNewShape(n1) || myDS->IsNewShape(n2)) {
248       continue;
249     }
250     aPK.SetIndices(n1, n2);
251     aMPK.Add(aPK);
252   }
253   //
254   // 3
255   BOPDS_VectorOfInterfVF& aVFs=myDS->InterfVF();
256   aNb=aVFs.Length();
257   for (i=0; i!=aNb; ++i) {
258     const BOPDS_InterfVF& aVF=aVFs(i);
259     aVF.Indices(n1, n2);
260     if (myDS->IsNewShape(n1) || myDS->IsNewShape(n2)) {
261       continue;
262     }
263     aPK.SetIndices(n1, n2);
264     aMPK.Add(aPK);
265   }
266   //
267   // 4
268   BOPDS_VectorOfInterfEF& aEFs=myDS->InterfEF();
269   aNb=aEFs.Length();
270   for (i=0; i!=aNb; ++i) {
271     const BOPDS_InterfEF& aEF=aEFs(i);
272     if (aEF.CommonPart().Type()==TopAbs_SHAPE) {
273       continue;
274     }
275     aEF.Indices(n1, n2);
276     if (myDS->IsNewShape(n1) || myDS->IsNewShape(n2)) {
277       continue;
278     }
279     aPK.SetIndices(n1, n2);
280     aMPK.Add(aPK);
281   }
282   //
283   // 5
284   BOPDS_VectorOfInterfFF& aFFs=myDS->InterfFF();
285   aNb=aFFs.Length();
286   for (i=0; i!=aNb; ++i) {
287     Standard_Boolean bTangentFaces, bFlag;
288     Standard_Integer aNbC, aNbP, j, iFound;
289     //
290     const BOPDS_InterfFF& aFF=aFFs(i);
291     aFF.Indices(n1, n2);
292     //
293     bTangentFaces=aFF.TangentFaces();
294     aNbP=aFF.Points().Length();
295     const BOPDS_VectorOfCurve& aVC=aFF.Curves();
296     aNbC=aVC.Length();
297     if (!aNbP && !aNbC && !bTangentFaces) {
298       continue;
299     }
300     //
301     iFound = (n1 == n2) ? 1 : 0;
302     //case of self-intersection inside one face
303     if (!iFound)
304     {
305       if (bTangentFaces) {
306         const TopoDS_Face& aF1=*((TopoDS_Face*)&myDS->Shape(n1));
307         const TopoDS_Face& aF2=*((TopoDS_Face*)&myDS->Shape(n2));
308         bFlag=BOPTools_AlgoTools::AreFacesSameDomain
309           (aF1, aF2, myContext, myFuzzyValue);
310         if (bFlag) {
311           ++iFound;
312         }
313       }
314       else {
315         for (j=0; j!=aNbC; ++j) {
316           const BOPDS_Curve& aNC=aVC(j);
317           const BOPDS_ListOfPaveBlock& aLPBC=aNC.PaveBlocks();
318           if (aLPBC.Extent()) {
319             ++iFound;
320             break;
321           }
322         }
323       }
324     }
325     //
326     if (!iFound) {
327       continue;
328     }
329     //
330     aPK.SetIndices(n1, n2);
331     aMPK.Add(aPK);
332   }
333   //
334   //
335   // 6
336   BOPDS_VectorOfInterfVZ& aVZs=myDS->InterfVZ();
337   aNb=aVZs.Length();
338   for (i=0; i!=aNb; ++i) {
339     //
340     const BOPDS_InterfVZ& aVZ=aVZs(i);
341     aVZ.Indices(n1, n2);
342     if (myDS->IsNewShape(n1) || myDS->IsNewShape(n2)) {
343       continue;
344     }
345     aPK.SetIndices(n1, n2);
346     aMPK.Add(aPK);
347   }
348   //
349   // 7
350   BOPDS_VectorOfInterfEZ& aEZs=myDS->InterfEZ();
351   aNb=aEZs.Length();
352   for (i=0; i!=aNb; ++i) {
353     //
354     const BOPDS_InterfEZ& aEZ=aEZs(i);
355     aEZ.Indices(n1, n2);
356     aPK.SetIndices(n1, n2);
357     aMPK.Add(aPK);
358   }
359   //
360   // 8
361   BOPDS_VectorOfInterfFZ& aFZs=myDS->InterfFZ();
362   aNb=aFZs.Length();
363   for (i=0; i!=aNb; ++i) {
364     //
365     const BOPDS_InterfFZ& aFZ=aFZs(i);
366     aFZ.Indices(n1, n2);
367     aPK.SetIndices(n1, n2);
368     aMPK.Add(aPK);
369   }
370   //
371   // 9
372   BOPDS_VectorOfInterfZZ& aZZs=myDS->InterfZZ();
373   aNb=aZZs.Length();
374   for (i=0; i!=aNb; ++i) {
375     //
376     const BOPDS_InterfZZ& aZZ=aZZs(i);
377     aZZ.Indices(n1, n2);
378     aPK.SetIndices(n1, n2);
379     aMPK.Add(aPK);
380   }
381 }
382
383 //=======================================================================
384 //function : CheckFaceSelfIntersection
385 //purpose  : 
386 //=======================================================================
387 void BOPAlgo_CheckerSI::CheckFaceSelfIntersection()
388 {
389   if (myLevelOfCheck < 5)
390     return;
391   
392   BOPDS_Pair aPK;
393
394   BOPDS_MapOfPair& aMPK=
395     *((BOPDS_MapOfPair*)&myDS->Interferences());
396   aMPK.Clear();
397   
398   BOPAlgo_VectorOfFaceSelfIntersect aVFace;
399   
400   Standard_Integer aNbS=myDS->NbSourceShapes();
401   
402   //
403   for (Standard_Integer i = 0; i < aNbS; i++)
404   {
405     const BOPDS_ShapeInfo& aSI = myDS->ShapeInfo(i);
406     if (aSI.ShapeType() != TopAbs_FACE)
407       continue;
408     //
409     const TopoDS_Face& aF = (*(TopoDS_Face*)(&aSI.Shape()));
410     BRepAdaptor_Surface BAsurf(aF, Standard_False);
411     GeomAbs_SurfaceType aSurfType = BAsurf.GetType();
412     if (aSurfType == GeomAbs_Plane ||
413         aSurfType == GeomAbs_Cylinder ||
414         aSurfType == GeomAbs_Cone ||
415         aSurfType == GeomAbs_Sphere)
416       continue;
417
418     if (aSurfType == GeomAbs_Torus)
419     {
420       gp_Torus aTorus = BAsurf.Torus();
421       Standard_Real aMajorRadius = aTorus.MajorRadius();
422       Standard_Real aMinorRadius = aTorus.MinorRadius();
423       if (aMajorRadius > aMinorRadius + Precision::Confusion())
424         continue;
425     }
426
427     Standard_Real aTolF = BRep_Tool::Tolerance(aF);
428     
429     BOPAlgo_FaceSelfIntersect& aFaceSelfIntersect = aVFace.Appended();
430     //
431     aFaceSelfIntersect.SetIndex(i);
432     aFaceSelfIntersect.SetFace(aF);
433     aFaceSelfIntersect.SetTolF(aTolF);
434     //
435     aFaceSelfIntersect.SetProgressIndicator(*myProgressScope);
436   }
437   
438   Standard_Integer aNbFace = aVFace.Length();
439   //======================================================
440   BOPTools_Parallel::Perform (myRunParallel, aVFace);
441   //======================================================
442   //
443   for (Standard_Integer k = 0; k < aNbFace; k++)
444   {
445     BOPAlgo_FaceSelfIntersect& aFaceSelfIntersect = aVFace(k);
446     //
447     Standard_Integer nF = aFaceSelfIntersect.IndexOfFace();
448
449     Standard_Boolean bIsDone = aFaceSelfIntersect.IsDone();
450     if (bIsDone)
451     {
452       const IntTools_SequenceOfCurves& aCvsX = aFaceSelfIntersect.Lines();
453       const IntTools_SequenceOfPntOn2Faces& aPntsX = aFaceSelfIntersect.Points();
454       //
455       Standard_Integer aNbCurves = aCvsX.Length();
456       Standard_Integer aNbPoints = aPntsX.Length();
457       //
458       if (aNbCurves || aNbPoints)
459       {
460         aPK.SetIndices(nF, nF);
461         aMPK.Add(aPK);
462       }
463     }
464   }
465 }