compare agaist DBLMIN rather than a somewhat arbitrary value near it

This commit is contained in:
Martin Kroeker 2023-03-15 21:12:26 +01:00 committed by GitHub
parent 8d6813ec41
commit 6052fd97a4
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 1 additions and 1 deletions

View File

@ -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) {