Update cscal.c

Author: martin-frbg

kernel/zarch/cscal.c

@@ -234,7 +234,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]))
+          if (isnan(x[i]) || isinf(x[i]) || isnan(x[i+1])) 
             temp0   = NAN;
 	  else
 	    temp0 = -da_i * x[i + 1];
@@ -243,7 +243,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]))
+	  if (isnan(x[i+inc_x]) || isinf(x[i+inc_x]) || isnan(x[i+1+inc_x]))
 	    temp1 = NAN;
 	  else
 	    temp1 = -da_i * x[i + 1 + inc_x];
@@ -258,7 +258,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]))
+	  if (isnan(x[i]) || isinf(x[i]) || isnan(x[i+1]))
 	    temp0 = NAN;
 	  else
             temp0 = -da_i * x[i + 1];