0023985: There is no section between attached faces.
[occt.git] / src / IntPatch / IntPatch_ImpImpIntersection_2.gxx
CommitLineData
b311480e 1// Created on: 1992-05-07
2// Created by: Jacques GOUSSARD
3// Copyright (c) 1992-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
c5c34473
J
22//=======================================================================
23//function : IntPatch_ImpImpIntersection
24//purpose :
25//=======================================================================
7fd59977 26IntPatch_ImpImpIntersection::IntPatch_ImpImpIntersection ():
27 done(Standard_False)
28{
29}
c5c34473
J
30//=======================================================================
31//function : IntPatch_ImpImpIntersection
32//purpose :
33//=======================================================================
7fd59977 34IntPatch_ImpImpIntersection::IntPatch_ImpImpIntersection
35 (const Handle(Adaptor3d_HSurface)& S1,
36 const Handle(Adaptor3d_TopolTool)& D1,
37 const Handle(Adaptor3d_HSurface)& S2,
38 const Handle(Adaptor3d_TopolTool)& D2,
39 const Standard_Real TolArc,
40 const Standard_Real TolTang)
41{
42 Perform(S1,D1,S2,D2,TolArc,TolTang);
43}
c5c34473
J
44//=======================================================================
45//function : Perform
46//purpose :
47//=======================================================================
7fd59977 48void IntPatch_ImpImpIntersection::Perform(const Handle(Adaptor3d_HSurface)& S1,
49 const Handle(Adaptor3d_TopolTool)& D1,
50 const Handle(Adaptor3d_HSurface)& S2,
51 const Handle(Adaptor3d_TopolTool)& D2,
52 const Standard_Real TolArc,
53 const Standard_Real TolTang) {
54 done = Standard_False;
55 spnt.Clear();
56 slin.Clear();
57
58 empt = Standard_True;
59 tgte = Standard_False;
60 oppo = Standard_False;
61
62 Standard_Boolean all1 = Standard_False;
63 Standard_Boolean all2 = Standard_False;
64 Standard_Boolean SameSurf = Standard_False;
65 Standard_Boolean multpoint = Standard_False;
66
67 Standard_Boolean nosolonS1 = Standard_False;
c5c34473 68 // indique s il y a des points sur restriction du carreau 1
7fd59977 69 Standard_Boolean nosolonS2 = Standard_False;
c5c34473 70 // indique s il y a des points sur restriction du carreau 2
7fd59977 71 Standard_Integer i, nbpt, nbseg;
7fd59977 72 IntPatch_SequenceOfSegmentOfTheSOnBounds edg1,edg2;
73 IntPatch_SequenceOfPathPointOfTheSOnBounds pnt1,pnt2;
c5c34473 74 //
7fd59977 75 // On commence par intersecter les supports des surfaces
7fd59977 76 IntSurf_Quadric quad1;
77 IntSurf_Quadric quad2;
78 IntPatch_ArcFunction AFunc;
79 Standard_Real Tolang = 1.e-8;
80 GeomAbs_SurfaceType typs1 = S1->GetType();
81 GeomAbs_SurfaceType typs2 = S2->GetType();
c5c34473 82 //
7fd59977 83 switch (typs1) {
84
85 case GeomAbs_Plane :
86 {
87 quad1.SetValue(S1->Plane());
88
89 switch (typs2) {
90
91 case GeomAbs_Plane:
92 {
93 quad2.SetValue(S2->Plane());
94 if (!IntPP(quad1,quad2,Tolang,TolTang,SameSurf,slin)) {
95 return;
96 }
97 }
98 break;
99
100
101 case GeomAbs_Cylinder:
102 {
103 quad2.SetValue(S2->Cylinder());
04cbc9d3 104 Standard_Real H = S1->LastVParameter() - S1->FirstVParameter();
105 if (!IntPCy(quad1,quad2,Tolang,TolTang,Standard_False,empt,slin,H)) {
7fd59977 106 return;
107 }
108 if (empt) {
109 done = Standard_True;
110 return;
111 }
112 }
113 break;
114
115 case GeomAbs_Sphere:
116 {
117 quad2.SetValue(S2->Sphere());
a9f7b6b5
P
118 //modified by NIZNHY-PKV Tue Sep 20 09:03:06 2011f
119 if (!IntPSp(quad1,quad2,Tolang,TolTang,Standard_False,empt,slin,spnt)) {
120 //if (!IntPSp(quad1,quad2,TolTang,Standard_False,empt,slin,spnt)) {
121 //modified by NIZNHY-PKV Tue Sep 20 09:03:10 2011t
7fd59977 122 return;
123 }
124 if (empt) {
125 done = Standard_True;
126 return;
127 }
128 }
129 break;
130
131 case GeomAbs_Cone:
132 {
133 quad2.SetValue(S2->Cone());
134 if (!IntPCo(quad1,quad2,Tolang,TolTang,Standard_False,
135 empt,multpoint,slin,spnt)) {
136 return;
137 }
138 if (empt) {
139 done = Standard_True;
140 return;
141 }
142 }
143 break;
144 default:
145 {
146 Standard_ConstructionError::Raise();
147 break;
148 }
149 }
150 }
151 break;
152
153 case GeomAbs_Cylinder:
154 {
155 quad1.SetValue(S1->Cylinder());
156 switch (typs2){
157
158 case GeomAbs_Plane:
159 {
160 quad2.SetValue(S2->Plane());
04cbc9d3 161 Standard_Real H = S1->LastVParameter() - S1->FirstVParameter();
162 if (!IntPCy(quad1,quad2,Tolang,TolTang,Standard_True,empt,slin,H)) {
7fd59977 163 return;
164 }
165 if (empt) {
166 done = Standard_True;
167 return;
168 }
169 }
170 break;
171
172 case GeomAbs_Cylinder:
173 {
174 quad2.SetValue(S2->Cylinder());
175 if (!IntCyCy(quad1,quad2,TolTang,empt,SameSurf,multpoint,slin,spnt)) {
176 return;
177 }
178 if (empt) {
179 done = Standard_True;
180 return;
181 }
182 }
183 break;
184
185 case GeomAbs_Sphere:
186 {
187 quad2.SetValue(S2->Sphere());
188 if (!IntCySp(quad1,quad2,TolTang,Standard_False,empt,multpoint,
189 slin,spnt)) {
190 return;
191 }
192 if (empt) {
193 done = Standard_True;
194
195 return;
196 }
197 }
198 break;
199
200 case GeomAbs_Cone:
201 {
202 quad2.SetValue(S2->Cone());
203 if (!IntCyCo(quad1,quad2,TolTang,Standard_False,empt,multpoint,
204 slin,spnt)) {
205 return;
206 }
207 if (empt) {
208 done = Standard_True;
209 return;
210 }
211 }
212 break;
213 default:
214 {
215 Standard_ConstructionError::Raise();
216 break;
217 }
218 }
219
220 }
221 break;
222
223 case GeomAbs_Sphere:
224 {
225 quad1.SetValue(S1->Sphere());
226
227 switch (typs2){
228
229 case GeomAbs_Plane:
230 {
231 quad2.SetValue(S2->Plane());
a9f7b6b5
P
232 //modified by NIZNHY-PKV Tue Sep 20 09:03:35 2011f
233 if (!IntPSp(quad1,quad2,Tolang,TolTang,Standard_True,empt,slin,spnt)) {
234 //if (!IntPSp(quad1,quad2,TolTang,Standard_True,empt,slin,spnt)) {
235 //modified by NIZNHY-PKV Tue Sep 20 09:03:38 2011t
7fd59977 236 return;
237 }
238 if (empt) {
239 done = Standard_True;
240 return;
241 }
242 }
243 break;
244
245 case GeomAbs_Cylinder:
246 {
247 quad2.SetValue(S2->Cylinder());
248 if (!IntCySp(quad1,quad2,TolTang,Standard_True,empt,multpoint,
249 slin,spnt)) {
250 return;
251 }
252 if (empt) {
253 done = Standard_True;
254 return;
255 }
256 }
257 break;
258
259 case GeomAbs_Sphere:
260 {
261 quad2.SetValue(S2->Sphere());
262 if (!IntSpSp(quad1,quad2,TolTang,empt,SameSurf,slin,spnt)) {
263 return;
264 }
265 if (empt) {
266 done = Standard_True;
267 return;
268 }
269 }
270 break;
271
272 case GeomAbs_Cone:
273 {
274 quad2.SetValue(S2->Cone());
275 if (!IntCoSp(quad1,quad2,TolTang,Standard_True,empt,multpoint,
276 slin,spnt)) {
277 return;
278 }
279 if (empt) {
280 done = Standard_True;
281 return;
282 }
283 }
284 break;
285 default:
286 {
287 Standard_ConstructionError::Raise();
288 break;
289 }
290 }
291
292 }
293 break;
294
295 case GeomAbs_Cone:
296 {
297 quad1.SetValue(S1->Cone());
298
299 switch (typs2){
300
301 case GeomAbs_Plane:
302 {
303 quad2.SetValue(S2->Plane());
304 if (!IntPCo(quad1,quad2,Tolang,TolTang,Standard_True,
305 empt,multpoint,slin,spnt)) {
306 return;
307 }
308 if (empt) {
309 done = Standard_True;
310 return;
311 }
312 }
313 break;
314
315 case GeomAbs_Cylinder:
316 {
317 quad2.SetValue(S2->Cylinder());
318 if (!IntCyCo(quad1,quad2,TolTang,Standard_True,empt,multpoint,
319 slin,spnt)) {
320 return;
321 }
322 if (empt) {
323 done = Standard_True;
324 return;
325 }
326 }
327 break;
328
329 case GeomAbs_Sphere:
330 {
331 quad2.SetValue(S2->Sphere());
332 if (!IntCoSp(quad1,quad2,TolTang,Standard_False,empt,multpoint,
333 slin,spnt)) {
334 return;
335 }
336 if (empt) {
337 done = Standard_True;
338 return;
339 }
340 }
341 break;
342
343 case GeomAbs_Cone:
344 {
345 quad2.SetValue(S2->Cone());
346 if (!IntCoCo(quad1,quad2,TolTang,empt,SameSurf,multpoint,
347 slin,spnt)) {
348 return;
349 }
350 if (empt) {
351 done = Standard_True;
352 return;
353 }
354 }
355 break;
356
357 default:
358 {
359 Standard_ConstructionError::Raise();
360 break;
361 }
362 }
363
364 }
365 break;
366 default:
367 {
368 Standard_ConstructionError::Raise();
369 break;
370 }
c5c34473
J
371 } //switch (typs1) {
372 //
7fd59977 373 if (!SameSurf) {
7fd59977 374 AFunc.SetQuadric(quad2);
375 AFunc.Set(S1);
376
377 solrst.Perform(AFunc, D1, TolArc, TolTang);
7fd59977 378 if (!solrst.IsDone()) {
379 return;
380 }
381
7fd59977 382 if (solrst.AllArcSolution() && typs1 == typs2) {
383 all1 = Standard_True;
384 }
385 nbpt = solrst.NbPoints();
386 nbseg= solrst.NbSegments();
387 for (i=1; i<= nbpt; i++) {
388 pnt1.Append(solrst.Point(i));
389 }
390 for (i=1; i<= nbseg; i++) {
391 edg1.Append(solrst.Segment(i));
392 }
393 nosolonS1 = (nbpt == 0) && (nbseg == 0);
394
395 if (nosolonS1 && all1) { // cas de face sans restrictions
396 all1 = Standard_False;
397 }
c5c34473 398 }//if (!SameSurf) {
7fd59977 399 else {
400 nosolonS1 = Standard_True;
401 }
402
403 if (!SameSurf) {
7fd59977 404 AFunc.SetQuadric(quad1);
405 AFunc.Set(S2);
406
407 solrst.Perform(AFunc, D2, TolArc, TolTang);
7fd59977 408 if (!solrst.IsDone()) {
409 return;
410 }
411
412 if (solrst.AllArcSolution() && typs1 == typs2) {
413 all2 = Standard_True;
414 }
415 nbpt = solrst.NbPoints();
416 nbseg= solrst.NbSegments();
417 for (i=1; i<= nbpt; i++) {
418 pnt2.Append(solrst.Point(i));
419 }
420
421 for (i=1; i<= nbseg; i++) {
422 edg2.Append(solrst.Segment(i));
423 }
424 nosolonS2 = (nbpt == 0) && (nbseg == 0);
425
426 if (nosolonS2 && all2) { // cas de face sans restrictions
427 all2 = Standard_False;
428 }
c5c34473 429 }// if (!SameSurf) {
7fd59977 430 else {
431 nosolonS2 = Standard_True;
432 }
c5c34473 433 //
7fd59977 434 if (SameSurf || (all1 && all2)) {
7fd59977 435 // faces "paralleles" parfaites
7fd59977 436 empt = Standard_False;
437 tgte = Standard_True;
438 slin.Clear();
439 spnt.Clear();
440
441 gp_Pnt Ptreference;
442
443 switch (typs1) {
c5c34473
J
444 case GeomAbs_Plane: {
445 Ptreference = (S1->Plane()).Location();
446 }
7fd59977 447 break;
c5c34473
J
448 case GeomAbs_Cylinder: {
449 Ptreference = ElSLib::Value(0.,0.,S1->Cylinder());
450 }
7fd59977 451 break;
c5c34473 452 case GeomAbs_Sphere: {
c6541a0c 453 Ptreference = ElSLib::Value(M_PI/4.,M_PI/4.,S1->Sphere());
c5c34473 454 }
7fd59977 455 break;
c5c34473
J
456 case GeomAbs_Cone: {
457 Ptreference = ElSLib::Value(0.,10.,S1->Cone());
458 }
7fd59977 459 break;
7fd59977 460 default:
461 break;
7fd59977 462 }
c5c34473 463 //
7fd59977 464 oppo = quad1.Normale(Ptreference).Dot(quad2.Normale(Ptreference)) < 0.0;
7fd59977 465 done = Standard_True;
466 return;
c5c34473 467 }// if (SameSurf || (all1 && all2)) {
7fd59977 468
469 if (!nosolonS1 || !nosolonS2) {
7fd59977 470 empt = Standard_False;
7fd59977 471 // C est la qu il faut commencer a bosser...
7fd59977 472 PutPointsOnLine(S1,S2,pnt1, slin, Standard_True, D1, quad1,quad2,
473 multpoint,TolArc);
474
475 PutPointsOnLine(S1,S2,pnt2, slin, Standard_False,D2, quad2,quad1,
476 multpoint,TolArc);
477
478 if (edg1.Length() != 0) {
479 ProcessSegments(edg1,slin,quad1,quad2,Standard_True,TolArc);
480 }
481
482 if (edg2.Length() != 0) {
483 ProcessSegments(edg2,slin,quad1,quad2,Standard_False,TolArc);
484 }
485
486 if (edg1.Length() !=0 || edg2.Length() !=0) {
c5c34473 487 // ProcessRLine(slin,S1,S2,TolArc);
7fd59977 488 ProcessRLine(slin,quad1,quad2,TolArc);
489 }
c5c34473 490 }//if (!nosolonS1 || !nosolonS2) {
7fd59977 491 else {
492 empt = ((slin.Length()==0) && (spnt.Length()==0));
493 }
c5c34473
J
494 //
495 Standard_Integer nblin, aNbPnt;
496 //
497 //modified by NIZNHY-PKV Tue Sep 06 10:03:35 2011f
498 aNbPnt=spnt.Length();
499 if (aNbPnt) {
500 IntPatch_SequenceOfPoint aSIP;
501 //
502 for(i=1; i<=aNbPnt; ++i) {
503 Standard_Real aU1, aV1, aU2, aV2;
504 gp_Pnt2d aP2D;
505 TopAbs_State aState1, aState2;
506 //
507 const IntPatch_Point& aIP=spnt(i);
508 aIP.Parameters(aU1, aV1, aU2, aV2);
509 //
510 aP2D.SetCoord(aU1, aV1);
511 aState1=D1->Classify(aP2D, TolArc);
512 //
513 aP2D.SetCoord(aU2, aV2);
514 aState2=D2->Classify(aP2D, TolArc);
515 //
516 if(aState1!=TopAbs_OUT && aState2!=TopAbs_OUT) {
517 aSIP.Append(aIP);
518 }
519 }
520 //
521 spnt.Clear();
522 //
523 aNbPnt=aSIP.Length();
524 for(i=1; i<=aNbPnt; ++i) {
525 const IntPatch_Point& aIP=aSIP(i);
526 spnt.Append(aIP);
527 }
528 //
529 }// if (aNbPnt) {
530 //modified by NIZNHY-PKV Tue Sep 06 10:18:20 2011t
531 //
532 nblin = slin.Length();
7fd59977 533 for(i=1; i<=nblin; i++) {
534 IntPatch_IType thetype = slin.Value(i)->ArcType();
535 if( (thetype == IntPatch_Ellipse)
536 ||(thetype == IntPatch_Circle)
537 ||(thetype == IntPatch_Lin)
538 ||(thetype == IntPatch_Parabola)
539 ||(thetype == IntPatch_Hyperbola)) {
540 Handle(IntPatch_GLine)& glin = *((Handle(IntPatch_GLine)*)&slin.Value(i));
7fd59977 541 glin->ComputeVertexParameters(TolArc);
542 }
543 else if(thetype == IntPatch_Analytic) {
544 Handle(IntPatch_ALine)& aligold = *((Handle(IntPatch_ALine)*)&slin.Value(i));
545 aligold->ComputeVertexParameters(TolArc);
546 }
547 else if(thetype == IntPatch_Restriction) {
548 Handle(IntPatch_RLine)& rlig = *((Handle(IntPatch_RLine)*)&slin.Value(i));
549 rlig->ComputeVertexParameters(TolArc);
550 }
551 }
c5c34473 552 //
7fd59977 553 //----------------------------------------------------------------
554 //-- On place 2 vertex sur les courbes de GLine qui n en
555 //-- contiennent pas.
7fd59977 556 for(i=1; i<=nblin; i++) {
557 gp_Pnt P;
558 IntPatch_Point point;
559 Standard_Real u1,v1,u2,v2;
7fd59977 560 Standard_Integer nbv;
561 if(slin.Value(i)->ArcType() == IntPatch_Circle) {
562 const Handle(IntPatch_GLine)& glin = *((Handle(IntPatch_GLine)*)&slin.Value(i));
563 nbv = glin->NbVertex();
564 if(glin->NbVertex() == 0) {
565 gp_Circ Circ = glin->Circle();
566 P=ElCLib::Value(0.0,Circ);
567 quad1.Parameters(P,u1,v1);
568 quad2.Parameters(P,u2,v2);
569 point.SetValue(P,TolArc,Standard_False);
570 point.SetParameters(u1,v1,u2,v2);
571 point.SetParameter(0.0);
572 glin->AddVertex(point);
573
574 P=ElCLib::Value(0.0,Circ);
575 quad1.Parameters(P,u1,v1);
576 quad2.Parameters(P,u2,v2);
577 point.SetValue(P,TolArc,Standard_False);
578 point.SetParameters(u1,v1,u2,v2);
c6541a0c 579 point.SetParameter(M_PI+M_PI);
7fd59977 580 glin->AddVertex(point);
581 }
582 }
583
584 else if(slin.Value(i)->ArcType() == IntPatch_Ellipse) {
585 const Handle(IntPatch_GLine)& glin = *((Handle(IntPatch_GLine)*)&slin.Value(i));
586 nbv = glin->NbVertex();
587 if(glin->NbVertex() == 0) {
588 gp_Elips Elips = glin->Ellipse();
589 P=ElCLib::Value(0.0,Elips);
590 quad1.Parameters(P,u1,v1);
591 quad2.Parameters(P,u2,v2);
592 point.SetValue(P,TolArc,Standard_False);
593 point.SetParameters(u1,v1,u2,v2);
594 point.SetParameter(0.0);
595 glin->AddVertex(point);
596
597 P=ElCLib::Value(0.0,Elips);
598 quad1.Parameters(P,u1,v1);
599 quad2.Parameters(P,u2,v2);
600 point.SetValue(P,TolArc,Standard_False);
601 point.SetParameters(u1,v1,u2,v2);
c6541a0c 602 point.SetParameter(M_PI+M_PI);
7fd59977 603 glin->AddVertex(point);
604 }
605 }
606 }
7fd59977 607 done = Standard_True;
608}
7fd59977 609