From 416594fec8f40cd05e97e0c789dc7c78f1afc11a Mon Sep 17 00:00:00 2001
From: vsr <vsr@opencascade.com>
Date: Thu, 17 May 2012 10:24:01 +0400
Subject: [PATCH] 0023158: ApproxInt_PrmPrmSvSurfaces raises FPE (division by
 zero) signal

Prevent division by zero in ApproxInt_PrmPrmSvSurfaces::Compute
---
 src/ApproxInt/ApproxInt_PrmPrmSvSurfaces.gxx | 11 +++++++++++
 1 file changed, 11 insertions(+)

diff --git a/src/ApproxInt/ApproxInt_PrmPrmSvSurfaces.gxx b/src/ApproxInt/ApproxInt_PrmPrmSvSurfaces.gxx
index 028e52db73..785069b232 100755
--- a/src/ApproxInt/ApproxInt_PrmPrmSvSurfaces.gxx
+++ b/src/ApproxInt/ApproxInt_PrmPrmSvSurfaces.gxx
@@ -24,6 +24,7 @@
 
 #include <TColStd_Array1OfReal.hxx>
 #include <math_FunctionSetRoot.hxx>
+#include <Precision.hxx>
 
 #define Debug(expr)  cout<<" expr :"<<expr;
 #define MySurf1 MyIntersectionOn2S.Function().AuxillarSurface1()
@@ -173,6 +174,11 @@ Standard_Boolean ApproxInt_PrmPrmSvSurfaces::Compute( Standard_Real& u1
   TgTU = Tg.Dot(TU);
   TgTV = Tg.Dot(TV);
   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 ; 
   DeltaV = (TgTV * TUTU - TgTU * TUTV ) / DIS ;
@@ -190,6 +196,11 @@ Standard_Boolean ApproxInt_PrmPrmSvSurfaces::Compute( Standard_Real& u1
   TgTU = Tg.Dot(TU);
   TgTV = Tg.Dot(TV);
   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 ; 
   DeltaV = (TgTV * TUTU - TgTU * TUTV ) / DIS ;