src/Adaptor3d/Adaptor3d_IsoCurve.cxx

   1 // Copyright (c) 1999-2014 OPEN CASCADE SAS
   2 //
   3 // This file is part of Open CASCADE Technology software library.
   4 //
   5 // This library is free software; you can redistribute it and/or modify it under
   6 // the terms of the GNU Lesser General Public License version 2.1 as published
   7 // by the Free Software Foundation, with special exception defined in the file
   8 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
   9 // distribution for complete text of the license and disclaimer of any warranty.
  10 //
  11 // Alternatively, this file may be used under the terms of Open CASCADE
  12 // commercial license or contractual agreement.
  13
  14
  15 #include <Adaptor3d_HCurve.hxx>
  16 #include <Adaptor3d_HIsoCurve.hxx>
  17 #include <Adaptor3d_HSurface.hxx>
  18 #include <Adaptor3d_IsoCurve.hxx>
  19 #include <BSplCLib.hxx>
  20 #include <BSplSLib.hxx>
  21 #include <ElCLib.hxx>
  22 #include <ElSLib.hxx>
  23 #include <Geom_BezierCurve.hxx>
  24 #include <Geom_BezierSurface.hxx>
  25 #include <Geom_BSplineCurve.hxx>
  26 #include <Geom_BSplineSurface.hxx>
  27 #include <GeomAbs_SurfaceType.hxx>
  28 #include <gp_Ax2.hxx>
  29 #include <gp_Circ.hxx>
  30 #include <gp_Elips.hxx>
  31 #include <gp_Hypr.hxx>
  32 #include <gp_Lin.hxx>
  33 #include <gp_Parab.hxx>
  34 #include <gp_Pnt.hxx>
  35 #include <gp_Vec.hxx>
  36 #include <Precision.hxx>
  37 #include <Standard_DomainError.hxx>
  38 #include <Standard_NoSuchObject.hxx>
  39 #include <Standard_NotImplemented.hxx>
  40 #include <Standard_OutOfRange.hxx>
  41 #include <TColgp_Array1OfPnt.hxx>
  42 #include <TColgp_Array2OfPnt.hxx>
  43
  44 //=======================================================================
  45 //function : Adaptor3d_IsoCurve
  46 //purpose  :
  47 //=======================================================================
  48 Adaptor3d_IsoCurve::Adaptor3d_IsoCurve()
  49 : myIso      (GeomAbs_NoneIso),
  50   myFirst    (0.0),
  51   myLast     (0.0),
  52   myParameter(0.0)
  53 {
  54 }
  55
  56 //=======================================================================
  57 //function : Adaptor3d_IsoCurve
  58 //purpose  :
  59 //=======================================================================
  60
  61 Adaptor3d_IsoCurve::Adaptor3d_IsoCurve(const Handle(Adaptor3d_HSurface)& S)
  62 : mySurface  (S),
  63   myIso      (GeomAbs_NoneIso),
  64   myFirst    (0.0),
  65   myLast     (0.0),
  66   myParameter(0.0)
  67 {
  68 }
  69
  70 //=======================================================================
  71 //function : Adaptor3d_IsoCurve
  72 //purpose  :
  73 //=======================================================================
  74
  75 Adaptor3d_IsoCurve::Adaptor3d_IsoCurve(const Handle(Adaptor3d_HSurface)& S,
  76                                        const GeomAbs_IsoType theIso,
  77                                        const Standard_Real theParam)
  78 : mySurface  (S),
  79   myIso      (GeomAbs_NoneIso),
  80   myFirst    (0.0),
  81   myLast     (0.0),
  82   myParameter(0.0)
  83 {
  84   Load(theIso, theParam);
  85 }
  86
  87 //=======================================================================
  88 //function : Adaptor3d_IsoCurve
  89 //purpose  :
  90 //=======================================================================
  91
  92 Adaptor3d_IsoCurve::Adaptor3d_IsoCurve(const Handle(Adaptor3d_HSurface)& theS,
  93                                        const GeomAbs_IsoType theIso,
  94                                        const Standard_Real theParam,
  95                                        const Standard_Real theWFirst,
  96                                        const Standard_Real theWLast)
  97 : mySurface  (theS),
  98   myIso      (theIso),
  99   myFirst    (theWFirst),
 100   myLast     (theWLast),
 101   myParameter(theParam)
 102 {
 103   Load(theIso, theParam, theWFirst, theWLast);
 104 }
 105
 106 //=======================================================================
 107 //function : Load
 108 //purpose  :
 109 //=======================================================================
 110
 111 void Adaptor3d_IsoCurve::Load(const Handle(Adaptor3d_HSurface)& S )
 112 {
 113   mySurface = S;
 114   myIso = GeomAbs_NoneIso;
 115 }
 116
 117 //=======================================================================
 118 //function : Load
 119 //purpose  :
 120 //=======================================================================
 121
 122 void Adaptor3d_IsoCurve::Load(const GeomAbs_IsoType Iso,
 123                             const Standard_Real Param)
 124 {
 125   switch (Iso) {
 126
 127   case GeomAbs_IsoU:
 128     Load(Iso,Param,
 129          mySurface->FirstVParameter(),
 130          mySurface->LastVParameter());
 131     break;
 132
 133   case GeomAbs_IsoV:
 134     Load(Iso,Param,
 135          mySurface->FirstUParameter(),
 136          mySurface->LastUParameter());
 137     break;
 138
 139   case GeomAbs_NoneIso:
 140     Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 141     break;
 142   }
 143 }
 144
 145 //=======================================================================
 146 //function : Load
 147 //purpose  :
 148 //=======================================================================
 149
 150 void Adaptor3d_IsoCurve::Load(const GeomAbs_IsoType Iso,
 151                             const Standard_Real Param,
 152                             const Standard_Real WFirst,
 153                             const Standard_Real WLast)
 154 {
 155   myIso =Iso;
 156   myParameter = Param;
 157   myFirst = WFirst;
 158   myLast = WLast;
 159
 160
 161   if (myIso == GeomAbs_IsoU) {
 162     myFirst = Max(myFirst, mySurface->FirstVParameter());
 163     myLast  = Min(myLast,  mySurface->LastVParameter());
 164   }
 165   else {
 166     myFirst = Max(myFirst, mySurface->FirstUParameter());
 167     myLast  = Min(myLast,  mySurface->LastUParameter());
 168   }
 169
 170 // Adjust the parameters on periodic surfaces
 171
 172   Standard_Real dummy = myParameter;
 173
 174   if (mySurface->IsUPeriodic()) {
 175
 176     if (myIso == GeomAbs_IsoU) {
 177       ElCLib::AdjustPeriodic
 178         (mySurface->FirstUParameter(),
 179          mySurface->FirstUParameter()+
 180          mySurface->UPeriod(),
 181          mySurface->UResolution(Precision::Confusion()),
 182          myParameter,dummy);
 183     }
 184     else {
 185       ElCLib::AdjustPeriodic
 186         (mySurface->FirstUParameter(),
 187          mySurface->FirstUParameter()+
 188          mySurface->UPeriod(),
 189          mySurface->UResolution(Precision::Confusion()),
 190          myFirst,myLast);
 191     }
 192   }
 193
 194   if (mySurface->IsVPeriodic()) {
 195
 196     if (myIso == GeomAbs_IsoV) {
 197       ElCLib::AdjustPeriodic
 198         (mySurface->FirstVParameter(),
 199          mySurface->FirstVParameter() +
 200          mySurface->VPeriod(),
 201          mySurface->VResolution(Precision::Confusion()),
 202          myParameter,dummy);
 203     }
 204     else {
 205       ElCLib::AdjustPeriodic
 206         (mySurface->FirstVParameter(),
 207          mySurface->FirstVParameter() +
 208          mySurface->VPeriod(),
 209          mySurface->VResolution(Precision::Confusion()),
 210          myFirst,myLast);
 211     }
 212   }
 213
 214 }
 215
 216 //=======================================================================
 217 //function : Continuity
 218 //purpose  :
 219 //=======================================================================
 220
 221 GeomAbs_Shape Adaptor3d_IsoCurve::Continuity() const
 222 {
 223   switch (myIso) {
 224   case GeomAbs_IsoU:
 225     return mySurface->VContinuity();
 226   case GeomAbs_IsoV:
 227     return mySurface->UContinuity();
 228   case GeomAbs_NoneIso:
 229   default:
 230     break;
 231   }
 232
 233   Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 234   return GeomAbs_C0;
 235 }
 236
 237 //=======================================================================
 238 //function : NbIntervals
 239 //purpose  :
 240 //=======================================================================
 241
 242 Standard_Integer Adaptor3d_IsoCurve::NbIntervals(const GeomAbs_Shape S) const
 243 {
 244   if (myIso == GeomAbs_NoneIso) Standard_NoSuchObject::Raise();
 245   Standard_Boolean UIso =  (myIso == GeomAbs_IsoU);
 246
 247   Standard_Integer nbInter = UIso ?
 248       mySurface->NbVIntervals(S) :
 249       mySurface->NbUIntervals(S);
 250
 251   TColStd_Array1OfReal T(1,nbInter+1);
 252
 253   if (UIso)
 254     mySurface->VIntervals(T,S);
 255   else
 256     mySurface->UIntervals(T,S);
 257
 258   if(nbInter == 1) return nbInter;
 259
 260   Standard_Integer first = 1;
 261   while (T(first) <= myFirst) first++;
 262   Standard_Integer last = nbInter+1;
 263   while (T(last) >= myLast) last--;
 264   return (last - first + 2);
 265 }
 266
 267 //=======================================================================
 268 //function : Intervals
 269 //purpose  :
 270 //=======================================================================
 271
 272 void Adaptor3d_IsoCurve::Intervals(TColStd_Array1OfReal& TI,
 273                                    const GeomAbs_Shape S) const
 274 {
 275   if (myIso == GeomAbs_NoneIso) Standard_NoSuchObject::Raise();
 276   Standard_Boolean UIso =  (myIso == GeomAbs_IsoU);
 277
 278   Standard_Integer nbInter = UIso ?
 279       mySurface->NbVIntervals(S) :
 280       mySurface->NbUIntervals(S);
 281
 282   TColStd_Array1OfReal T(1,nbInter+1);
 283
 284   if (UIso)
 285     mySurface->VIntervals(T,S);
 286   else
 287     mySurface->UIntervals(T,S);
 288
 289   if(nbInter == 1) {
 290     TI(TI.Lower()) = myFirst ;
 291     TI(TI.Lower() + 1) = myLast ;
 292     return;
 293   }
 294
 295   Standard_Integer first = 1;
 296   while (T(first) <= myFirst) first++;
 297   Standard_Integer last = nbInter+1;
 298   while (T(last) >= myLast) last--;
 299
 300   Standard_Integer i = TI.Lower(), j;
 301   for (j = first-1; j <= last+1; j++) {
 302     TI(i) = T(j);
 303     i++;
 304   }
 305   TI(TI.Lower()) = myFirst ;
 306   TI(TI.Lower() + last-first + 2) = myLast ;
 307 }
 308
 309 //=======================================================================
 310 //function : Trim
 311 //purpose  :
 312 //=======================================================================
 313
 314 Handle(Adaptor3d_HCurve) Adaptor3d_IsoCurve::Trim
 315  (const Standard_Real First,
 316   const Standard_Real Last,
 317   const Standard_Real) const
 318 {
 319   Handle(Adaptor3d_HIsoCurve) HI = new Adaptor3d_HIsoCurve(*this);
 320   ((Adaptor3d_IsoCurve *)&(HI->Curve()))->Load(myIso,myParameter,First,Last);
 321   return HI;
 322 }
 323
 324 //=======================================================================
 325 //function : IsClosed
 326 //purpose  :
 327 //=======================================================================
 328
 329 Standard_Boolean Adaptor3d_IsoCurve::IsClosed() const
 330 {
 331   switch (myIso) {
 332   case GeomAbs_IsoU:
 333     return mySurface->IsVClosed();
 334   case GeomAbs_IsoV:
 335     return mySurface->IsUClosed();
 336   case GeomAbs_NoneIso:
 337   default:
 338     break;
 339   }
 340
 341   Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 342   return Standard_False;
 343 }
 344
 345 //=======================================================================
 346 //function : IsPeriodic
 347 //purpose  :
 348 //=======================================================================
 349
 350 Standard_Boolean Adaptor3d_IsoCurve::IsPeriodic() const
 351 {
 352   switch (myIso) {
 353   case GeomAbs_IsoU:
 354     return mySurface->IsVPeriodic();
 355   case GeomAbs_IsoV:
 356     return mySurface->IsUPeriodic();
 357   case GeomAbs_NoneIso:
 358   default:
 359     break;
 360   }
 361
 362   Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 363   return Standard_False;
 364 }
 365
 366 //=======================================================================
 367 //function : Period
 368 //purpose  :
 369 //=======================================================================
 370
 371 Standard_Real Adaptor3d_IsoCurve::Period() const
 372 {
 373   switch (myIso) {
 374   case GeomAbs_IsoU:
 375     return mySurface->VPeriod();
 376   case GeomAbs_IsoV:
 377     return mySurface->UPeriod();
 378   case GeomAbs_NoneIso:
 379   default:
 380     break;
 381   }
 382
 383   Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 384   return 0.;
 385 }
 386
 387 //=======================================================================
 388 //function : Value
 389 //purpose  :
 390 //=======================================================================
 391
 392 gp_Pnt Adaptor3d_IsoCurve::Value(const Standard_Real T) const
 393 {
 394   switch (myIso) {
 395
 396   case GeomAbs_IsoU:
 397     return mySurface->Value(myParameter,T);
 398
 399   case GeomAbs_IsoV:
 400     return mySurface->Value(T,myParameter);
 401
 402   case GeomAbs_NoneIso:
 403     {
 404       Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 405       break;
 406     }
 407   }
 408   // portage WNT
 409   return gp_Pnt();
 410 }
 411
 412
 413 //=======================================================================
 414 //function : D0
 415 //purpose  :
 416 //=======================================================================
 417
 418 void Adaptor3d_IsoCurve::D0(const Standard_Real T, gp_Pnt& P) const
 419 {
 420   switch (myIso) {
 421
 422   case GeomAbs_IsoU:
 423     mySurface->D0(myParameter,T,P);
 424     break;
 425
 426   case GeomAbs_IsoV:
 427     mySurface->D0(T,myParameter,P);
 428     break;
 429
 430   case GeomAbs_NoneIso:
 431     Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 432     break;
 433   }
 434 }
 435
 436 //=======================================================================
 437 //function : D1
 438 //purpose  :
 439 //=======================================================================
 440
 441 void Adaptor3d_IsoCurve::D1(const Standard_Real T, gp_Pnt& P, gp_Vec& V) const
 442 {
 443   gp_Vec dummy;
 444   switch (myIso) {
 445
 446   case GeomAbs_IsoU:
 447     mySurface->D1(myParameter,T,P,dummy,V);
 448     break;
 449
 450   case GeomAbs_IsoV:
 451     mySurface->D1(T,myParameter,P,V,dummy);
 452     break;
 453
 454   case GeomAbs_NoneIso:
 455     Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 456     break;
 457   }
 458 }
 459
 460 //=======================================================================
 461 //function : D2
 462 //purpose  :
 463 //=======================================================================
 464
 465 void Adaptor3d_IsoCurve::D2(const Standard_Real T, gp_Pnt& P,
 466                           gp_Vec& V1, gp_Vec& V2) const
 467 {
 468   gp_Vec dummy1,dummy2,dummy3;
 469   switch (myIso) {
 470
 471   case GeomAbs_IsoU:
 472     mySurface->D2(myParameter,T,P,
 473                   dummy1,V1,dummy2,V2,dummy3);
 474     break;
 475   case GeomAbs_IsoV:
 476     mySurface->D2(T,myParameter,
 477                   P,V1,dummy1,V2,dummy2,dummy3);
 478     break;
 479   case GeomAbs_NoneIso:
 480     Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 481     break;
 482   }
 483 }
 484
 485 //=======================================================================
 486 //function : D3
 487 //purpose  :
 488 //=======================================================================
 489
 490 void Adaptor3d_IsoCurve::D3(const Standard_Real T, gp_Pnt& P,
 491                           gp_Vec& V1, gp_Vec& V2, gp_Vec& V3) const
 492 {
 493   gp_Vec dummy[6];
 494   switch (myIso) {
 495
 496   case GeomAbs_IsoU:
 497     mySurface->D3(myParameter,T,P,dummy[0],V1,dummy[1],
 498                   V2,dummy[2],dummy[3],V3,dummy[4],dummy[5]);
 499     break;
 500
 501   case GeomAbs_IsoV:
 502     mySurface->D3(T,myParameter,P,V1,dummy[0],V2,dummy[1],
 503                   dummy[2],V3,dummy[3],dummy[4],dummy[5]);
 504     break;
 505
 506   case GeomAbs_NoneIso:
 507     Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 508     break;
 509   }
 510 }
 511
 512 //=======================================================================
 513 //function : DN
 514 //purpose  :
 515 //=======================================================================
 516
 517 gp_Vec Adaptor3d_IsoCurve::DN(const Standard_Real T,
 518                             const Standard_Integer N) const
 519 {
 520   switch (myIso) {
 521
 522   case GeomAbs_IsoU:
 523     return mySurface->DN(myParameter,T,0,N);
 524   case GeomAbs_IsoV:
 525     return mySurface->DN(T,myParameter,N,0);
 526   case GeomAbs_NoneIso:
 527     {
 528       Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 529       break;
 530     }
 531   }
 532
 533   // portage WNT
 534   return gp_Vec();
 535 }
 536
 537
 538 //=======================================================================
 539 //function : Resolution
 540 //purpose  :
 541 //=======================================================================
 542
 543 Standard_Real Adaptor3d_IsoCurve::Resolution(const Standard_Real R3D) const
 544 {
 545   // Peut-on faire mieux ??
 546   return Precision::Parametric(R3D);
 547 }
 548
 549
 550
 551 //=======================================================================
 552 //function : GetType
 553 //purpose  :
 554 //=======================================================================
 555
 556 GeomAbs_CurveType Adaptor3d_IsoCurve::GetType() const {
 557
 558   switch (mySurface->GetType()) {
 559
 560   case GeomAbs_Plane:
 561     return GeomAbs_Line;
 562
 563   case GeomAbs_Cylinder:
 564   case GeomAbs_Cone:
 565     {
 566       switch (myIso) {
 567       case GeomAbs_IsoU:
 568         return GeomAbs_Line;
 569
 570       case GeomAbs_IsoV:
 571         return GeomAbs_Circle;
 572
 573       case GeomAbs_NoneIso:
 574         {
 575           Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 576           // portage WNT
 577           return GeomAbs_OtherCurve;
 578         }
 579       }
 580       break;
 581     }
 582
 583   case GeomAbs_Sphere:
 584   case GeomAbs_Torus:
 585     return GeomAbs_Circle;
 586
 587   case GeomAbs_BezierSurface:
 588     return GeomAbs_BezierCurve;
 589
 590   case GeomAbs_BSplineSurface:
 591     return GeomAbs_BSplineCurve;
 592
 593   case GeomAbs_SurfaceOfRevolution:
 594     {
 595       switch (myIso) {
 596       case GeomAbs_IsoU:
 597         return mySurface->BasisCurve()->GetType();
 598
 599       case GeomAbs_IsoV:
 600         return GeomAbs_Circle;
 601
 602       case GeomAbs_NoneIso:
 603         Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 604         // portage WNT
 605         return GeomAbs_OtherCurve;
 606         break;
 607       }
 608       break;
 609     }
 610
 611   case GeomAbs_SurfaceOfExtrusion:
 612     {
 613       switch (myIso) {
 614       case GeomAbs_IsoU:
 615         return GeomAbs_Line;
 616
 617       case GeomAbs_IsoV:
 618         return mySurface->BasisCurve()->GetType();
 619
 620       case GeomAbs_NoneIso:
 621         Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 622         // portage WNT
 623         return GeomAbs_OtherCurve;
 624         break;
 625       }
 626       break;
 627     }
 628   default:
 629     return GeomAbs_OtherCurve;
 630   }
 631
 632   // portage WNT
 633   return GeomAbs_OtherCurve;
 634 }
 635
 636 //=======================================================================
 637 //function : Line
 638 //purpose  :
 639 //=======================================================================
 640
 641 gp_Lin Adaptor3d_IsoCurve::Line() const
 642 {
 643   gp_Pnt P;
 644   gp_Vec V;
 645   D1(0,P,V);
 646   return gp_Lin(P,V);
 647 }
 648
 649 //=======================================================================
 650 //function : computeHR
 651 //purpose  :
 652 //=======================================================================
 653
 654 static void computeHR(const gp_Ax3&        axes,
 655                       const gp_Pnt&        P,
 656                             Standard_Real& h,
 657                             Standard_Real& radius)
 658 {
 659   gp_Vec V(axes.Location(),P);
 660   h = V * axes.Direction();
 661   radius = V * axes.XDirection();
 662 }
 663
 664 //=======================================================================
 665 //function : Circle
 666 //purpose  :
 667 //=======================================================================
 668
 669 gp_Circ Adaptor3d_IsoCurve::Circle() const
 670 {
 671   gp_Ax3 axes;
 672   Standard_Real radius,h = 0.;
 673
 674   switch (mySurface->GetType()) {
 675
 676   case GeomAbs_Cylinder:
 677     {
 678       gp_Cylinder cyl = mySurface->Cylinder();
 679
 680       switch (myIso) {
 681
 682       case GeomAbs_IsoU:
 683         {
 684           Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:UIso");
 685           return gp_Circ();
 686         }
 687       case GeomAbs_IsoV:
 688         {
 689           return ElSLib::CylinderVIso(cyl.Position(),cyl.Radius(),myParameter);
 690         }
 691       case GeomAbs_NoneIso:
 692         {
 693           Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 694           return gp_Circ();
 695         }
 696       }
 697       break;
 698     }
 699
 700   case GeomAbs_Cone:
 701     {
 702       gp_Cone cone = mySurface->Cone();
 703
 704       switch (myIso) {
 705
 706       case GeomAbs_IsoU:
 707         {
 708           Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:UIso");
 709           return gp_Circ();
 710         }
 711       case GeomAbs_IsoV:
 712         {
 713           return ElSLib::ConeVIso(cone.Position(),cone.RefRadius(),
 714                                   cone.SemiAngle(),myParameter);
 715         }
 716       case GeomAbs_NoneIso:
 717         {
 718           Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 719           return gp_Circ();
 720         }
 721       }
 722       break;
 723     }
 724
 725   case GeomAbs_Sphere:
 726     {
 727       gp_Sphere sph = mySurface->Sphere();
 728
 729       switch (myIso) {
 730
 731       case GeomAbs_IsoU:
 732         {
 733           return ElSLib::SphereUIso(sph.Position(),sph.Radius(),myParameter);
 734         }
 735
 736       case GeomAbs_IsoV:
 737         {
 738           return ElSLib::SphereVIso(sph.Position(),sph.Radius(),myParameter);
 739         }
 740
 741       case GeomAbs_NoneIso:
 742         {
 743           Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 744           return gp_Circ();
 745         }
 746       }
 747       break;
 748     }
 749
 750   case GeomAbs_Torus: {
 751     gp_Torus tor = mySurface->Torus();
 752
 753     switch (myIso) {
 754
 755     case GeomAbs_IsoU:
 756       {
 757         return ElSLib::TorusUIso(tor.Position(),tor.MajorRadius(),
 758                                  tor.MinorRadius(),myParameter);
 759       }
 760
 761     case GeomAbs_IsoV:
 762       {
 763         return ElSLib::TorusVIso(tor.Position(),tor.MajorRadius(),
 764                                  tor.MinorRadius(),myParameter);
 765       }
 766
 767     case GeomAbs_NoneIso:
 768       {
 769         Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 770         return gp_Circ();
 771       }
 772     }
 773     break;
 774   }
 775
 776   case GeomAbs_SurfaceOfRevolution:
 777     {
 778       if (myIso == GeomAbs_IsoV) {
 779         const gp_Pnt aVal0 = Value (0.0);
 780         gp_Ax1 Ax1 = mySurface->AxeOfRevolution();
 781         if (gp_Lin (Ax1).Contains (aVal0, Precision::Confusion())) {
 782           return gp_Circ(gp_Ax2(aVal0, Ax1.Direction()),0);
 783         }
 784         else {
 785           gp_Vec DX(Ax1.Location(), aVal0);
 786           axes = gp_Ax3(Ax1.Location(), Ax1.Direction(), DX);
 787           computeHR(axes,aVal0,h,radius);
 788           gp_Vec VT = axes.Direction();
 789           axes.Translate(VT * h);
 790           return gp_Circ(axes.Ax2(),radius);
 791         }
 792       }
 793       else {
 794         return mySurface->BasisCurve()->Circle().Rotated
 795           (mySurface->AxeOfRevolution(),myParameter);
 796       }
 797     }
 798
 799   case GeomAbs_SurfaceOfExtrusion: {
 800     return  mySurface->BasisCurve()->Circle().Translated
 801       (myParameter * gp_Vec(mySurface->Direction()));
 802   }
 803   default:
 804     {
 805       Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:Circle");
 806       return gp_Circ();
 807     }
 808
 809   }
 810   // portage WNT
 811   return gp_Circ();
 812 }
 813
 814 //=======================================================================
 815 //function : Ellipse
 816 //purpose  :
 817 //=======================================================================
 818
 819 gp_Elips Adaptor3d_IsoCurve::Ellipse() const
 820 {
 821   switch (mySurface->GetType()) {
 822
 823   case GeomAbs_SurfaceOfExtrusion: {
 824     return  mySurface->BasisCurve()->Ellipse().Translated
 825       (myParameter * gp_Vec(mySurface->Direction()));
 826   }
 827   default:
 828     {
 829       Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:Ellipse");
 830     }
 831   }
 832   return gp_Elips();
 833 }
 834
 835 //=======================================================================
 836 //function : Hyperbola
 837 //purpose  :
 838 //=======================================================================
 839
 840 gp_Hypr Adaptor3d_IsoCurve::Hyperbola() const
 841 {
 842   Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:Hyperbola");
 843   return gp_Hypr();
 844 }
 845
 846 //=======================================================================
 847 //function : Parabola
 848 //purpose  :
 849 //=======================================================================
 850
 851 gp_Parab Adaptor3d_IsoCurve::Parabola() const
 852 {
 853   Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:Parabola");
 854   return gp_Parab();
 855 }
 856
 857 //=======================================================================
 858 //function : Degree
 859 //purpose  :
 860 //=======================================================================
 861
 862 Standard_Integer Adaptor3d_IsoCurve::Degree() const
 863 {
 864   Standard_Integer degree = 0 ;
 865   GeomAbs_SurfaceType type = mySurface->GetType() ;
 866   switch(type) {
 867   case GeomAbs_BezierSurface:
 868   case GeomAbs_BSplineSurface:
 869     {
 870       switch (myIso) {
 871       case GeomAbs_IsoU:
 872         degree = mySurface->VDegree() ;
 873         break ;
 874       case GeomAbs_IsoV:
 875         degree = mySurface->UDegree() ;
 876         break ;
 877
 878       case GeomAbs_NoneIso:
 879       default:
 880         Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 881       }
 882     }
 883     break ;
 884   case GeomAbs_SurfaceOfRevolution:
 885     {
 886       switch (myIso) {
 887       case GeomAbs_IsoU:
 888         degree = mySurface->BasisCurve()->Degree();
 889         break;
 890       default:
 891         Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 892       }
 893     }
 894     break;
 895   case GeomAbs_SurfaceOfExtrusion:
 896     {
 897       switch (myIso) {
 898       case GeomAbs_IsoV:
 899         degree = mySurface->BasisCurve()->Degree();
 900         break;
 901       default:
 902         Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 903       }
 904     }
 905     break;
 906   default:
 907     Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 908     break ;
 909   }
 910   return degree ;
 911 }
 912
 913 //=======================================================================
 914 //function : IsRational
 915 //purpose  :
 916 //=======================================================================
 917
 918 Standard_Boolean Adaptor3d_IsoCurve::IsRational() const
 919 {
 920   Standard_Integer is_rational = Standard_False ;
 921   GeomAbs_SurfaceType type = mySurface->GetType() ;
 922   switch(type) {
 923   case GeomAbs_BezierSurface:
 924   case GeomAbs_BSplineSurface:
 925     {
 926       switch (myIso) {
 927       case GeomAbs_IsoU:
 928         is_rational = mySurface->IsVRational() ;
 929         break ;
 930       case GeomAbs_IsoV:
 931         is_rational = mySurface->IsURational() ;
 932         break ;
 933
 934       case GeomAbs_NoneIso:
 935       default:
 936         Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 937       }
 938     }
 939     break ;
 940   case GeomAbs_SurfaceOfRevolution:
 941     {
 942       switch (myIso) {
 943       case GeomAbs_IsoU:
 944         is_rational = mySurface->BasisCurve()->IsRational();
 945         break;
 946       default:
 947         Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 948       }
 949     }
 950     break;
 951   case GeomAbs_SurfaceOfExtrusion:
 952     {
 953       switch (myIso) {
 954       case GeomAbs_IsoV:
 955         is_rational = mySurface->BasisCurve()->IsRational();
 956         break;
 957       default:
 958         Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 959       }
 960     }
 961     break;
 962   default:
 963     Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 964   }
 965   return is_rational;
 966 }
 967
 968
 969 //=======================================================================
 970 //function : NbPoles
 971 //purpose  :
 972 //=======================================================================
 973
 974 Standard_Integer Adaptor3d_IsoCurve::NbPoles() const
 975 {
 976   Standard_Integer nb_poles = 0 ;
 977   GeomAbs_SurfaceType type = mySurface->GetType() ;
 978   switch(type) {
 979   case GeomAbs_BezierSurface:
 980   case GeomAbs_BSplineSurface:
 981     switch (myIso) {
 982     case GeomAbs_IsoU:
 983       nb_poles = mySurface->NbVPoles() ;
 984       break ;
 985     case GeomAbs_IsoV:
 986       nb_poles = mySurface->NbUPoles() ;
 987       break ;
 988
 989     case GeomAbs_NoneIso:
 990     default:
 991       Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
 992     }
 993     break ;
 994   case GeomAbs_SurfaceOfRevolution:
 995     {
 996       switch( myIso) {
 997       case GeomAbs_IsoU:
 998         {
 999           nb_poles = mySurface->BasisCurve()->NbPoles();
1000         }
1001         break;
1002       default:
1003         Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
1004       }
1005     }
1006     break;
1007   case GeomAbs_SurfaceOfExtrusion:
1008     {
1009       switch( myIso) {
1010       case GeomAbs_IsoV:
1011         {
1012           nb_poles = mySurface->BasisCurve()->NbPoles();
1013         }
1014         break;
1015       default:
1016         Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
1017       }
1018     }
1019     break;
1020
1021   default:
1022     Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
1023     break ;
1024   }
1025  return nb_poles ;
1026  }
1027
1028 //=======================================================================
1029 //function : NbKnots
1030 //purpose  :
1031 //=======================================================================
1032
1033 Standard_Integer Adaptor3d_IsoCurve::NbKnots() const
1034 {
1035   Standard_Integer nb_knots = 0 ;
1036   GeomAbs_SurfaceType type = mySurface->GetType() ;
1037   switch(type) {
1038   case GeomAbs_BSplineSurface:
1039     {
1040       switch (myIso) {
1041       case GeomAbs_IsoU:
1042         nb_knots = mySurface->NbVKnots() ;
1043         break ;
1044       case GeomAbs_IsoV:
1045         nb_knots = mySurface->NbUKnots() ;
1046         break ;
1047
1048       case GeomAbs_NoneIso:
1049       default:
1050         Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
1051       }
1052     }
1053     break ;
1054   case GeomAbs_SurfaceOfRevolution:
1055     {
1056       switch (myIso) {
1057       case GeomAbs_IsoU:
1058         {
1059           nb_knots = mySurface->BasisCurve()->NbKnots() ;
1060           break ;
1061         }
1062       default:
1063         Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
1064       }
1065     }
1066     break ;
1067   case GeomAbs_SurfaceOfExtrusion:
1068     {
1069       switch (myIso) {
1070       case GeomAbs_IsoV:
1071         {
1072           nb_knots = mySurface->BasisCurve()->NbKnots() ;
1073           break ;
1074         }
1075       default:
1076         Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
1077       }
1078     }
1079     break ;
1080   default:
1081     Standard_NoSuchObject::Raise("Adaptor3d_IsoCurve:NoneIso");
1082     break ;
1083   }
1084   return nb_knots ;
1085  }
1086
1087 //=======================================================================
1088 //function : Bezier
1089 //purpose  :
1090 //=======================================================================
1091
1092 Handle(Geom_BezierCurve) Adaptor3d_IsoCurve::Bezier() const
1093 {
1094   Handle(Geom_BezierCurve) C;
1095   if (mySurface->GetType() == GeomAbs_SurfaceOfRevolution) {
1096     C = mySurface->BasisCurve()->Bezier();
1097     C = Handle(Geom_BezierCurve)::DownCast(C->Copy());
1098     C->Rotate(mySurface->AxeOfRevolution(),myParameter);
1099   }
1100   else if (mySurface->GetType() == GeomAbs_SurfaceOfExtrusion) {
1101     C = mySurface->BasisCurve()->Bezier();
1102     C = Handle(Geom_BezierCurve)::DownCast(C->Copy());
1103     C->Translate(myParameter * gp_Vec(mySurface->Direction()));
1104   }
1105   else if (myIso == GeomAbs_IsoU) {
1106     C = Handle(Geom_BezierCurve)::DownCast
1107       (mySurface->Bezier()->UIso(myParameter)) ;
1108   }
1109   else {
1110     C = Handle(Geom_BezierCurve)::DownCast
1111       (mySurface->Bezier()->VIso(myParameter));
1112   }
1113 //  C->Segment(myFirst,myLast);
1114   return C;
1115 }
1116
1117 //=======================================================================
1118 //function : BSpline
1119 //purpose  :
1120 //=======================================================================
1121
1122 Handle(Geom_BSplineCurve) Adaptor3d_IsoCurve::BSpline() const
1123 {
1124   Handle(Geom_BSplineCurve) C;
1125   if (mySurface->GetType() == GeomAbs_SurfaceOfRevolution) {
1126     C = mySurface->BasisCurve()->BSpline();
1127     C = Handle(Geom_BSplineCurve)::DownCast(C->Copy());
1128     C->Rotate(mySurface->AxeOfRevolution(),myParameter);
1129   }
1130   else if (mySurface->GetType() == GeomAbs_SurfaceOfExtrusion) {
1131     C = mySurface->BasisCurve()->BSpline();
1132     C = Handle(Geom_BSplineCurve)::DownCast(C->Copy());
1133     C->Translate(myParameter * gp_Vec(mySurface->Direction()));
1134   }
1135   else if (myIso == GeomAbs_IsoU) {
1136     C = Handle(Geom_BSplineCurve)::DownCast
1137       (mySurface->BSpline()->UIso(myParameter)) ;
1138   }
1139   else {
1140     C = Handle(Geom_BSplineCurve)::DownCast
1141       (mySurface->BSpline()->VIso(myParameter));
1142   }
1143 //  C->Segment(myFirst,myLast);
1144   return C;
1145 }
1146