X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;f=src%2FGeomAPI%2FGeomAPI_Edge.cpp;h=9364605b9b1b95c6cc6862e8145285e6024ff76b;hb=a211165b993cec59888b300199b44b525c0daf69;hp=7966720d34217bff89d855912dc5560cc3e332d4;hpb=a36dc3bdc0503158c453605807b2c5205f03d52c;p=modules%2Fshaper.git diff --git a/src/GeomAPI/GeomAPI_Edge.cpp b/src/GeomAPI/GeomAPI_Edge.cpp index 7966720d3..9364605b9 100644 --- a/src/GeomAPI/GeomAPI_Edge.cpp +++ b/src/GeomAPI/GeomAPI_Edge.cpp @@ -1,4 +1,4 @@ -// Copyright (C) 2014-2017 CEA/DEN, EDF R&D +// Copyright (C) 2014-2019 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 @@ -12,10 +12,9 @@ // // 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 +// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA // -// See http://www.salome-platform.org/ or -// email : webmaster.salome@opencascade.com +// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com // #include @@ -26,6 +25,7 @@ #include #include #include +#include #include @@ -36,12 +36,15 @@ #include #include #include +#include #include #include #include #include #include #include +#include +#include #include #include #include @@ -68,6 +71,48 @@ GeomAPI_Edge::GeomAPI_Edge(const std::shared_ptr& theShape) } } +void GeomAPI_Edge::vertices(std::shared_ptr& theStartVertex, + std::shared_ptr& theEndVertex) const +{ + const TopoDS_Edge& anEdge = impl(); + TopoDS_Vertex aStart, aEnd; + TopExp::Vertices(anEdge, aStart, aEnd); + theStartVertex.reset(new GeomAPI_Vertex); + theStartVertex->setImpl(new TopoDS_Vertex(aStart)); + theEndVertex.reset(new GeomAPI_Vertex); + theEndVertex->setImpl(new TopoDS_Vertex(aEnd)); +} + +static Handle(Geom_Curve) baseCurve(const TopoDS_Edge& theEdge) +{ + double aFirst, aLast; + Handle(Geom_Curve) aCurve = BRep_Tool::Curve(theEdge, aFirst, aLast); + while (aCurve->IsKind(STANDARD_TYPE(Geom_TrimmedCurve))) { + Handle(Geom_TrimmedCurve) tc = Handle(Geom_TrimmedCurve)::DownCast(aCurve); + aCurve = tc->BasisCurve(); + } + return aCurve; +} + +bool GeomAPI_Edge::isSameGeometry(const std::shared_ptr theShape) const +{ + if (!theShape->isEdge()) + return false; + if (isSame(theShape)) + return true; + + TopoDS_Edge anOwnEdge = TopoDS::Edge(impl()); + TopoDS_Edge anOtherEdge = TopoDS::Edge(theShape->impl()); + + Handle(Geom_Curve) anOwnCurve = baseCurve(anOwnEdge); + Handle(Geom_Curve) anOtherCurve = baseCurve(anOtherEdge); + GeomAPI_ExtremaCurveCurve anExtrema(anOwnCurve, anOtherCurve); + + bool isSame = anExtrema.Extrema().IsParallel() && + anExtrema.TotalLowerDistance() < Precision::Confusion(); + return isSame; +} + bool GeomAPI_Edge::isLine() const { const TopoDS_Shape& aShape = const_cast(this)->impl(); @@ -135,6 +180,18 @@ bool GeomAPI_Edge::isEllipse() const return false; } +bool GeomAPI_Edge::isBSpline() 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()) // degenerative edge + return false; + while (aCurve->IsKind(STANDARD_TYPE(Geom_TrimmedCurve))) + aCurve = Handle(Geom_TrimmedCurve)::DownCast(aCurve)->BasisCurve(); + return aCurve->IsKind(STANDARD_TYPE(Geom_BSplineCurve)); +} + std::shared_ptr GeomAPI_Edge::firstPoint() { const TopoDS_Shape& aShape = const_cast(this)->impl(); @@ -250,12 +307,26 @@ bool GeomAPI_Edge::isEqual(const std::shared_ptr theEdge) const return true; } +void GeomAPI_Edge::setRange(const double& theFirst, const double& theLast) +{ + TopoDS_Edge anEdge = impl(); + double aFirst, aLast; + Handle(Geom_Curve) aCurve = BRep_Tool::Curve(anEdge, aFirst, aLast); + double aTolerance = BRep_Tool::Tolerance(anEdge); + if (aCurve->DynamicType() == STANDARD_TYPE(Geom_TrimmedCurve)) { + aCurve = Handle(Geom_TrimmedCurve)::DownCast(aCurve)->BasisCurve(); + BRep_Builder().UpdateEdge(anEdge, aCurve, aTolerance); + } + BRep_Builder().Range(anEdge, theFirst, theLast); +} + 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); } +// LCOV_EXCL_START bool GeomAPI_Edge::isInPlane(std::shared_ptr thePlane) const { double aFirst, aLast; @@ -294,6 +365,7 @@ bool GeomAPI_Edge::isInPlane(std::shared_ptr thePlane) const } return inPlane; } +// LCOV_EXCL_STOP void GeomAPI_Edge::intersectWithPlane(const std::shared_ptr thePlane, std::list>& theResult) const @@ -305,8 +377,9 @@ void GeomAPI_Edge::intersectWithPlane(const std::shared_ptr thePlan double A, B, C, D; thePlane->coefficients(A, B, C, D); gp_Pln aPln(A, B, C, D); - Handle(Geom_Plane) aPlane = new Geom_Plane(aPln); + + // intersect the plane with the curve GeomAPI_IntCS aIntersect; aIntersect.Perform(aCurve, aPlane); if (aIntersect.IsDone() && (aIntersect.NbPoints() > 0)) { @@ -325,6 +398,20 @@ void GeomAPI_Edge::intersectWithPlane(const std::shared_ptr thePlan theResult.push_back(aPntPtr); } } + else { + // find minimal distance between the plane and the curve + GeomAPI_ExtremaCurveSurface anExtrema(aCurve, aPlane); + double aTolerance = BRep_Tool::Tolerance(TopoDS::Edge(aShape)); + if (anExtrema.NbExtrema() > 0 && + !anExtrema.Extrema().IsParallel() && + anExtrema.LowerDistance() < aTolerance) { + // distance is lower than tolerance => tangent case + gp_Pnt aPntC, aPntS; + anExtrema.NearestPoints(aPntC, aPntS); + std::shared_ptr aPntPtr(new GeomAPI_Pnt(aPntS.X(), aPntS.Y(), aPntS.Z())); + theResult.push_back(aPntPtr); + } + } } }