b311480e |
1 | // Created on: 1995-06-13 |
2 | // Created by: Bruno DUMORTIER |
3 | // Copyright (c) 1995-1999 Matra Datavision |
973c2be1 |
4 | // Copyright (c) 1999-2014 OPEN CASCADE SAS |
b311480e |
5 | // |
973c2be1 |
6 | // This file is part of Open CASCADE Technology software library. |
b311480e |
7 | // |
d5f74e42 |
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 |
973c2be1 |
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. |
b311480e |
13 | // |
973c2be1 |
14 | // Alternatively, this file may be used under the terms of Open CASCADE |
15 | // commercial license or contractual agreement. |
7fd59977 |
16 | |
42cf5bc1 |
17 | |
18 | #include <AdvApprox_ApproxAFunction.hxx> |
19 | #include <Geom2d_BSplineCurve.hxx> |
20 | #include <Geom_BSplineCurve.hxx> |
21 | #include <GeomLib_MakeCurvefromApprox.hxx> |
7fd59977 |
22 | #include <gp_Pnt2d.hxx> |
42cf5bc1 |
23 | #include <Standard_OutOfRange.hxx> |
24 | #include <StdFail_NotDone.hxx> |
7fd59977 |
25 | #include <TColgp_Array1OfPnt.hxx> |
26 | #include <TColgp_Array1OfPnt2d.hxx> |
42cf5bc1 |
27 | #include <TColStd_Array1OfInteger.hxx> |
28 | #include <TColStd_Array1OfReal.hxx> |
29 | #include <TColStd_HArray1OfInteger.hxx> |
30 | #include <TColStd_HArray1OfReal.hxx> |
7fd59977 |
31 | |
32 | //======================================================================= |
33 | //function : GeomLib_MakeCurvefromApprox |
34 | //purpose : |
35 | //======================================================================= |
7fd59977 |
36 | GeomLib_MakeCurvefromApprox::GeomLib_MakeCurvefromApprox |
37 | (const AdvApprox_ApproxAFunction& Approx) |
38 | :myApprox(Approx) |
39 | { |
40 | } |
41 | |
42 | |
43 | //======================================================================= |
44 | //function : Nb1DSpaces |
45 | //purpose : |
46 | //======================================================================= |
47 | |
48 | Standard_Integer GeomLib_MakeCurvefromApprox::Nb1DSpaces() const |
49 | { |
50 | return myApprox.NumSubSpaces(1); |
51 | } |
52 | |
53 | |
54 | //======================================================================= |
55 | //function : Nb2DSpaces |
56 | //purpose : |
57 | //======================================================================= |
58 | |
59 | Standard_Integer GeomLib_MakeCurvefromApprox::Nb2DSpaces() const |
60 | { |
61 | return myApprox.NumSubSpaces(2); |
62 | } |
63 | |
64 | |
65 | //======================================================================= |
66 | //function : Nb3DSpaces |
67 | //purpose : |
68 | //======================================================================= |
69 | |
70 | Standard_Integer GeomLib_MakeCurvefromApprox::Nb3DSpaces() const |
71 | { |
72 | return myApprox.NumSubSpaces(3); |
73 | } |
74 | |
75 | |
76 | //======================================================================= |
77 | //function : Curve2d |
78 | //purpose : |
79 | //======================================================================= |
80 | |
81 | Handle(Geom2d_BSplineCurve) GeomLib_MakeCurvefromApprox::Curve2d |
82 | (const Standard_Integer Index2d) const |
83 | { |
84 | Standard_OutOfRange_Raise_if |
85 | ( Index2d < 0 || Index2d > Nb2DSpaces(), |
86 | " GeomLib_MakeCurvefromApprox : Curve2d"); |
87 | StdFail_NotDone_Raise_if |
88 | ( !IsDone(), |
89 | " GeomLib_MakeCurvefromApprox : Curve2d"); |
90 | |
91 | TColgp_Array1OfPnt2d Poles( 1, myApprox.NbPoles()); |
92 | TColStd_Array1OfReal Knots( 1, myApprox.NbKnots()); |
93 | TColStd_Array1OfInteger Mults( 1, myApprox.NbKnots()); |
94 | |
95 | myApprox.Poles2d(Index2d, Poles); |
96 | Knots = myApprox.Knots()->Array1(); |
97 | Mults = myApprox.Multiplicities()->Array1(); |
98 | |
99 | Handle(Geom2d_BSplineCurve) C = |
100 | new Geom2d_BSplineCurve( Poles, Knots, Mults, myApprox.Degree()); |
101 | |
102 | return C; |
103 | } |
104 | |
105 | |
106 | //======================================================================= |
107 | //function : Curve2d |
108 | //purpose : |
109 | //======================================================================= |
110 | |
111 | Handle(Geom2d_BSplineCurve) GeomLib_MakeCurvefromApprox::Curve2d |
112 | (const Standard_Integer Index1d, |
113 | const Standard_Integer Index2d) const |
114 | { |
115 | Standard_OutOfRange_Raise_if |
116 | ( Index1d < 0 || Index1d > Nb1DSpaces() || |
117 | Index2d < 0 || Index2d > Nb2DSpaces(), |
118 | " GeomLib_MakeCurvefromApprox : Curve2d"); |
119 | StdFail_NotDone_Raise_if |
120 | ( !IsDone() && ! myApprox.HasResult(), |
121 | " GeomLib_MakeCurvefromApprox : Curve2d"); |
122 | |
123 | TColgp_Array1OfPnt2d Poles ( 1, myApprox.NbPoles()); |
124 | TColStd_Array1OfReal Weigths( 1, myApprox.NbPoles()); |
125 | TColStd_Array1OfReal Knots ( 1, myApprox.NbKnots()); |
126 | TColStd_Array1OfInteger Mults ( 1, myApprox.NbKnots()); |
127 | |
128 | myApprox.Poles2d(Index2d, Poles); |
129 | myApprox.Poles1d(Index1d, Weigths); |
130 | Knots = myApprox.Knots()->Array1(); |
131 | Mults = myApprox.Multiplicities()->Array1(); |
132 | |
133 | Standard_Real X,Y,W; |
134 | for ( Standard_Integer i = 1; i <= myApprox.NbPoles(); i++) { |
135 | W = Weigths(i); |
136 | Poles(i).Coord(X,Y); |
137 | Poles(i).SetCoord(X/W, Y/W); |
138 | } |
139 | |
140 | Handle(Geom2d_BSplineCurve) C = |
141 | new Geom2d_BSplineCurve( Poles, Knots, Mults, myApprox.Degree()); |
142 | |
143 | return C; |
144 | } |
145 | //======================================================================= |
146 | //function : Curve2dFromTwo1d |
147 | //purpose : |
148 | //======================================================================= |
149 | |
150 | Handle(Geom2d_BSplineCurve) GeomLib_MakeCurvefromApprox::Curve2dFromTwo1d |
151 | (const Standard_Integer Index1d, |
152 | const Standard_Integer Index2d) const |
153 | { |
154 | Standard_OutOfRange_Raise_if |
155 | ( Index1d < 0 || Index1d > Nb1DSpaces() || |
156 | Index2d < 0 || Index2d > Nb1DSpaces(), |
157 | " GeomLib_MakeCurvefromApprox : Curve2d"); |
158 | StdFail_NotDone_Raise_if |
159 | ( !IsDone() && ! myApprox.HasResult(), |
160 | " GeomLib_MakeCurvefromApprox : Curve2d"); |
161 | |
162 | TColgp_Array1OfPnt2d Poles ( 1, myApprox.NbPoles()); |
163 | TColStd_Array1OfReal Poles1d1( 1, myApprox.NbPoles()); |
164 | TColStd_Array1OfReal Poles1d2( 1, myApprox.NbPoles()); |
165 | TColStd_Array1OfReal Knots ( 1, myApprox.NbKnots()); |
166 | TColStd_Array1OfInteger Mults ( 1, myApprox.NbKnots()); |
167 | |
168 | myApprox.Poles1d(Index2d, Poles1d2); |
169 | myApprox.Poles1d(Index1d, Poles1d1); |
170 | |
171 | Knots = myApprox.Knots()->Array1(); |
172 | Mults = myApprox.Multiplicities()->Array1(); |
173 | |
174 | // Standard_Real X,Y,W; |
175 | for ( Standard_Integer i = 1; i <= myApprox.NbPoles(); i++) { |
176 | Poles(i).SetCoord(Poles1d1.Value(i),Poles1d2.Value(i)); |
177 | } |
178 | |
179 | Handle(Geom2d_BSplineCurve) C = |
180 | new Geom2d_BSplineCurve( Poles, Knots, Mults, myApprox.Degree()); |
181 | |
182 | return C; |
183 | } |
184 | |
185 | |
186 | |
187 | //======================================================================= |
188 | //function : Curve |
189 | //purpose : |
190 | //======================================================================= |
191 | |
192 | Handle(Geom_BSplineCurve) GeomLib_MakeCurvefromApprox::Curve |
193 | (const Standard_Integer Index3d) const |
194 | { |
195 | Standard_OutOfRange_Raise_if |
196 | ( Index3d < 0 || Index3d > Nb3DSpaces(), |
197 | " GeomLib_MakeCurvefromApprox : Curve"); |
198 | StdFail_NotDone_Raise_if |
199 | ( !IsDone() && ! myApprox.HasResult(), |
200 | " GeomLib_MakeCurvefromApprox : Curve"); |
201 | |
202 | TColgp_Array1OfPnt Poles( 1, myApprox.NbPoles()); |
203 | TColStd_Array1OfReal Knots( 1, myApprox.NbKnots()); |
204 | TColStd_Array1OfInteger Mults( 1, myApprox.NbKnots()); |
205 | |
206 | myApprox.Poles(Index3d, Poles); |
207 | Knots = myApprox.Knots()->Array1(); |
208 | Mults = myApprox.Multiplicities()->Array1(); |
209 | |
210 | Handle(Geom_BSplineCurve) C = |
211 | new Geom_BSplineCurve( Poles, Knots, Mults, myApprox.Degree()); |
212 | |
213 | return C; |
214 | } |
215 | |
216 | |
217 | //======================================================================= |
218 | //function : Curve |
219 | //purpose : |
220 | //======================================================================= |
221 | |
222 | Handle(Geom_BSplineCurve) GeomLib_MakeCurvefromApprox::Curve |
223 | (const Standard_Integer Index1d, |
224 | const Standard_Integer Index3d) const |
225 | { |
226 | Standard_OutOfRange_Raise_if |
227 | ( Index1d < 0 || Index1d > Nb1DSpaces() || |
228 | Index3d < 0 || Index3d > Nb3DSpaces(), |
229 | " GeomLib_MakeCurvefromApprox : Curve3d"); |
230 | StdFail_NotDone_Raise_if |
231 | ( !IsDone(), |
232 | " GeomLib_MakeCurvefromApprox : Curve3d"); |
233 | |
234 | TColgp_Array1OfPnt Poles ( 1, myApprox.NbPoles()); |
235 | TColStd_Array1OfReal Weigths( 1, myApprox.NbPoles()); |
236 | TColStd_Array1OfReal Knots ( 1, myApprox.NbKnots()); |
237 | TColStd_Array1OfInteger Mults ( 1, myApprox.NbKnots()); |
238 | |
239 | myApprox.Poles (Index3d, Poles); |
240 | myApprox.Poles1d(Index1d, Weigths); |
241 | Knots = myApprox.Knots()->Array1(); |
242 | Mults = myApprox.Multiplicities()->Array1(); |
243 | |
244 | Standard_Real X,Y,Z,W; |
245 | for ( Standard_Integer i = 1; i <= myApprox.NbPoles(); i++) { |
246 | W = Weigths(i); |
247 | Poles(i).Coord(X,Y,Z); |
248 | Poles(i).SetCoord(X/W, Y/W, Z/W); |
249 | } |
250 | |
251 | Handle(Geom_BSplineCurve) C = |
252 | new Geom_BSplineCurve( Poles, Knots, Mults, myApprox.Degree()); |
253 | |
254 | return C; |
255 | } |