OpenCascade/occt
1
0
Fork 0
mirror of https://git.dev.opencascade.org/repos/occt.git synced 2025-04-03 17:56:21 +03:00
Code Packages Releases Activity

0031616: Modeling algorithm - Section between two shells returns wire with gaps (720)

Browse Source 
Adaptor3d_TopolTool.cxx, IntTools_TopolTool.cxx:

Anisotropy of BSpline surface along U and V direction is taken in account for calculation of numbers of sample points;

bug31616: test case added
This commit is contained in:
ifv 2020-06-17 17:10:12 +03:00 committed by bugmaster
parent 832a6f4412
commit 79e9ba31d2
3 changed files with 62 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

45
src/Adaptor3d/Adaptor3d_TopolTool.cxx
View File

@ -654,7 +654,10 @@ static void Analyse(const TColgp_Array2OfPnt& array2,
}
void Adaptor3d_TopolTool::ComputeSamplePoints() {
void Adaptor3d_TopolTool::ComputeSamplePoints()
{
const Standard_Integer aMaxNbSample = 50;
Standard_Real uinf,usup,vinf,vsup;
uinf = myS->FirstUParameter(); usup = myS->LastUParameter();
vinf = myS->FirstVParameter(); vsup = myS->LastVParameter();
@ -692,31 +695,49 @@ void Adaptor3d_TopolTool::ComputeSamplePoints() {
//--
if(nbsu<6) nbsu=6;
if(nbsv<6) nbsv=6;
myNbSamplesU = nbsu;
myNbSamplesV = nbsv;
if(nbsv<6) nbsv=6;
if(nbsu>8 || nbsv>8) {
if(typS == GeomAbs_BSplineSurface) {
if (typS == GeomAbs_BSplineSurface) {
if (nbsu > 8 || nbsv>8) {
const Handle(Geom_BSplineSurface)& Bspl = myS->BSpline();
Standard_Integer nbup = Bspl->NbUPoles();
Standard_Integer nbvp = Bspl->NbVPoles();
TColgp_Array2OfPnt array2(1,nbup,1,nbvp);
TColgp_Array2OfPnt array2(1, nbup, 1, nbvp);
Bspl->Poles(array2);
Analyse(array2,nbup,nbvp,myNbSamplesU,myNbSamplesV);
Analyse(array2, nbup, nbvp, nbsu, nbsv);
}
else if(typS == GeomAbs_BezierSurface) {
// Check anisotropy
Standard_Real anULen = (usup - uinf) / myS->UResolution(1.);
Standard_Real anVLen = (vsup - vinf) / myS->VResolution(1.);
Standard_Real aRatio = anULen / anVLen;
if (aRatio >= 10.)
{
nbsu *= 2;
nbsu = Min(nbsu, aMaxNbSample);
}
else if (aRatio <= 0.1)
{
nbsv *= 2;
nbsv = Min(nbsv, aMaxNbSample);
}
}
else if (typS == GeomAbs_BezierSurface) {
if (nbsu > 8 || nbsv > 8) {
const Handle(Geom_BezierSurface)& Bez = myS->Bezier();
Standard_Integer nbup = Bez->NbUPoles();
Standard_Integer nbvp = Bez->NbVPoles();
TColgp_Array2OfPnt array2(1,nbup,1,nbvp);
TColgp_Array2OfPnt array2(1, nbup, 1, nbvp);
Bez->Poles(array2);
Analyse(array2,nbup,nbvp,myNbSamplesU,myNbSamplesV);
Analyse(array2, nbup, nbvp, nbsu, nbsv);
}
}
myNbSamplesU = nbsu;
myNbSamplesV = nbsv;
}
Standard_Integer Adaptor3d_TopolTool::NbSamplesU()
{
if(myNbSamplesU <0) {

14
src/IntTools/IntTools_TopolTool.cxx
View File

@ -246,6 +246,20 @@ void IntTools_TopolTool::ComputeSamplePoints()
}
if(nbsu < 10) nbsu = 10;
if(nbsv < 10) nbsv = 10;
// Check anisotropy
Standard_Real anULen = (usup - uinf) / myS->UResolution(1.);
Standard_Real anVLen = (vsup - vinf) / myS->VResolution(1.);
Standard_Real aRatio = anULen / anVLen;
if (aRatio >= 10.)
{
nbsu *= 2;
nbsu = Min(nbsu, aMaxNbSample);
}
else if (aRatio <= 0.1 )
{
nbsv *= 2;
nbsv = Min(nbsv, aMaxNbSample);
}
}
break;
case GeomAbs_SurfaceOfExtrusion: {

15
tests/bugs/modalg_7/bug31616 Normal file
View File

@ -0,0 +1,15 @@
puts "============================================================================================="
puts "OCC31616: Modeling algorithm - Section between two shells returns wire with gaps (720)"
puts "============================================================================================="
puts ""
restore [locate_data_file bug31616.brep] s
explode s
bsection sec s_1 s_2
checksection sec -r 2
checknbshapes sec -vertex 69 -edge 68
checkprops sec -l 20.7581