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