b7fe2406a64b6ad3620ab8218f27052407631e69
[occt.git] / src / Geom / Geom_ConicalSurface.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_ConicalSurface.ixx>
18
19 #include <Standard_ConstructionError.hxx>
20 #include <Standard_RangeError.hxx>
21
22 #include <Precision.hxx>
23 #include <gp.hxx>
24 #include <gp_Lin.hxx>
25 #include <gp_Circ.hxx>
26 #include <gp_Ax2d.hxx>
27 #include <gp_XYZ.hxx>
28 #include <gp_Dir.hxx>
29 #include <ElSLib.hxx>
30 #include <GeomAbs_UVSense.hxx>
31 #include <Geom_Line.hxx>
32 #include <Geom_Circle.hxx>
33
34 typedef Geom_ConicalSurface         ConicalSurface;
35 typedef Handle(Geom_ConicalSurface) Handle(ConicalSurface);
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 //=======================================================================
51 //function : Copy
52 //purpose  : 
53 //=======================================================================
54
55 Handle(Geom_Geometry) Geom_ConicalSurface::Copy () const {
56  
57    Handle(ConicalSurface) Cs;
58    Cs = new ConicalSurface (pos, semiAngle, radius);
59    return Cs;
60 }
61
62 //=======================================================================
63 //function : Geom_ConicalSurface
64 //purpose  : 
65 //=======================================================================
66
67 Geom_ConicalSurface::Geom_ConicalSurface ( const Ax3& A3 , 
68                                            const Standard_Real Ang, 
69                                            const Standard_Real R) :
70        radius(R), semiAngle (Ang) 
71 {
72
73   if (R < 0.0 || Abs(Ang) <= gp::Resolution() || Abs(Ang) >= M_PI/2.0 - gp::Resolution()) 
74     Standard_ConstructionError::Raise();
75   
76   pos = A3;
77 }
78
79
80 //=======================================================================
81 //function : Geom_ConicalSurface
82 //purpose  : 
83 //=======================================================================
84
85 Geom_ConicalSurface::Geom_ConicalSurface ( const gp_Cone& C ) 
86 : radius (C.RefRadius()), semiAngle (C.SemiAngle()) 
87 {
88    pos = C.Position();
89 }
90
91
92 //=======================================================================
93 //function : UReversedParameter
94 //purpose  : 
95 //=======================================================================
96
97 Standard_Real Geom_ConicalSurface::UReversedParameter( const Standard_Real U) const
98 {
99   return ( 2.*M_PI - U);
100 }
101
102
103 //=======================================================================
104 //function : VReversedParameter
105 //purpose  : 
106 //=======================================================================
107
108 Standard_Real Geom_ConicalSurface::VReversedParameter( const Standard_Real V) const
109 {
110   return ( -V);
111 }
112
113 //=======================================================================
114 //function : VReverse
115 //purpose  : 
116 //=======================================================================
117
118 void Geom_ConicalSurface::VReverse()
119 {
120   semiAngle = -semiAngle;
121   pos.ZReverse();
122 }
123
124 //=======================================================================
125 //function : RefRadius
126 //purpose  : 
127 //=======================================================================
128
129 Standard_Real Geom_ConicalSurface::RefRadius () const               
130 { return radius; }
131
132 //=======================================================================
133 //function : SemiAngle
134 //purpose  : 
135 //=======================================================================
136
137 Standard_Real Geom_ConicalSurface::SemiAngle () const               
138 { return semiAngle;}
139
140 //=======================================================================
141 //function : IsUClosed
142 //purpose  : 
143 //=======================================================================
144
145 Standard_Boolean Geom_ConicalSurface::IsUClosed () const            
146 { return Standard_True; }
147
148 //=======================================================================
149 //function : IsVClosed
150 //purpose  : 
151 //=======================================================================
152
153 Standard_Boolean Geom_ConicalSurface::IsVClosed () const            
154 { return Standard_False; }
155
156 //=======================================================================
157 //function : IsUPeriodic
158 //purpose  : 
159 //=======================================================================
160
161 Standard_Boolean Geom_ConicalSurface::IsUPeriodic () const          
162 { return Standard_True; }
163
164 //=======================================================================
165 //function : IsVPeriodic
166 //purpose  : 
167 //=======================================================================
168
169 Standard_Boolean Geom_ConicalSurface::IsVPeriodic () const          
170 { return Standard_False; }
171
172 //=======================================================================
173 //function : Cone
174 //purpose  : 
175 //=======================================================================
176
177 gp_Cone Geom_ConicalSurface::Cone () const {
178
179   return gp_Cone (pos, semiAngle, radius);
180 }
181
182
183 //=======================================================================
184 //function : SetCone
185 //purpose  : 
186 //=======================================================================
187
188 void Geom_ConicalSurface::SetCone (const gp_Cone& C) {
189
190   radius = C.RefRadius ();
191   semiAngle   = C.SemiAngle ();
192   pos    = C.Position  ();
193 }
194
195
196 //=======================================================================
197 //function : SetRadius
198 //purpose  : 
199 //=======================================================================
200
201 void Geom_ConicalSurface::SetRadius (const Standard_Real R) {
202
203   if (R < 0.0)  Standard_ConstructionError::Raise();
204   radius = R;
205 }
206
207
208 //=======================================================================
209 //function : SetSemiAngle
210 //purpose  : 
211 //=======================================================================
212
213 void Geom_ConicalSurface::SetSemiAngle (const Standard_Real Ang) {
214
215   if (Abs(Ang) <= gp::Resolution() || Abs(Ang) >= M_PI/2.0 - gp::Resolution()) {
216     Standard_ConstructionError::Raise();
217   }
218   semiAngle = Ang;
219 }
220
221
222 //=======================================================================
223 //function : Apex
224 //purpose  : 
225 //=======================================================================
226
227 Pnt Geom_ConicalSurface::Apex () const 
228 {
229
230    XYZ Coord = Position().Direction().XYZ();
231    Coord.Multiply (-radius / Tan (semiAngle));
232    Coord.Add      (Position().Location().XYZ());
233    return Pnt     (Coord);
234 }
235
236
237 //=======================================================================
238 //function : Bounds
239 //purpose  : 
240 //=======================================================================
241
242 void Geom_ConicalSurface::Bounds (Standard_Real& U1, Standard_Real& U2, 
243                                   Standard_Real& V1, Standard_Real& V2) const {
244
245    U1 = 0.0;  U2 = 2.0 * M_PI;  
246    V1 = -Precision::Infinite();  V2 = Precision::Infinite();
247 }
248
249
250 //=======================================================================
251 //function : Coefficients
252 //purpose  : 
253 //=======================================================================
254
255 void Geom_ConicalSurface::Coefficients (Standard_Real& A1, Standard_Real& A2, Standard_Real& A3,
256                                         Standard_Real& B1, Standard_Real& B2, Standard_Real& B3,
257                                         Standard_Real& C1, Standard_Real& C2, Standard_Real& C3, 
258                                         Standard_Real& D)  const 
259 {
260    // Dans le repere du cone :
261    // X**2 + Y**2 - (Myradius - Z * Tan(semiAngle))**2 = 0.0
262
263       Trsf T;
264       T.SetTransformation (pos);
265       Standard_Real KAng = Tan (semiAngle);
266       Standard_Real T11 = T.Value (1, 1);
267       Standard_Real T12 = T.Value (1, 2);
268       Standard_Real T13 = T.Value (1, 3);
269       Standard_Real T14 = T.Value (1, 4);
270       Standard_Real T21 = T.Value (2, 1);
271       Standard_Real T22 = T.Value (2, 2);
272       Standard_Real T23 = T.Value (2, 3);
273       Standard_Real T24 = T.Value (2, 4);
274       Standard_Real T31 = T.Value (3, 1) * KAng;
275       Standard_Real T32 = T.Value (3, 2) * KAng;
276       Standard_Real T33 = T.Value (3, 3) * KAng;
277       Standard_Real T34 = T.Value (3, 4) * KAng;
278       A1 = T11 * T11 + T21 * T21 - T31 * T31;
279       A2 = T12 * T12 + T22 * T22 - T32 * T32;
280       A3 = T13 * T13 + T23 * T23 - T33 * T33;
281       B1 = T11 * T12 + T21 * T22 - T31 * T32;
282       B2 = T11 * T13 + T21 * T23 - T31 * T33;
283       B3 = T12 * T13 + T22 * T23 - T32 * T33;
284       C1 = T11 * T14 + T21 * T24 + radius * T31;
285       C2 = T12 * T14 + T22 * T24 + radius * T32;
286       C3 = T13 * T14 + T23 * T24 + radius * T33;
287       D = T14 * T14 + T24 * T24 - radius * radius - T34 * T34 +
288           2.0 * radius * T34;  
289 }
290
291
292
293 //=======================================================================
294 //function : D0
295 //purpose  : 
296 //=======================================================================
297
298 void Geom_ConicalSurface::D0 (const Standard_Real U, const Standard_Real V, Pnt& P) const 
299 {
300
301   P = ElSLib::ConeValue (U, V, pos, radius, semiAngle);
302 }
303
304
305 //=======================================================================
306 //function : D1
307 //purpose  : 
308 //=======================================================================
309
310 void Geom_ConicalSurface::D1 (const Standard_Real U  , const Standard_Real V, 
311                                     Pnt& P  , 
312                                     Vec& D1U, Vec& D1V     ) const 
313 {
314   ElSLib::ConeD1 (U, V, pos, radius, semiAngle, P, D1U, D1V);
315 }
316
317
318 //=======================================================================
319 //function : D2
320 //purpose  : 
321 //=======================================================================
322
323 void Geom_ConicalSurface::D2 ( const Standard_Real U  , const Standard_Real V, 
324                                      Pnt& P  ,  
325                                      Vec& D1U, Vec& D1V, 
326                                      Vec& D2U, Vec& D2V, Vec& D2UV) const 
327 {
328   ElSLib::ConeD2 (U, V, pos, radius, semiAngle, P, D1U, D1V, 
329                   D2U, D2V, D2UV);
330 }
331
332
333 //=======================================================================
334 //function : D3
335 //purpose  : 
336 //=======================================================================
337
338 void Geom_ConicalSurface::D3 (const Standard_Real U, const Standard_Real V, 
339                               Pnt& P      , 
340                               Vec& D1U    , Vec& D1V,
341                               Vec& D2U    , Vec& D2V, Vec& D2UV,
342                               Vec& D3U    , Vec& D3V, Vec& D3UUV, Vec& D3UVV
343                               ) const
344 {
345   ElSLib::ConeD3 (U, V, pos, radius, semiAngle, P, D1U, D1V, D2U, D2V,
346                   D2UV, D3U, D3V, D3UUV, D3UVV);
347 }
348
349
350 //=======================================================================
351 //function : DN
352 //purpose  : 
353 //=======================================================================
354
355 Vec Geom_ConicalSurface::DN (const Standard_Real    U , const Standard_Real     V, 
356                              const Standard_Integer Nu, const Standard_Integer Nv ) const
357 {
358   Standard_RangeError_Raise_if (Nu + Nv < 1 || Nu < 0 || Nv < 0, " ");
359   if (Nv > 1) { return Vec (0.0, 0.0, 0.0); }
360   else {
361     return ElSLib::ConeDN (U, V, pos, radius, semiAngle, Nu, Nv);
362   }
363 }
364
365
366 //=======================================================================
367 //function : UIso
368 //purpose  : 
369 //=======================================================================
370
371 Handle(Geom_Curve) Geom_ConicalSurface::UIso (const Standard_Real U) const 
372 {
373   Handle(Geom_Line) 
374     GL = new Geom_Line(ElSLib::ConeUIso(pos,radius,semiAngle,U));
375   return GL;
376 }
377
378
379 //=======================================================================
380 //function : VIso
381 //purpose  : 
382 //=======================================================================
383
384 Handle(Geom_Curve) Geom_ConicalSurface::VIso (const Standard_Real V) const 
385 {
386   Handle(Geom_Circle) 
387     GC = new Geom_Circle(ElSLib::ConeVIso(pos,radius,semiAngle,V));
388   return GC;
389 }
390
391
392 //=======================================================================
393 //function : Transform
394 //purpose  : 
395 //=======================================================================
396
397 void Geom_ConicalSurface::Transform (const Trsf& T) 
398 {
399   radius = radius * Abs(T.ScaleFactor());
400   pos.Transform (T);
401 }
402
403 //=======================================================================
404 //function : TransformParameters
405 //purpose  : 
406 //=======================================================================
407
408 void Geom_ConicalSurface::TransformParameters(Standard_Real& ,
409                                               Standard_Real& V,
410                                               const gp_Trsf& T) 
411 const
412 {
413   if (!Precision::IsInfinite(V)) V *= Abs(T.ScaleFactor());
414 }
415
416 //=======================================================================
417 //function : ParametricTransformation
418 //purpose  : 
419 //=======================================================================
420
421 gp_GTrsf2d Geom_ConicalSurface::ParametricTransformation(const gp_Trsf& T)
422 const
423 {
424   gp_GTrsf2d T2;
425   gp_Ax2d Axis(gp::Origin2d(),gp::DX2d());
426   T2.SetAffinity(Axis, Abs(T.ScaleFactor()));
427   return T2;
428 }
429