fix handling of INF or NAN

This commit is contained in:
Martin Kroeker 2024-07-04 17:47:19 +02:00 committed by GitHub
parent b422742899
commit 536200bc9e
No known key found for this signature in database
GPG Key ID: B5690EEEBB952194
1 changed files with 4 additions and 0 deletions

View File

@ -48,7 +48,11 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da, FLOAT *x, BLAS
{
if ( da == 0.0 )
if (!isnan(x[i] && !isinf(x[i])) {
x[i]=0.0;
} else {
x[i]=NAN;
}
else
x[i] = da * x[i] ;