From dca6e5062075498badfe0f9074aed26e47c43343 Mon Sep 17 00:00:00 2001 From: gka Date: Mon, 24 Apr 2006 07:32:40 +0000 Subject: [PATCH] Addition method to get tangent on curve by parameter --- src/GEOMImpl/GEOMImpl_VectorDriver.cxx | 31 ++++++++++++++++++++++++-- 1 file changed, 29 insertions(+), 2 deletions(-) diff --git a/src/GEOMImpl/GEOMImpl_VectorDriver.cxx b/src/GEOMImpl/GEOMImpl_VectorDriver.cxx index 531077599..cb0b81c2b 100644 --- a/src/GEOMImpl/GEOMImpl_VectorDriver.cxx +++ b/src/GEOMImpl/GEOMImpl_VectorDriver.cxx @@ -34,6 +34,8 @@ #include #include +#include +#include #include #include @@ -67,7 +69,7 @@ Standard_Integer GEOMImpl_VectorDriver::Execute(TFunction_Logbook& log) const GEOMImpl_IVector aPI (aFunction); Standard_Integer aType = aFunction->GetType(); - if (aType != VECTOR_DX_DY_DZ && aType != VECTOR_TWO_PNT) return 0; + if (aType != VECTOR_DX_DY_DZ && aType != VECTOR_TWO_PNT && aType != VECTOR_TANGENT_CURVE_PAR) return 0; TopoDS_Shape aShape; @@ -98,7 +100,32 @@ Standard_Integer GEOMImpl_VectorDriver::Execute(TFunction_Logbook& log) const Standard_ConstructionError::Raise("The end points are too close"); } aShape = BRepBuilderAPI_MakeEdge(V1, V2).Shape(); - } else { + } + else if(aType == VECTOR_TANGENT_CURVE_PAR) { + Handle(GEOM_Function) aRefCurve = aPI.GetCurve(); + TopoDS_Shape aRefShape = aRefCurve->GetValue(); + if (aRefShape.ShapeType() != TopAbs_EDGE) { + Standard_TypeMismatch::Raise + ("Tangent On Curve creation aborted : curve shape is not an edge"); + } + Standard_Real aFParam =0., aLParam =0., aParam =0.; + Handle(Geom_Curve) aCurve = BRep_Tool::Curve(TopoDS::Edge(aRefShape), aFParam, aLParam); + if(aCurve.IsNull()) { + Standard_TypeMismatch::Raise + ("Tangent On Curve creation aborted : curve is null"); + } + + aParam = aFParam + (aLParam - aFParam) * aPI.GetParameter(); + gp_Pnt aPoint1,aPoint2; + gp_Vec aVec; + aCurve->D1(aParam,aPoint1,aVec); + if(aVec.Magnitude() < gp::Resolution()) + Standard_TypeMismatch::Raise + ("Tangent On Curve creation aborted : invalid value of tangent"); + aPoint2.SetXYZ(aPoint1.XYZ() + aVec.XYZ()); + BRepBuilderAPI_MakeEdge aBuilder(aPoint1,aPoint2); + if(aBuilder.IsDone()) + aShape = aBuilder.Shape(); } if (aShape.IsNull()) return 0; -- 2.39.2