OpenCascade/occt
1
0
Fork 0
mirror of https://git.dev.opencascade.org/repos/occt.git synced 2025-06-20 11:54:07 +03:00
Code

Modeling - Periodic BSpline curve bounding

Browse Source 
Enhance periodic curve handling in BndLib_Add3dCurve::Add method
This commit is contained in:
dpasukhi 2025-04-11 16:56:17 +01:00
parent 7c8a432c8e
commit 4f3be4a125
1 changed files with 24 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
src/ModelingData/TKGeomBase/BndLib/BndLib_Add3dCurve.cxx
View File

@ -215,12 +215,32 @@ void BndLib_Add3dCurve::Add(const Adaptor3d_Curve& C,
u2 = Bsaux->LastParameter();
// modified by NIZHNY-EAP Fri Dec 3 14:29:18 1999 ___END___
}
const Standard_Real anU2Periodic =
Bsaux->IsPeriodic() ? (u2 - (Bsaux->LastParameter() - Bsaux->FirstParameter())) : u2;
Standard_Real aSegmentTol = 2. * Precision::PConfusion();
if (Abs((anU2Periodic - u1) < aSegmentTol))
// For periodic curves, check if parameters are close in either direction
if (Bsaux->IsPeriodic())
{
aSegmentTol = Abs(anU2Periodic - u1) * 0.01;
const Standard_Real aPeriod = Bsaux->LastParameter() - Bsaux->FirstParameter();
// Check direct distance between parameters
const Standard_Real aDirectDiff = Abs(u2 - u1);
// Check distances across period boundary (in both directions)
const Standard_Real aCrossPeriodDiff1 = Abs(u2 - aPeriod - u1);
const Standard_Real aCrossPeriodDiff2 = Abs(u1 - aPeriod - u2);
// Find the minimum difference (closest approach)
const Standard_Real aMinDiff = Min(aDirectDiff, Min(aCrossPeriodDiff1, aCrossPeriodDiff2));
if (aMinDiff < aSegmentTol)
{
aSegmentTol = aMinDiff * 0.01;
}
}
// For non-periodic curves, just check direct parameter difference
else if (Abs(u2 - u1) < aSegmentTol)
{
aSegmentTol = Abs(u2 - u1) * 0.01;
}
Bsaux->Segment(u1, u2, aSegmentTol);
Bs = Bsaux;