additional fixes for handling INF arguments
This commit is contained in:
parent
ce130f11d2
commit
9ff4e9714e
|
@ -61,7 +61,9 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da_r,FLOAT da_i, F
|
||||||
{
|
{
|
||||||
temp = - da_i * x[ip+1] ;
|
temp = - da_i * x[ip+1] ;
|
||||||
if (isnan(x[ip]) || isinf(x[ip])) temp = NAN;
|
if (isnan(x[ip]) || isinf(x[ip])) temp = NAN;
|
||||||
x[ip+1] = da_i * x[ip] ;
|
if (!isinf(x[ip+1]))
|
||||||
|
x[ip+1] = da_i * x[ip] ;
|
||||||
|
else x[ip+1] = NAN;
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
else
|
else
|
||||||
|
|
Loading…
Reference in New Issue