1 // Created on: 1993-03-24
3 // Copyright (c) 1993-1999 Matra Datavision
4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
6 // This file is part of Open CASCADE Technology software library.
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.
14 // Alternatively, this file may be used under the terms of Open CASCADE
15 // commercial license or contractual agreement.
19 #include <Geom2d_Ellipse.hxx>
20 #include <Geom2d_Geometry.hxx>
22 #include <gp_Ax2d.hxx>
23 #include <gp_Ax22d.hxx>
24 #include <gp_Dir2d.hxx>
25 #include <gp_Elips2d.hxx>
26 #include <gp_Pnt2d.hxx>
27 #include <gp_Trsf2d.hxx>
28 #include <gp_Vec2d.hxx>
30 #include <Standard_ConstructionError.hxx>
31 #include <Standard_RangeError.hxx>
32 #include <Standard_Type.hxx>
34 typedef Geom2d_Ellipse Ellipse;
36 typedef gp_Dir2d Dir2d;
37 typedef gp_Pnt2d Pnt2d;
38 typedef gp_Vec2d Vec2d;
39 typedef gp_Trsf2d Trsf2d;
42 //=======================================================================
45 //=======================================================================
47 Handle(Geom2d_Geometry) Geom2d_Ellipse::Copy() const
49 Handle(Geom2d_Ellipse) E;
50 E = new Ellipse (pos, majorRadius, minorRadius);
55 //=======================================================================
56 //function : Geom2d_Ellipse
58 //=======================================================================
60 Geom2d_Ellipse::Geom2d_Ellipse (const gp_Elips2d& E) {
62 majorRadius = E.MajorRadius();
63 minorRadius = E.MinorRadius();
69 //=======================================================================
70 //function : Geom2d_Ellipse
72 //=======================================================================
74 Geom2d_Ellipse::Geom2d_Ellipse (const Ax2d& MajorAxis,
75 const Standard_Real MajorRadius,
76 const Standard_Real MinorRadius,
77 const Standard_Boolean Sense)
78 :majorRadius (MajorRadius), minorRadius (MinorRadius){
80 if (MajorRadius < MinorRadius || MinorRadius < 0.0 ) {
81 Standard_ConstructionError::Raise();
83 pos = gp_Ax22d(MajorAxis, Sense);
87 //=======================================================================
88 //function : Geom2d_Ellipse
90 //=======================================================================
92 Geom2d_Ellipse::Geom2d_Ellipse (const gp_Ax22d& Axis,
93 const Standard_Real MajorRadius,
94 const Standard_Real MinorRadius)
95 : majorRadius (MajorRadius), minorRadius (MinorRadius)
97 if (MajorRadius < MinorRadius || MinorRadius < 0.0 ) {
98 Standard_ConstructionError::Raise();
104 //=======================================================================
105 //function : SetElips2d
107 //=======================================================================
109 void Geom2d_Ellipse::SetElips2d (const gp_Elips2d& E)
111 majorRadius = E.MajorRadius();
112 minorRadius = E.MinorRadius();
117 //=======================================================================
118 //function : SetMajorRadius
120 //=======================================================================
122 void Geom2d_Ellipse::SetMajorRadius (const Standard_Real MajorRadius)
124 if (MajorRadius < minorRadius)
125 Standard_ConstructionError::Raise();
127 majorRadius = MajorRadius;
131 //=======================================================================
132 //function : SetMinorRadius
134 //=======================================================================
136 void Geom2d_Ellipse::SetMinorRadius (const Standard_Real MinorRadius)
138 if (MinorRadius < 0 || majorRadius < MinorRadius)
139 { Standard_ConstructionError::Raise(); }
141 { minorRadius = MinorRadius; }
145 //=======================================================================
148 //=======================================================================
150 gp_Elips2d Geom2d_Ellipse::Elips2d () const
152 return gp_Elips2d (pos, majorRadius, minorRadius);
155 //=======================================================================
156 //function : ReversedParameter
158 //=======================================================================
160 Standard_Real Geom2d_Ellipse::ReversedParameter( const Standard_Real U) const
162 return (2. * M_PI - U);
165 //=======================================================================
166 //function : Directrix1
168 //=======================================================================
170 Ax2d Geom2d_Ellipse::Directrix1 () const
172 gp_Elips2d Ev (pos, majorRadius, minorRadius);
173 return Ev.Directrix1();
177 //=======================================================================
178 //function : Directrix2
180 //=======================================================================
182 Ax2d Geom2d_Ellipse::Directrix2 () const
184 gp_Elips2d Ev (pos, majorRadius, minorRadius);
185 return Ev.Directrix2();
189 //=======================================================================
190 //function : Eccentricity
192 //=======================================================================
194 Standard_Real Geom2d_Ellipse::Eccentricity () const
196 if (majorRadius == 0.0)
200 (Sqrt(majorRadius*majorRadius-minorRadius*minorRadius))/majorRadius;
205 //=======================================================================
208 //=======================================================================
210 Standard_Real Geom2d_Ellipse::Focal () const
212 return 2.0 * Sqrt(majorRadius * majorRadius - minorRadius * minorRadius);
216 //=======================================================================
219 //=======================================================================
221 Pnt2d Geom2d_Ellipse::Focus1 () const
223 Standard_Real C = Sqrt(majorRadius * majorRadius - minorRadius * minorRadius);
224 return Pnt2d (pos.Location().X() + C * pos.XDirection().X(),
225 pos.Location().Y() + C * pos.XDirection().Y());
229 //=======================================================================
232 //=======================================================================
234 Pnt2d Geom2d_Ellipse::Focus2 () const
236 Standard_Real C = Sqrt(majorRadius * majorRadius - minorRadius * minorRadius);
237 return Pnt2d (pos.Location().X() - C * pos.XDirection().X(),
238 pos.Location().Y() - C * pos.XDirection().Y());
241 //=======================================================================
242 //function : MajorRadius
244 //=======================================================================
246 Standard_Real Geom2d_Ellipse::MajorRadius () const
251 //=======================================================================
252 //function : MinorRadius
254 //=======================================================================
256 Standard_Real Geom2d_Ellipse::MinorRadius () const
261 //=======================================================================
262 //function : Parameter
264 //=======================================================================
266 Standard_Real Geom2d_Ellipse::Parameter () const
268 if (majorRadius == 0.0)
271 return (minorRadius * minorRadius)/majorRadius;
275 //=======================================================================
276 //function : FirstParameter
278 //=======================================================================
280 Standard_Real Geom2d_Ellipse::FirstParameter () const
285 //=======================================================================
286 //function : LastParameter
288 //=======================================================================
290 Standard_Real Geom2d_Ellipse::LastParameter () const
295 //=======================================================================
296 //function : IsClosed
298 //=======================================================================
300 Standard_Boolean Geom2d_Ellipse::IsClosed () const
302 return Standard_True;
305 //=======================================================================
306 //function : IsPeriodic
308 //=======================================================================
310 Standard_Boolean Geom2d_Ellipse::IsPeriodic () const
312 return Standard_True;
315 //=======================================================================
318 //=======================================================================
320 void Geom2d_Ellipse::D0 (const Standard_Real U, Pnt2d& P) const
322 P = ElCLib::EllipseValue (U, pos, majorRadius, minorRadius);
325 //=======================================================================
328 //=======================================================================
330 void Geom2d_Ellipse::D1 (const Standard_Real U, Pnt2d& P, Vec2d& V1) const
332 ElCLib::EllipseD1 (U, pos, majorRadius, minorRadius, P, V1);
336 //=======================================================================
339 //=======================================================================
341 void Geom2d_Ellipse::D2 (const Standard_Real U,
343 Vec2d& V1, Vec2d& V2) const
345 ElCLib::EllipseD2 (U, pos, majorRadius, minorRadius, P, V1, V2);
349 //=======================================================================
352 //=======================================================================
354 void Geom2d_Ellipse::D3 (const Standard_Real U,
356 Vec2d& V1, Vec2d& V2, Vec2d& V3) const
358 ElCLib::EllipseD3 (U, pos, majorRadius, minorRadius, P, V1, V2, V3);
362 //=======================================================================
365 //=======================================================================
367 Vec2d Geom2d_Ellipse::DN (const Standard_Real U, const Standard_Integer N) const
369 Standard_RangeError_Raise_if (N < 1, " ");
370 return ElCLib::EllipseDN (U, pos, majorRadius, minorRadius, N);
373 //=======================================================================
374 //function : Transform
376 //=======================================================================
378 void Geom2d_Ellipse::Transform (const Trsf2d& T)
380 majorRadius = majorRadius * Abs(T.ScaleFactor());
381 minorRadius = minorRadius * Abs(T.ScaleFactor());