From 3d27cbd9a3c9d36cb920f61a3c2ccce42ca47232 Mon Sep 17 00:00:00 2001
From: Martin Kroeker <martin@ruby.chemie.uni-freiburg.de>
Date: Sun, 26 Feb 2023 23:44:14 +0100
Subject: [PATCH] avoid overflow in division

---
 lapack/getf2/getf2_k.c | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/lapack/getf2/getf2_k.c b/lapack/getf2/getf2_k.c
index 073a0251d..5d3b4fdc2 100644
--- a/lapack/getf2/getf2_k.c
+++ b/lapack/getf2/getf2_k.c
@@ -99,7 +99,8 @@ blasint CNAME(blas_arg_t *args, BLASLONG *range_m, BLASLONG *range_n, FLOAT *sa,
       jp--;
       temp1 = *(b + jp);
 
-      if (temp1 != ZERO) {
+      //if (temp1 != ZERO) {
+	if (fabs(temp1) > 1.e-305) {    
 	temp1 = dp1 / temp1;
 
 	if (jp != j) {