-// 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
//
// 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<mailto:webmaster.salome@opencascade.com>
+// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
//
#include<GeomAPI_Edge.h>
#include <Geom_Curve.hxx>
#include <Geom_Line.hxx>
#include <Geom_Circle.hxx>
+#include <Geom_TrimmedCurve.hxx>
#include <Geom_Ellipse.hxx>
#include <Geom_Plane.hxx>
#include <GeomAPI_IntCS.hxx>
return false;
}
+/// extracts a circle curve from the arbitrary curve, returns null is it is different type
+static Handle(Geom_Circle) circ(const Handle(Geom_Curve) theCurve)
+{
+ Handle(Geom_Circle) aResult = Handle(Geom_Circle)::DownCast(theCurve);
+ if (!aResult.IsNull())
+ return aResult;
+ // check this may be a trimmed curve that contains circle inside
+ Handle(Geom_TrimmedCurve) aTrimmed = Handle(Geom_TrimmedCurve)::DownCast(theCurve);
+ while(!aTrimmed.IsNull()) {
+ aResult = Handle(Geom_Circle)::DownCast(aTrimmed->BasisCurve());
+ if (!aResult.IsNull())
+ return aResult;
+ aTrimmed = Handle(Geom_TrimmedCurve)::DownCast(aTrimmed->BasisCurve());
+ }
+ return aResult; // null, not circle
+}
+
bool GeomAPI_Edge::isCircle() const
{
const TopoDS_Shape& aShape = const_cast<GeomAPI_Edge*>(this)->impl<TopoDS_Shape>();
double aFirst, aLast;
Handle(Geom_Curve) aCurve = BRep_Tool::Curve((const TopoDS_Edge&)aShape, aFirst, aLast);
- if (aCurve.IsNull()) // degenerative edge
- return false;
- if (aCurve->IsKind(STANDARD_TYPE(Geom_Circle)))
- {
+ if (!circ(aCurve).IsNull()) {
// Check the difference of first and last parameters to be equal to the curve period
if (Abs(aLast - aFirst - aCurve->Period()) < Precision::PConfusion())
return true;
const TopoDS_Shape& aShape = const_cast<GeomAPI_Edge*>(this)->impl<TopoDS_Shape>();
double aFirst, aLast;
Handle(Geom_Curve) aCurve = BRep_Tool::Curve((const TopoDS_Edge&)aShape, aFirst, aLast);
- if (aCurve.IsNull()) // degenerative edge
- return false;
- if (aCurve->IsKind(STANDARD_TYPE(Geom_Circle)))
- {
+ if (!circ(aCurve).IsNull()) {
// Check the difference of first and last parameters is not equal the curve period
if (Abs(aLast - aFirst - aCurve->Period()) >= Precision::PConfusion())
return true;
const TopoDS_Shape& aShape = const_cast<GeomAPI_Edge*>(this)->impl<TopoDS_Shape>();
double aFirst, aLast;
Handle(Geom_Curve) aCurve = BRep_Tool::Curve((const TopoDS_Edge&)aShape, aFirst, aLast);
- if (!aCurve.IsNull()) {
- Handle(Geom_Circle) aCirc = Handle(Geom_Circle)::DownCast(aCurve);
- if (!aCirc.IsNull()) {
- gp_Pnt aLoc = aCirc->Location();
- std::shared_ptr<GeomAPI_Pnt> aCenter(new GeomAPI_Pnt(aLoc.X(), aLoc.Y(), aLoc.Z()));
- gp_Dir anAxis = aCirc->Axis().Direction();
- std::shared_ptr<GeomAPI_Dir> aDir(new GeomAPI_Dir(anAxis.X(), anAxis.Y(), anAxis.Z()));
- return std::shared_ptr<GeomAPI_Circ>(new GeomAPI_Circ(aCenter, aDir, aCirc->Radius()));
- }
+ Handle(Geom_Circle) aCirc = circ(aCurve);
+ if (!aCirc.IsNull()) {
+ gp_Pnt aLoc = aCirc->Location();
+ std::shared_ptr<GeomAPI_Pnt> aCenter(new GeomAPI_Pnt(aLoc.X(), aLoc.Y(), aLoc.Z()));
+ gp_Dir anAxis = aCirc->Axis().Direction();
+ std::shared_ptr<GeomAPI_Dir> aDir(new GeomAPI_Dir(anAxis.X(), anAxis.Y(), anAxis.Z()));
+ return std::shared_ptr<GeomAPI_Circ>(new GeomAPI_Circ(aCenter, aDir, aCirc->Radius()));
}
return std::shared_ptr<GeomAPI_Circ>(); // not circle
}
return aEllipse;
}
}
- return std::shared_ptr<GeomAPI_Ellipse>(); // not elipse
+ return std::shared_ptr<GeomAPI_Ellipse>(); // not ellipse
}
std::shared_ptr<GeomAPI_Lin> GeomAPI_Edge::line() const
double aInStart, aInEnd;
Handle(Geom_Curve) aInCurve = BRep_Tool::Curve(TopoDS::Edge(aInShape), aInStart, aInEnd);
+ // Check that end point parameters are the same
+ if ((aMyStart != aInStart) || (aMyEnd != aInEnd))
+ return false;
+
// Check that curves a the same type
GeomAdaptor_Curve aMyAdaptor(aMyCurve);
GeomAdaptor_Curve aInAdaptor(aInCurve);
if (aMyAdaptor.GetType() != aInAdaptor.GetType())
return false;
- // Check that end point parameters are the same
- if ((aMyStart != aInStart) || (aMyEnd != aInEnd))
- return false;
-
// Check that end points are equal
gp_Pnt aMyPnt1 = aMyAdaptor.Value(aMyStart);
gp_Pnt aMyPnt2 = aMyAdaptor.Value(aMyEnd);
return true;
}
+// LCOV_EXCL_START
void GeomAPI_Edge::getRange(double& theFirst, double& theLast) const
{
const TopoDS_Shape& aShape = const_cast<GeomAPI_Edge*>(this)->impl<TopoDS_Shape>();
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();
+ Handle(Geom_Circle) aCirc = circ(aCurve);
+ if (!aCirc.IsNull()) {
+ 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;
}
+// LCOV_EXCL_STOP
void GeomAPI_Edge::intersectWithPlane(const std::shared_ptr<GeomAPI_Pln> thePlane,
std::list<std::shared_ptr<GeomAPI_Pnt>>& theResult) const
