-// Copyright (C) 2014-20xx CEA/DEN, EDF R&D
-
-// File: GeomAlgoAPI_EdgeBuilder.cpp
-// Created: 23 Apr 2014
-// Author: Mikhail PONIKAROV
+// 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<mailto:webmaster.salome@opencascade.com>
+//
#include <GeomAlgoAPI_EdgeBuilder.h>
#include <gp_Pln.hxx>
#include <gp_Ax2.hxx>
#include <gp_Circ.hxx>
+#include <gp_Elips.hxx>
#include <Bnd_Box.hxx>
#include <BRepBndLib.hxx>
aRes->setImpl(new TopoDS_Shape(anEdge));
return aRes;
}
+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;
+}
+
+std::shared_ptr<GeomAPI_Edge> GeomAlgoAPI_EdgeBuilder::line(
+ const std::shared_ptr<GeomAPI_Lin> theLin)
+{
+ if(!theLin.get()) {
+ return std::shared_ptr<GeomAPI_Edge>();
+ }
+
+ 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;
+}
std::shared_ptr<GeomAPI_Edge> GeomAlgoAPI_EdgeBuilder::cylinderAxis(
std::shared_ptr<GeomAPI_Shape> theCylindricalFace)
Handle(Geom_Surface) aSurf = BRep_Tool::Surface(aFace, aLoc);
if (aSurf.IsNull())
return aResult;
- Handle(Geom_RectangularTrimmedSurface) aTrimmed =
+ Handle(Geom_RectangularTrimmedSurface) aTrimmed =
Handle(Geom_RectangularTrimmedSurface)::DownCast(aSurf);
if (!aTrimmed.IsNull())
aSurf = aTrimmed->BasisSurface();
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(),
+ 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());
gp_Pnt anEnd(anAxis.Location().XYZ() - anAxis.Direction().XYZ() * 100.);
anEnd.Transform(aLoc.Transformation());
*/
-
+
BRepBuilderAPI_MakeEdge anEdgeBuilder(aStart, anEnd);
std::shared_ptr<GeomAPI_Edge> aRes(new GeomAPI_Edge);
TopoDS_Edge anEdge = anEdgeBuilder.Edge();
return aRes;
}
+std::shared_ptr<GeomAPI_Edge> GeomAlgoAPI_EdgeBuilder::lineCircle(
+ std::shared_ptr<GeomAPI_Circ> theCircle)
+{
+ if(!theCircle.get()) {
+ return std::shared_ptr<GeomAPI_Edge>();
+ }
+
+ 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;
+}
+
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)
/// 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);
gp_Circ aCircle(gp_Ax2(aCenter, aDir), aRadius);
const gp_Pnt& aStart = theStartPoint->impl<gp_Pnt>();
- const gp_Pnt& anEnd = theEndPoint->impl<gp_Pnt>();
+ const gp_Pnt& anEndInter = theEndPoint->impl<gp_Pnt>();
+
+ // project end point to a circle
+ gp_XYZ aEndDir = anEndInter.XYZ() - aCenter.XYZ();
+ gp_Pnt anEnd(aCenter.XYZ() + aEndDir.Normalized() * aRadius);
BRepBuilderAPI_MakeEdge anEdgeBuilder;
- if (aStart.IsEqual(anEnd, Precision::Confusion()))
- anEdgeBuilder = BRepBuilderAPI_MakeEdge(aCircle);
- else
- anEdgeBuilder = BRepBuilderAPI_MakeEdge(aCircle, aStart, anEnd);
+ anEdgeBuilder = BRepBuilderAPI_MakeEdge(aCircle, aStart, anEnd);
anEdgeBuilder.Build();
}
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;
+}