-// File: HLRBRep.cxx
-// Created: Thu Aug 27 12:33:14 1992
-// Author: Christophe MARION
-// <cma@sdsun2>
+// Created on: 1992-08-27
+// Created by: Christophe MARION
+// Copyright (c) 1992-1999 Matra Datavision
+// Copyright (c) 1999-2014 OPEN CASCADE SAS
+//
+// This file is part of Open CASCADE Technology software library.
+//
+// This library is free software; you can redistribute it and/or modify it under
+// the terms of the GNU Lesser General Public License version 2.1 as published
+// by the Free Software Foundation, with special exception defined in the file
+// OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
+// distribution for complete text of the license and disclaimer of any warranty.
+//
+// Alternatively, this file may be used under the terms of Open CASCADE
+// commercial license or contractual agreement.
#include <HLRBRep.ixx>
#include <BRepLib_MakeEdge2d.hxx>
#include <Geom2d_BezierCurve.hxx>
#include <Geom2d_BSplineCurve.hxx>
+#include <Geom_BSplineCurve.hxx>
+#include <Geom_TrimmedCurve.hxx>
#include <TColStd_Array1OfInteger.hxx>
#include <TColStd_Array1OfReal.hxx>
#include <TColgp_Array1OfPnt2d.hxx>
+#include <BRep_Tool.hxx>
+#include <TopoDS.hxx>
+#include <TopExp.hxx>
+#include <BRepLib_MakeVertex.hxx>
+#include <BRep_Builder.hxx>
+
//=======================================================================
//function : MakeEdge
//=======================================================================
TopoDS_Edge HLRBRep::MakeEdge (const HLRBRep_Curve& ec,
- const Standard_Real U1,
- const Standard_Real U2)
+ const Standard_Real U1,
+ const Standard_Real U2)
{
TopoDS_Edge Edg;
- //gp_Pnt2d P,P1,P2;
- Standard_Real sta3d = U1;
- Standard_Real end3d = U2;
- Standard_Real sta = ec.Parameter2d(U1);
- Standard_Real end = ec.Parameter2d(U2);
-
- if (ec.GetType() == GeomAbs_Line) {
+ const Standard_Real sta = ec.Parameter2d(U1);
+ const Standard_Real end = ec.Parameter2d(U2);
+
+ switch (ec.GetType())
+ {
+ case GeomAbs_Line:
Edg = BRepLib_MakeEdge2d(ec.Line(),sta,end);
- }
- else if (ec.GetType() == GeomAbs_Circle) {
+ break;
+
+ case GeomAbs_Circle:
Edg = BRepLib_MakeEdge2d(ec.Circle(),sta,end);
- }
- else if (ec.GetType() == GeomAbs_Ellipse) {
+ break;
+
+ case GeomAbs_Ellipse:
Edg = BRepLib_MakeEdge2d(ec.Ellipse(),sta,end);
- }
- else if (ec.GetType() == GeomAbs_Hyperbola) {
+ break;
+
+ case GeomAbs_Hyperbola:
Edg = BRepLib_MakeEdge2d(ec.Hyperbola(),sta,end);
- }
- else if (ec.GetType() == GeomAbs_Parabola) {
+ break;
+
+ case GeomAbs_Parabola:
Edg = BRepLib_MakeEdge2d(ec.Parabola(),sta,end);
- }
- else if (ec.GetType() == GeomAbs_BezierCurve) {
+ break;
+
+ case GeomAbs_BezierCurve: {
TColgp_Array1OfPnt2d Poles(1,ec.NbPoles());
+ Handle(Geom2d_BezierCurve) ec2d;
if (ec.IsRational()) {
TColStd_Array1OfReal Weights(1,ec.NbPoles());
ec.PolesAndWeights(Poles,Weights);
- Edg = BRepLib_MakeEdge2d(new Geom2d_BezierCurve(Poles,Weights),sta,end);
+ ec2d = new Geom2d_BezierCurve(Poles,Weights);
}
else {
ec.Poles(Poles);
- Edg = BRepLib_MakeEdge2d(new Geom2d_BezierCurve(Poles),sta,end);
+ ec2d = new Geom2d_BezierCurve(Poles);
}
+ BRepLib_MakeEdge2d mke2d(ec2d,sta,end);
+ if (mke2d.IsDone())
+ Edg = mke2d.Edge();
+ break;
}
- else if (ec.GetType() == GeomAbs_BSplineCurve) {
- TColgp_Array1OfPnt2d Poles(1,ec.NbPoles());
- TColStd_Array1OfReal knots(1,ec.NbKnots());
- TColStd_Array1OfInteger mults(1,ec.NbKnots());
- //-- ec.KnotsAndMultiplicities(knots,mults);
- ec.Knots(knots);
- ec.Multiplicities(mults);
- if (ec.IsRational()) {
- TColStd_Array1OfReal Weights(1,ec.NbPoles());
- ec.PolesAndWeights(Poles,Weights);
- Edg = BRepLib_MakeEdge2d
- (new Geom2d_BSplineCurve
- (Poles,Weights,knots,mults,ec.Degree()),sta,end);
+
+ case GeomAbs_BSplineCurve: {
+ Handle(Geom2d_BSplineCurve) ec2d;
+ GeomAdaptor_Curve GAcurve = ec.GetCurve().Curve();
+ TopoDS_Edge anEdge = ec.GetCurve().Edge();
+ Standard_Real fpar, lpar;
+ Handle(Geom_Curve) aCurve = BRep_Tool::Curve(anEdge, fpar, lpar);
+ if (aCurve->DynamicType() == STANDARD_TYPE(Geom_TrimmedCurve))
+ aCurve = (Handle(Geom_TrimmedCurve)::DownCast(aCurve))->BasisCurve();
+ const Handle(Geom_BSplineCurve)& BSplCurve = Handle(Geom_BSplineCurve)::DownCast(aCurve);
+ Handle(Geom_BSplineCurve) theCurve = Handle(Geom_BSplineCurve)::DownCast(BSplCurve->Copy());
+ if (theCurve->IsPeriodic() && !GAcurve.IsClosed())
+ {
+ theCurve->Segment(sta, end);
+ TColgp_Array1OfPnt2d Poles(1, theCurve->NbPoles());
+ TColStd_Array1OfReal knots(1, theCurve->NbKnots());
+ TColStd_Array1OfInteger mults(1, theCurve->NbKnots());
+ //-- ec.KnotsAndMultiplicities(knots,mults);
+ theCurve->Knots(knots);
+ theCurve->Multiplicities(mults);
+ if (theCurve->IsRational()) {
+ TColStd_Array1OfReal Weights(1, theCurve->NbPoles());
+ ec.PolesAndWeights(theCurve, Poles, Weights);
+ ec2d = new Geom2d_BSplineCurve(Poles, Weights, knots, mults,
+ theCurve->Degree(), theCurve->IsPeriodic());
+ }
+ else {
+ ec.Poles(theCurve, Poles);
+ ec2d = new Geom2d_BSplineCurve(Poles, knots, mults,
+ theCurve->Degree(), theCurve->IsPeriodic());
+ }
}
- else {
- ec.Poles(Poles);
- Edg = BRepLib_MakeEdge2d
- (new Geom2d_BSplineCurve(Poles,knots,mults,ec.Degree()),sta,end);
+ else
+ {
+ TColgp_Array1OfPnt2d Poles(1,ec.NbPoles());
+ TColStd_Array1OfReal knots(1,ec.NbKnots());
+ TColStd_Array1OfInteger mults(1,ec.NbKnots());
+ //-- ec.KnotsAndMultiplicities(knots,mults);
+ ec.Knots(knots);
+ ec.Multiplicities(mults);
+ if (ec.IsRational()) {
+ TColStd_Array1OfReal Weights(1,ec.NbPoles());
+ ec.PolesAndWeights(Poles,Weights);
+ ec2d = new Geom2d_BSplineCurve(Poles,Weights,knots,mults,ec.Degree(),ec.IsPeriodic());
+ }
+ else {
+ ec.Poles(Poles);
+ ec2d = new Geom2d_BSplineCurve(Poles,knots,mults,ec.Degree(),ec.IsPeriodic());
+ }
}
+ BRepLib_MakeEdge2d mke2d(ec2d, sta, end);
+ if (mke2d.IsDone())
+ Edg = mke2d.Edge();
+ break;
}
- else {
- Standard_Integer nbPnt = 15;
+ default: {
+ const Standard_Integer nbPnt = 15;
TColgp_Array1OfPnt2d Poles(1,nbPnt);
TColStd_Array1OfReal knots(1,nbPnt);
TColStd_Array1OfInteger mults(1,nbPnt);
mults.Init(1);
mults(1 ) = 2;
mults(nbPnt) = 2;
- Standard_Real step = (end3d-sta3d)/(nbPnt-1);
-
- for (Standard_Integer i = 1; i <= nbPnt; i++) {
- Poles(i) = ec.Value(sta3d);
- knots(i) = sta3d;
- sta3d += step;
+ const Standard_Real step = (U2-U1)/(nbPnt-1);
+ Standard_Real par3d = U1;
+ for (Standard_Integer i = 1; i < nbPnt; i++) {
+ Poles(i) = ec.Value(par3d);
+ knots(i) = par3d;
+ par3d += step;
}
- Edg = BRepLib_MakeEdge2d
- (new Geom2d_BSplineCurve(Poles,knots,mults,1),sta,end);
+ Poles(nbPnt) = ec.Value(U2);
+ knots(nbPnt) = U2;
+
+ Handle(Geom2d_BSplineCurve) ec2d = new Geom2d_BSplineCurve(Poles,knots,mults,1);
+ BRepLib_MakeEdge2d mke2d(ec2d,sta,end);
+ if (mke2d.IsDone())
+ Edg = mke2d.Edge();
+ }
}
return Edg;
}
+//=======================================================================
+//function : MakeEdge3d
+//purpose :
+//=======================================================================
+
+TopoDS_Edge HLRBRep::MakeEdge3d(const HLRBRep_Curve& ec,
+ const Standard_Real U1,
+ const Standard_Real U2)
+{
+ TopoDS_Edge Edg;
+ //const Standard_Real sta = ec.Parameter2d(U1);
+ //const Standard_Real end = ec.Parameter2d(U2);
+
+ TopoDS_Edge anEdge = ec.GetCurve().Edge();
+ Standard_Real fpar, lpar;
+ //BRep_Tool::Range(anEdge, fpar, lpar);
+ //Handle(Geom_Curve) aCurve = BRep_Tool::Curve(anEdge, fpar, lpar);
+ BRepAdaptor_Curve BAcurve(anEdge);
+ fpar = BAcurve.FirstParameter();
+ lpar = BAcurve.LastParameter();
+
+ Edg = TopoDS::Edge(anEdge.EmptyCopied());
+ Edg.Orientation(TopAbs_FORWARD);
+ BRep_Builder BB;
+ BB.Range(Edg, U1, U2);
+
+ //Share vertices if possible
+ TopoDS_Vertex V1, V2, V1new, V2new;
+ TopExp::Vertices(anEdge, V1, V2);
+
+ Standard_Real Tol = Precision::PConfusion();
+ if (Abs(fpar - U1) <= Tol)
+ V1new = V1;
+ else
+ {
+ gp_Pnt aPnt = BAcurve.Value(U1);
+ V1new = BRepLib_MakeVertex(aPnt);
+ }
+ if (Abs(lpar - U2) <= Tol)
+ V2new = V2;
+ else
+ {
+ gp_Pnt aPnt = BAcurve.Value(U2);
+ V2new = BRepLib_MakeVertex(aPnt);
+ }
+
+ V1new.Orientation(TopAbs_FORWARD);
+ V2new.Orientation(TopAbs_REVERSED);
+ BB.Add(Edg, V1new);
+ BB.Add(Edg, V2new);
+ return Edg;
+}
+
//=======================================================================
//function : PolyHLRAngleAndDeflection
//purpose :
Standard_Real& OutAngl,
Standard_Real& OutDefl)
{
- static Standard_Real HAngMin = 1*PI/180;
- static Standard_Real HAngLim = 5*PI/180;
- static Standard_Real HAngMax = 35*PI/180;
+ static Standard_Real HAngMin = 1*M_PI/180;
+ static Standard_Real HAngLim = 5*M_PI/180;
+ static Standard_Real HAngMax = 35*M_PI/180;
OutAngl = InAngl;
if (OutAngl < HAngMin) OutAngl = HAngMin;
(HAngMax - HAngLim) / (HAngMax - HAngMin));
OutDefl = OutAngl * OutAngl * 0.5;
}
-