1 // Created on: 1992-04-06
2 // Created by: Jacques GOUSSARD
3 // Copyright (c) 1992-1999 Matra Datavision
4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
6 // This file is part of Open CASCADE Technology software library.
8 // This library is free software; you can redistribute it and/or modify it under
9 // the terms of the GNU Lesser General Public License version 2.1 as published
10 // by the Free Software Foundation, with special exception defined in the file
11 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
12 // distribution for complete text of the license and disclaimer of any warranty.
14 // Alternatively, this file may be used under the terms of Open CASCADE
15 // commercial license or contractual agreement.
20 #include <IntAna_Curve.hxx>
21 #include <IntPatch_ALine.hxx>
22 #include <IntPatch_Point.hxx>
23 #include <Standard_DomainError.hxx>
24 #include <Standard_OutOfRange.hxx>
25 #include <Standard_Type.hxx>
27 IMPLEMENT_STANDARD_RTTIEXT(IntPatch_ALine,IntPatch_Line)
31 IntPatch_ALine::IntPatch_ALine (const IntAna_Curve& C,
32 const Standard_Boolean Tang,
33 const IntSurf_TypeTrans Trans1,
34 const IntSurf_TypeTrans Trans2) :
35 IntPatch_Line(Tang,Trans1,Trans2),
36 fipt(Standard_False),lapt(Standard_False)
39 typ = IntPatch_Analytic;
44 IntPatch_ALine::IntPatch_ALine (const IntAna_Curve& C,
45 const Standard_Boolean Tang,
46 const IntSurf_Situation Situ1,
47 const IntSurf_Situation Situ2) :
48 IntPatch_Line(Tang,Situ1,Situ2),
49 fipt(Standard_False),lapt(Standard_False)
52 typ = IntPatch_Analytic;
57 IntPatch_ALine::IntPatch_ALine (const IntAna_Curve& C,
58 const Standard_Boolean Tang) :
60 fipt(Standard_False),lapt(Standard_False)
63 typ = IntPatch_Analytic;
67 const IntAna_Curve& IntPatch_ALine::Curve() const {
72 #define PCONFUSION 0.00001
75 void IntPatch_ALine::AddVertex (const IntPatch_Point& VTXj) {
77 Standard_Integer n = NbVertex();
79 Standard_Real par = VTXj.ParameterOnLine();
81 for(int i=1; i<=n ;i++) {
82 const IntPatch_Point& VTXi = svtx.Value(i);
83 if((VTXj.IsOnDomS1()==Standard_False) && (VTXj.IsOnDomS2()==Standard_False)) {
84 if((VTXi.IsOnDomS1()==Standard_False) && (VTXi.IsOnDomS2()==Standard_False)) {
85 if(Abs(par-VTXi.ParameterOnLine())<=PCONFUSION) {
87 cout<<" Rejet IntPatch_ALine::AddVertex (0) "<<endl;
95 const IntPatch_Point& VTXi = svtx.Value(i);
96 if((VTXj.IsOnDomS1()==Standard_True) && (VTXj.IsOnDomS2()==Standard_False)) {
97 if((VTXi.IsOnDomS1()==Standard_True) && (VTXi.IsOnDomS2()==Standard_False)) {
98 if(Abs(VTXi.ParameterOnArc1()-VTXj.ParameterOnArc1())<=PCONFUSION) {
99 if(VTXi.ArcOnS1() == VTXj.ArcOnS1()) {
101 cout<<" Rejet IntPatch_ALine::AddVertex (1) "<<endl;
109 for(i=1; i<=n ;i++) {
110 const IntPatch_Point& VTXi = svtx.Value(i);
111 if((VTXj.IsOnDomS2()==Standard_True) && (VTXj.IsOnDomS1()==Standard_False)) {
112 if((VTXi.IsOnDomS2()==Standard_True) && (VTXi.IsOnDomS1()==Standard_False)) {
113 if(Abs(VTXi.ParameterOnArc2()-VTXj.ParameterOnArc2())<=PCONFUSION) {
114 if(VTXi.ArcOnS2() == VTXj.ArcOnS2()) {
116 cout<<" Rejet IntPatch_ALine::AddVertex (2) "<<endl;
124 for(i=1; i<=n ;i++) {
125 const IntPatch_Point& VTXi = svtx.Value(i);
126 if((VTXj.IsOnDomS2()==Standard_True) && (VTXj.IsOnDomS1()==Standard_True)) {
127 if((VTXi.IsOnDomS2()==Standard_True) && (VTXi.IsOnDomS1()==Standard_True)) {
128 if(Abs(VTXi.ParameterOnArc2()-VTXj.ParameterOnArc2())<=PCONFUSION) {
129 if(Abs(VTXi.ParameterOnArc1()-VTXj.ParameterOnArc1())<=PCONFUSION) {
130 if(VTXi.ArcOnS2() == VTXj.ArcOnS2()) {
131 if(VTXi.ArcOnS1() == VTXj.ArcOnS1()) {
133 cout<<" Rejet IntPatch_ALine::AddVertex (3) "<<endl;
144 //-- Est ce que VTXj present sur 1 et 2 remplace un point VTXi present sur 1
145 for(i=1; i<=n ;i++) {
146 const IntPatch_Point& VTXi = svtx.Value(i);
147 if((VTXj.IsOnDomS2()==Standard_True) && (VTXj.IsOnDomS1()==Standard_True)) {
148 if((VTXi.IsOnDomS2()==Standard_False) && (VTXi.IsOnDomS1()==Standard_True)) {
149 Standard_Real p = Abs(VTXi.ParameterOnArc1()-VTXj.ParameterOnArc1());
151 cout<<" Est ce que VTXj present sur 1 et 2 remplace un point VTXi present sur 1 : "<<p<<endl;
154 if(VTXi.ArcOnS1() == VTXj.ArcOnS1()) {
156 cout<<" Replace IntPatch_ALine::AddVertex (1) "<<endl;
165 //-- Est ce que VTXj present sur 1 et 2 remplace un point VTXi present sur 2
166 for(i=1; i<=n ;i++) {
167 const IntPatch_Point& VTXi = svtx.Value(i);
168 if((VTXj.IsOnDomS2()==Standard_True) && (VTXj.IsOnDomS1()==Standard_True)) {
169 if((VTXi.IsOnDomS2()==Standard_True) && (VTXi.IsOnDomS1()==Standard_False)) {
170 Standard_Real p = Abs(VTXi.ParameterOnArc2()-VTXj.ParameterOnArc2());
172 cout<<" Est ce que VTXj present sur 1 et 2 remplace un point VTXi present sur 2 : "<<p<<endl;
175 if(VTXi.ArcOnS2() == VTXj.ArcOnS2()) {
177 cout<<" Replace IntPatch_ALine::AddVertex (2) "<<endl;
188 //-- Est ce que VTXi deja present sur 1 et 2 et un point VTXj present sur 1
189 for(i=1; i<=n ;i++) {
190 const IntPatch_Point& VTXi = svtx.Value(i);
191 if((VTXi.IsOnDomS2()==Standard_True) && (VTXi.IsOnDomS1()==Standard_True)) {
192 if((VTXj.IsOnDomS2()==Standard_False) && (VTXj.IsOnDomS1()==Standard_True)) {
193 Standard_Real p = Abs(VTXi.ParameterOnArc1()-VTXj.ParameterOnArc1());
195 if(VTXi.ArcOnS1() == VTXj.ArcOnS1()) {
197 cout<<" Replace IntPatch_ALine::AddVertex (1) -> RIEN "<<endl;
205 //-- Est ce que VTXj present sur 1 et 2 remplace un point VTXi present sur 2
206 for(i=1; i<=n ;i++) {
207 const IntPatch_Point& VTXi = svtx.Value(i);
208 if((VTXi.IsOnDomS2()==Standard_True) && (VTXi.IsOnDomS1()==Standard_True)) {
209 if((VTXj.IsOnDomS2()==Standard_True) && (VTXj.IsOnDomS1()==Standard_False)) {
210 Standard_Real p = Abs(VTXi.ParameterOnArc2()-VTXj.ParameterOnArc2());
212 if(VTXi.ArcOnS2() == VTXj.ArcOnS2()) {
234 void IntPatch_ALine::ComputeVertexParameters(const Standard_Real Tol) {
235 Standard_Boolean SortIsOK,APointDeleted;
236 Standard_Boolean SortAgain = Standard_True;
237 Standard_Integer nbvtx,i,j;
238 Standard_Real ParamMinOnLine,ParamMaxOnLine;
239 Standard_Boolean OpenFirst,OpenLast;
241 ParamMinOnLine = FirstParameter(OpenFirst);
242 ParamMaxOnLine = LastParameter(OpenLast);
245 //----------------------------------------------------------
246 //-- F i l t r e s u r r e s t r i c t i o n s --
247 //----------------------------------------------------------
248 //-- deux vertex sur la meme restriction et seulement
249 //-- sur celle ci ne doivent pas avoir le meme parametre
253 //-- Le tri est necessaire si suppression du first ou du last point
257 //-- On verifie qu un vertex a bien toute ses representations :
258 //-- Cas tres rare : point de tangence sur un debut de ligne
259 //-- et la ligne fait 2 * 2 PI de parametrage.
261 for(i=1; i<=nbvtx; i++) {
262 IntPatch_Point& VTX = svtx.ChangeValue(i);
263 Standard_Real p=VTX.ParameterOnLine();
264 Standard_Real pmpimpi=p-M_PI-M_PI;
265 if(pmpimpi >= ParamMinOnLine) {
266 gp_Pnt P1 = Value(pmpimpi);
267 Standard_Real d1 = P1.Distance(VTX.Value());
269 IntPatch_Point OVTX(VTX);
270 OVTX.SetParameter(pmpimpi);
275 if(pmpimpi <= ParamMaxOnLine) {
276 gp_Pnt P1 = Value(pmpimpi);
277 Standard_Real d1 = P1.Distance(VTX.Value());
279 IntPatch_Point OVTX(VTX);
280 OVTX.SetParameter(pmpimpi);
293 SortIsOK = Standard_True;
294 for(i=2; i<=nbvtx; i++) {
295 if(svtx.Value(i-1).ParameterOnLine() > svtx.Value(i).ParameterOnLine()) {
296 SortIsOK = Standard_False;
297 svtx.Exchange(i-1,i);
298 if(fipt) { if(indf==i) indf=i-1; else if(indf==i-1) indf=i; }
299 if(lapt) { if(indl==i) indl=i-1; else if(indl==i-1) indl=i; }
306 APointDeleted = Standard_False;
307 for(i=1; (i<=nbvtx) && (APointDeleted==Standard_False) ;i++) {
308 const IntPatch_Point& VTXi = svtx.Value(i);
309 if((VTXi.IsOnDomS1()==Standard_True) && (VTXi.IsOnDomS2()==Standard_False)) {
310 for(j=1; (j<=nbvtx) && (APointDeleted==Standard_False) ;j++) {
312 const IntPatch_Point& VTXj = svtx.Value(j);
313 if((VTXj.IsOnDomS1()==Standard_True) && (VTXj.IsOnDomS2()==Standard_False)) {
314 if(Abs(VTXi.ParameterOnArc1()-VTXj.ParameterOnArc1())<=PCONFUSION) {
315 if(Abs(VTXi.ParameterOnLine()-VTXj.ParameterOnLine())<=PCONFUSION) {
316 if(VTXi.ArcOnS1() == VTXj.ArcOnS1()) {
319 if(lapt) { if(indl>=j) indl--; }
320 if(fipt) { if(indf>=j) indf--; }
321 APointDeleted = Standard_True;
331 while(APointDeleted == Standard_True);
334 APointDeleted = Standard_False;
335 for(i=1; (i<=nbvtx) && (APointDeleted==Standard_False) ;i++) {
336 const IntPatch_Point& VTXi = svtx.Value(i);
337 if((VTXi.IsOnDomS2()==Standard_True) && (VTXi.IsOnDomS1()==Standard_False)) {
338 for(j=1; (j<=nbvtx) && (APointDeleted==Standard_False) ;j++) {
340 const IntPatch_Point& VTXj = svtx.Value(j);
341 if((VTXj.IsOnDomS2()==Standard_True) && (VTXj.IsOnDomS1()==Standard_False)) {
342 if(Abs(VTXi.ParameterOnArc2()-VTXj.ParameterOnArc2())<=PCONFUSION) {
343 if(Abs(VTXi.ParameterOnLine()-VTXj.ParameterOnLine())<=PCONFUSION) {
344 if(VTXi.ArcOnS2() == VTXj.ArcOnS2()) {
347 if(lapt) { if(indl>=j) indl--; }
348 if(fipt) { if(indf>=j) indf--; }
349 APointDeleted = Standard_True;
359 while(APointDeleted == Standard_True);
362 //----------------------------------------------------------
363 //-- Tri des vertex et suppression des Vtx superflus
369 SortIsOK = Standard_True;
370 for(i=2; i<=nbvtx; i++) {
371 if(svtx.Value(i-1).ParameterOnLine() > svtx.Value(i).ParameterOnLine()) {
372 SortIsOK = Standard_False;
373 svtx.Exchange(i-1,i);
374 if(fipt) { if(indf==i) indf=i-1; else if(indf==i-1) indf=i; }
375 if(lapt) { if(indl==i) indl=i-1; else if(indl==i-1) indl=i; }
381 SortAgain = Standard_False;
382 SortIsOK = Standard_True;
383 for(i=2; i<=nbvtx && SortIsOK; i++) {
384 IntPatch_Point& VTX = svtx.ChangeValue(i);
385 for(j=1; j<i && SortIsOK; j++) {
386 IntPatch_Point& VTXM1 = svtx.ChangeValue(j);
387 Standard_Boolean kill = Standard_False;
388 Standard_Boolean killm1 = Standard_False;
389 if(Abs(VTXM1.ParameterOnLine()-VTX.ParameterOnLine())<PCONFUSION) {
390 if(VTXM1.IsOnDomS1() && VTX.IsOnDomS1()) { //-- OnS1 OnS1
391 if(VTXM1.ArcOnS1() == VTX.ArcOnS1()) { //-- OnS1 == OnS1
392 if(VTXM1.IsOnDomS2()) { //-- OnS1 == OnS1 OnS2
393 if(VTX.IsOnDomS2()==Standard_False) { //-- OnS1 == OnS1 OnS2 PasOnS2
397 if(VTXM1.ArcOnS2() == VTX.ArcOnS2()) { //-- OnS1 == OnS1 OnS2 == OnS2
402 else { //-- OnS1 == OnS1 PasOnS2
403 if(VTX.IsOnDomS2()) { //-- OnS1 == OnS1 PasOnS2 OnS2
404 killm1=Standard_True;
409 else { //-- Pas OnS1 et OnS1
410 if(VTXM1.IsOnDomS2()==Standard_False && VTX.IsOnDomS2()==Standard_False) {
411 if(VTXM1.IsOnDomS1() && VTX.IsOnDomS1()==Standard_False) {
414 else if(VTX.IsOnDomS1() && VTXM1.IsOnDomS1()==Standard_False) {
415 killm1=Standard_True;
420 if(!(kill || killm1)) {
421 if(VTXM1.IsOnDomS2() && VTX.IsOnDomS2()) { //-- OnS2 OnS2
422 if(VTXM1.ArcOnS2() == VTX.ArcOnS2()) { //-- OnS2 == OnS2
423 if(VTXM1.IsOnDomS1()) { //-- OnS2 == OnS2 OnS1
424 if(VTX.IsOnDomS1()==Standard_False) { //-- OnS2 == OnS2 OnS1 PasOnS1
428 if(VTXM1.ArcOnS1() == VTX.ArcOnS1()) { //-- OnS2 == OnS2 OnS1 == OnS1
433 else { //-- OnS2 == OnS2 PasOnS1
434 if(VTX.IsOnDomS1()) { //-- OnS2 == OnS2 PasOnS1 OnS1
435 killm1=Standard_True;
440 else { //-- Pas OnS2 et OnS2
441 if(VTXM1.IsOnDomS1()==Standard_False && VTX.IsOnDomS1()==Standard_False) {
442 if(VTXM1.IsOnDomS2() && VTX.IsOnDomS2()==Standard_False) {
445 else if(VTX.IsOnDomS2() && VTXM1.IsOnDomS2()==Standard_False) {
446 killm1=Standard_True;
453 SortIsOK = Standard_False;
454 if(lapt) { if(indl>i) indl--; else if(indl==i) indl=j; }
455 if(fipt) { if(indf>i) indf--; else if(indf==i) indf=j; }
460 SortIsOK = Standard_False;
461 if(lapt) { if(indl>j) indl--; else if(indl==j) indl=i-1;}
462 if(fipt) { if(indf>j) indf--; else if(indf==j) indf=i-1;}
473 //----------------------------------------------------------
474 //-- Traitement des lignes periodiques --
475 //----------------------------------------------------------
476 if(OpenFirst == Standard_False && OpenLast == Standard_False) {
479 IntPatch_Point& VTX0 = svtx.ChangeValue(1);
480 IntPatch_Point& VTXN = svtx.ChangeValue(nbvtx);
481 if(VTX0.ParameterOnLine() == ParamMinOnLine) {
482 if(VTXN.ParameterOnLine() !=ParamMaxOnLine) {
483 gp_Pnt PN=Value(ParamMaxOnLine);
484 Standard_Real d = PN.Distance(VTX0.Value());
486 IntPatch_Point OVTX(VTX0);
487 OVTX.SetParameter(ParamMaxOnLine);
492 if(VTXN.ParameterOnLine() == ParamMaxOnLine) {
493 if(VTX0.ParameterOnLine() !=ParamMinOnLine) {
494 gp_Pnt P0=Value(ParamMinOnLine);
495 Standard_Real d = P0.Distance(VTX0.Value());
497 IntPatch_Point OVTX(VTXN);
498 OVTX.SetParameter(ParamMinOnLine);
506 //---------------------------------------------------------
507 //-- Faut il supprimer le premier et le dernier point
511 IntPatch_Point& VTX0 = svtx.ChangeValue(1);
512 if( (VTX0.IsOnDomS1() == Standard_False)
513 && (VTX0.IsOnDomS2() == Standard_False)) {
522 IntPatch_Point& VTX0 = svtx.ChangeValue(nbvtx);
523 if( (VTX0.IsOnDomS1() == Standard_False)
524 && (VTX0.IsOnDomS2() == Standard_False)) {
532 //-- Si 2 vertex ont le meme parametre on identifie le p3d
535 SortIsOK = Standard_True;
536 for(i=2; i<=nbvtx; i++) {
537 IntPatch_Point& VTX = svtx.ChangeValue(i);
538 IntPatch_Point& VTXm1 = svtx.ChangeValue(i-1);
539 if(Abs(VTX.ParameterOnLine()-VTXm1.ParameterOnLine())<PCONFUSION) {
540 if(VTX.IsOnDomS1() && VTXm1.IsOnDomS1()==Standard_False) {
541 VTXm1.SetArc(Standard_True,
543 VTX.ParameterOnArc1(),
544 VTX.TransitionLineArc1(),
545 VTX.TransitionOnS1());
547 else if(VTXm1.IsOnDomS1() && VTX.IsOnDomS1()==Standard_False) {
548 VTX.SetArc(Standard_True,
550 VTXm1.ParameterOnArc1(),
551 VTXm1.TransitionLineArc1(),
552 VTXm1.TransitionOnS1());
554 if(VTX.IsVertexOnS1() && VTXm1.IsVertexOnS1()==Standard_False) {
555 VTXm1.SetVertex(Standard_True, VTX.VertexOnS1());
556 VTXm1.SetArc(Standard_True,
558 VTX.ParameterOnArc1(),
559 VTX.TransitionLineArc1(),
560 VTX.TransitionOnS1());
562 else if(VTXm1.IsVertexOnS1() && VTX.IsVertexOnS1()==Standard_False) {
563 VTX.SetVertex(Standard_True, VTXm1.VertexOnS1());
564 VTX.SetArc(Standard_True,
566 VTXm1.ParameterOnArc1(),
567 VTXm1.TransitionLineArc1(),
568 VTXm1.TransitionOnS1());
571 if(VTX.IsOnDomS2() && VTXm1.IsOnDomS2()==Standard_False) {
572 VTXm1.SetArc(Standard_False,
574 VTX.ParameterOnArc2(),
575 VTX.TransitionLineArc2(),
576 VTX.TransitionOnS2());
578 else if(VTXm1.IsOnDomS2() && VTX.IsOnDomS2()==Standard_False) {
579 VTX.SetArc(Standard_False,
581 VTXm1.ParameterOnArc2(),
582 VTXm1.TransitionLineArc2(),
583 VTXm1.TransitionOnS2());
585 if(VTX.IsVertexOnS2() && VTXm1.IsVertexOnS2()==Standard_False) {
586 VTXm1.SetVertex(Standard_False, VTX.VertexOnS2());
587 VTXm1.SetArc(Standard_False,
589 VTX.ParameterOnArc2(),
590 VTX.TransitionLineArc2(),
591 VTX.TransitionOnS2());
593 else if(VTXm1.IsVertexOnS2() && VTX.IsVertexOnS2()==Standard_False) {
594 VTX.SetVertex(Standard_False, VTXm1.VertexOnS2());
595 VTX.SetArc(Standard_False,
597 VTXm1.ParameterOnArc2(),
598 VTXm1.TransitionLineArc2(),
599 VTXm1.TransitionOnS2());
602 if(VTX.Value().SquareDistance(VTXm1.Value()) > 1e-12) {
603 IntPatch_Point CopyVtx = VTXm1;
604 VTXm1.SetParameter(VTX.ParameterOnLine());
605 VTXm1.SetValue(VTX.Value(),VTX.Tolerance(),VTX.IsTangencyPoint());
606 Standard_Real u1,v1,u2,v2;
607 VTX.Parameters(u1,v1,u2,v2);
608 VTXm1.SetParameters(u1,v1,u2,v2);
609 if(CopyVtx.IsOnDomS1()) {
610 VTXm1.SetArc(Standard_True,
612 CopyVtx.ParameterOnArc1(),
613 CopyVtx.TransitionLineArc1(),
614 CopyVtx.TransitionOnS1());
616 if(CopyVtx.IsOnDomS2()) {
617 VTXm1.SetArc(Standard_False,
619 CopyVtx.ParameterOnArc2(),
620 CopyVtx.TransitionLineArc2(),
621 CopyVtx.TransitionOnS2());
623 if(CopyVtx.IsVertexOnS1()) {
624 VTXm1.SetVertex(Standard_True,CopyVtx.VertexOnS1());
625 VTXm1.SetArc(Standard_True,
627 CopyVtx.ParameterOnArc1(),
628 CopyVtx.TransitionLineArc1(),
629 CopyVtx.TransitionOnS1());
631 if(CopyVtx.IsVertexOnS2()) {
632 VTXm1.SetVertex(Standard_False,CopyVtx.VertexOnS2());
633 VTXm1.SetArc(Standard_False,
635 CopyVtx.ParameterOnArc2(),
636 CopyVtx.TransitionLineArc2(),
637 CopyVtx.TransitionOnS2());
641 SortIsOK=Standard_False;
642 //-- cout<<" IntPatch_ALine : ComputeVertexParameters : Ajust "<<endl;
649 /*nbvtx = NbVertex();
650 for(Standard_Integer opopo = 1; opopo<=nbvtx; opopo++) {
651 svtx.Value(opopo).Dump();