-// Copyright (C) 2014-2017 CEA/DEN, EDF R&D
+// Copyright (C) 2014-2024 CEA, EDF
//
// 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 <GeomAlgoAPI_EdgeBuilder.h>
+
+#include <GeomAPI_Ax2.h>
+#include <GeomAPI_Ax3.h>
+#include <GeomAPI_BSpline2d.h>
+#include <GeomAPI_Ellipse.h>
+#include <GeomAPI_Pnt2d.h>
+
#include <gp_Pln.hxx>
#include <BRepBuilderAPI_MakeEdge.hxx>
#include <TopoDS_Edge.hxx>
#include <TopoDS_Face.hxx>
#include <TopoDS.hxx>
#include <BRep_Tool.hxx>
+#include <Geom2d_BSplineCurve.hxx>
#include <Geom_Plane.hxx>
#include <Geom_CylindricalSurface.hxx>
#include <Geom_RectangularTrimmedSurface.hxx>
+#include <GeomLib.hxx>
#include <gp_Ax2.hxx>
#include <gp_Circ.hxx>
+#include <gp_Elips.hxx>
#include <Bnd_Box.hxx>
#include <BRepBndLib.hxx>
+static GeomEdgePtr createLine(const gp_Pnt& theStart, const gp_Pnt& theEnd)
+{
+ GeomEdgePtr aRes;
+ if (!theStart.IsEqual(theEnd, Precision::Confusion()) &&
+ Abs(theStart.SquareDistance(theEnd)) < 1.e+100) {
+ BRepBuilderAPI_MakeEdge anEdgeBuilder(theStart, theEnd);
+ TopoDS_Edge anEdge = anEdgeBuilder.Edge();
+ aRes = GeomEdgePtr(new GeomAPI_Edge);
+ aRes->setImpl(new TopoDS_Shape(anEdge));
+ }
+ return aRes;
+}
+
std::shared_ptr<GeomAPI_Edge> GeomAlgoAPI_EdgeBuilder::line(
std::shared_ptr<GeomAPI_Pnt> theStart, std::shared_ptr<GeomAPI_Pnt> theEnd)
{
const gp_Pnt& aStart = theStart->impl<gp_Pnt>();
const gp_Pnt& anEnd = theEnd->impl<gp_Pnt>();
-
- if (aStart.IsEqual(anEnd, Precision::Confusion()))
- return std::shared_ptr<GeomAPI_Edge>();
- if (Abs(aStart.SquareDistance(anEnd)) > 1.e+100)
- return std::shared_ptr<GeomAPI_Edge>();
- BRepBuilderAPI_MakeEdge anEdgeBuilder(aStart, anEnd);
- std::shared_ptr<GeomAPI_Edge> aRes(new GeomAPI_Edge);
- TopoDS_Edge anEdge = anEdgeBuilder.Edge();
- aRes->setImpl(new TopoDS_Shape(anEdge));
- return aRes;
+ return createLine(aStart, anEnd);
}
std::shared_ptr<GeomAPI_Edge> GeomAlgoAPI_EdgeBuilder::line(
double theDX, double theDY, double theDZ)
{
-
const gp_Pnt& aStart = gp_Pnt(0, 0, 0);
const gp_Pnt& anEnd = gp_Pnt(theDX, theDY, theDZ);
-
- if (aStart.IsEqual(anEnd, Precision::Confusion()))
- return std::shared_ptr<GeomAPI_Edge>();
- if (Abs(aStart.SquareDistance(anEnd)) > 1.e+100)
- return std::shared_ptr<GeomAPI_Edge>();
- BRepBuilderAPI_MakeEdge anEdgeBuilder(aStart, anEnd);
- std::shared_ptr<GeomAPI_Edge> aRes(new GeomAPI_Edge);
- TopoDS_Edge anEdge = anEdgeBuilder.Edge();
- aRes->setImpl(new TopoDS_Shape(anEdge));
- return aRes;
+ return createLine(aStart, anEnd);
}
std::shared_ptr<GeomAPI_Edge> GeomAlgoAPI_EdgeBuilder::line(
const std::shared_ptr<GeomAPI_Lin> theLin)
{
- if(!theLin.get()) {
- return std::shared_ptr<GeomAPI_Edge>();
+ GeomEdgePtr aRes;
+ if (theLin.get()) {
+ const gp_Lin& aLin = theLin->impl<gp_Lin>();
+ BRepBuilderAPI_MakeEdge anEdgeBuilder(aLin);
+ TopoDS_Edge anEdge = anEdgeBuilder.Edge();
+ aRes = GeomEdgePtr(new GeomAPI_Edge);
+ aRes->setImpl(new TopoDS_Shape(anEdge));
}
-
- const gp_Lin& aLin = theLin->impl<gp_Lin>();
- BRepBuilderAPI_MakeEdge anEdgeBuilder(aLin);
- std::shared_ptr<GeomAPI_Edge> aRes(new GeomAPI_Edge());
- TopoDS_Edge anEdge = anEdgeBuilder.Edge();
- aRes->setImpl(new TopoDS_Shape(anEdge));
return aRes;
}
aParamMax += aDelta * 0.1;
gp_Pnt aStart(aParamMin * anAxis.Direction().XYZ() + anAxis.Location().XYZ());
+ aStart.Transform(aLoc.Transformation());
gp_Pnt anEnd(aParamMax * anAxis.Direction().XYZ() + anAxis.Location().XYZ());
+ anEnd.Transform(aLoc.Transformation());
/*
gp_Pnt aStart(anAxis.Location().Transformed(aLoc.Transformation()));
// edge length is 100, "-" because cylinder of extrusion has negative direction with the cylinder
std::shared_ptr<GeomAPI_Edge> GeomAlgoAPI_EdgeBuilder::lineCircle(
std::shared_ptr<GeomAPI_Circ> theCircle)
{
- if(!theCircle.get()) {
- return std::shared_ptr<GeomAPI_Edge>();
+ GeomEdgePtr aRes;
+ if (theCircle.get()) {
+ const gp_Circ& aCirc = theCircle->impl<gp_Circ>();
+ BRepBuilderAPI_MakeEdge anEdgeBuilder(aCirc);
+ TopoDS_Edge anEdge = anEdgeBuilder.Edge();
+ aRes = GeomEdgePtr(new GeomAPI_Edge);
+ aRes->setImpl(new TopoDS_Shape(anEdge));
}
-
- const gp_Circ& aCirc = theCircle->impl<gp_Circ>();
- BRepBuilderAPI_MakeEdge anEdgeBuilder(aCirc);
- std::shared_ptr<GeomAPI_Edge> aRes(new GeomAPI_Edge());
- TopoDS_Edge anEdge = anEdgeBuilder.Edge();
- aRes->setImpl(new TopoDS_Shape(anEdge));
return aRes;
}
/// OCCT creates an edge on a circle with empty radius, but visualization
/// is not able to process it
- if (theCenter->isEqual(theStartPoint))
+ if (theCenter->isEqual(theStartPoint) || theCenter->isEqual(theEndPoint))
return aRes;
double aRadius = theCenter->distance(theStartPoint);
}
return aRes;
}
+
+std::shared_ptr<GeomAPI_Edge> GeomAlgoAPI_EdgeBuilder::ellipse(
+ const std::shared_ptr<GeomAPI_Pnt>& theCenter,
+ const std::shared_ptr<GeomAPI_Dir>& theNormal,
+ const std::shared_ptr<GeomAPI_Dir>& theMajorAxis,
+ const double theMajorRadius,
+ const double theMinorRadius)
+{
+ const gp_Pnt& aCenter = theCenter->impl<gp_Pnt>();
+ const gp_Dir& aNormal = theNormal->impl<gp_Dir>();
+ const gp_Dir& aMajorAxis = theMajorAxis->impl<gp_Dir>();
+
+ gp_Elips anEllipse(gp_Ax2(aCenter, aNormal, aMajorAxis), theMajorRadius, theMinorRadius);
+
+ BRepBuilderAPI_MakeEdge anEdgeBuilder(anEllipse);
+ std::shared_ptr<GeomAPI_Edge> aRes(new GeomAPI_Edge);
+ TopoDS_Edge anEdge = anEdgeBuilder.Edge();
+ aRes->setImpl(new TopoDS_Shape(anEdge));
+ return aRes;
+}
+
+std::shared_ptr<GeomAPI_Edge> GeomAlgoAPI_EdgeBuilder::ellipticArc(
+ const std::shared_ptr<GeomAPI_Pnt>& theCenter,
+ const std::shared_ptr<GeomAPI_Dir>& theNormal,
+ const std::shared_ptr<GeomAPI_Dir>& theMajorAxis,
+ const double theMajorRadius,
+ const double theMinorRadius,
+ const std::shared_ptr<GeomAPI_Pnt>& theStart,
+ const std::shared_ptr<GeomAPI_Pnt>& theEnd)
+{
+ std::shared_ptr<GeomAPI_Ax2> anAx2(new GeomAPI_Ax2(theCenter, theNormal, theMajorAxis));
+ GeomAPI_Ellipse anEllipse(anAx2, theMajorRadius, theMinorRadius);
+
+ GeomPointPtr aStartPnt = anEllipse.project(theStart);
+ GeomPointPtr aEndPnt = anEllipse.project(theEnd);
+
+ double aStartParam, aEndParam;
+ anEllipse.parameter(aStartPnt, Precision::Confusion(), aStartParam);
+ anEllipse.parameter(aEndPnt, Precision::Confusion(), aEndParam);
+
+ BRepBuilderAPI_MakeEdge anEdgeBuilder(anEllipse.impl<gp_Elips>(), aStartParam, aEndParam);
+ GeomEdgePtr aRes(new GeomAPI_Edge);
+ TopoDS_Edge anEdge = anEdgeBuilder.Edge();
+ aRes->setImpl(new TopoDS_Shape(anEdge));
+ return aRes;
+}
+
+GeomEdgePtr GeomAlgoAPI_EdgeBuilder::bsplineOnPlane(
+ const std::shared_ptr<GeomAPI_Ax3>& thePlane,
+ const std::list<GeomPnt2dPtr>& thePoles,
+ const std::list<double>& theWeights,
+ const std::list<double>& theKnots,
+ const std::list<int>& theMults,
+ const int theDegree,
+ const bool thePeriodic)
+{
+ std::shared_ptr<GeomAPI_BSpline2d> aBSplineCurve(
+ new GeomAPI_BSpline2d(theDegree, thePoles, theWeights, theKnots, theMults, thePeriodic));
+ return bsplineOnPlane(thePlane, aBSplineCurve);
+}
+
+GeomEdgePtr GeomAlgoAPI_EdgeBuilder::bsplineOnPlane(
+ const std::shared_ptr<GeomAPI_Ax3>& thePlane,
+ const std::shared_ptr<GeomAPI_BSpline2d>& theCurve)
+{
+ Handle(Geom_Curve) aCurve3D = GeomLib::To3d(thePlane->impl<gp_Ax3>().Ax2(),
+ theCurve->impl<Handle_Geom2d_BSplineCurve>());
+
+ BRepBuilderAPI_MakeEdge anEdgeBuilder(aCurve3D);
+ GeomEdgePtr aRes(new GeomAPI_Edge);
+ TopoDS_Edge anEdge = anEdgeBuilder.Edge();
+ aRes->setImpl(new TopoDS_Shape(anEdge));
+ return aRes;
+}