disable the fast path for inc=1, alpha=0 as it does not handle x=NaN or Inf

martin-frbg · web-flow · commit 541e1b695979 · 2024-06-23T10:37:55.000+02:00
diff --git a/kernel/mips/sscal_msa.c b/kernel/mips/sscal_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)
         {
             v4f32 zero_v = {0.0, 0.0, 0.0, 0.0};
 
@@ -259,7 +259,10 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da, FLOAT *x,
         {
             for (i = n; i--;)
             {
-                *x = 0;
+                if (isfinite(*x)
+                          *x = 0;
+                else
+                          *x = NAN;
                 x += inc_x;
             }
         }

Original file line number	Diff line number	Diff line change
`@@ -42,7 +42,7 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da, FLOAT *x,`
`42`	`42`
`43`	`43`	`if (1 == inc_x)`
`44`	`44`	`{`
`45`		`- if (0.0 == da)`
	`45`	`+ if (0) // if (0.0 == da)`
`46`	`46`	`{`
`47`	`47`	`v4f32 zero_v = {0.0, 0.0, 0.0, 0.0};`
`48`	`48`
`@@ -259,7 +259,10 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da, FLOAT *x,`
`259`	`259`	`{`
`260`	`260`	`for (i = n; i--;)`
`261`	`261`	`{`
`262`		`- *x = 0;`
	`262`	`+ if (isfinite(*x)`
	`263`	`+ *x = 0;`
	`264`	`+ else`
	`265`	`+ *x = NAN;`
`263`	`266`	`x += inc_x;`
`264`	`267`	`}`
`265`	`268`	`}`