0030675: Visualization - remove redundant proxy classes in hierarchy of PrsMgr_Presen...
[occt.git] / src / Geom2dConvert / Geom2dConvert_BSplineCurveToBezierCurve.cxx
1 // Created on: 1996-03-12
2 // Created by: Bruno DUMORTIER
3 // Copyright (c) 1995-1999 Matra Datavision
4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
5 //
6 // This file is part of Open CASCADE Technology software library.
7 //
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.
13 //
14 // Alternatively, this file may be used under the terms of Open CASCADE
15 // commercial license or contractual agreement.
16
17
18 #include <Geom2d_BezierCurve.hxx>
19 #include <Geom2d_BSplineCurve.hxx>
20 #include <Geom2dConvert_BSplineCurveToBezierCurve.hxx>
21 #include <Standard_DimensionError.hxx>
22 #include <Standard_DomainError.hxx>
23 #include <Standard_OutOfRange.hxx>
24 #include <TColgp_Array1OfPnt2d.hxx>
25 #include <TColStd_Array1OfReal.hxx>
26
27 //=======================================================================
28 //function : Geom2dConvert_BSplineCurveToBezierCurve
29 //purpose  : 
30 //=======================================================================
31 Geom2dConvert_BSplineCurveToBezierCurve::Geom2dConvert_BSplineCurveToBezierCurve (const Handle(Geom2d_BSplineCurve)& BasisCurve)
32 {
33   myCurve = Handle(Geom2d_BSplineCurve)::DownCast(BasisCurve->Copy());
34   // periodic curve can't be converted correctly by two main reasons:
35   // last pole (equal to first one) is missing;
36   // poles recomputation using default boor scheme is fails.
37   if(myCurve->IsPeriodic()) myCurve->SetNotPeriodic();
38   Standard_Real Uf = myCurve->FirstParameter();
39   Standard_Real Ul = myCurve->LastParameter();
40   myCurve->Segment(Uf,Ul);
41   myCurve->IncreaseMultiplicity(myCurve->FirstUKnotIndex(),
42                                 myCurve->LastUKnotIndex(),
43                                 myCurve->Degree());
44 }
45
46
47 //=======================================================================
48 //function : Geom2dConvert_BSplineCurveToBezierCurve
49 //purpose  : 
50 //=======================================================================
51
52 Geom2dConvert_BSplineCurveToBezierCurve::Geom2dConvert_BSplineCurveToBezierCurve
53 (const Handle(Geom2d_BSplineCurve)& BasisCurve, 
54  const Standard_Real U1, 
55  const Standard_Real U2,
56  const Standard_Real ParametricTolerance)
57 {
58  if (U2 - U1 < ParametricTolerance)
59       throw Standard_DomainError("GeomConvert_BSplineCurveToBezierSurface");
60
61   Standard_Real Uf = U1, Ul = U2; 
62   Standard_Real PTol = ParametricTolerance/2 ;
63
64   Standard_Integer I1, I2;
65   myCurve = Handle(Geom2d_BSplineCurve)::DownCast(BasisCurve->Copy());
66   if(myCurve->IsPeriodic()) myCurve->SetNotPeriodic();
67
68   myCurve->LocateU(U1, PTol, I1, I2);
69   if (I1==I2) { // On est sur le noeud
70     if ( myCurve->Knot(I1) > U1) Uf = myCurve->Knot(I1);
71   }
72
73   myCurve->LocateU(U2, PTol, I1, I2);
74   if (I1==I2) { // On est sur le noeud
75     if ( myCurve->Knot(I1) < U2) Ul = myCurve->Knot(I1);
76   }
77
78
79   myCurve->Segment(Uf, Ul);
80   myCurve->IncreaseMultiplicity(myCurve->FirstUKnotIndex(),
81                                 myCurve->LastUKnotIndex(),
82                                 myCurve->Degree());
83 }
84
85
86 //=======================================================================
87 //function : Arc
88 //purpose  : 
89 //=======================================================================
90
91 Handle(Geom2d_BezierCurve) Geom2dConvert_BSplineCurveToBezierCurve::Arc
92 (const Standard_Integer Index)
93 {
94   if ( Index < 1 || Index > myCurve->NbKnots()-1) {
95     throw Standard_OutOfRange("Geom2dConvert_BSplineCurveToBezierCurve");
96   }
97   Standard_Integer Deg = myCurve->Degree();
98
99   TColgp_Array1OfPnt2d Poles(1,Deg+1);
100
101   Handle(Geom2d_BezierCurve) C;
102   if ( myCurve->IsRational()) {
103     TColStd_Array1OfReal Weights(1,Deg+1);
104     for ( Standard_Integer i = 1; i <= Deg+1; i++) {
105       Poles(i)   = myCurve->Pole(i+Deg*(Index-1));
106       Weights(i) = myCurve->Weight(i+Deg*(Index-1));
107     }
108     C = new Geom2d_BezierCurve(Poles,Weights);
109   }
110   else {
111     for ( Standard_Integer i = 1; i <= Deg+1; i++) {
112       Poles(i)   = myCurve->Pole(i+Deg*(Index-1));
113     }
114     C = new Geom2d_BezierCurve(Poles);
115   }
116   return C;
117 }
118
119
120 //=======================================================================
121 //function : Arcs
122 //purpose  : 
123 //=======================================================================
124
125 void Geom2dConvert_BSplineCurveToBezierCurve::Arcs
126 (TColGeom2d_Array1OfBezierCurve& Curves)
127 {
128   Standard_Integer n = NbArcs();
129   for ( Standard_Integer i = 1; i <= n; i++) {
130     Curves(i) = Arc(i);
131   } 
132 }
133
134 //=======================================================================
135 //function : Knots
136 //purpose  : 
137 //=======================================================================
138
139 void Geom2dConvert_BSplineCurveToBezierCurve::Knots
140      (TColStd_Array1OfReal& TKnots) const
141 {
142  Standard_Integer ii, kk;
143   for (ii = 1, kk = TKnots.Lower();
144        ii <=myCurve->NbKnots(); ii++, kk++)
145     TKnots(kk) = myCurve->Knot(ii);
146 }
147
148 //=======================================================================
149 //function : NbArcs
150 //purpose  : 
151 //=======================================================================
152
153 Standard_Integer Geom2dConvert_BSplineCurveToBezierCurve::NbArcs() const 
154 {
155   return (myCurve->NbKnots()-1);
156 }