From 6052fd97a43b5ff950aa744f2edadafca6de02f9 Mon Sep 17 00:00:00 2001 From: Martin Kroeker Date: Wed, 15 Mar 2023 21:12:26 +0100 Subject: [PATCH] compare agaist DBLMIN rather than a somewhat arbitrary value near it --- lapack/getf2/getf2_k.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/lapack/getf2/getf2_k.c b/lapack/getf2/getf2_k.c index 5d3b4fdc2..11dea4439 100644 --- a/lapack/getf2/getf2_k.c +++ b/lapack/getf2/getf2_k.c @@ -100,7 +100,7 @@ blasint CNAME(blas_arg_t *args, BLASLONG *range_m, BLASLONG *range_n, FLOAT *sa, temp1 = *(b + jp); //if (temp1 != ZERO) { - if (fabs(temp1) > 1.e-305) { + if (fabs(temp1) >= DBL_MIN temp1 = dp1 / temp1; if (jp != j) {