Handle NAN and INF
This commit is contained in:
parent
a782103b9c
commit
0d2e486edf
|
@ -60,6 +60,7 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da_r,FLOAT da_i, F
|
||||||
else
|
else
|
||||||
{
|
{
|
||||||
temp = - da_i * x[ip+1] ;
|
temp = - da_i * x[ip+1] ;
|
||||||
|
if (isnan(x[ip]) || isinf(x[ip])) temp = NAN;
|
||||||
x[ip+1] = da_i * x[ip] ;
|
x[ip+1] = da_i * x[ip] ;
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
Loading…
Reference in New Issue