| @@ -260,7 +260,7 @@ 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+1])) | if (!isinf(x[i+1])) | ||||
| x[i+1] = da_i * x[i]; | x[i+1] = da_i * x[i]; | ||||
| else x[i+1] = NaN; | |||||
| else x[i+1] = NAN; | |||||
| x[i] = temp0; | x[i] = temp0; | ||||
| if (isnan(x[i+inc_x]) || isinf(x[i+inc_x])) | if (isnan(x[i+inc_x]) || isinf(x[i+inc_x])) | ||||
| temp1 = NAN; | temp1 = NAN; | ||||
| @@ -268,7 +268,7 @@ int CNAME(BLASLONG n, BLASLONG dummy0, BLASLONG dummy1, FLOAT da_r, FLOAT da_i, | |||||
| temp1 = -da_i * x[i+1+inc_x]; | temp1 = -da_i * x[i+1+inc_x]; | ||||
| if (!isinf(x[i+1+inc_x])) | if (!isinf(x[i+1+inc_x])) | ||||
| x[i+1+inc_x] = da_i * x[i+inc_x]; | x[i+1+inc_x] = da_i * x[i+inc_x]; | ||||
| else x[i+1+inc_x] = NaN; | |||||
| else x[i+1+inc_x] = NAN; | |||||
| x[i+inc_x] = temp1; | x[i+inc_x] = temp1; | ||||
| i += 2*inc_x ; | i += 2*inc_x ; | ||||
| j+=2; | j+=2; | ||||