+// Copyright (C) 2014-20xx CEA/DEN, EDF R&D
+
// File: GeomAlgoAPI_EdgeBuilder.cpp
// Created: 23 Apr 2014
// Author: Mikhail PONIKAROV
#include <gp_Pln.hxx>
#include <BRepBuilderAPI_MakeEdge.hxx>
#include <TopoDS_Edge.hxx>
+#include <TopoDS_Face.hxx>
#include <TopoDS.hxx>
#include <BRep_Tool.hxx>
#include <Geom_Plane.hxx>
+#include <Geom_CylindricalSurface.hxx>
+#include <Geom_RectangularTrimmedSurface.hxx>
+
+#include <gp_Ax2.hxx>
+#include <gp_Circ.hxx>
+#include <Bnd_Box.hxx>
+#include <BRepBndLib.hxx>
-boost::shared_ptr<GeomAPI_Shape> GeomAlgoAPI_EdgeBuilder::line(
- boost::shared_ptr<GeomAPI_Pnt> theStart, boost::shared_ptr<GeomAPI_Pnt> theEnd)
+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 boost::shared_ptr<GeomAPI_Shape>();
+ 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;
+}
+
+std::shared_ptr<GeomAPI_Edge> GeomAlgoAPI_EdgeBuilder::cylinderAxis(
+ std::shared_ptr<GeomAPI_Shape> theCylindricalFace)
+{
+ std::shared_ptr<GeomAPI_Edge> aResult;
+ const TopoDS_Shape& aShape = theCylindricalFace->impl<TopoDS_Shape>();
+ if (aShape.IsNull())
+ return aResult;
+ TopoDS_Face aFace = TopoDS::Face(aShape);
+ if (aFace.IsNull())
+ return aResult;
+ TopLoc_Location aLoc;
+ Handle(Geom_Surface) aSurf = BRep_Tool::Surface(aFace, aLoc);
+ if (aSurf.IsNull())
+ return aResult;
+ Handle(Geom_RectangularTrimmedSurface) aTrimmed =
+ Handle(Geom_RectangularTrimmedSurface)::DownCast(aSurf);
+ if (!aTrimmed.IsNull())
+ aSurf = aTrimmed->BasisSurface();
+ Handle(Geom_CylindricalSurface) aCyl = Handle(Geom_CylindricalSurface)::DownCast(aSurf);
+ if (aCyl.IsNull())
+ return aResult;
+ gp_Ax1 anAxis = aCyl->Axis();
+ // compute the start and the end points of the resulting edge by the bounding box of the face
+ // (vertices projected to the axis) plus 10%
+ Bnd_Box aFaceBnd;
+ BRepBndLib::Add(aFace, aFaceBnd);
+ gp_Pnt aBoxMin(aFaceBnd.CornerMin()), aBoxMax(aFaceBnd.CornerMax());
+ bool isFirst = true;
+ double aParamMin = 0, aParamMax = 0;
+ for(int aX = 0; aX < 2; aX++) {
+ for(int aY = 0; aY < 2; aY++) {
+ for(int aZ = 0; aZ < 2; aZ++) {
+ gp_XYZ aBoxVertex(aX == 0 ? aBoxMin.X() : aBoxMax.X(),
+ aY == 0 ? aBoxMin.Y() : aBoxMax.Y(), aZ == 0 ? aBoxMin.Z() : aBoxMax.Z());
+ gp_XYZ aVec(aBoxVertex - anAxis.Location().XYZ());
+ double aProjParam = aVec.Dot(anAxis.Direction().XYZ());
+ if (isFirst) {
+ isFirst = false;
+ aParamMin = aProjParam;
+ aParamMax = aProjParam;
+ } else {
+ if (aParamMin > aProjParam)
+ aParamMin = aProjParam;
+ else if (aParamMax < aProjParam)
+ aParamMax = aProjParam;
+ }
+ }
+ }
+ }
+ // add 10%
+ double aDelta = aParamMax - aParamMin;
+ if (aDelta < 1.e-4) aDelta = 1.e-4;
+ aParamMin -= aDelta * 0.1;
+ aParamMax += aDelta * 0.1;
+
+ gp_Pnt aStart(aParamMin * anAxis.Direction().XYZ() + anAxis.Location().XYZ());
+ gp_Pnt anEnd(aParamMax * anAxis.Direction().XYZ() + anAxis.Location().XYZ());
+ /*
+ gp_Pnt aStart(anAxis.Location().Transformed(aLoc.Transformation()));
+ // edge length is 100, "-" because cylinder of extrusion has negative direction with the cylinder
+ gp_Pnt anEnd(anAxis.Location().XYZ() - anAxis.Direction().XYZ() * 100.);
+ anEnd.Transform(aLoc.Transformation());
+ */
+
BRepBuilderAPI_MakeEdge anEdgeBuilder(aStart, anEnd);
- boost::shared_ptr<GeomAPI_Shape> aRes(new GeomAPI_Shape);
+ std::shared_ptr<GeomAPI_Edge> aRes(new GeomAPI_Edge);
+ TopoDS_Edge anEdge = anEdgeBuilder.Edge();
+ // an axis is an infinite object
+ anEdge.Infinite(Standard_True);
+ aRes->setImpl(new TopoDS_Shape(anEdge));
+ return aRes;
+}
+
+std::shared_ptr<GeomAPI_Edge> GeomAlgoAPI_EdgeBuilder::lineCircle(
+ std::shared_ptr<GeomAPI_Pnt> theCenter, std::shared_ptr<GeomAPI_Dir> theNormal,
+ double theRadius)
+{
+ const gp_Pnt& aCenter = theCenter->impl<gp_Pnt>();
+ const gp_Dir& aDir = theNormal->impl<gp_Dir>();
+
+ gp_Circ aCircle(gp_Ax2(aCenter, aDir), theRadius);
+
+ BRepBuilderAPI_MakeEdge anEdgeBuilder(aCircle);
+ 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::lineCircleArc(
+ std::shared_ptr<GeomAPI_Pnt> theCenter, std::shared_ptr<GeomAPI_Pnt> theStartPoint,
+ std::shared_ptr<GeomAPI_Pnt> theEndPoint, std::shared_ptr<GeomAPI_Dir> theNormal)
+{
+ const gp_Pnt& aCenter = theCenter->impl<gp_Pnt>();
+ const gp_Dir& aDir = theNormal->impl<gp_Dir>();
+
+ double aRadius = theCenter->distance(theStartPoint);
+ gp_Circ aCircle(gp_Ax2(aCenter, aDir), aRadius);
+
+ const gp_Pnt& aStart = theStartPoint->impl<gp_Pnt>();
+ const gp_Pnt& anEnd = theEndPoint->impl<gp_Pnt>();
+
+ BRepBuilderAPI_MakeEdge anEdgeBuilder;
+ if (aStart.IsEqual(anEnd, Precision::Confusion()))
+ anEdgeBuilder = BRepBuilderAPI_MakeEdge(aCircle);
+ else
+ anEdgeBuilder = BRepBuilderAPI_MakeEdge(aCircle, aStart, anEnd);
+
+ std::shared_ptr<GeomAPI_Edge> aRes(new GeomAPI_Edge);
+ anEdgeBuilder.Build();
+
+ if (anEdgeBuilder.IsDone())
+ aRes->setImpl(new TopoDS_Shape(anEdgeBuilder.Edge()));
+ else
+ aRes = std::shared_ptr<GeomAPI_Edge>();
+ return aRes;
+}