floraachy
/
OpenBLAS
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

fix another corner case involving infinity

This commit is contained in:
Martin Kroeker 2024-05-31 01:06:58 +02:00 committed by GitHub
parent b9a1c9a06c
commit ad2b5c67c8
No known key found for this signature in database
GPG Key ID: B5690EEEBB952194
1 changed files with 6 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
kernel/x86_64/zscal.c
View File

@ -278,7 +278,9 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da_r, FLOAT da_i,
temp0 = NAN; temp0 = NAN;
else else
temp0 = -da_i * x[i+1]; temp0 = -da_i * x[i+1];
x[i+1] = da_i * x[i]; if (!isinf(x[i+1]))
x[i+1] = da_i * x[i];
else x[i+1] = NAN;
x[i] = temp0; x[i] = temp0;
i += inc_x ; i += inc_x ;
j++; j++;
@ -412,7 +414,9 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da_r, FLOAT da_i,
temp0 = -da_i * x[i+1]; temp0 = -da_i * x[i+1];
if (isinf(x[i])) if (isinf(x[i]))
temp0 = NAN; temp0 = NAN;
x[i+1] = da_i * x[i]; if (!isinf(x[i+1]))
x[i+1] = da_i * x[i];
else x[i+1] = NAN;
if ( x[i] == x[i]) //preserve NaN if ( x[i] == x[i]) //preserve NaN
x[i] = temp0; x[i] = temp0;
i += 2 ; i += 2 ;