X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;f=src%2FGeomAlgoAPI%2FGeomAlgoAPI_EdgeBuilder.cpp;h=3c177a1787e65cce0425633c6ba3540bf2cf2c2e;hb=50c26f8bb70e02bfe1f70e0f3545611378433e77;hp=0e622ba36958ba6581158c8babd8b1d65c87ee25;hpb=3ca758d1ecab7b2e601b36425ea4c20e9e857412;p=modules%2Fshaper.git diff --git a/src/GeomAlgoAPI/GeomAlgoAPI_EdgeBuilder.cpp b/src/GeomAlgoAPI/GeomAlgoAPI_EdgeBuilder.cpp index 0e622ba36..3c177a178 100644 --- a/src/GeomAlgoAPI/GeomAlgoAPI_EdgeBuilder.cpp +++ b/src/GeomAlgoAPI/GeomAlgoAPI_EdgeBuilder.cpp @@ -1,3 +1,5 @@ +// Copyright (C) 2014-20xx CEA/DEN, EDF R&D + // File: GeomAlgoAPI_EdgeBuilder.cpp // Created: 23 Apr 2014 // Author: Mikhail PONIKAROV @@ -6,33 +8,111 @@ #include #include #include +#include #include #include #include +#include +#include #include #include +#include +#include -boost::shared_ptr GeomAlgoAPI_EdgeBuilder::line( - boost::shared_ptr theStart, boost::shared_ptr theEnd) +std::shared_ptr GeomAlgoAPI_EdgeBuilder::line( + std::shared_ptr theStart, std::shared_ptr theEnd) { const gp_Pnt& aStart = theStart->impl(); const gp_Pnt& anEnd = theEnd->impl(); if (aStart.IsEqual(anEnd, Precision::Confusion())) - return boost::shared_ptr(); + return std::shared_ptr(); if (Abs(aStart.SquareDistance(anEnd)) > 1.e+100) - return boost::shared_ptr(); + return std::shared_ptr(); + BRepBuilderAPI_MakeEdge anEdgeBuilder(aStart, anEnd); + std::shared_ptr aRes(new GeomAPI_Edge); + TopoDS_Edge anEdge = anEdgeBuilder.Edge(); + aRes->setImpl(new TopoDS_Shape(anEdge)); + return aRes; +} + +std::shared_ptr GeomAlgoAPI_EdgeBuilder::cylinderAxis( + std::shared_ptr theCylindricalFace) +{ + std::shared_ptr aResult; + const TopoDS_Shape& aShape = theCylindricalFace->impl(); + 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 aRes(new GeomAPI_Edge); + std::shared_ptr 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; } -boost::shared_ptr GeomAlgoAPI_EdgeBuilder::lineCircle( - boost::shared_ptr theCenter, - boost::shared_ptr theNormal, double theRadius) +std::shared_ptr GeomAlgoAPI_EdgeBuilder::lineCircle( + std::shared_ptr theCenter, std::shared_ptr theNormal, + double theRadius) { const gp_Pnt& aCenter = theCenter->impl(); const gp_Dir& aDir = theNormal->impl(); @@ -40,17 +120,15 @@ boost::shared_ptr GeomAlgoAPI_EdgeBuilder::lineCircle( gp_Circ aCircle(gp_Ax2(aCenter, aDir), theRadius); BRepBuilderAPI_MakeEdge anEdgeBuilder(aCircle); - boost::shared_ptr aRes(new GeomAPI_Edge); + std::shared_ptr aRes(new GeomAPI_Edge); TopoDS_Edge anEdge = anEdgeBuilder.Edge(); aRes->setImpl(new TopoDS_Shape(anEdge)); return aRes; } -boost::shared_ptr GeomAlgoAPI_EdgeBuilder::lineCircleArc( - boost::shared_ptr theCenter, - boost::shared_ptr theStartPoint, - boost::shared_ptr theEndPoint, - boost::shared_ptr theNormal) +std::shared_ptr GeomAlgoAPI_EdgeBuilder::lineCircleArc( + std::shared_ptr theCenter, std::shared_ptr theStartPoint, + std::shared_ptr theEndPoint, std::shared_ptr theNormal) { const gp_Pnt& aCenter = theCenter->impl(); const gp_Dir& aDir = theNormal->impl(); @@ -62,18 +140,18 @@ boost::shared_ptr GeomAlgoAPI_EdgeBuilder::lineCircleArc( const gp_Pnt& anEnd = theEndPoint->impl(); BRepBuilderAPI_MakeEdge anEdgeBuilder; - if (aStart.IsEqual(anEnd, Precision::Confusion()) || - gp_Pnt(0, 0, 0).IsEqual(anEnd, Precision::Confusion())) + if (aStart.IsEqual(anEnd, Precision::Confusion()) + || gp_Pnt(0, 0, 0).IsEqual(anEnd, Precision::Confusion())) anEdgeBuilder = BRepBuilderAPI_MakeEdge(aCircle); else anEdgeBuilder = BRepBuilderAPI_MakeEdge(aCircle, aStart, anEnd); - boost::shared_ptr aRes(new GeomAPI_Edge); + std::shared_ptr aRes(new GeomAPI_Edge); anEdgeBuilder.Build(); if (anEdgeBuilder.IsDone()) aRes->setImpl(new TopoDS_Shape(anEdgeBuilder.Edge())); else - aRes = boost::shared_ptr(); + aRes = std::shared_ptr(); return aRes; }