OpenCascade/occt
1
0
Fork 0
mirror of https://git.dev.opencascade.org/repos/occt.git synced 2025-04-05 18:16:23 +03:00
Code Packages Releases Activity

0023076: Empty result done by intersection algorithm for a curve and a surface

Browse Source
This commit is contained in:
pkv 2012-04-06 11:16:24 +04:00
parent 664cae74a9
commit 0316739bfe
1 changed files with 349 additions and 354 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

703
src/IntCurveSurface/IntCurveSurface_Inter.gxx
View File

@ -18,10 +18,10 @@
// purpose or non-infringement. Please see the License for the specific terms // purpose or non-infringement. Please see the License for the specific terms
// and conditions governing the rights and limitations under the License. // and conditions governing the rights and limitations under the License.
#ifndef DEB //#ifndef DEB
#define No_Standard_RangeError //#define No_Standard_RangeError
#define No_Standard_OutOfRange //#define No_Standard_OutOfRange
#endif //#endif
#define TOLTANGENCY 0.00000001 #define TOLTANGENCY 0.00000001
@ -60,11 +60,7 @@
#include <Extrema_ExtElC.hxx> #include <Extrema_ExtElC.hxx>
#include <Extrema_POnCurv.hxx> #include <Extrema_POnCurv.hxx>
#if 0
//-- jgv patch (from)
#include <Extrema_ExtPS.hxx>
//-- jgv patch (to)
#endif
#include <ProjLib_Plane.hxx> #include <ProjLib_Plane.hxx>
#include <IntAna2d_AnaIntersection.hxx> #include <IntAna2d_AnaIntersection.hxx>
@ -77,146 +73,120 @@
#include <Adaptor3d_HCurve.hxx> #include <Adaptor3d_HCurve.hxx>
#include <Adaptor3d_HSurface.hxx> #include <Adaptor3d_HSurface.hxx>
#if 0
//-- jgv patch (from)
#include <BndLib_Add3dCurve.hxx>
#include <GeomAdaptor_Curve.hxx>
#include <Geom_Line.hxx>
#include <Geom_Circle.hxx>
#include <Geom_Ellipse.hxx>
#include <Geom_Hyperbola.hxx>
#include <Geom_Parabola.hxx>
#include <Geom_TrimmedCurve.hxx>
#include <GeomAdaptor_Surface.hxx>
#include <Geom_Plane.hxx>
#include <Geom_CylindricalSurface.hxx>
#include <Geom_ConicalSurface.hxx>
#include <Geom_ToroidalSurface.hxx>
#include <Geom_SphericalSurface.hxx>
#include <Geom_SurfaceOfRevolution.hxx>
#include <Geom_SurfaceOfLinearExtrusion.hxx>
#include <Geom_OffsetSurface.hxx>
#include <Geom_RectangularTrimmedSurface.hxx>
//-- jgv patch (to)
#endif
#include <TColgp_Array2OfPnt.hxx> #include <TColgp_Array2OfPnt.hxx>
#include <TColStd_HArray1OfReal.hxx> #include <TColStd_HArray1OfReal.hxx>
#include <Adaptor3d_TopolTool.hxx> #include <Adaptor3d_TopolTool.hxx>
#include <ElCLib.hxx> #include <ElCLib.hxx>
//================================================================================ static
static void EstLimForInfExtr(const gp_Lin& Line, void EstLimForInfExtr(const gp_Lin& Line,
const TheSurface& surface, const TheSurface& surface,
const Standard_Boolean IsOffSurf, const Standard_Boolean IsOffSurf,
const Standard_Integer nbsu, const Standard_Integer nbsu,
const Standard_Boolean U1inf, const Standard_Boolean U1inf,
const Standard_Boolean U2inf, const Standard_Boolean U2inf,
const Standard_Boolean V1inf, const Standard_Boolean V1inf,
const Standard_Boolean V2inf, const Standard_Boolean V2inf,
Standard_Real& U1new, Standard_Real& U1new,
Standard_Real& U2new, Standard_Real& U2new,
Standard_Real& V1new, Standard_Real& V1new,
Standard_Real& V2new, Standard_Real& V2new,
Standard_Boolean& NoIntersection); Standard_Boolean& NoIntersection);
//================================================================================ static
static void EstLimForInfRevl(const gp_Lin& Line, void EstLimForInfRevl(const gp_Lin& Line,
const TheSurface& surface, const TheSurface& surface,
const Standard_Boolean U1inf, const Standard_Boolean U1inf,
const Standard_Boolean U2inf, const Standard_Boolean U2inf,
const Standard_Boolean V1inf, const Standard_Boolean V1inf,
const Standard_Boolean V2inf, const Standard_Boolean V2inf,
Standard_Real& U1new, Standard_Real& U1new,
Standard_Real& U2new, Standard_Real& U2new,
Standard_Real& V1new, Standard_Real& V1new,
Standard_Real& V2new, Standard_Real& V2new,
Standard_Boolean& NoIntersection); Standard_Boolean& NoIntersection);
//================================================================================ static
static void EstLimForInfOffs(const gp_Lin& Line, void EstLimForInfOffs(const gp_Lin& Line,
const TheSurface& surface, const TheSurface& surface,
const Standard_Integer nbsu, const Standard_Integer nbsu,
const Standard_Boolean U1inf, const Standard_Boolean U1inf,
const Standard_Boolean U2inf, const Standard_Boolean U2inf,
const Standard_Boolean V1inf, const Standard_Boolean V1inf,
const Standard_Boolean V2inf, const Standard_Boolean V2inf,
Standard_Real& U1new, Standard_Real& U1new,
Standard_Real& U2new, Standard_Real& U2new,
Standard_Real& V1new, Standard_Real& V1new,
Standard_Real& V2new, Standard_Real& V2new,
Standard_Boolean& NoIntersection); Standard_Boolean& NoIntersection);
//================================================================================ static
static void EstLimForInfSurf(Standard_Real& U1new, void EstLimForInfSurf(Standard_Real& U1new,
Standard_Real& U2new, Standard_Real& U2new,
Standard_Real& V1new, Standard_Real& V1new,
Standard_Real& V2new); Standard_Real& V2new);
//================================================================================ static
static void SectionPointToParameters(const Intf_SectionPoint& Sp, void SectionPointToParameters(const Intf_SectionPoint& Sp,
const IntCurveSurface_ThePolyhedron& Surf, const IntCurveSurface_ThePolyhedron& Surf,
const IntCurveSurface_ThePolygon& Curv, const IntCurveSurface_ThePolygon& Curv,
Standard_Real& u, Standard_Real& u,
Standard_Real& v, Standard_Real& v,
Standard_Real& w); Standard_Real& w);
//================================================================================
static void IntCurveSurface_ComputeTransitions(const TheCurve& curve, static
const Standard_Real w, void IntCurveSurface_ComputeTransitions(const TheCurve& curve,
IntCurveSurface_TransitionOnCurve& TransOnCurve, const Standard_Real w,
const TheSurface& surface, IntCurveSurface_TransitionOnCurve& TransOnCurve,
const Standard_Real u, const TheSurface& surface,
const Standard_Real v); const Standard_Real u,
//================================================================================ const Standard_Real v);
static void IntCurveSurface_ComputeParamsOnQuadric(const TheSurface& surface,
const gp_Pnt& P, static
Standard_Real& u, void IntCurveSurface_ComputeParamsOnQuadric(const TheSurface& surface,
Standard_Real& v); const gp_Pnt& P,
//================================================================================ Standard_Real& u,
#if 0 Standard_Real& v);
static Handle(Geom_Curve) GetCurve(const Handle(Adaptor3d_HCurve) AdCurve);
//================================================================================ static
static Handle(Geom_Surface) GetSurface(const Handle(Adaptor3d_HSurface) AdSurface); Standard_Boolean DoTrim(const TheCurve& curve,
#endif const TheSurface& surface);
//================================================================================ static
//== void DoCurve(const TheCurve& curve,
IntCurveSurface_Inter::IntCurveSurface_Inter() { Bnd_Box& boxCurve);
static
void DoCommon(TColStd_Array1OfReal& Coords,
Standard_Real& Cmin,
Standard_Real& Cmax);
static
void DoCommonBox(const Bnd_Box& boxSurface,
const Bnd_Box& boxCurve,
TColStd_Array1OfReal& X,
TColStd_Array1OfReal& Y,
TColStd_Array1OfReal& Z);
static
void ProjectIntersectAndEstLim(const gp_Lin& theLine,
const gp_Pln& thePln,
const ProjLib_Plane& theBasCurvProj,
Standard_Real& theVmin,
Standard_Real& theVmax,
Standard_Boolean& theNoIntersection);
//=======================================================================
//function : IntCurveSurface_Inter
//purpose :
//=======================================================================
IntCurveSurface_Inter::IntCurveSurface_Inter()
{
} }
//=======================================================================
static Standard_Boolean DoTrim(const TheCurve& curve, //function : DoSurface
const TheSurface& surface) //purpose :
{ //=======================================================================
Standard_Boolean isAnaCurve = Standard_False, isAnaSurface = Standard_False;
GeomAbs_CurveType CurveType = TheCurveTool::GetType(curve);
switch (CurveType) {
case GeomAbs_Line:
case GeomAbs_Circle:
case GeomAbs_Ellipse:
case GeomAbs_Hyperbola:
case GeomAbs_Parabola: isAnaCurve = Standard_True; break;
default: break;
}
GeomAbs_SurfaceType SurfaceType = TheSurfaceTool::GetType(surface);
switch (SurfaceType) {
case GeomAbs_Plane:
case GeomAbs_Cylinder:
case GeomAbs_Cone:
case GeomAbs_Sphere:
case GeomAbs_Torus: isAnaSurface = Standard_True; break;
default: break;
}
Standard_Boolean result = (isAnaCurve && isAnaSurface) ? Standard_False : Standard_True;
if(result) {
Standard_Boolean isUClosed = (TheSurfaceTool::IsUClosed(surface) || TheSurfaceTool::IsUPeriodic(surface));
Standard_Boolean isVClosed = (TheSurfaceTool::IsVClosed(surface) || TheSurfaceTool::IsVPeriodic(surface));
if(isUClosed && isVClosed)
result = Standard_False;
}
return result;
}
// modified by NIZHNY-MKK Tue Jul 26 14:41:59 2005
// static void DoSurface(const TheSurface& surface,
void IntCurveSurface_Inter::DoSurface(const TheSurface& surface, void IntCurveSurface_Inter::DoSurface(const TheSurface& surface,
const Standard_Real u0, const Standard_Real u0,
const Standard_Real u1, const Standard_Real u1,
@ -261,8 +231,12 @@ void IntCurveSurface_Inter::DoSurface(const TheSurface& surface,
gap = Max(Ures,Vres); gap = Max(Ures,Vres);
} }
static void DoCurve(const TheCurve& curve, //=======================================================================
Bnd_Box& boxCurve) //function : DoCurve
//purpose :
//=======================================================================
void DoCurve(const TheCurve& curve,
Bnd_Box& boxCurve)
{ {
Standard_Real CF = TheCurveTool::FirstParameter(curve); Standard_Real CF = TheCurveTool::FirstParameter(curve);
Standard_Real CL = TheCurveTool::LastParameter(curve); Standard_Real CL = TheCurveTool::LastParameter(curve);
@ -282,9 +256,14 @@ static void DoCurve(const TheCurve& curve,
boxCurve.Add(aPnt); boxCurve.Add(aPnt);
} }
} }
static void DoCommon(TColStd_Array1OfReal& Coords,
Standard_Real& Cmin, //=======================================================================
Standard_Real& Cmax) //function : DoCommon
//purpose :
//=======================================================================
void DoCommon(TColStd_Array1OfReal& Coords,
Standard_Real& Cmin,
Standard_Real& Cmax)
{ {
Standard_Integer i = 0, j = 0; Standard_Integer i = 0, j = 0;
for(i = 1; i <= 4; i++) { for(i = 1; i <= 4; i++) {
@ -300,6 +279,10 @@ static void DoCommon(TColStd_Array1OfReal& Coords,
Cmin = Coords(3); Cmin = Coords(3);
} }
//=======================================================================
//function : DoCommonBox
//purpose :
//=======================================================================
static void DoCommonBox(const Bnd_Box& boxSurface, static void DoCommonBox(const Bnd_Box& boxSurface,
const Bnd_Box& boxCurve, const Bnd_Box& boxCurve,
TColStd_Array1OfReal& X, TColStd_Array1OfReal& X,
@ -333,9 +316,11 @@ static void DoCommonBox(const Bnd_Box& boxSurface,
Z.SetValue(1,CZmin); Z.SetValue(2,CZmax); Z.SetValue(1,CZmin); Z.SetValue(2,CZmax);
} }
// modified by NIZHNY-MKK Tue Jul 26 14:41:42 2005 //=======================================================================
//function : DoNewBounds
//purpose :
//=======================================================================
void IntCurveSurface_Inter::DoNewBounds( void IntCurveSurface_Inter::DoNewBounds(
// static void DoNewBounds(
const TheSurface& surface, const TheSurface& surface,
const Standard_Real u0, const Standard_Real u0,
const Standard_Real u1, const Standard_Real u1,
@ -424,22 +409,17 @@ void IntCurveSurface_Inter::DoNewBounds(
} }
} }
//================================================================================ //=======================================================================
//== P e r f o r m g e n e r a l //function : Perform
//== //purpose : Decompose la surface si besoin est
//== Decompose la surface si besoin est //=======================================================================
//================================================================================
void IntCurveSurface_Inter::Perform(const TheCurve& curve, void IntCurveSurface_Inter::Perform(const TheCurve& curve,
const TheSurface& surface) { const TheSurface& surface) {
ResetFields(); ResetFields();
done = Standard_True; done = Standard_True;
Standard_Integer NbUOnS = TheSurfaceTool::NbUIntervals(surface,GeomAbs_C2); Standard_Integer NbUOnS = TheSurfaceTool::NbUIntervals(surface,GeomAbs_C2);
Standard_Integer NbVOnS = TheSurfaceTool::NbVIntervals(surface,GeomAbs_C2); Standard_Integer NbVOnS = TheSurfaceTool::NbVIntervals(surface,GeomAbs_C2);
#ifdef DEB
Standard_Integer NbOnC = TheCurveTool::NbIntervals(curve,GeomAbs_C2); Standard_Integer NbOnC = TheCurveTool::NbIntervals(curve,GeomAbs_C2);
#else
TheCurveTool::NbIntervals(curve,GeomAbs_C2);
#endif
Standard_Real U0,U1,V0,V1; Standard_Real U0,U1,V0,V1;
if(NbUOnS > 1) { if(NbUOnS > 1) {
@ -653,8 +633,10 @@ void IntCurveSurface_Inter::Perform(const TheCurve& curve,
Perform(curve,surface,U0,V0,U1,V1); Perform(curve,surface,U0,V0,U1,V1);
} }
} }
//================================================================================ //=======================================================================
//================================================================================ //function : Perform
//purpose :
//=======================================================================
void IntCurveSurface_Inter::Perform(const TheCurve& curve, void IntCurveSurface_Inter::Perform(const TheCurve& curve,
const TheSurface& surface, const TheSurface& surface,
const Standard_Real U1,const Standard_Real V1, const Standard_Real U1,const Standard_Real V1,
@ -736,7 +718,10 @@ void IntCurveSurface_Inter::Perform(const TheCurve& curve,
} }
} }
} }
//================================================================================ //=======================================================================
//function : Perform
//purpose :
//=======================================================================
void IntCurveSurface_Inter::Perform(const TheCurve& curve, void IntCurveSurface_Inter::Perform(const TheCurve& curve,
const IntCurveSurface_ThePolygon& polygon, const IntCurveSurface_ThePolygon& polygon,
const TheSurface& surface) { const TheSurface& surface) {
@ -755,7 +740,10 @@ void IntCurveSurface_Inter::Perform(const TheCurve& curve,
IntCurveSurface_ThePolyhedron polyhedron(surface,nbsu,nbsv,u1,v1,u2,v2); IntCurveSurface_ThePolyhedron polyhedron(surface,nbsu,nbsv,u1,v1,u2,v2);
Perform(curve,polygon,surface,polyhedron); Perform(curve,polygon,surface,polyhedron);
} }
//================================================================================ //=======================================================================
//function : Perform
//purpose :
//=======================================================================
void IntCurveSurface_Inter::Perform(const TheCurve& curve, void IntCurveSurface_Inter::Perform(const TheCurve& curve,
const TheSurface& surface, const TheSurface& surface,
const IntCurveSurface_ThePolyhedron& polyhedron) { const IntCurveSurface_ThePolyhedron& polyhedron) {
@ -766,7 +754,10 @@ void IntCurveSurface_Inter::Perform(const TheCurve& curve,
IntCurveSurface_ThePolygon polygon(curve,TheCurveTool::NbSamples(curve,u1,u2)); IntCurveSurface_ThePolygon polygon(curve,TheCurveTool::NbSamples(curve,u1,u2));
Perform(curve,polygon,surface,polyhedron); Perform(curve,polygon,surface,polyhedron);
} }
//================================================================================ //=======================================================================
//function : Perform
//purpose :
//=======================================================================
void IntCurveSurface_Inter::Perform(const TheCurve& curve, void IntCurveSurface_Inter::Perform(const TheCurve& curve,
const IntCurveSurface_ThePolygon& polygon, const IntCurveSurface_ThePolygon& polygon,
const TheSurface& surface, const TheSurface& surface,
@ -781,6 +772,10 @@ void IntCurveSurface_Inter::Perform(const TheCurve& curve,
InternalPerform(curve,polygon,surface,polyhedron,u1,v1,u2,v2); InternalPerform(curve,polygon,surface,polyhedron,u1,v1,u2,v2);
} }
//=======================================================================
//function : Perform
//purpose :
//=======================================================================
void IntCurveSurface_Inter::Perform(const TheCurve& curve, void IntCurveSurface_Inter::Perform(const TheCurve& curve,
const IntCurveSurface_ThePolygon& polygon, const IntCurveSurface_ThePolygon& polygon,
const TheSurface& surface, const TheSurface& surface,
@ -795,10 +790,11 @@ void IntCurveSurface_Inter::Perform(const TheCurve& curve,
v2 = TheSurfaceTool::LastVParameter(surface); v2 = TheSurfaceTool::LastVParameter(surface);
InternalPerform(curve,polygon,surface,polyhedron,u1,v1,u2,v2,BndBSB); InternalPerform(curve,polygon,surface,polyhedron,u1,v1,u2,v2,BndBSB);
} }
//================================================================================ //=======================================================================
//== C a l c u l d u p o i n t a p p r o c h e == //function : InternalPerform
//== p u i s d u p o i n t E x a c t == //purpose : C a l c u l d u p o i n t a p p r o c h e
//================================================================================ //== p u i s d u p o i n t E x a c t
//=======================================================================
void IntCurveSurface_Inter::InternalPerform(const TheCurve& curve, void IntCurveSurface_Inter::InternalPerform(const TheCurve& curve,
const IntCurveSurface_ThePolygon& polygon, const IntCurveSurface_ThePolygon& polygon,
const TheSurface& surface, const TheSurface& surface,
@ -940,6 +936,10 @@ void IntCurveSurface_Inter::InternalPerform(const TheCurve& curve,
} }
} }
//=======================================================================
//function : InternalPerform
//purpose :
//=======================================================================
void IntCurveSurface_Inter::InternalPerform(const TheCurve& curve, void IntCurveSurface_Inter::InternalPerform(const TheCurve& curve,
const IntCurveSurface_ThePolygon& polygon, const IntCurveSurface_ThePolygon& polygon,
const TheSurface& surface, const TheSurface& surface,
@ -1079,7 +1079,10 @@ void IntCurveSurface_Inter::InternalPerform(const TheCurve& curve,
delete [] TabU; delete [] TabU;
} }
} }
//================================================================================ //=======================================================================
//function : InternalPerformCurveQuadric
//purpose :
//=======================================================================
void IntCurveSurface_Inter::InternalPerformCurveQuadric(const TheCurve& curve, void IntCurveSurface_Inter::InternalPerformCurveQuadric(const TheCurve& curve,
const TheSurface& surface) { const TheSurface& surface) {
IntCurveSurface_TheQuadCurvExactInter QuadCurv(surface,curve); IntCurveSurface_TheQuadCurvExactInter QuadCurv(surface,curve);
@ -1095,7 +1098,10 @@ void IntCurveSurface_Inter::InternalPerformCurveQuadric(const TheCurve&
} }
} }
//================================================================================ //=======================================================================
//function : InternalPerform
//purpose :
//=======================================================================
void IntCurveSurface_Inter::InternalPerform(const TheCurve& curve, void IntCurveSurface_Inter::InternalPerform(const TheCurve& curve,
const IntCurveSurface_ThePolygon& polygon, const IntCurveSurface_ThePolygon& polygon,
const TheSurface& surface, const TheSurface& surface,
@ -1148,7 +1154,10 @@ void IntCurveSurface_Inter::InternalPerform(const TheCurve& curve,
} }
} //-- Fin : la Surface est une quadrique } //-- Fin : la Surface est une quadrique
} }
//================================================================================ //=======================================================================
//function : PerformConicSurf
//purpose :
//=======================================================================
void IntCurveSurface_Inter::PerformConicSurf(const gp_Lin& Line, void IntCurveSurface_Inter::PerformConicSurf(const gp_Lin& Line,
const TheCurve& curve, const TheCurve& curve,
const TheSurface& surface, const TheSurface& surface,
@ -1275,7 +1284,10 @@ void IntCurveSurface_Inter::PerformConicSurf(const gp_Lin& Line,
} }
} }
} }
//================================================================================ //=======================================================================
//function : PerformConicSurf
//purpose :
//=======================================================================
void IntCurveSurface_Inter::PerformConicSurf(const gp_Circ& Circle, void IntCurveSurface_Inter::PerformConicSurf(const gp_Circ& Circle,
const TheCurve& curve, const TheCurve& curve,
const TheSurface& surface, const TheSurface& surface,
@ -1317,7 +1329,10 @@ void IntCurveSurface_Inter::PerformConicSurf(const gp_Circ& Circle,
} }
} }
} }
//================================================================================ //=======================================================================
//function : PerformConicSurf
//purpose :
//=======================================================================
void IntCurveSurface_Inter::PerformConicSurf(const gp_Elips& Ellipse, void IntCurveSurface_Inter::PerformConicSurf(const gp_Elips& Ellipse,
const TheCurve& curve, const TheCurve& curve,
const TheSurface& surface, const TheSurface& surface,
@ -1359,7 +1374,10 @@ void IntCurveSurface_Inter::PerformConicSurf(const gp_Elips& Ellipse,
} }
} }
} }
//================================================================================ //=======================================================================
//function : PerformConicSurf
//purpose :
//=======================================================================
void IntCurveSurface_Inter::PerformConicSurf(const gp_Parab& Parab, void IntCurveSurface_Inter::PerformConicSurf(const gp_Parab& Parab,
const TheCurve& curve, const TheCurve& curve,
const TheSurface& surface, const TheSurface& surface,
@ -1415,7 +1433,10 @@ void IntCurveSurface_Inter::PerformConicSurf(const gp_Parab& Parab,
} }
} }
} }
//================================================================================ //=======================================================================
//function : PerformConicSurf
//purpose :
//=======================================================================
void IntCurveSurface_Inter::PerformConicSurf(const gp_Hypr& Hypr, void IntCurveSurface_Inter::PerformConicSurf(const gp_Hypr& Hypr,
const TheCurve& curve, const TheCurve& curve,
const TheSurface& surface, const TheSurface& surface,
@ -1471,7 +1492,10 @@ void IntCurveSurface_Inter::PerformConicSurf(const gp_Hypr& Hypr,
} }
} }
} }
//================================================================================ //=======================================================================
//function : AppendIntAna
//purpose :
//=======================================================================
void IntCurveSurface_Inter::AppendIntAna(const TheCurve& curve, void IntCurveSurface_Inter::AppendIntAna(const TheCurve& curve,
const TheSurface& surface, const TheSurface& surface,
const IntAna_IntConicQuad& intana_ConicQuad) { const IntAna_IntConicQuad& intana_ConicQuad) {
@ -1497,7 +1521,10 @@ void IntCurveSurface_Inter::AppendIntAna(const TheCurve& curve,
//-- cout<<" IntAna Conic Quad Not Done "<<endl; //-- cout<<" IntAna Conic Quad Not Done "<<endl;
} }
} }
//================================================================================ //=======================================================================
//function : AppendPoint
//purpose :
//=======================================================================
void IntCurveSurface_Inter::AppendPoint(const TheCurve& curve, void IntCurveSurface_Inter::AppendPoint(const TheCurve& curve,
const Standard_Real lw, const Standard_Real lw,
const TheSurface& surface, const TheSurface& surface,
@ -1547,17 +1574,23 @@ void IntCurveSurface_Inter::AppendPoint(const TheCurve& curve,
} }
//================================================================================ //=======================================================================
//function : AppendSegment
//purpose :
//=======================================================================
void IntCurveSurface_Inter::AppendSegment(const TheCurve& , void IntCurveSurface_Inter::AppendSegment(const TheCurve& ,
const Standard_Real , const Standard_Real ,
const Standard_Real , const Standard_Real ,
const TheSurface& ) { const TheSurface& ) {
cout<<" !!! Not Yet Implemented IntCurveSurface_Inter::Append(const IntCurveSurf ...)"<<endl; //cout<<" !!! Not Yet Implemented
//IntCurveSurface_Inter::Append(const IntCurveSurf ...)"<<endl;
} }
//================================================================================ //=======================================================================
//== P o i n t d i n t e r f e r e n c e - - > U , V e t W == //function : SectionPointToParameters
//================================================================================ //purpose : P o i n t d i n t e r f e r e n c e - - >
// U , V e t W
//=======================================================================
void SectionPointToParameters(const Intf_SectionPoint& Sp, void SectionPointToParameters(const Intf_SectionPoint& Sp,
const IntCurveSurface_ThePolyhedron& Polyhedron, const IntCurveSurface_ThePolyhedron& Polyhedron,
const IntCurveSurface_ThePolygon& Polygon, const IntCurveSurface_ThePolygon& Polygon,
@ -1633,7 +1666,10 @@ void SectionPointToParameters(const Intf_SectionPoint& Sp,
U = u1; U = u1;
V = v1; V = v1;
} }
//================================================================================ //=======================================================================
//function : IntCurveSurface_ComputeTransitions
//purpose :
//=======================================================================
void IntCurveSurface_ComputeTransitions(const TheCurve& curve, void IntCurveSurface_ComputeTransitions(const TheCurve& curve,
const Standard_Real w, const Standard_Real w,
IntCurveSurface_TransitionOnCurve& TransOnCurve, IntCurveSurface_TransitionOnCurve& TransOnCurve,
@ -1670,7 +1706,10 @@ void IntCurveSurface_ComputeTransitions(const TheCurve& curve,
TransOnCurve = IntCurveSurface_Tangent; TransOnCurve = IntCurveSurface_Tangent;
} }
} }
//================================================================================ //=======================================================================
//function : IntCurveSurface_ComputeParamsOnQuadric
//purpose :
//=======================================================================
void IntCurveSurface_ComputeParamsOnQuadric(const TheSurface& surface, void IntCurveSurface_ComputeParamsOnQuadric(const TheSurface& surface,
const gp_Pnt& P, const gp_Pnt& P,
Standard_Real& u, Standard_Real& u,
@ -1702,42 +1741,23 @@ void IntCurveSurface_ComputeParamsOnQuadric(const TheSurface& surface,
#endif #endif
} }
} }
//================================================================================ //=======================================================================
//function : EstLimForInfExtr
//purpose : Estimation of limits for infinite surfaces
//=======================================================================
void EstLimForInfExtr(const gp_Lin& Line,
/* const TheSurface& surface,
const Standard_Boolean IsOffSurf,
Standard_Real u1,v1,u2,v2; const Standard_Integer nbsu,
u1 = TheSurfaceTool::FirstUParameter(surface); const Standard_Boolean U1inf,
v1 = TheSurfaceTool::FirstVParameter(surface); const Standard_Boolean U2inf,
u2 = TheSurfaceTool::LastUParameter(surface); const Standard_Boolean V1inf,
v2 = TheSurfaceTool::LastVParameter(surface); const Standard_Boolean V2inf,
Standard_Integer nbsu,nbsv; Standard_Real& U1new,
nbsu = TheSurfaceTool::NbSamplesU(surface,u1,u2); Standard_Real& U2new,
nbsv = TheSurfaceTool::NbSamplesV(surface,v1,v2); Standard_Real& V1new,
IntCurveSurface_ThePolyhedron polyhedron(surface,nbsu,nbsv,u1,v1,u2,v2); Standard_Real& V2new,
Standard_Boolean& NoIntersection)
*/
//====================================================================================
// Estimation of limits for infinite surfaces
//====================================================================================
//================================================================================
static void EstLimForInfExtr(const gp_Lin& Line,
const TheSurface& surface,
const Standard_Boolean IsOffSurf,
const Standard_Integer nbsu,
const Standard_Boolean U1inf,
const Standard_Boolean U2inf,
const Standard_Boolean V1inf,
const Standard_Boolean V2inf,
Standard_Real& U1new,
Standard_Real& U2new,
Standard_Real& V1new,
Standard_Real& V2new,
Standard_Boolean& NoIntersection)
{ {
@ -1904,8 +1924,6 @@ static void EstLimForInfExtr(const gp_Lin& Line,
return; return;
} }
//================================================================================
//======================================================================= //=======================================================================
//function : ProjectIntersectAndEstLim //function : ProjectIntersectAndEstLim
//purpose : project <theLine> and it's X-axe symmetric line to <thePln> and //purpose : project <theLine> and it's X-axe symmetric line to <thePln> and
@ -1914,20 +1932,19 @@ static void EstLimForInfExtr(const gp_Lin& Line,
// <theBasCurvProj>. // <theBasCurvProj>.
// Is called from EstLimForInfRevl() // Is called from EstLimForInfRevl()
//======================================================================= //=======================================================================
void ProjectIntersectAndEstLim(const gp_Lin& theLine,
static void ProjectIntersectAndEstLim(const gp_Lin& theLine, const gp_Pln& thePln,
const gp_Pln& thePln, const ProjLib_Plane& theBasCurvProj,
const ProjLib_Plane& theBasCurvProj, Standard_Real& theVmin,
Standard_Real& theVmin, Standard_Real& theVmax,
Standard_Real& theVmax, Standard_Boolean& theNoIntersection)
Standard_Boolean& theNoIntersection)
{ {
ProjLib_Plane aLineProj( thePln, theLine ); ProjLib_Plane aLineProj( thePln, theLine );
if (!aLineProj.IsDone()) { if (!aLineProj.IsDone()) {
#ifdef DEB #ifdef DEB
cout cout
<< "Info: IntCurveSurface_Inter::ProjectIntersectAndEstLim(), !aLineProj.IsDone()" << "Info: IntCurveSurface_Inter::ProjectIntersectAndEstLim(), !aLineProj.IsDone()"
<< endl; << endl;
#endif #endif
return; return;
} }
@ -1995,18 +2012,17 @@ static void ProjectIntersectAndEstLim(const gp_Lin& theLine,
// Algo: intersect projections of Line and basis curve on the // Algo: intersect projections of Line and basis curve on the
// plane passing through revolution axe // plane passing through revolution axe
//======================================================================= //=======================================================================
void EstLimForInfRevl(const gp_Lin& Line,
static void EstLimForInfRevl(const gp_Lin& Line, const TheSurface& surface,
const TheSurface& surface, const Standard_Boolean U1inf,
const Standard_Boolean U1inf, const Standard_Boolean U2inf,
const Standard_Boolean U2inf, const Standard_Boolean V1inf,
const Standard_Boolean V1inf, const Standard_Boolean V2inf,
const Standard_Boolean V2inf, Standard_Real& U1new,
Standard_Real& U1new, Standard_Real& U2new,
Standard_Real& U2new, Standard_Real& V1new,
Standard_Real& V1new, Standard_Real& V2new,
Standard_Real& V2new, Standard_Boolean& NoIntersection)
Standard_Boolean& NoIntersection)
{ {
NoIntersection = Standard_False; NoIntersection = Standard_False;
@ -2096,19 +2112,22 @@ static void EstLimForInfRevl(const gp_Lin& Line,
return; return;
} }
//================================================================================ //=======================================================================
static void EstLimForInfOffs(const gp_Lin& Line, //function : EstLimForInfOffs
const TheSurface& surface, //purpose :
const Standard_Integer nbsu, //=======================================================================
const Standard_Boolean U1inf, void EstLimForInfOffs(const gp_Lin& Line,
const Standard_Boolean U2inf, const TheSurface& surface,
const Standard_Boolean V1inf, const Standard_Integer nbsu,
const Standard_Boolean V2inf, const Standard_Boolean U1inf,
Standard_Real& U1new, const Standard_Boolean U2inf,
Standard_Real& U2new, const Standard_Boolean V1inf,
Standard_Real& V1new, const Standard_Boolean V2inf,
Standard_Real& V2new, Standard_Real& U1new,
Standard_Boolean& NoIntersection) Standard_Real& U2new,
Standard_Real& V1new,
Standard_Real& V2new,
Standard_Boolean& NoIntersection)
{ {
NoIntersection = Standard_False; NoIntersection = Standard_False;
@ -2311,117 +2330,93 @@ static void EstLimForInfOffs(const gp_Lin& Line,
V1new = Max(V1new, -1.e10); V1new = Max(V1new, -1.e10);
V2new = Min(V2new, 1.e10); V2new = Min(V2new, 1.e10);
} }
} }
//=======================================================================
//================================================================================ //function : EstLimForInfSurf
static void EstLimForInfSurf(Standard_Real& U1new, //purpose :
Standard_Real& U2new, //=======================================================================
Standard_Real& V1new, void EstLimForInfSurf(Standard_Real& U1new,
Standard_Real& V2new) Standard_Real& U2new,
Standard_Real& V1new,
Standard_Real& V2new)
{ {
U1new = Max(U1new, -1.e10); U1new = Max(U1new, -1.e10);
U2new = Min(U2new, 1.e10); U2new = Min(U2new, 1.e10);
V1new = Max(V1new, -1.e10); V1new = Max(V1new, -1.e10);
V2new = Min(V2new, 1.e10); V2new = Min(V2new, 1.e10);
} }
//=======================================================================
#if 0 //function : DoTrim
//-- jgv patch (from) //purpose :
static Handle(Geom_Curve) GetCurve(const Handle(Adaptor3d_HCurve) AdCurve) //=======================================================================
Standard_Boolean DoTrim(const TheCurve& curve,
const TheSurface& surface)
{ {
Handle(Geom_Curve) theCurve; Standard_Boolean isAnaCurve, isAnaSurface, result;
GeomAbs_CurveType CurveType = AdCurve->GetType(); GeomAbs_CurveType aCT;
switch (CurveType) GeomAbs_SurfaceType aST;
{
case GeomAbs_Line: isAnaCurve = Standard_False;
theCurve = new Geom_Line( AdCurve->Line() ); isAnaSurface = Standard_False;
break; //
case GeomAbs_Circle: aCT=TheCurveTool::GetType(curve);
theCurve = new Geom_Circle( AdCurve->Circle() ); //
break; isAnaCurve= (aCT==GeomAbs_Line ||
case GeomAbs_Ellipse: aCT==GeomAbs_Circle ||
theCurve = new Geom_Ellipse( AdCurve->Ellipse() ); aCT==GeomAbs_Ellipse ||
break; aCT==GeomAbs_Hyperbola ||
case GeomAbs_Hyperbola: aCT==GeomAbs_Parabola);
theCurve = new Geom_Hyperbola( AdCurve->Hyperbola() ); //
break; aST=TheSurfaceTool::GetType(surface);
case GeomAbs_Parabola: isAnaSurface=(aST==GeomAbs_Plane ||
theCurve = new Geom_Parabola( AdCurve->Parabola() ); aST==GeomAbs_Cylinder ||
break; aST==GeomAbs_Cone ||
case GeomAbs_BezierCurve: aST==GeomAbs_Sphere ||
theCurve = AdCurve->Bezier(); aST==GeomAbs_Torus);
break; //
case GeomAbs_BSplineCurve: result = !(isAnaCurve && isAnaSurface);
theCurve = AdCurve->BSpline(); //
break; if(result) {
Standard_Boolean isUClosed, isVClosed;
isUClosed=(TheSurfaceTool::IsUClosed(surface) ||
TheSurfaceTool::IsUPeriodic(surface));
isVClosed=(TheSurfaceTool::IsVClosed(surface)
|| TheSurfaceTool::IsVPeriodic(surface));
if(isUClosed && isVClosed) {
result =!result;
} }
if (!theCurve.IsNull()) }
{ //modified by NIZNHY-PKV Wed Apr 04 12:11:01 2012f
Standard_Real f = AdCurve->FirstParameter(); if (result) {
Standard_Real l = AdCurve->LastParameter(); Standard_Real aU0,aU1,aV0,aV1, dU, dV, dA, dB, aR, aTresh;
theCurve = new Geom_TrimmedCurve( theCurve, f, l ); //
aTresh=100.;
aV0 = TheSurfaceTool::FirstVParameter(surface);
aV1 = TheSurfaceTool::LastVParameter(surface);
aU0 = TheSurfaceTool::FirstUParameter(surface);
aU1 = TheSurfaceTool::LastUParameter(surface);
dU=aU1-aU0;
dV=aV1-aV0;
//
dA=dU;
dB=dV;
if (dV>dU) {
dA=dV;
dB=dU;
} }
return theCurve; //
aR=dA/dB;
if (dB<Precision::PConfusion()) {
result=!result;
return result;
}
if (aR>aTresh) {
result=!result;
}
}
//modified by NIZNHY-PKV Wed Apr 04 12:11:03 2012t
return result;
} }
static Handle(Geom_Surface) GetSurface(const Handle(Adaptor3d_HSurface) AdSurface)
{
Handle(Geom_Surface) theSurface;
GeomAbs_SurfaceType SurfaceType = AdSurface->GetType();
switch (SurfaceType)
{
case GeomAbs_Plane:
theSurface = new Geom_Plane( AdSurface->Plane() );
break;
case GeomAbs_Cylinder:
theSurface = new Geom_CylindricalSurface( AdSurface->Cylinder() );
break;
case GeomAbs_Cone:
theSurface = new Geom_ConicalSurface( AdSurface->Cone() );
break;
case GeomAbs_Torus:
theSurface = new Geom_ToroidalSurface( AdSurface->Torus() );
break;
case GeomAbs_Sphere:
theSurface = new Geom_SphericalSurface( AdSurface->Sphere() );
break;
case GeomAbs_BezierSurface:
theSurface = AdSurface->Bezier();
break;
case GeomAbs_BSplineSurface:
theSurface = AdSurface->BSpline();
break;
case GeomAbs_SurfaceOfRevolution:
{
gp_Ax1 Axis = AdSurface->AxeOfRevolution();
Handle(Adaptor3d_HCurve) AdBC = AdSurface->BasisCurve();
Handle(Geom_Curve) BC = GetCurve(AdBC);
if (!BC.IsNull())
theSurface = new Geom_SurfaceOfRevolution( BC, Axis );
break;
}
case GeomAbs_SurfaceOfExtrusion:
{
gp_Dir Direction = AdSurface->Direction();
Handle(Adaptor3d_HCurve) AdBC = AdSurface->BasisCurve();
Handle(Geom_Curve) BC = GetCurve(AdBC);
if (!BC.IsNull())
theSurface = new Geom_SurfaceOfLinearExtrusion( BC, Direction );
break;
}
}
if (!theSurface.IsNull())
{
Standard_Real uf = AdSurface->FirstUParameter();
Standard_Real ul = AdSurface->LastUParameter();
Standard_Real vf = AdSurface->FirstVParameter();
Standard_Real vl = AdSurface->LastVParameter();
theSurface = new Geom_RectangularTrimmedSurface( theSurface, uf, ul, vf, vl );
}
return theSurface;
}
//-- jgv patch (to)
#endif