55e97b7208b1560bf27ba031710e076911514b91
[occt.git] / src / Geom / Geom_CylindricalSurface.cxx
1 // Created on: 1993-03-10
2 // Created by: JCV
3 // Copyright (c) 1993-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 #include <Geom_CylindricalSurface.ixx>
18
19 #include <Precision.hxx>
20 #include <GeomAbs_UVSense.hxx>
21 #include <Geom_Line.hxx>
22 #include <Geom_Circle.hxx>
23 #include <Standard_ConstructionError.hxx>
24 #include <Standard_RangeError.hxx>
25
26 #include <gp_Lin.hxx>
27 #include <gp_Circ.hxx>
28 #include <gp_XYZ.hxx>
29 #include <gp_Ax2d.hxx>
30 #include <gp.hxx>
31 #include <ElSLib.hxx>
32
33 typedef Geom_CylindricalSurface         CylindricalSurface;
34 typedef Handle(Geom_CylindricalSurface) Handle(CylindricalSurface);
35
36 typedef gp_Ax1  Ax1;
37 typedef gp_Ax2  Ax2;
38 typedef gp_Ax3  Ax3;
39 typedef gp_Circ Circ;
40 typedef gp_Dir  Dir;
41 typedef gp_Lin  Lin;
42 typedef gp_Pnt  Pnt;
43 typedef gp_Trsf Trsf;
44 typedef gp_Vec  Vec;
45 typedef gp_XYZ  XYZ;
46
47
48
49 //=======================================================================
50 //function : Copy
51 //purpose  : 
52 //=======================================================================
53
54 Handle(Geom_Geometry) Geom_CylindricalSurface::Copy () const {
55
56    Handle(CylindricalSurface) Cs;
57    Cs = new CylindricalSurface (pos, radius);
58    return Cs;
59 }
60
61
62 //=======================================================================
63 //function : Geom_CylindricalSurface
64 //purpose  : 
65 //=======================================================================
66
67 Geom_CylindricalSurface::Geom_CylindricalSurface (const gp_Cylinder& C) 
68 : radius (C.Radius()) {
69
70   pos = C.Position (); 
71 }
72
73
74 //=======================================================================
75 //function : Geom_CylindricalSurface
76 //purpose  : 
77 //=======================================================================
78
79 Geom_CylindricalSurface::Geom_CylindricalSurface ( const Ax3& A3, 
80                                                    const Standard_Real R) 
81 : radius (R) {
82
83   if (R < 0.0) Standard_ConstructionError::Raise();
84   pos = A3;
85 }
86
87
88 //=======================================================================
89 //function : UReversedParameter
90 //purpose  : 
91 //=======================================================================
92
93 Standard_Real Geom_CylindricalSurface::UReversedParameter( const Standard_Real U) const
94 {
95   return (2.*M_PI - U);
96 }
97
98 //=======================================================================
99 //function : VReversedParameter
100 //purpose  : 
101 //=======================================================================
102
103 Standard_Real Geom_CylindricalSurface::VReversedParameter( const Standard_Real V) const 
104 {
105   return (-V);
106 }
107
108 //=======================================================================
109 //function : Radius
110 //purpose  : 
111 //=======================================================================
112
113 Standard_Real Geom_CylindricalSurface::Radius () const             { return radius; }
114
115 //=======================================================================
116 //function : IsUClosed
117 //purpose  : 
118 //=======================================================================
119
120 Standard_Boolean Geom_CylindricalSurface::IsUClosed () const       { return Standard_True; }
121
122 //=======================================================================
123 //function : IsVClosed
124 //purpose  : 
125 //=======================================================================
126
127 Standard_Boolean Geom_CylindricalSurface::IsVClosed () const       { return Standard_False; }
128
129 //=======================================================================
130 //function : IsUPeriodic
131 //purpose  : 
132 //=======================================================================
133
134 Standard_Boolean Geom_CylindricalSurface::IsUPeriodic () const     { return Standard_True; }
135
136 //=======================================================================
137 //function : IsVPeriodic
138 //purpose  : 
139 //=======================================================================
140
141 Standard_Boolean Geom_CylindricalSurface::IsVPeriodic () const     { return Standard_False; }
142
143 //=======================================================================
144 //function : SetCylinder
145 //purpose  : 
146 //=======================================================================
147
148 void Geom_CylindricalSurface::SetCylinder (const gp_Cylinder& C) {
149
150    radius = C.Radius();
151    pos = C.Position();
152 }
153
154
155 //=======================================================================
156 //function : SetRadius
157 //purpose  : 
158 //=======================================================================
159
160 void Geom_CylindricalSurface::SetRadius (const Standard_Real R) {
161
162   if (R < 0.0) { Standard_ConstructionError::Raise(); }
163   radius = R;
164 }
165
166
167 //=======================================================================
168 //function : Bounds
169 //purpose  : 
170 //=======================================================================
171
172 void Geom_CylindricalSurface::Bounds (Standard_Real& U1, Standard_Real& U2, 
173                                       Standard_Real& V1, Standard_Real& V2) const {
174
175    U1 = 0.0;  U2 = 2.0 * M_PI;
176    V1 = - Precision::Infinite();   V2 = Precision::Infinite();
177 }
178
179
180 //=======================================================================
181 //function : Coefficients
182 //purpose  : 
183 //=======================================================================
184
185 void Geom_CylindricalSurface::Coefficients (Standard_Real& A1, Standard_Real& A2, Standard_Real& A3, 
186                                             Standard_Real& B1, Standard_Real& B2, Standard_Real& B3, 
187                                             Standard_Real& C1, Standard_Real& C2, Standard_Real& C3, 
188                                             Standard_Real& D) const {
189 // Dans le repere local du cylindre :
190 // X**2 + Y**2 - radius = 0.0
191
192   Trsf T;
193   T.SetTransformation (pos);
194   Standard_Real T11 = T.Value (1, 1);
195   Standard_Real T12 = T.Value (1, 2);
196   Standard_Real T13 = T.Value (1, 3);
197   Standard_Real T14 = T.Value (1, 4);
198   Standard_Real T21 = T.Value (2, 1);
199   Standard_Real T22 = T.Value (2, 2);
200   Standard_Real T23 = T.Value (2, 3);
201   Standard_Real T24 = T.Value (2, 4);
202   A1 = T11 * T11 + T21 * T21;
203   A2 = T12 * T12 + T22 * T22;
204   A3 = T13 * T13 + T23 * T23;
205   B1 = T11 * T12 + T21 * T22;
206   B2 = T11 * T13 + T21 * T23;
207   B3 = T12 * T13 + T22 * T23;
208   C1 = T11 * T14 + T21 * T24;
209   C2 = T12 * T14 + T22 * T24;
210   C3 = T13 * T14 + T23 * T24;
211   D = T14 * T14 + T24 * T24 - radius * radius;  
212 }
213
214
215
216 //=======================================================================
217 //function : Cylinder
218 //purpose  : 
219 //=======================================================================
220
221 gp_Cylinder Geom_CylindricalSurface::Cylinder () const {
222   
223   return gp_Cylinder (pos, radius);
224 }
225
226
227 //=======================================================================
228 //function : D0
229 //purpose  : 
230 //=======================================================================
231
232 void Geom_CylindricalSurface::D0 (const Standard_Real U, 
233                                   const Standard_Real V,
234                                         Pnt& P ) const 
235 {
236   ElSLib::CylinderD0 (U, V, pos, radius, P);
237 }
238
239
240 //=======================================================================
241 //function : D1
242 //purpose  : 
243 //=======================================================================
244
245 void Geom_CylindricalSurface::D1 (const Standard_Real U, const Standard_Real V , 
246                                         Pnt& P, 
247                                         Vec& D1U,     Vec& D1V) const 
248 {
249   ElSLib::CylinderD1 (U, V, pos, radius, P, D1U, D1V);
250 }
251
252
253
254 //=======================================================================
255 //function : D2
256 //purpose  : 
257 //=======================================================================
258
259 void Geom_CylindricalSurface::D2 (const Standard_Real U, const Standard_Real V, 
260                                   Pnt& P,
261                                   Vec& D1U, Vec& D1V,
262                                   Vec& D2U, Vec& D2V, Vec& D2UV) const 
263 {
264   ElSLib::CylinderD2 (U, V, pos, radius, P, D1U, D1V, D2U, D2V, D2UV);
265 }
266
267
268
269 //=======================================================================
270 //function : D3
271 //purpose  : 
272 //=======================================================================
273
274 void Geom_CylindricalSurface::D3 (const Standard_Real U, const Standard_Real V,
275                                   Pnt& P, 
276                                   Vec& D1U, Vec& D1V, 
277                                   Vec& D2U, Vec& D2V, Vec& D2UV,
278                                   Vec& D3U, Vec& D3V, Vec& D3UUV, Vec& D3UVV) 
279                                   const 
280 {
281   ElSLib::CylinderD3 (U, V, pos, radius, P, D1U, D1V, D2U, D2V,
282                       D2UV, D3U, D3V, D3UUV, D3UVV);
283 }
284
285
286 //=======================================================================
287 //function : DN
288 //purpose  : 
289 //=======================================================================
290
291 Vec Geom_CylindricalSurface::DN (const Standard_Real    U , const Standard_Real    V , 
292                                  const Standard_Integer Nu, const Standard_Integer Nv ) const 
293 {
294   Standard_RangeError_Raise_if (Nu + Nv < 1 || Nu < 0 || Nv <0, " ");
295
296   if (Nv > 1) { return Vec (0.0, 0.0, 0.0); }
297   else {
298     return ElSLib::CylinderDN (U, V, pos, radius, Nu, Nv);
299   }
300 }
301
302
303
304
305 //=======================================================================
306 //function : UIso
307 //purpose  : 
308 //=======================================================================
309
310 Handle(Geom_Curve) Geom_CylindricalSurface::UIso (const Standard_Real U) const 
311 {
312   Handle(Geom_Line) GL = new Geom_Line(ElSLib::CylinderUIso(pos,radius,U));
313   return GL;
314 }
315
316
317
318 //=======================================================================
319 //function : VIso
320 //purpose  : 
321 //=======================================================================
322
323 Handle(Geom_Curve) Geom_CylindricalSurface::VIso (const Standard_Real V) const 
324 {
325   Handle(Geom_Circle) GC = new Geom_Circle(ElSLib::CylinderVIso(pos,radius,V));
326   return GC;
327 }
328
329
330
331 //=======================================================================
332 //function : Transform
333 //purpose  : 
334 //=======================================================================
335
336 void Geom_CylindricalSurface::Transform (const Trsf& T) {
337
338    radius = radius * Abs(T.ScaleFactor());
339    pos.Transform (T);
340 }
341
342 //=======================================================================
343 //function : TransformParameters
344 //purpose  : 
345 //=======================================================================
346
347 void Geom_CylindricalSurface::TransformParameters(Standard_Real& ,
348                                                   Standard_Real& V,
349                                                   const gp_Trsf& T) 
350 const
351 {
352   if (!Precision::IsInfinite(V)) V *= Abs(T.ScaleFactor());
353 }
354
355 //=======================================================================
356 //function : ParametricTransformation
357 //purpose  : 
358 //=======================================================================
359
360 gp_GTrsf2d Geom_CylindricalSurface::ParametricTransformation(const gp_Trsf& T)
361 const
362 {
363   gp_GTrsf2d T2;
364   gp_Ax2d Axis(gp::Origin2d(),gp::DX2d());
365   T2.SetAffinity(Axis, Abs(T.ScaleFactor()));
366   return T2;
367 }