0028599: Replacement of old Boolean operations with new ones in BRepProj_Projection...
[occt.git] / src / IntPolyh / IntPolyh_Intersection_1.cxx
1 // Created on: 2005-09-26
2 // Created by: Igor FEOKTISTOV
3 // Copyright (c) 2005-2014 OPEN CASCADE SAS
4 //
5 // This file is part of Open CASCADE Technology software library.
6 //
7 // This library is free software; you can redistribute it and/or modify it under
8 // the terms of the GNU Lesser General Public License version 2.1 as published
9 // by the Free Software Foundation, with special exception defined in the file
10 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
11 // distribution for complete text of the license and disclaimer of any warranty.
12 //
13 // Alternatively, this file may be used under the terms of Open CASCADE
14 // commercial license or contractual agreement.
15
16
17 #include <Adaptor3d_HSurface.hxx>
18 #include <IntPolyh_ListOfCouples.hxx>
19 #include <IntPolyh_Couple.hxx>
20 #include <IntPolyh_Intersection.hxx>
21 #include <IntPolyh_MaillageAffinage.hxx>
22 #include <IntPolyh_PMaillageAffinage.hxx>
23 #include <IntPolyh_Triangle.hxx>
24
25 Standard_Integer MYPRINT1   = 0;
26
27
28 //=======================================================================
29 //function : IntPolyh_Intersection
30 //purpose  : 
31 //=======================================================================
32 IntPolyh_Intersection::IntPolyh_Intersection(const Handle(Adaptor3d_HSurface)& S1,
33                                              const TColStd_Array1OfReal& Upars1,
34                                              const TColStd_Array1OfReal& Vpars1,
35                                              const Handle(Adaptor3d_HSurface)& S2,
36                                              const TColStd_Array1OfReal& Upars2,
37                                              const TColStd_Array1OfReal& Vpars2)
38 {
39   myNbSU1 = Upars1.Length();
40   myNbSV1 = Vpars1.Length();
41   myNbSU2 = Upars2.Length(); 
42   myNbSV2 = Vpars2.Length(); 
43   mySurf1 = S1;
44   mySurf2 = S2;
45   done = Standard_False;
46   TSectionLines.Init(1000);
47   TTangentZones.Init(10000);
48   Perform(Upars1, Vpars1, Upars2, Vpars2);
49 }
50
51 //=======================================================================
52 //function : Perform
53 //purpose  : 
54 //=======================================================================
55
56 void IntPolyh_Intersection::Perform(const TColStd_Array1OfReal& Upars1,
57                                     const TColStd_Array1OfReal& Vpars1,
58                                     const TColStd_Array1OfReal& Upars2,
59                                     const TColStd_Array1OfReal& Vpars2) { 
60
61   done = Standard_True;
62
63   Standard_Boolean isStdDone = Standard_False;
64   Standard_Boolean isAdvDone = Standard_False;
65   Standard_Integer nbCouplesStd = 0;
66   Standard_Integer nbCouplesAdv = 0;
67   
68   
69   IntPolyh_PMaillageAffinage aPMaillageStd = 0;
70   IntPolyh_PMaillageAffinage aPMaillageFF = 0;
71   IntPolyh_PMaillageAffinage aPMaillageFR = 0;
72   IntPolyh_PMaillageAffinage aPMaillageRF = 0;
73   IntPolyh_PMaillageAffinage aPMaillageRR = 0;
74
75   isStdDone = PerformStd( Upars1, Vpars1, Upars2, Vpars2, 
76                           aPMaillageStd,nbCouplesStd);
77
78   // default interference done well, use it
79   if(isStdDone && nbCouplesStd > 10) {
80     aPMaillageStd->StartPointsChain(TSectionLines, TTangentZones);
81   }
82   // default interference done, but too few interferences foud;
83   // use advanced interference
84   else if(isStdDone && nbCouplesStd <= 10) {
85     isAdvDone = PerformAdv( Upars1, Vpars1, Upars2, Vpars2,
86                             aPMaillageFF,aPMaillageFR,aPMaillageRF,
87                             aPMaillageRR,nbCouplesAdv);
88       
89     // advanced interference found
90     if(isAdvDone && nbCouplesAdv > 0) {
91       aPMaillageFF->StartPointsChain(TSectionLines,TTangentZones);
92       aPMaillageFR->StartPointsChain(TSectionLines,TTangentZones);
93       aPMaillageRF->StartPointsChain(TSectionLines,TTangentZones);
94       aPMaillageRR->StartPointsChain(TSectionLines,TTangentZones);
95     }
96     else {
97       // use result of default
98       if(nbCouplesStd > 0)
99         aPMaillageStd->StartPointsChain(TSectionLines, TTangentZones);
100     }
101   }
102   // default interference faild, use advanced
103   else {
104 //       isAdvDone = PerformAdv(aPMaillageFF,aPMaillageFR,aPMaillageRF,aPMaillageRR,nbCouplesAdv);
105       
106 //       if(isAdvDone && nbCouplesAdv > 0) {cout << "4adv done, nbc: " << nbCouplesAdv << endl;
107 //      aPMaillageFF->StartPointsChain(TSectionLines,TTangentZones);
108 //      aPMaillageFR->StartPointsChain(TSectionLines,TTangentZones);
109 //      aPMaillageRF->StartPointsChain(TSectionLines,TTangentZones);
110 //      aPMaillageRR->StartPointsChain(TSectionLines,TTangentZones);
111 //       }
112   }
113
114   // accept result
115   nbsectionlines = TSectionLines.NbItems();
116   nbtangentzones = TTangentZones.NbItems();
117
118   // clean up
119   if(aPMaillageStd) delete aPMaillageStd;
120   if(aPMaillageFF) delete aPMaillageFF;
121   if(aPMaillageFR) delete aPMaillageFR;
122   if(aPMaillageRF) delete aPMaillageRF;
123   if(aPMaillageRR) delete aPMaillageRR;
124
125   // verify
126   if(!isStdDone && !isAdvDone)
127     done = Standard_False;
128 }
129
130 //=======================================================================
131 //function : PerformMaillage
132 //purpose  : Computes MaillageAffinage
133 //=======================================================================
134 Standard_Boolean IntPolyh_Intersection::PerformMaillage(const Standard_Boolean isFirstFwd,
135                                                         const Standard_Boolean isSecondFwd,
136                                                         const TColStd_Array1OfReal& Upars1,
137                                                         const TColStd_Array1OfReal& Vpars1,
138                                                         const TColStd_Array1OfReal& Upars2,
139                                                         const TColStd_Array1OfReal& Vpars2,
140                                                         IntPolyh_PMaillageAffinage &theMaillageS)
141 {
142   theMaillageS = new IntPolyh_MaillageAffinage(mySurf1, Upars1.Length(), Vpars1.Length(),
143                                                mySurf2, Upars2.Length(), Vpars2.Length(),
144                                                MYPRINT1);
145   
146   
147   theMaillageS->FillArrayOfPnt(1, isFirstFwd, Upars1, Vpars1);
148   theMaillageS->FillArrayOfPnt(2, isSecondFwd, Upars2, Vpars2);
149   
150   
151   
152   Standard_Real xx0,yy0,zz0,xx1,yy1,zz1;
153   theMaillageS->CommonBox(theMaillageS->GetBox(1), theMaillageS->GetBox(2),
154                           xx0, yy0, zz0, xx1, yy1, zz1);
155   
156   theMaillageS->FillArrayOfTriangles(1);
157   theMaillageS->FillArrayOfTriangles(2);
158
159   theMaillageS->FillArrayOfEdges(1);
160   theMaillageS->FillArrayOfEdges(2);
161
162   theMaillageS->TrianglesDeflectionsRefinementBSB();
163
164   Standard_Integer FinTTC = theMaillageS->TriangleCompare();
165
166   // if too many intersections, consider surfaces parallel (eap)
167   // test for parallel surf
168   if(FinTTC > 200) {
169     const Standard_Real eps = .996; //~ cos of 5deg.
170     IntPolyh_ListOfCouples& Couples = theMaillageS->GetCouples();
171
172     Standard_Integer npara = 0;
173     IntPolyh_ListIteratorOfListOfCouples aIt(Couples);
174     for(; aIt.More(); aIt.Next()) {
175       Standard_Real cosa = Abs(aIt.Value().Angle());
176       if(cosa > eps) ++npara;
177     }
178     if(npara >= theMaillageS->GetArrayOfTriangles(1).NbItems() ||
179        npara >= theMaillageS->GetArrayOfTriangles(2).NbItems() ) {
180       return Standard_False;
181     }
182   }
183
184   return Standard_True;
185 }
186
187 //=======================================================================
188 //function : PerformMaillage
189 //purpose  : Computes MaillageAffinage
190 //=======================================================================
191 Standard_Boolean IntPolyh_Intersection::PerformMaillage(const TColStd_Array1OfReal& Upars1,
192                                                         const TColStd_Array1OfReal& Vpars1,
193                                                         const TColStd_Array1OfReal& Upars2,
194                                                         const TColStd_Array1OfReal& Vpars2,
195                                                         IntPolyh_PMaillageAffinage &theMaillageS)
196 {
197     
198   theMaillageS = new IntPolyh_MaillageAffinage(mySurf1, Upars1.Length(), Vpars1.Length(),
199                                                mySurf2, Upars2.Length(), Vpars2.Length(),
200                                                MYPRINT1);
201   
202   theMaillageS->FillArrayOfPnt(1, Upars1, Vpars1);
203   theMaillageS->FillArrayOfPnt(2, Upars2, Vpars2);
204   
205   
206   Standard_Real xx0,yy0,zz0,xx1,yy1,zz1;
207   theMaillageS->CommonBox(theMaillageS->GetBox(1), theMaillageS->GetBox(2),
208                           xx0, yy0, zz0, xx1, yy1, zz1);
209
210   theMaillageS->FillArrayOfTriangles(1);
211   theMaillageS->FillArrayOfTriangles(2);
212
213   theMaillageS->FillArrayOfEdges(1);
214   theMaillageS->FillArrayOfEdges(2);
215
216   theMaillageS->TrianglesDeflectionsRefinementBSB();
217
218   Standard_Integer FinTTC = theMaillageS->TriangleCompare();
219
220   if( FinTTC == 0 ) {
221     Standard_Boolean myZone = Standard_True;
222     theMaillageS->SetEnlargeZone( myZone );
223     theMaillageS->FillArrayOfPnt(1);
224     theMaillageS->FillArrayOfPnt(2);
225     theMaillageS->CommonBox(theMaillageS->GetBox(1), theMaillageS->GetBox(2),
226                             xx0, yy0, zz0, xx1, yy1, zz1);
227     theMaillageS->FillArrayOfTriangles(1);
228     theMaillageS->FillArrayOfTriangles(2);
229     theMaillageS->FillArrayOfEdges(1);
230     theMaillageS->FillArrayOfEdges(2);
231     theMaillageS->TrianglesDeflectionsRefinementBSB();
232     FinTTC = theMaillageS->TriangleCompare();
233     myZone = Standard_False;
234     theMaillageS->SetEnlargeZone( myZone );
235   }
236
237   // if too many intersections, consider surfaces parallel (eap)
238 //IFV test for parallel surf
239   if(FinTTC > 200) {
240     const Standard_Real eps = .996; //~ cos of 5deg.
241     IntPolyh_ListOfCouples& Couples = theMaillageS->GetCouples();
242
243     Standard_Integer npara = 0;
244     IntPolyh_ListIteratorOfListOfCouples aIt(Couples);
245     for(; aIt.More(); aIt.Next()) {
246       Standard_Real cosa = Abs(aIt.Value().Angle());
247       if(cosa > eps) ++npara;
248     }
249     if(npara >= theMaillageS->GetArrayOfTriangles(1).NbItems() ||
250        npara >= theMaillageS->GetArrayOfTriangles(2).NbItems() ) {
251       return Standard_False;
252     }
253   }
254   
255   return Standard_True;
256 }
257
258 //=======================================================================
259 //function : PerformAdv
260 //purpose  : 
261 //=======================================================================
262 Standard_Boolean IntPolyh_Intersection::PerformAdv(const TColStd_Array1OfReal& Upars1,
263                                                    const TColStd_Array1OfReal& Vpars1,
264                                                    const TColStd_Array1OfReal& Upars2,
265                                                    const TColStd_Array1OfReal& Vpars2,
266                                                    IntPolyh_PMaillageAffinage& MaillageFF,
267                                                    IntPolyh_PMaillageAffinage& MaillageFR,
268                                                    IntPolyh_PMaillageAffinage& MaillageRF,
269                                                    IntPolyh_PMaillageAffinage& MaillageRR,
270                                                    Standard_Integer&           NbCouples)
271 {
272   Standard_Boolean isdone = Standard_True;
273   NbCouples = 0;
274
275   if(!PerformMaillage(Standard_True,Standard_False,
276                       Upars1, Vpars1, Upars2, Vpars2, 
277                       MaillageFR) ||
278      !PerformMaillage(Standard_False,Standard_True,
279                       Upars1, Vpars1, Upars2, Vpars2, 
280                       MaillageRF) ||
281      !PerformMaillage(Standard_True,Standard_True,
282                       Upars1, Vpars1, Upars2, Vpars2, 
283                       MaillageFF)  ||
284      !PerformMaillage(Standard_False,Standard_False,
285                       Upars1, Vpars1, Upars2, Vpars2, 
286                       MaillageRR) )
287     isdone = Standard_False; 
288
289   if(isdone) {
290     NbCouples = MaillageFF->GetCouples().Extent() +
291       MaillageFR->GetCouples().Extent() +
292         MaillageRF->GetCouples().Extent() +
293           MaillageRR->GetCouples().Extent();
294
295     if(NbCouples > 0)
296       MergeCouples(MaillageFF->GetCouples(),MaillageFR->GetCouples(),
297                    MaillageRF->GetCouples(),MaillageRR->GetCouples());
298   }
299   return isdone;
300 }
301 //=======================================================================
302 //function : PerformStd
303 //purpose  : 
304 //=======================================================================
305
306 Standard_Boolean IntPolyh_Intersection::PerformStd(const TColStd_Array1OfReal& Upars1,
307                                                    const TColStd_Array1OfReal& Vpars1,
308                                                    const TColStd_Array1OfReal& Upars2,
309                                                    const TColStd_Array1OfReal& Vpars2,
310                                                    IntPolyh_PMaillageAffinage& MaillageS,
311                                                    Standard_Integer&           NbCouples)
312 {
313   Standard_Boolean isdone = PerformMaillage(Upars1, Vpars1, Upars2, Vpars2, 
314                                             MaillageS);
315   NbCouples = (isdone) ? (MaillageS->GetCouples().Extent()) : 0;
316   return isdone;
317 }