1 // Copyright (C) 2017-2019 CEA/DEN, EDF R&D
3 // This library is free software; you can redistribute it and/or
4 // modify it under the terms of the GNU Lesser General Public
5 // License as published by the Free Software Foundation; either
6 // version 2.1 of the License, or (at your option) any later version.
8 // This library is distributed in the hope that it will be useful,
9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
11 // Lesser General Public License for more details.
13 // You should have received a copy of the GNU Lesser General Public
14 // License along with this library; if not, write to the Free Software
15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
20 // File: GeomAPI_Ellipse.cpp
21 // Created: 25 April 2017
22 // Author: Vitaly Smetannikov
24 #include "GeomAPI_Ellipse.h"
25 #include "GeomAPI_Ax2.h"
26 #include "GeomAPI_Curve.h"
27 #include "GeomAPI_Pnt.h"
29 #include <Geom_Ellipse.hxx>
30 #include <GeomAPI_ProjectPointOnCurve.hxx>
31 #include <gp_Elips.hxx>
33 #define MY_ELIPS implPtr<gp_Elips>()
35 GeomAPI_Ellipse::GeomAPI_Ellipse(const std::shared_ptr<GeomAPI_Ax2>& theAx2,
36 double theMajorRadius, double theMinorRadius)
37 : GeomAPI_Interface(new gp_Elips(theAx2->impl<gp_Ax2>(), theMajorRadius, theMinorRadius))
41 GeomAPI_Ellipse::GeomAPI_Ellipse(GeomCurvePtr theCurve)
43 Handle(Geom_Curve) aCurve = theCurve->impl<Handle(Geom_Curve)>();
44 Handle(Geom_Ellipse) anEllipse = Handle(Geom_Ellipse)::DownCast(aCurve);
45 if (anEllipse.IsNull())
46 throw Standard_ConstructionError("GeomAPI_Ellipse: Curve is not an ellipse");
47 setImpl(new gp_Elips(anEllipse->Elips()));
50 std::shared_ptr<GeomAPI_Pnt> GeomAPI_Ellipse::center() const
52 const gp_Pnt& aCenter = MY_ELIPS->Location();
53 return std::shared_ptr<GeomAPI_Pnt>(new GeomAPI_Pnt(aCenter.X(), aCenter.Y(), aCenter.Z()));
56 GeomPointPtr GeomAPI_Ellipse::firstFocus() const
58 const gp_Pnt& aFirst = MY_ELIPS->Focus1();
59 return std::shared_ptr<GeomAPI_Pnt>(new GeomAPI_Pnt(aFirst.X(), aFirst.Y(), aFirst.Z()));
62 GeomPointPtr GeomAPI_Ellipse::secondFocus() const
64 const gp_Pnt& aSecond = MY_ELIPS->Focus2();
65 return std::shared_ptr<GeomAPI_Pnt>(new GeomAPI_Pnt(aSecond.X(), aSecond.Y(), aSecond.Z()));
68 std::shared_ptr<GeomAPI_Dir> GeomAPI_Ellipse::normal() const
70 const gp_Ax1& anAxis = MY_ELIPS->Axis();
71 const gp_Dir& aDir = anAxis.Direction();
72 return std::shared_ptr<GeomAPI_Dir>(new GeomAPI_Dir(aDir.X(), aDir.Y(), aDir.Z()));
75 double GeomAPI_Ellipse::minorRadius() const
77 return MY_ELIPS->MinorRadius();
80 double GeomAPI_Ellipse::majorRadius() const
82 return MY_ELIPS->MajorRadius();
85 const std::shared_ptr<GeomAPI_Pnt> GeomAPI_Ellipse::project(
86 const std::shared_ptr<GeomAPI_Pnt>& thePoint) const
88 std::shared_ptr<GeomAPI_Pnt> aResult;
92 Handle(Geom_Ellipse) aEllipse = new Geom_Ellipse(*MY_ELIPS);
94 const gp_Pnt& aPoint = thePoint->impl<gp_Pnt>();
96 GeomAPI_ProjectPointOnCurve aProj(aPoint, aEllipse);
97 Standard_Integer aNbPoint = aProj.NbPoints();
99 gp_Pnt aNearest = aProj.NearestPoint();
100 aResult = GeomPointPtr(new GeomAPI_Pnt(aNearest.X(), aNearest.Y(), aNearest.Z()));