1 // Created on: 1993-03-24
3 // Copyright (c) 1993-1999 Matra Datavision
4 // Copyright (c) 1999-2012 OPEN CASCADE SAS
6 // The content of this file is subject to the Open CASCADE Technology Public
7 // License Version 6.5 (the "License"). You may not use the content of this file
8 // except in compliance with the License. Please obtain a copy of the License
9 // at http://www.opencascade.org and read it completely before using this file.
11 // The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
12 // main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
14 // The Original Code and all software distributed under the License is
15 // distributed on an "AS IS" basis, without warranty of any kind, and the
16 // Initial Developer hereby disclaims all such warranties, including without
17 // limitation, any warranties of merchantability, fitness for a particular
18 // purpose or non-infringement. Please see the License for the specific terms
19 // and conditions governing the rights and limitations under the License.
24 #include <Geom2d_Ellipse.ixx>
26 #include <gp_Dir2d.hxx>
29 #include <Standard_RangeError.hxx>
30 #include <Standard_ConstructionError.hxx>
33 typedef Geom2d_Ellipse Ellipse;
34 typedef Handle(Geom2d_Ellipse) Handle(Ellipse);
36 typedef gp_Dir2d Dir2d;
37 typedef gp_Pnt2d Pnt2d;
38 typedef gp_Vec2d Vec2d;
39 typedef gp_Trsf2d Trsf2d;
46 //=======================================================================
49 //=======================================================================
51 Handle(Geom2d_Geometry) Geom2d_Ellipse::Copy() const
54 E = new Ellipse (pos, majorRadius, minorRadius);
59 //=======================================================================
60 //function : Geom2d_Ellipse
62 //=======================================================================
64 Geom2d_Ellipse::Geom2d_Ellipse (const gp_Elips2d& E) {
66 majorRadius = E.MajorRadius();
67 minorRadius = E.MinorRadius();
73 //=======================================================================
74 //function : Geom2d_Ellipse
76 //=======================================================================
78 Geom2d_Ellipse::Geom2d_Ellipse (const Ax2d& MajorAxis,
79 const Standard_Real MajorRadius,
80 const Standard_Real MinorRadius,
81 const Standard_Boolean Sense)
82 :majorRadius (MajorRadius), minorRadius (MinorRadius){
84 if (MajorRadius < MinorRadius || MinorRadius < 0.0 ) {
85 Standard_ConstructionError::Raise();
87 pos = gp_Ax22d(MajorAxis, Sense);
91 //=======================================================================
92 //function : Geom2d_Ellipse
94 //=======================================================================
96 Geom2d_Ellipse::Geom2d_Ellipse (const gp_Ax22d& Axis,
97 const Standard_Real MajorRadius,
98 const Standard_Real MinorRadius)
99 : majorRadius (MajorRadius), minorRadius (MinorRadius)
101 if (MajorRadius < MinorRadius || MinorRadius < 0.0 ) {
102 Standard_ConstructionError::Raise();
108 //=======================================================================
109 //function : SetElips2d
111 //=======================================================================
113 void Geom2d_Ellipse::SetElips2d (const gp_Elips2d& E)
115 majorRadius = E.MajorRadius();
116 minorRadius = E.MinorRadius();
121 //=======================================================================
122 //function : SetMajorRadius
124 //=======================================================================
126 void Geom2d_Ellipse::SetMajorRadius (const Standard_Real MajorRadius)
128 if (MajorRadius < minorRadius)
129 Standard_ConstructionError::Raise();
131 majorRadius = MajorRadius;
135 //=======================================================================
136 //function : SetMinorRadius
138 //=======================================================================
140 void Geom2d_Ellipse::SetMinorRadius (const Standard_Real MinorRadius)
142 if (MinorRadius < 0 || majorRadius < MinorRadius)
143 { Standard_ConstructionError::Raise(); }
145 { minorRadius = MinorRadius; }
149 //=======================================================================
152 //=======================================================================
154 gp_Elips2d Geom2d_Ellipse::Elips2d () const
156 return gp_Elips2d (pos, majorRadius, minorRadius);
159 //=======================================================================
160 //function : ReversedParameter
162 //=======================================================================
164 Standard_Real Geom2d_Ellipse::ReversedParameter( const Standard_Real U) const
166 return (2. * M_PI - U);
169 //=======================================================================
170 //function : Directrix1
172 //=======================================================================
174 Ax2d Geom2d_Ellipse::Directrix1 () const
176 gp_Elips2d Ev (pos, majorRadius, minorRadius);
177 return Ev.Directrix1();
181 //=======================================================================
182 //function : Directrix2
184 //=======================================================================
186 Ax2d Geom2d_Ellipse::Directrix2 () const
188 gp_Elips2d Ev (pos, majorRadius, minorRadius);
189 return Ev.Directrix2();
193 //=======================================================================
194 //function : Eccentricity
196 //=======================================================================
198 Standard_Real Geom2d_Ellipse::Eccentricity () const
200 if (majorRadius == 0.0)
204 (Sqrt(majorRadius*majorRadius-minorRadius*minorRadius))/majorRadius;
209 //=======================================================================
212 //=======================================================================
214 Standard_Real Geom2d_Ellipse::Focal () const
216 return 2.0 * Sqrt(majorRadius * majorRadius - minorRadius * minorRadius);
220 //=======================================================================
223 //=======================================================================
225 Pnt2d Geom2d_Ellipse::Focus1 () const
227 Standard_Real C = Sqrt(majorRadius * majorRadius - minorRadius * minorRadius);
228 return Pnt2d (pos.Location().X() + C * pos.XDirection().X(),
229 pos.Location().Y() + C * pos.XDirection().Y());
233 //=======================================================================
236 //=======================================================================
238 Pnt2d Geom2d_Ellipse::Focus2 () const
240 Standard_Real C = Sqrt(majorRadius * majorRadius - minorRadius * minorRadius);
241 return Pnt2d (pos.Location().X() - C * pos.XDirection().X(),
242 pos.Location().Y() - C * pos.XDirection().Y());
245 //=======================================================================
246 //function : MajorRadius
248 //=======================================================================
250 Standard_Real Geom2d_Ellipse::MajorRadius () const
255 //=======================================================================
256 //function : MinorRadius
258 //=======================================================================
260 Standard_Real Geom2d_Ellipse::MinorRadius () const
265 //=======================================================================
266 //function : Parameter
268 //=======================================================================
270 Standard_Real Geom2d_Ellipse::Parameter () const
272 if (majorRadius == 0.0)
275 return (minorRadius * minorRadius)/majorRadius;
279 //=======================================================================
280 //function : FirstParameter
282 //=======================================================================
284 Standard_Real Geom2d_Ellipse::FirstParameter () const
289 //=======================================================================
290 //function : LastParameter
292 //=======================================================================
294 Standard_Real Geom2d_Ellipse::LastParameter () const
299 //=======================================================================
300 //function : IsClosed
302 //=======================================================================
304 Standard_Boolean Geom2d_Ellipse::IsClosed () const
306 return Standard_True;
309 //=======================================================================
310 //function : IsPeriodic
312 //=======================================================================
314 Standard_Boolean Geom2d_Ellipse::IsPeriodic () const
316 return Standard_True;
319 //=======================================================================
322 //=======================================================================
324 void Geom2d_Ellipse::D0 (const Standard_Real U, Pnt2d& P) const
326 P = ElCLib::EllipseValue (U, pos, majorRadius, minorRadius);
329 //=======================================================================
332 //=======================================================================
334 void Geom2d_Ellipse::D1 (const Standard_Real U, Pnt2d& P, Vec2d& V1) const
336 ElCLib::EllipseD1 (U, pos, majorRadius, minorRadius, P, V1);
340 //=======================================================================
343 //=======================================================================
345 void Geom2d_Ellipse::D2 (const Standard_Real U,
347 Vec2d& V1, Vec2d& V2) const
349 ElCLib::EllipseD2 (U, pos, majorRadius, minorRadius, P, V1, V2);
353 //=======================================================================
356 //=======================================================================
358 void Geom2d_Ellipse::D3 (const Standard_Real U,
360 Vec2d& V1, Vec2d& V2, Vec2d& V3) const
362 ElCLib::EllipseD3 (U, pos, majorRadius, minorRadius, P, V1, V2, V3);
366 //=======================================================================
369 //=======================================================================
371 Vec2d Geom2d_Ellipse::DN (const Standard_Real U, const Standard_Integer N) const
373 Standard_RangeError_Raise_if (N < 1, " ");
374 return ElCLib::EllipseDN (U, pos, majorRadius, minorRadius, N);
377 //=======================================================================
378 //function : Transform
380 //=======================================================================
382 void Geom2d_Ellipse::Transform (const Trsf2d& T)
384 majorRadius = majorRadius * Abs(T.ScaleFactor());
385 minorRadius = minorRadius * Abs(T.ScaleFactor());