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