]> SALOME platform Git repositories - modules/shaper.git/blob - src/GeomAPI/GeomAPI_Ellipse.cpp
Salome HOME
Improve code coverage for elliptic arcs.
[modules/shaper.git] / src / GeomAPI / GeomAPI_Ellipse.cpp
1 // Copyright (C) 2017-2019  CEA/DEN, EDF R&D
2 //
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.
7 //
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.
12 //
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
16 //
17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
18 //
19
20 // File:        GeomAPI_Ellipse.cpp
21 // Created:     25 April 2017
22 // Author:      Vitaly Smetannikov
23
24 #include "GeomAPI_Ellipse.h"
25 #include "GeomAPI_Ax2.h"
26 #include "GeomAPI_Curve.h"
27 #include "GeomAPI_Pnt.h"
28
29 #include <Geom_Ellipse.hxx>
30 #include <GeomAPI_ProjectPointOnCurve.hxx>
31 #include <GeomLib_Tool.hxx>
32 #include <gp_Elips.hxx>
33
34 #define MY_ELIPS implPtr<gp_Elips>()
35
36 GeomAPI_Ellipse::GeomAPI_Ellipse(const std::shared_ptr<GeomAPI_Ax2>& theAx2,
37                                  double theMajorRadius, double theMinorRadius)
38 : GeomAPI_Interface(new gp_Elips(theAx2->impl<gp_Ax2>(), theMajorRadius, theMinorRadius))
39 {
40 }
41
42 GeomAPI_Ellipse::GeomAPI_Ellipse(GeomCurvePtr theCurve)
43 {
44   GeomCurvePtr anUntrimmedCurve = theCurve->basisCurve();
45   Handle(Geom_Curve) aCurve = anUntrimmedCurve->impl<Handle(Geom_Curve)>();
46   Handle(Geom_Ellipse) anEllipse = Handle(Geom_Ellipse)::DownCast(aCurve);
47   if (anEllipse.IsNull())
48     throw Standard_ConstructionError("GeomAPI_Ellipse: Curve is not an ellipse");
49   setImpl(new gp_Elips(anEllipse->Elips()));
50 }
51
52 std::shared_ptr<GeomAPI_Pnt> GeomAPI_Ellipse::center() const
53 {
54   const gp_Pnt& aCenter = MY_ELIPS->Location();
55   return std::shared_ptr<GeomAPI_Pnt>(new GeomAPI_Pnt(aCenter.X(), aCenter.Y(), aCenter.Z()));
56 }
57
58 GeomPointPtr GeomAPI_Ellipse::firstFocus() const
59 {
60   const gp_Pnt& aFirst = MY_ELIPS->Focus1();
61   return std::shared_ptr<GeomAPI_Pnt>(new GeomAPI_Pnt(aFirst.X(), aFirst.Y(), aFirst.Z()));
62 }
63
64 GeomPointPtr GeomAPI_Ellipse::secondFocus() const
65 {
66   const gp_Pnt& aSecond = MY_ELIPS->Focus2();
67   return std::shared_ptr<GeomAPI_Pnt>(new GeomAPI_Pnt(aSecond.X(), aSecond.Y(), aSecond.Z()));
68 }
69
70 std::shared_ptr<GeomAPI_Dir> GeomAPI_Ellipse::normal() const
71 {
72   const gp_Ax1& anAxis = MY_ELIPS->Axis();
73   const gp_Dir& aDir = anAxis.Direction();
74   return std::shared_ptr<GeomAPI_Dir>(new GeomAPI_Dir(aDir.X(), aDir.Y(), aDir.Z()));
75 }
76
77 double GeomAPI_Ellipse::minorRadius() const
78 {
79   return MY_ELIPS->MinorRadius();
80 }
81
82 double GeomAPI_Ellipse::majorRadius() const
83 {
84   return MY_ELIPS->MajorRadius();
85 }
86
87 const std::shared_ptr<GeomAPI_Pnt> GeomAPI_Ellipse::project(
88     const std::shared_ptr<GeomAPI_Pnt>& thePoint) const
89 {
90   std::shared_ptr<GeomAPI_Pnt> aResult;
91   if (!MY_ELIPS)
92     return aResult;
93
94   Handle(Geom_Ellipse) aEllipse = new Geom_Ellipse(*MY_ELIPS);
95
96   const gp_Pnt& aPoint = thePoint->impl<gp_Pnt>();
97
98   GeomAPI_ProjectPointOnCurve aProj(aPoint, aEllipse);
99   Standard_Integer aNbPoint = aProj.NbPoints();
100   if (aNbPoint > 0) {
101     gp_Pnt aNearest = aProj.NearestPoint();
102     aResult = GeomPointPtr(new GeomAPI_Pnt(aNearest.X(), aNearest.Y(), aNearest.Z()));
103   }
104   return aResult;
105 }
106
107 const bool GeomAPI_Ellipse::parameter(const std::shared_ptr<GeomAPI_Pnt> thePoint,
108                                       const double theTolerance,
109                                       double& theParameter) const
110 {
111   Handle(Geom_Ellipse) aCurve = new Geom_Ellipse(*MY_ELIPS);
112   return GeomLib_Tool::Parameter(aCurve, thePoint->impl<gp_Pnt>(),
113                                  theTolerance, theParameter) == Standard_True;
114 }