From a39ced05515a084fd9b2aeddf7e0b2eed18192b3 Mon Sep 17 00:00:00 2001 From: Martin Kroeker Date: Sun, 26 Feb 2023 23:42:20 +0100 Subject: [PATCH 1/2] 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); From 3d27cbd9a3c9d36cb920f61a3c2ccce42ca47232 Mon Sep 17 00:00:00 2001 From: Martin Kroeker Date: Sun, 26 Feb 2023 23:44:14 +0100 Subject: [PATCH 2/2] 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) {