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

---
 lapack/getf2/zgetf2_k.c | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/lapack/getf2/zgetf2_k.c b/lapack/getf2/zgetf2_k.c
index 512adf8ba..466b66e26 100644
--- a/lapack/getf2/zgetf2_k.c
+++ b/lapack/getf2/zgetf2_k.c
@@ -105,8 +105,9 @@ blasint CNAME(blas_arg_t *args, BLASLONG *range_m, BLASLONG *range_n, FLOAT *sa,
       temp1 = *(b + jp * 2 + 0);
       temp2 = *(b + jp * 2 + 1);
 
-      if ((temp1 != ZERO) || (temp2 != ZERO)) {
-
+  //    if ((temp1 != ZERO) || (temp2 != ZERO)) {
+	if ((fabs(temp1) > 1.e-305) || (fabs(temp2) > 1.e-305)) {
+		
 	if (jp != j) {
 	  SWAP_K(j + 1, 0, 0, ZERO, ZERO, a + j * 2, lda,
 		 a + jp * 2, lda, NULL, 0);