src/Geom2d/Geom2d_Ellipse.cxx

   1 // Created on: 1993-03-24
   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 <Geom2d_Ellipse.ixx>
  18 #include <gp.hxx>
  19 #include <gp_Dir2d.hxx>
  20 #include <gp_XYZ.hxx>
  21 #include <ElCLib.hxx>
  22 #include <Standard_RangeError.hxx>
  23 #include <Standard_ConstructionError.hxx>
  24
  25 typedef Geom2d_Ellipse         Ellipse;
  26 typedef gp_Ax2d   Ax2d;
  27 typedef gp_Dir2d  Dir2d;
  28 typedef gp_Pnt2d  Pnt2d;
  29 typedef gp_Vec2d  Vec2d;
  30 typedef gp_Trsf2d Trsf2d;
  31 typedef gp_XY     XY;
  32
  33 //=======================================================================
  34 //function : Copy
  35 //purpose  :
  36 //=======================================================================
  37
  38 Handle(Geom2d_Geometry) Geom2d_Ellipse::Copy() const
  39 {
  40   Handle(Geom2d_Ellipse) E;
  41   E = new Ellipse (pos, majorRadius, minorRadius);
  42   return E;
  43 }
  44
  45
  46 //=======================================================================
  47 //function : Geom2d_Ellipse
  48 //purpose  :
  49 //=======================================================================
  50
  51 Geom2d_Ellipse::Geom2d_Ellipse (const gp_Elips2d& E) {
  52
  53   majorRadius = E.MajorRadius();
  54   minorRadius = E.MinorRadius();
  55   pos = E.Axis();
  56 }
  57
  58
  59
  60 //=======================================================================
  61 //function : Geom2d_Ellipse
  62 //purpose  :
  63 //=======================================================================
  64
  65 Geom2d_Ellipse::Geom2d_Ellipse (const Ax2d& MajorAxis,
  66                                 const Standard_Real MajorRadius,
  67                                 const Standard_Real MinorRadius,
  68                                 const Standard_Boolean Sense)
  69 :majorRadius (MajorRadius), minorRadius (MinorRadius){
  70
  71   if (MajorRadius < MinorRadius || MinorRadius < 0.0 ) {
  72     Standard_ConstructionError::Raise();
  73   }
  74   pos = gp_Ax22d(MajorAxis, Sense);
  75 }
  76
  77
  78 //=======================================================================
  79 //function : Geom2d_Ellipse
  80 //purpose  :
  81 //=======================================================================
  82
  83 Geom2d_Ellipse::Geom2d_Ellipse (const gp_Ax22d& Axis,
  84                                 const Standard_Real MajorRadius,
  85                                 const Standard_Real MinorRadius)
  86 : majorRadius (MajorRadius), minorRadius (MinorRadius)
  87 {
  88   if (MajorRadius < MinorRadius || MinorRadius < 0.0 ) {
  89     Standard_ConstructionError::Raise();
  90   }
  91   pos = Axis;
  92 }
  93
  94
  95 //=======================================================================
  96 //function : SetElips2d
  97 //purpose  :
  98 //=======================================================================
  99
 100 void Geom2d_Ellipse::SetElips2d (const gp_Elips2d& E)
 101 {
 102   majorRadius = E.MajorRadius();
 103   minorRadius = E.MinorRadius();
 104   pos = E.Axis();
 105 }
 106
 107
 108 //=======================================================================
 109 //function : SetMajorRadius
 110 //purpose  :
 111 //=======================================================================
 112
 113 void Geom2d_Ellipse::SetMajorRadius (const Standard_Real MajorRadius)
 114 {
 115   if (MajorRadius < minorRadius)
 116     Standard_ConstructionError::Raise();
 117   else
 118     majorRadius = MajorRadius;
 119 }
 120
 121
 122 //=======================================================================
 123 //function : SetMinorRadius
 124 //purpose  :
 125 //=======================================================================
 126
 127 void Geom2d_Ellipse::SetMinorRadius (const Standard_Real MinorRadius)
 128 {
 129    if (MinorRadius < 0 || majorRadius < MinorRadius)
 130      { Standard_ConstructionError::Raise(); }
 131    else
 132      { minorRadius = MinorRadius; }
 133 }
 134
 135
 136 //=======================================================================
 137 //function : Elips2d
 138 //purpose  :
 139 //=======================================================================
 140
 141 gp_Elips2d Geom2d_Ellipse::Elips2d () const
 142 {
 143   return gp_Elips2d (pos, majorRadius, minorRadius);
 144 }
 145
 146 //=======================================================================
 147 //function : ReversedParameter
 148 //purpose  :
 149 //=======================================================================
 150
 151 Standard_Real Geom2d_Ellipse::ReversedParameter( const Standard_Real U) const
 152 {
 153   return (2. * M_PI - U);
 154 }
 155
 156 //=======================================================================
 157 //function : Directrix1
 158 //purpose  :
 159 //=======================================================================
 160
 161 Ax2d Geom2d_Ellipse::Directrix1 () const
 162 {
 163   gp_Elips2d Ev (pos, majorRadius, minorRadius);
 164   return Ev.Directrix1();
 165 }
 166
 167
 168 //=======================================================================
 169 //function : Directrix2
 170 //purpose  :
 171 //=======================================================================
 172
 173 Ax2d Geom2d_Ellipse::Directrix2 () const
 174 {
 175   gp_Elips2d Ev (pos, majorRadius, minorRadius);
 176   return Ev.Directrix2();
 177 }
 178
 179
 180 //=======================================================================
 181 //function : Eccentricity
 182 //purpose  :
 183 //=======================================================================
 184
 185 Standard_Real Geom2d_Ellipse::Eccentricity () const
 186 {
 187   if (majorRadius == 0.0)
 188     { return 0.0; }
 189   else
 190     { return
 191         (Sqrt(majorRadius*majorRadius-minorRadius*minorRadius))/majorRadius;
 192     }
 193 }
 194
 195
 196 //=======================================================================
 197 //function : Focal
 198 //purpose  :
 199 //=======================================================================
 200
 201 Standard_Real Geom2d_Ellipse::Focal () const
 202 {
 203   return 2.0 * Sqrt(majorRadius * majorRadius - minorRadius * minorRadius);
 204 }
 205
 206
 207 //=======================================================================
 208 //function : Focus1
 209 //purpose  :
 210 //=======================================================================
 211
 212 Pnt2d Geom2d_Ellipse::Focus1 () const
 213 {
 214   Standard_Real C = Sqrt(majorRadius * majorRadius - minorRadius * minorRadius);
 215   return Pnt2d (pos.Location().X() + C * pos.XDirection().X(),
 216                 pos.Location().Y() + C * pos.XDirection().Y());
 217 }
 218
 219
 220 //=======================================================================
 221 //function : Focus2
 222 //purpose  :
 223 //=======================================================================
 224
 225 Pnt2d Geom2d_Ellipse::Focus2 () const
 226 {
 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());
 230 }
 231
 232 //=======================================================================
 233 //function : MajorRadius
 234 //purpose  :
 235 //=======================================================================
 236
 237 Standard_Real Geom2d_Ellipse::MajorRadius () const
 238 {
 239   return majorRadius;
 240 }
 241
 242 //=======================================================================
 243 //function : MinorRadius
 244 //purpose  :
 245 //=======================================================================
 246
 247 Standard_Real Geom2d_Ellipse::MinorRadius () const
 248 {
 249   return minorRadius;
 250 }
 251
 252 //=======================================================================
 253 //function : Parameter
 254 //purpose  :
 255 //=======================================================================
 256
 257 Standard_Real Geom2d_Ellipse::Parameter () const
 258 {
 259   if (majorRadius == 0.0)
 260     return 0.0;
 261   else
 262     return (minorRadius * minorRadius)/majorRadius;
 263 }
 264
 265
 266 //=======================================================================
 267 //function : FirstParameter
 268 //purpose  :
 269 //=======================================================================
 270
 271 Standard_Real Geom2d_Ellipse::FirstParameter () const
 272 {
 273   return 0.0;
 274 }
 275
 276 //=======================================================================
 277 //function : LastParameter
 278 //purpose  :
 279 //=======================================================================
 280
 281 Standard_Real Geom2d_Ellipse::LastParameter () const
 282 {
 283   return 2.0 * M_PI;
 284 }
 285
 286 //=======================================================================
 287 //function : IsClosed
 288 //purpose  :
 289 //=======================================================================
 290
 291 Standard_Boolean Geom2d_Ellipse::IsClosed () const
 292 {
 293   return Standard_True;
 294 }
 295
 296 //=======================================================================
 297 //function : IsPeriodic
 298 //purpose  :
 299 //=======================================================================
 300
 301 Standard_Boolean Geom2d_Ellipse::IsPeriodic () const
 302 {
 303   return Standard_True;
 304 }
 305
 306 //=======================================================================
 307 //function : D0
 308 //purpose  :
 309 //=======================================================================
 310
 311 void Geom2d_Ellipse::D0 (const Standard_Real U, Pnt2d& P) const
 312 {
 313   P = ElCLib::EllipseValue (U, pos, majorRadius, minorRadius);
 314 }
 315
 316 //=======================================================================
 317 //function : D1
 318 //purpose  :
 319 //=======================================================================
 320
 321 void Geom2d_Ellipse::D1 (const Standard_Real U, Pnt2d& P, Vec2d& V1) const
 322 {
 323   ElCLib::EllipseD1 (U, pos, majorRadius, minorRadius, P, V1);
 324 }
 325
 326
 327 //=======================================================================
 328 //function : D2
 329 //purpose  :
 330 //=======================================================================
 331
 332 void Geom2d_Ellipse::D2 (const Standard_Real U,
 333                                Pnt2d& P,
 334                                Vec2d& V1, Vec2d& V2) const
 335 {
 336   ElCLib::EllipseD2 (U, pos, majorRadius, minorRadius, P, V1, V2);
 337 }
 338
 339
 340 //=======================================================================
 341 //function : D3
 342 //purpose  :
 343 //=======================================================================
 344
 345 void Geom2d_Ellipse::D3 (const Standard_Real U,
 346                                Pnt2d& P,
 347                                Vec2d& V1, Vec2d& V2, Vec2d& V3) const
 348 {
 349   ElCLib::EllipseD3 (U, pos, majorRadius, minorRadius, P, V1, V2, V3);
 350 }
 351
 352
 353 //=======================================================================
 354 //function : DN
 355 //purpose  :
 356 //=======================================================================
 357
 358 Vec2d Geom2d_Ellipse::DN (const Standard_Real U, const Standard_Integer N) const
 359 {
 360   Standard_RangeError_Raise_if (N < 1, " ");
 361   return ElCLib::EllipseDN (U, pos, majorRadius, minorRadius, N);
 362 }
 363
 364 //=======================================================================
 365 //function : Transform
 366 //purpose  :
 367 //=======================================================================
 368
 369 void Geom2d_Ellipse::Transform (const Trsf2d& T)
 370 {
 371   majorRadius = majorRadius * Abs(T.ScaleFactor());
 372   minorRadius = minorRadius * Abs(T.ScaleFactor());
 373   pos.Transform(T);
 374 }