src/BRepFill/BRepFill_TrimEdgeTool.cxx

   1 // Created on: 1995-04-24
   2 // Created by: Bruno DUMORTIER
   3 // Copyright (c) 1995-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
  21
  22
  23 #include <BRepFill_TrimEdgeTool.ixx>
  24 #include <BRep_Tool.hxx>
  25 #include <Bisector_BisecAna.hxx>
  26 #include <Geom2d_Curve.hxx>
  27 #include <Geom2d_TrimmedCurve.hxx>
  28 #include <Geom2d_CartesianPoint.hxx>
  29 #include <Geom_Curve.hxx>
  30 #include <GeomProjLib.hxx>
  31 #include <Geom_TrimmedCurve.hxx>
  32 #include <Geom2dAPI_ProjectPointOnCurve.hxx>
  33 #include <Geom2dInt_GInter.hxx>
  34 #include <gp_Pnt.hxx>
  35 #include <TopLoc_Location.hxx>
  36 #include <TopoDS.hxx>
  37 #include <Precision.hxx>
  38 #include <IntRes2d_IntersectionPoint.hxx>
  39 #include <IntRes2d_IntersectionSegment.hxx>
  40
  41 #include <StdFail_NotDone.hxx>
  42
  43 #ifdef DRAW
  44 #include <DrawTrSurf.hxx>
  45 #include <DBRep.hxx>
  46 static Standard_Boolean Affich       = Standard_False;
  47 #endif
  48
  49
  50 //=======================================================================
  51 //function : SimpleExpression
  52 //purpose  :
  53 //=======================================================================
  54
  55 static void SimpleExpression (const Bisector_Bisec&        B,
  56                                     Handle(Geom2d_Curve)&  Bis)
  57 {
  58   Bis = B.Value();
  59
  60   Handle(Standard_Type) BT = Bis->DynamicType();
  61   if (BT == STANDARD_TYPE(Geom2d_TrimmedCurve)) {
  62     Handle(Geom2d_TrimmedCurve) TrBis
  63       = Handle(Geom2d_TrimmedCurve)::DownCast(Bis);
  64     Handle(Geom2d_Curve) BasBis = TrBis->BasisCurve();
  65     BT = BasBis->DynamicType();
  66     if (BT == STANDARD_TYPE(Bisector_BisecAna)) {
  67       Bis = Handle(Bisector_BisecAna)::DownCast(BasBis)->Geom2dCurve();
  68       Bis = new Geom2d_TrimmedCurve (Bis,
  69                                      TrBis->FirstParameter(),
  70                                      TrBis->LastParameter());
  71     }
  72   }
  73 }
  74
  75
  76 //=======================================================================
  77 //function : BRepFill_TrimEdgeTool
  78 //purpose  :
  79 //=======================================================================
  80
  81 BRepFill_TrimEdgeTool::BRepFill_TrimEdgeTool()
  82 {
  83 }
  84
  85
  86 //=======================================================================
  87 //function : BRepFill_TrimEdgeTool
  88 //purpose  :
  89 //=======================================================================
  90
  91 BRepFill_TrimEdgeTool::BRepFill_TrimEdgeTool
  92 (const Bisector_Bisec& Bisec,
  93  const Handle(Geom2d_Geometry)& S1,
  94  const Handle(Geom2d_Geometry)& S2,
  95  const Standard_Real   Offset) :
  96 myOffset(Offset),
  97 myBisec(Bisec)
  98 {
  99   isPoint1 = (S1->DynamicType() == STANDARD_TYPE(Geom2d_CartesianPoint));
 100   isPoint2 = (S2->DynamicType() == STANDARD_TYPE(Geom2d_CartesianPoint));
 101
 102 // return geometries of shapes.
 103 //  Standard_Real f,l;
 104   if (isPoint1) {
 105     myP1 = Handle(Geom2d_Point)::DownCast(S1)->Pnt2d();
 106   }
 107   else {
 108     myC1 = Handle(Geom2d_Curve)::DownCast(S1);
 109 #ifdef DRAW
 110     if ( Affich) {
 111 //POP pour NT
 112       char* myC1name = "myC1";
 113       DrawTrSurf::Set(myC1name,myC1);
 114 //      DrawTrSurf::Set("myC1",myC1);
 115     }
 116 #endif
 117   }
 118   if (isPoint2) {
 119     myP2 = Handle(Geom2d_Point)::DownCast(S2)->Pnt2d();
 120   }
 121   else {
 122     myC2 = Handle(Geom2d_Curve)::DownCast(S2);
 123 #ifdef DRAW
 124     if ( Affich) {
 125       char* myC2name = "myC2";
 126       DrawTrSurf::Set(myC2name,myC2);
 127 //      DrawTrSurf::Set("myC2",myC2);
 128     }
 129 #endif
 130   }
 131   // return the simple expression of the bissectrice
 132   Handle(Geom2d_Curve) Bis;
 133   SimpleExpression(myBisec, Bis);
 134   myBis = Geom2dAdaptor_Curve(Bis);
 135 #ifdef DRAW
 136   if ( Affich) {
 137       char* myBisname = "myBis";
 138       DrawTrSurf::Set(myBisname,Bis);
 139   }
 140 #endif
 141
 142 }
 143
 144 //=======================================================================
 145 //function : Bubble
 146 //purpose  : Order the sequence of points by increasing x.
 147 //=======================================================================
 148
 149 static void Bubble(TColgp_SequenceOfPnt& Seq)
 150 {
 151   Standard_Boolean Invert = Standard_True;
 152   Standard_Integer NbPoints = Seq.Length();
 153   while (Invert) {
 154     Invert = Standard_False;
 155     for ( Standard_Integer i = 1; i < NbPoints; i++) {
 156       gp_Pnt P1 = Seq.Value(i);
 157       gp_Pnt P2 = Seq.Value(i+1);
 158       if (P2.X()<P1.X())  {
 159         Seq.Exchange(i,i+1);
 160         Invert = Standard_True;
 161       }
 162     }
 163   }
 164 }
 165
 166
 167 //=======================================================================
 168 //function : EvalParameters
 169 //purpose  :
 170 //=======================================================================
 171
 172 static void EvalParameters(const Geom2dAdaptor_Curve& Bis,
 173                            const Geom2dAdaptor_Curve& AC,
 174                                  TColgp_SequenceOfPnt& Params)
 175 {
 176   Geom2dInt_GInter Intersector;
 177   Standard_Real Tol = Precision::Confusion();
 178 //  Standard_Real TolC = 1.e-9;
 179
 180   Geom2dAdaptor_Curve CBis(Bis);
 181   Geom2dAdaptor_Curve CAC (AC);
 182
 183   //Intersector = Geom2dInt_GInter(CBis, CAC, TolC, Tol);
 184   Intersector = Geom2dInt_GInter(CAC, CBis, Tol, Tol);
 185
 186   Standard_Integer NbPoints, NbSegments;
 187   Standard_Real U1, U2;
 188   gp_Pnt P;
 189
 190   if ( !Intersector.IsDone()) {
 191     StdFail_NotDone::Raise("BRepFill_TrimSurfaceTool::IntersectWith");
 192   }
 193
 194   NbPoints = Intersector.NbPoints();
 195
 196   if (NbPoints > 0) {
 197     for ( Standard_Integer i = 1; i <= NbPoints; i++) {
 198       U1 = Intersector.Point(i).ParamOnSecond();
 199       U2 = Intersector.Point(i).ParamOnFirst();
 200       P = gp_Pnt(U1,U2,0.);
 201       Params.Append(P);
 202     }
 203
 204   }
 205
 206   NbSegments = Intersector.NbSegments();
 207
 208   if (NbSegments > 0) {
 209     IntRes2d_IntersectionSegment Seg;
 210     for ( Standard_Integer i = 1; i <= NbSegments; i++) {
 211       Seg = Intersector.Segment(i);
 212       U1  = Seg.FirstPoint().ParamOnSecond();
 213       Standard_Real Ulast = Seg.LastPoint().ParamOnSecond();
 214       if ( Abs(U1    - CBis.FirstParameter()) <= Tol &&
 215            Abs(Ulast - CBis.LastParameter())  <= Tol    ) {
 216         P = gp_Pnt(U1,Seg.FirstPoint().ParamOnFirst(),0.);
 217         Params.Append(P);
 218         P = gp_Pnt(Ulast,Seg.LastPoint().ParamOnFirst(),0.);
 219         Params.Append(P);
 220       }
 221       else {
 222         U1 += Seg.LastPoint().ParamOnSecond();
 223         U1 /= 2.;
 224         U2  = Seg.FirstPoint().ParamOnFirst();
 225         U2 += Seg.LastPoint().ParamOnFirst();
 226         U2 /= 2.;
 227         P = gp_Pnt(U1,U2,0.);
 228         Params.Append(P);
 229       }
 230     }
 231   }
 232
 233   // Order the sequence by growing parameter on the bissectrice.
 234   Bubble( Params);
 235 }
 236
 237 static void EvalParametersBis(const Geom2dAdaptor_Curve& Bis,
 238                               const Geom2dAdaptor_Curve& AC,
 239                               TColgp_SequenceOfPnt& Params,
 240                               const Standard_Real Tol)
 241 {
 242   Geom2dInt_GInter Intersector;
 243   Standard_Real TolC = Tol;
 244
 245   Geom2dAdaptor_Curve CBis(Bis);
 246   Geom2dAdaptor_Curve CAC (AC);
 247
 248   Intersector = Geom2dInt_GInter(CAC, CBis, TolC, Tol);
 249
 250   Standard_Integer NbPoints, NbSegments;
 251   Standard_Real U1, U2;
 252   gp_Pnt P;
 253
 254   if ( !Intersector.IsDone()) {
 255     StdFail_NotDone::Raise("BRepFill_TrimSurfaceTool::IntersectWith");
 256   }
 257
 258   NbPoints = Intersector.NbPoints();
 259
 260   if (NbPoints > 0) {
 261     for ( Standard_Integer i = 1; i <= NbPoints; i++) {
 262       U1 = Intersector.Point(i).ParamOnSecond();
 263       U2 = Intersector.Point(i).ParamOnFirst();
 264       P = gp_Pnt(U1,U2,0.);
 265       Params.Append(P);
 266     }
 267
 268   }
 269
 270   NbSegments = Intersector.NbSegments();
 271
 272   if (NbSegments > 0) {
 273     IntRes2d_IntersectionSegment Seg;
 274     for ( Standard_Integer i = 1; i <= NbSegments; i++) {
 275       Seg = Intersector.Segment(i);
 276       U1  = Seg.FirstPoint().ParamOnSecond();
 277       Standard_Real Ulast = Seg.LastPoint().ParamOnSecond();
 278       if ( Abs(U1    - CBis.FirstParameter()) <= Tol &&
 279            Abs(Ulast - CBis.LastParameter())  <= Tol    ) {
 280         P = gp_Pnt(U1,Seg.FirstPoint().ParamOnFirst(),0.);
 281         Params.Append(P);
 282         P = gp_Pnt(Ulast,Seg.LastPoint().ParamOnFirst(),0.);
 283         Params.Append(P);
 284       }
 285       else {
 286         U1 += Seg.LastPoint().ParamOnSecond();
 287         U1 /= 2.;
 288         U2  = Seg.FirstPoint().ParamOnFirst();
 289         U2 += Seg.LastPoint().ParamOnFirst();
 290         U2 /= 2.;
 291         P = gp_Pnt(U1,U2,0.);
 292         Params.Append(P);
 293       }
 294     }
 295   }
 296
 297   // Order the sequence by parameter growing on the bissectrice.
 298   Bubble( Params);
 299 }
 300
 301
 302 //=======================================================================
 303 //function : IntersectWith
 304 //purpose  :
 305 //=======================================================================
 306
 307 void BRepFill_TrimEdgeTool::IntersectWith(const TopoDS_Edge& Edge1,
 308                                           const TopoDS_Edge& Edge2,
 309                                                 TColgp_SequenceOfPnt& Params)
 310 {
 311   Params.Clear();
 312
 313   // return curves associated to edges.
 314   TopLoc_Location L;
 315   Standard_Real   f,l;
 316   Handle(Geom_Surface) Surf;
 317
 318   Handle(Geom2d_Curve) C1;
 319   BRep_Tool::CurveOnSurface(Edge1,C1,Surf,L,f,l);
 320   Geom2dAdaptor_Curve AC1(C1,f,l);
 321
 322   Handle(Geom2d_Curve) C2;
 323   BRep_Tool::CurveOnSurface(Edge2,C2,Surf,L,f,l);
 324   Geom2dAdaptor_Curve AC2(C2,f,l);
 325
 326 #ifdef DRAW
 327   if ( Affich) {
 328     f = AC1.FirstParameter();
 329     l = AC1.LastParameter();
 330     char* CURVE1name = "CURVE1";
 331     DrawTrSurf::Set(CURVE1name, new Geom2d_TrimmedCurve(C1,f,l));
 332     f = AC2.FirstParameter();
 333     l = AC2.LastParameter();
 334     char* CURVE2name = "CURVE2";
 335     DrawTrSurf::Set(CURVE2name, new Geom2d_TrimmedCurve(C2,f,l));
 336     f = myBis.FirstParameter();
 337     l = myBis.LastParameter();
 338     char* bisname = "BIS";
 339     DrawTrSurf::Set(bisname, new Geom2d_TrimmedCurve(myBis.Curve(),f,l));
 340     char* Edge1name = "E1";
 341     DBRep::Set(Edge1name, Edge1);
 342     char* Edge2name = "E2";
 343     DBRep::Set(Edge2name, Edge2);
 344
 345   }
 346 #endif
 347
 348   // Calculate intersection
 349   TColgp_SequenceOfPnt Points2;
 350   gp_Pnt PSeq;
 351
 352   EvalParameters (myBis,AC1,Params);
 353   EvalParameters (myBis,AC2,Points2);
 354
 355
 356   Standard_Integer SeanceDeRattrapage=0;
 357   Standard_Real TolInit= 1.e-9;
 358   Standard_Integer nn = 7;
 359
 360   if(AC1.GetType() != GeomAbs_Circle ||
 361      AC1.GetType() != GeomAbs_Line   ||
 362      AC2.GetType() != GeomAbs_Circle ||
 363      AC2.GetType() != GeomAbs_Line     ) {
 364
 365     TolInit = 1.e-8;
 366     nn = 6;
 367   }
 368
 369   while (     SeanceDeRattrapage < nn // TolInit <= 0.01
 370          && ( Points2.Length() != Params.Length() ||
 371              (Points2.Length() == 0 && Params.Length() == 0) ) ) {
 372
 373 #ifdef DEB
 374     cout << "BRepFill_TrimEdgeTool: incoherent intersection. Try with a greater tolerance" << endl;
 375 #endif
 376
 377     Params.Clear();
 378     Points2.Clear();
 379
 380     TolInit*=10.0;
 381
 382     EvalParametersBis(myBis,AC1,Params,TolInit);
 383     EvalParametersBis(myBis,AC2,Points2,TolInit);
 384     SeanceDeRattrapage++;
 385   }
 386
 387 #ifdef DEB
 388   if(SeanceDeRattrapage != 0) cout << "SeanceDeRattrapage = " << SeanceDeRattrapage << endl;
 389   if(SeanceDeRattrapage == nn) {
 390     cout << "BRepFill_TrimEdgeTool: incoherent intersection" << endl;
 391   }
 392 #endif
 393
 394
 395   if(Params.Length() == 0 && Points2.Length() == 1) {
 396
 397     //cout << "Params.Length() == 0 && Points2.Length() == 1" << endl;
 398     Standard_Real dmin;
 399     Standard_Real tBis = Points2(1).X();
 400     gp_Pnt2d PBis = myBis.Value(tBis);
 401
 402     Standard_Real t = AC1.FirstParameter();
 403     gp_Pnt2d PC = AC1.Value(t);
 404     dmin = PC.SquareDistance(PBis);
 405     gp_Pnt P(tBis, t, 0.);
 406     Params.Append(P);
 407
 408     t = AC1.LastParameter();
 409     PC = AC1.Value(t);
 410     if(dmin > PC.SquareDistance(PBis)) {
 411       P.SetY(t);
 412       Params.SetValue(1,P);
 413     }
 414   }
 415   else if(Params.Length() == 1 && Points2.Length() == 0) {
 416
 417     //cout << "Params.Length() == 1 && Points2.Length() == 0" << endl;
 418     Standard_Real dmin;
 419     Standard_Real tBis = Params(1).X();
 420     gp_Pnt2d PBis = myBis.Value(tBis);
 421
 422     Standard_Real t = AC2.FirstParameter();
 423     gp_Pnt2d PC = AC2.Value(t);
 424     dmin = PC.SquareDistance(PBis);
 425     gp_Pnt P(tBis, t, 0.);
 426     Points2.Append(P);
 427
 428     t = AC2.LastParameter();
 429     PC = AC2.Value(t);
 430     if(dmin > PC.SquareDistance(PBis)) {
 431       P.SetY(t);
 432       Points2.SetValue(1,P);
 433     }
 434   }
 435
 436   // small manipulation to remove incorrect intersections:
 437   // return only common intersections (same parameter
 438   // on the bissectrice.).
 439   // The tolerance can be eventually changed.
 440
 441   gp_Pnt P1,P2;
 442   Standard_Real Tol = 4 * 100 * Precision::PConfusion();
 443   Standard_Integer i = 1;
 444   Standard_Integer NbPoints = Params.Length();
 445
 446   if(NbPoints == 1 && Points2.Length() == 1) {
 447     //cout << "NbPoints == 1 && Points2.Length() == 1" << endl;
 448     for ( i = 1; i <= NbPoints; i++) {
 449       PSeq = Params(i);
 450       PSeq.SetZ((Points2.Value(i)).Y());
 451       Params.SetValue(i,PSeq);
 452     }
 453     return;
 454   }
 455
 456   i = 1;
 457   while ( i <= Min( Params.Length(), Points2.Length())) {
 458     P1 = Params(i);
 459     P2 = Points2(i);
 460     Standard_Real P1xP2x=Abs( P1.X() - P2.X());
 461
 462     if ( P1xP2x > Tol ) {
 463 #ifdef DEB
 464       cout << "BRepFill_TrimEdgeTool: no same parameter on the bissectrice" << endl;
 465 #endif
 466       if(P1xP2x>TolInit) {
 467 #ifdef DEB
 468       cout << "BRepFill_TrimEdgeTool: Continue somehow" << endl;
 469 #endif
 470       i++;
 471       }
 472       else {
 473         if ( P1.X() < P2.X()) Params.Remove(i);
 474         else                  Points2.Remove(i);
 475       }
 476     }
 477     else i++;
 478   }
 479
 480   if ( Params.Length() > Points2.Length()) {
 481     Params.Remove(Points2.Length()+1, Params.Length());
 482   }
 483   else if ( Params.Length() < Points2.Length()) {
 484     Points2.Remove(Params.Length()+1, Points2.Length());
 485   }
 486
 487   NbPoints = Params.Length();
 488   for ( i = 1; i <= NbPoints; i++) {
 489     PSeq = Params(i);
 490     PSeq.SetZ((Points2.Value(i)).Y());
 491     Params.SetValue(i,PSeq);
 492   }
 493 }
 494
 495 //=======================================================================
 496 //function : AddOrConfuse
 497 //purpose  : the first or the last point of the bissectrice is on the
 498 //           parallel if it was not found in the intersections,
 499 //           it is projected on parallel lines and added in the parameters
 500 //=======================================================================
 501
 502 void BRepFill_TrimEdgeTool::AddOrConfuse(const Standard_Boolean  Start,
 503                                          const TopoDS_Edge&      Edge1,
 504                                          const TopoDS_Edge&      Edge2,
 505                                          TColgp_SequenceOfPnt&   Params)
 506 const
 507 {
 508   Standard_Boolean  ToProj = Standard_True;
 509   gp_Pnt2d          PBis;
 510   Standard_Real     Tol = 10*Precision::Confusion();
 511
 512   // return curves associated to edges.
 513   TopLoc_Location L;
 514   Standard_Real   f,l;
 515   Handle(Geom_Surface) Surf;
 516
 517   Handle(Geom2d_Curve) C1;
 518   BRep_Tool::CurveOnSurface(Edge1,C1,Surf,L,f,l);
 519   Geom2dAdaptor_Curve AC1(C1,f,l);
 520
 521
 522   if (Start) PBis = myBis.Value(myBis.FirstParameter());
 523   else       PBis = myBis.Value(myBis.LastParameter ());
 524
 525   // Test if the end of the bissectrice is in the set of intersection points.
 526   if (!Params.IsEmpty()) {
 527     gp_Pnt2d P;
 528     if (Start) P = AC1.Value(Params.First().Y());
 529     else       P = AC1.Value(Params.Last ().Y());
 530     ToProj     = !PBis.IsEqual(P,Tol);
 531   }
 532
 533   if (ToProj) {
 534 #ifdef DEB
 535     cout << " project extremity bissectrice on parallel."<<endl;
 536 #endif
 537
 538     // Project point on parallels and add in Params
 539
 540     Standard_Real f2,l2;
 541     Handle(Geom2d_Curve) C2;
 542     BRep_Tool::CurveOnSurface(Edge2,C2,Surf,L,f2,l2);
 543
 544     Geom2dAPI_ProjectPointOnCurve Projector1(PBis,C1,f,l);
 545     Geom2dAPI_ProjectPointOnCurve Projector2(PBis,C2,f2,l2);
 546
 547     if (Projector1.NbPoints() == 0) {
 548 #ifdef DEB
 549       cout << "Failed projection in BRepFill_TrimEdgeTool::AddOrConfuse"<<endl;
 550 #endif
 551       return;
 552     }
 553     if (!Projector1.NearestPoint().IsEqual(PBis,Tol)) {
 554 #ifdef DEB
 555       cout <<"Incorrect solution in BRepFill_TrimEdgeTool::AddOrConfuse"<<endl;
 556 #endif
 557       return;
 558     }
 559     if (Projector2.NbPoints() == 0) {
 560 #ifdef DEB
 561       cout << "Failed projection in BRepFill_TrimEdgeTool::AddOrConfuse"<<endl;
 562 #endif
 563       return;
 564     }
 565     if (!Projector2.NearestPoint().IsEqual(PBis,Tol)) {
 566 #ifdef DEB
 567       cout <<" Mauvaisesolution dans BRepFill_TrimEdgeTool::AddOrConfuse"<<endl;
 568 #endif
 569       return;
 570     }
 571     gp_Pnt PInt (0,
 572                  Projector1.LowerDistanceParameter(),
 573                  Projector2.LowerDistanceParameter());
 574     if (Start) {
 575       PInt.SetX (myBis.FirstParameter());
 576       Params.Prepend(PInt);
 577     }
 578     else {
 579       PInt.SetX (myBis.LastParameter());
 580       Params.Append(PInt);
 581     }
 582   }
 583 }
 584
 585 //=======================================================================
 586 //function : IsInside
 587 //purpose  :
 588 //=======================================================================
 589
 590 Standard_Boolean BRepFill_TrimEdgeTool::IsInside(const gp_Pnt2d& P) const
 591 {
 592 //  Modified by Sergey KHROMOV - Fri Sep 27 11:43:12 2002 Begin
 593 //   Standard_Real Dist;
 594   Standard_Real Dist = RealLast();
 595 //  Modified by Sergey KHROMOV - Fri Sep 27 11:43:12 2002 End
 596   if (isPoint1)
 597     Dist = P.Distance(myP1);
 598   else if (isPoint2)
 599     Dist = P.Distance(myP2);
 600   else {
 601     Geom2dAPI_ProjectPointOnCurve Projector(P,myC1);
 602     if (Projector.NbPoints() > 0) {
 603       Dist = Projector.LowerDistance();
 604     }
 605 //  Modified by Sergey KHROMOV - Fri Sep 27 11:43:43 2002 Begin
 606 //     else {
 607 //       gp_Pnt2d PF = myC1->Value(myC1->FirstParameter());
 608 //       gp_Pnt2d PL = myC1->Value(myC1->LastParameter());
 609 //       Dist = Min (P.Distance(PF),P.Distance(PL));
 610 //     }
 611
 612 // Check of distances between P and first and last point of the first curve
 613 // should be performed in any case, despite of the results of projection.
 614     gp_Pnt2d      PF       = myC1->Value(myC1->FirstParameter());
 615     gp_Pnt2d      PL       = myC1->Value(myC1->LastParameter());
 616     Standard_Real aDistMin = Min (P.Distance(PF),P.Distance(PL));
 617
 618     if (Dist > aDistMin)
 619       Dist = aDistMin;
 620 //  Modified by Sergey KHROMOV - Fri Sep 27 11:43:44 2002 End
 621   }
 622
 623 //  return (Dist < Abs(myOffset);
 624 // return (Dist < Abs(myOffset) + Precision::Confusion());
 625   return (Dist < Abs(myOffset) - Precision::Confusion());
 626 }
 627