Fix handling of NAN and INF

martin-frbg · web-flow · commit d3d99c34f24a · 2024-01-07T17:56:51.000+01:00
diff --git a/kernel/mips/zscal.c b/kernel/mips/zscal.c
@@ -63,6 +63,7 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da_r,FLOAT da_i, F
 				x[ip+1] = da_r * x[ip+1] + da_i * x[ip]   ;
 			}
 		}
+		if (!isnan(x_ip] && !isinf(x[ip]))
 		x[ip]   = temp;
 
 		ip += inc_x2;

Original file line number	Diff line number	Diff line change
`@@ -63,6 +63,7 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da_r,FLOAT da_i, F`
`63`	`63`	`x[ip+1] = da_r * x[ip+1] + da_i * x[ip] ;`
`64`	`64`	`}`
`65`	`65`	`}`
	`66`	`+ if (!isnan(x_ip] && !isinf(x[ip]))`
`66`	`67`	`x[ip] = temp;`
`67`	`68`
`68`	`69`	`ip += inc_x2;`