fix special cases of x= NAN or INF

martin-frbg · web-flow · commit c08113c2798b · 2024-06-23T01:12:33.000+02:00
diff --git 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;
             }
         }

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`	`v2f64 zero_v = {0.0, 0.0};`
`48`	`48`
`@@ -243,9 +243,11 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da, FLOAT *x,`
`243`	`243`	`if (da == 0.0)`
`244`	`244`	`{`
`245`	`245`	`for (i = n; i--;)`
`246`		`- {`
`247`		`- *x = 0.0;`
`248`		`-`
	`246`	`+ {`
	`247`	`+ if (isfinite(*x))`
	`248`	`+ *x = 0.0;`
	`249`	`+ else`
	`250`	`+ *x = NAN;`
`249`	`251`	`x += inc_x;`
`250`	`252`	`}`
`251`	`253`	`}`