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

0023158: ApproxInt_PrmPrmSvSurfaces raises FPE (division by zero) signal

Browse Source 
Prevent division by zero in ApproxInt_PrmPrmSvSurfaces::Compute
This commit is contained in:
vsr 2012-05-17 10:24:01 +04:00
parent ef8da89204
commit 416594fec8
1 changed files with 11 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
src/ApproxInt/ApproxInt_PrmPrmSvSurfaces.gxx
View File

@ -24,6 +24,7 @@
#include <TColStd_Array1OfReal.hxx> #include <TColStd_Array1OfReal.hxx>
#include <math_FunctionSetRoot.hxx> #include <math_FunctionSetRoot.hxx>
#include <Precision.hxx>
#define Debug(expr) cout<<" expr :"<<expr; #define Debug(expr) cout<<" expr :"<<expr;
#define MySurf1 MyIntersectionOn2S.Function().AuxillarSurface1() #define MySurf1 MyIntersectionOn2S.Function().AuxillarSurface1()
@ -173,6 +174,11 @@ Standard_Boolean ApproxInt_PrmPrmSvSurfaces::Compute( Standard_Real& u1
TgTU = Tg.Dot(TU); TgTU = Tg.Dot(TU);
TgTV = Tg.Dot(TV); TgTV = Tg.Dot(TV);
DIS = TUTU * TVTV - TUTV * TUTV; DIS = TUTU * TVTV - TUTV * TUTV;
if(fabs(DIS)<Precision::Angular()) {
MyIsTangent=Standard_False;
MyHasBeenComputed = MyHasBeenComputedbis = Standard_False;
return(Standard_False);
}
DeltaU = (TgTU * TVTV - TgTV * TUTV ) / DIS ; DeltaU = (TgTU * TVTV - TgTV * TUTV ) / DIS ;
DeltaV = (TgTV * TUTU - TgTU * TUTV ) / DIS ; DeltaV = (TgTV * TUTU - TgTU * TUTV ) / DIS ;
@ -190,6 +196,11 @@ Standard_Boolean ApproxInt_PrmPrmSvSurfaces::Compute( Standard_Real& u1
TgTU = Tg.Dot(TU); TgTU = Tg.Dot(TU);
TgTV = Tg.Dot(TV); TgTV = Tg.Dot(TV);
DIS = TUTU * TVTV - TUTV * TUTV; DIS = TUTU * TVTV - TUTV * TUTV;
if(fabs(DIS)<Precision::Angular()) {
MyIsTangent=Standard_False;
MyHasBeenComputed = MyHasBeenComputedbis = Standard_False;
return(Standard_False);
}
DeltaU = (TgTU * TVTV - TgTV * TUTV ) / DIS ; DeltaU = (TgTU * TVTV - TgTV * TUTV ) / DIS ;
DeltaV = (TgTV * TUTU - TgTU * TUTV ) / DIS ; DeltaV = (TgTV * TUTU - TgTU * TUTV ) / DIS ;