X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;f=src%2FGeomAPI%2FGeomAPI_Edge.cpp;h=9927c2ae9cc3bb899fcb53611073f05d3d67b436;hb=4df4bd61da1ea5d357671c819a8ced6ec9ba77ac;hp=b5736f4b618d70b52ad2e8c786c8db1ba4d20ff6;hpb=a24b7e6f4d112d5e7889fd76f030298fc428cd01;p=modules%2Fshaper.git diff --git a/src/GeomAPI/GeomAPI_Edge.cpp b/src/GeomAPI/GeomAPI_Edge.cpp index b5736f4b6..9927c2ae9 100644 --- a/src/GeomAPI/GeomAPI_Edge.cpp +++ b/src/GeomAPI/GeomAPI_Edge.cpp @@ -1,28 +1,58 @@ -// Copyright (C) 2014-20xx CEA/DEN, EDF R&D - -// File: GeomAPI_Edge.cpp -// Created: 24 Jul 2014 -// Author: Artem ZHIDKOV +// Copyright (C) 2014-2017 CEA/DEN, EDF R&D +// +// This library is free software; you can redistribute it and/or +// modify it under the terms of the GNU Lesser General Public +// License as published by the Free Software Foundation; either +// version 2.1 of the License, or (at your option) any later version. +// +// This library is distributed in the hope that it will be useful, +// but WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +// Lesser General Public License for more details. +// +// You should have received a copy of the GNU Lesser General Public +// License along with this library; if not, write to the Free Software +// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA +// +// See http://www.salome-platform.org/ or +// email : webmaster.salome@opencascade.com +// #include +#include #include #include #include #include +#include +#include + +#include #include #include #include +#include #include #include #include #include +#include #include #include +#include +#include + +#include GeomAPI_Edge::GeomAPI_Edge() - : GeomAPI_Shape() { + TopoDS_Edge* anEdge = new TopoDS_Edge; + + BRep_Builder aBuilder; + aBuilder.MakeEdge(*anEdge); + + setImpl(anEdge); } GeomAPI_Edge::GeomAPI_Edge(const std::shared_ptr& theShape) @@ -70,6 +100,16 @@ bool GeomAPI_Edge::isArc() const return false; } +bool GeomAPI_Edge::isEllipse() const +{ + const TopoDS_Shape& aShape = const_cast(this)->impl(); + double aFirst, aLast; + Handle(Geom_Curve) aCurve = BRep_Tool::Curve((const TopoDS_Edge&)aShape, aFirst, aLast); + if (aCurve->IsKind(STANDARD_TYPE(Geom_Ellipse))) + return true; + return false; +} + std::shared_ptr GeomAPI_Edge::firstPoint() { const TopoDS_Shape& aShape = const_cast(this)->impl(); @@ -90,14 +130,14 @@ std::shared_ptr GeomAPI_Edge::lastPoint() return std::shared_ptr(new GeomAPI_Pnt(aPoint.X(), aPoint.Y(), aPoint.Z())); } -std::shared_ptr GeomAPI_Edge::circle() +std::shared_ptr GeomAPI_Edge::circle() const { const TopoDS_Shape& aShape = const_cast(this)->impl(); double aFirst, aLast; Handle(Geom_Curve) aCurve = BRep_Tool::Curve((const TopoDS_Edge&)aShape, aFirst, aLast); - if (aCurve) { + if (!aCurve.IsNull()) { Handle(Geom_Circle) aCirc = Handle(Geom_Circle)::DownCast(aCurve); - if (aCirc) { + if (!aCirc.IsNull()) { gp_Pnt aLoc = aCirc->Location(); std::shared_ptr aCenter(new GeomAPI_Pnt(aLoc.X(), aLoc.Y(), aLoc.Z())); gp_Dir anAxis = aCirc->Axis().Direction(); @@ -108,7 +148,24 @@ std::shared_ptr GeomAPI_Edge::circle() return std::shared_ptr(); // not circle } -std::shared_ptr GeomAPI_Edge::line() +std::shared_ptr GeomAPI_Edge::ellipse() const +{ + const TopoDS_Shape& aShape = const_cast(this)->impl(); + double aFirst, aLast; + Handle(Geom_Curve) aCurve = BRep_Tool::Curve((const TopoDS_Edge&)aShape, aFirst, aLast); + if (!aCurve.IsNull()) { + Handle(Geom_Ellipse) aElips = Handle(Geom_Ellipse)::DownCast(aCurve); + if (!aElips.IsNull()) { + gp_Elips aGpElips = aElips->Elips(); + std::shared_ptr aEllipse(new GeomAPI_Ellipse()); + aEllipse->setImpl(new gp_Elips(aGpElips)); + return aEllipse; + } + } + return std::shared_ptr(); // not elipse +} + +std::shared_ptr GeomAPI_Edge::line() const { const TopoDS_Shape& aShape = const_cast(this)->impl(); double aFirst, aLast; @@ -136,6 +193,9 @@ bool GeomAPI_Edge::isEqual(const std::shared_ptr theEdge) const const TopoDS_Shape& aMyShape = const_cast(this)->impl(); const TopoDS_Shape& aInShape = theEdge->impl(); + if (aMyShape.IsNull() || aInShape.IsNull()) + return false; + if (aMyShape.ShapeType() != aInShape.ShapeType()) return false; @@ -160,9 +220,61 @@ bool GeomAPI_Edge::isEqual(const std::shared_ptr theEdge) const gp_Pnt aInPnt1 = aInAdaptor.Value(aInStart); gp_Pnt aInPnt2 = aInAdaptor.Value(aInEnd); - if ((!aMyPnt1.IsEqual(aInPnt1, Precision::Confusion())) || + if ((!aMyPnt1.IsEqual(aInPnt1, Precision::Confusion())) || (!aMyPnt2.IsEqual(aInPnt2, Precision::Confusion()))) return false; return true; } + +void GeomAPI_Edge::getRange(double& theFirst, double& theLast) const +{ + const TopoDS_Shape& aShape = const_cast(this)->impl(); + Handle(Geom_Curve) aCurve = BRep_Tool::Curve((const TopoDS_Edge&)aShape, theFirst, theLast); +} + +bool GeomAPI_Edge::isInPlane(std::shared_ptr thePlane) const +{ + double aFirst, aLast; + const TopoDS_Shape& aShape = const_cast(this)->impl(); + Handle(Geom_Curve) aCurve = BRep_Tool::Curve((const TopoDS_Edge&)aShape, aFirst, aLast); + if (aCurve.IsNull()) + return false; + + double A, B, C, D; + thePlane->coefficients(A, B, C, D); + gp_Pln aPlane(A, B, C, D); + + bool inPlane = false; + if (aCurve->IsKind(STANDARD_TYPE(Geom_Line))) { + // check start and end points on the plane + gp_Pnt aFirstPnt = aCurve->Value(aFirst); + gp_Pnt aLastPnt = aCurve->Value(aLast); + inPlane = aPlane.SquareDistance(aFirstPnt) < Precision::SquareConfusion() && + aPlane.SquareDistance(aLastPnt) < Precision::SquareConfusion(); + } else if (aCurve->IsKind(STANDARD_TYPE(Geom_Circle))) { + // check the center on the plane and normals are collinear + Handle(Geom_Circle) aCirc = Handle(Geom_Circle)::DownCast(aCurve); + gp_Pnt aCenter = aCirc->Location(); + Standard_Real aDot = aPlane.Axis().Direction().Dot(aCirc->Axis().Direction()); + inPlane = aPlane.SquareDistance(aCenter) < Precision::SquareConfusion() && + Abs(Abs(aDot) - 1.0) < Precision::Confusion(); + } else { + // three points checking + gp_Pnt aFirstPnt = aCurve->Value(aFirst); + gp_Pnt aMidPnt = aCurve->Value((aFirst + aLast) / 2.); + gp_Pnt aLastPnt = aCurve->Value(aLast); + inPlane = aPlane.SquareDistance(aFirstPnt) < Precision::SquareConfusion() && + aPlane.SquareDistance(aMidPnt) < Precision::SquareConfusion() && + aPlane.SquareDistance(aLastPnt) < Precision::SquareConfusion(); + } + return inPlane; +} + +double GeomAPI_Edge::length() const +{ + const TopoDS_Edge& anEdge = TopoDS::Edge(impl()); + BRepAdaptor_Curve aBRepAdaptor = BRepAdaptor_Curve(anEdge); + Adaptor3d_Curve* anAdaptor3d = &aBRepAdaptor; + return GCPnts_AbscissaPoint::Length(*anAdaptor3d); +}