From c08113c2798b907dce8efbea41b5f6545c71d913 Mon Sep 17 00:00:00 2001 From: Martin Kroeker Date: Sun, 23 Jun 2024 01:12:33 +0200 Subject: [PATCH] fix special cases of x= NAN or INF --- kernel/mips/dscal_msa.c | 10 ++++++---- 1 file changed, 6 insertions(+), 4 deletions(-) diff --git a/kernel/mips/dscal_msa.c b/kernel/mips/dscal_msa.c index 2e41d8bef..e95f0a655 100644 --- a/kernel/mips/dscal_msa.c +++ b/kernel/mips/dscal_msa.c @@ -42,7 +42,7 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da, FLOAT *x, if (1 == inc_x) { - if (0.0 == da) + if (0) //if (0.0 == da ) { v2f64 zero_v = {0.0, 0.0}; @@ -243,9 +243,11 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da, FLOAT *x, if (da == 0.0) { for (i = n; i--;) - { - *x = 0.0; - + { + if (isfinite(*x)) + *x = 0.0; + else + *x = NAN; x += inc_x; } }