Handle NAN in input

martin-frbg · web-flow · commit acf17a825dad · 2024-01-07T20:26:16.000+01:00
diff --git a/kernel/mips/zscal.c b/kernel/mips/zscal.c
@@ -63,8 +63,11 @@ 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]   ;
 			}
 		}
-		x[ip]   = temp;
-
+		if ( da_r != da_r ) 
+			x[ip] = da_r;
+		else
+			x[ip]   = temp;
+		
 		ip += inc_x2;
 	}
 

Original file line number	Diff line number	Diff line change
`@@ -63,8 +63,11 @@ 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`		`- x[ip] = temp;`
`67`		`-`
	`66`	`+ if ( da_r != da_r )`
	`67`	`+ x[ip] = da_r;`
	`68`	`+ else`
	`69`	`+ x[ip] = temp;`
	`70`	`+`
`68`	`71`	`ip += inc_x2;`
`69`	`72`	`}`
`70`	`73`