|
|
|
@@ -266,7 +266,7 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da_r, FLOAT da_i, |
|
|
|
} else { |
|
|
|
|
|
|
|
while (j < n1) { |
|
|
|
if (isnan(x[i]) || isinf(x[i]) || isnan(x[i+1])) |
|
|
|
if (isnan(x[i]) || isinf(x[i])) |
|
|
|
temp0 = NAN; |
|
|
|
else |
|
|
|
temp0 = -da_i * x[i + 1]; |
|
|
|
@@ -275,7 +275,7 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da_r, FLOAT da_i, |
|
|
|
else |
|
|
|
x[i + 1] = NAN; |
|
|
|
x[i] = temp0; |
|
|
|
if (isnan(x[i+inc_x]) || isinf(x[i+inc_x]) || isnan(x[i+1+inc_x])) |
|
|
|
if (isnan(x[i+inc_x]) || isinf(x[i+inc_x])) |
|
|
|
temp1 = NAN; |
|
|
|
else |
|
|
|
temp1 = -da_i * x[i + 1 + inc_x]; |
|
|
|
@@ -290,7 +290,7 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da_r, FLOAT da_i, |
|
|
|
} |
|
|
|
|
|
|
|
while (j < n) { |
|
|
|
if (isnan(x[i]) || isinf(x[i]) || isnan(x[i+1])) |
|
|
|
if (isnan(x[i]) || isinf(x[i])) |
|
|
|
temp0 = NAN; |
|
|
|
else |
|
|
|
temp0 = -da_i * x[i + 1]; |
|
|
|
|