eigen_solution_of_triangular_matrix: hanged size_t to long

src/Eof3d.cc

@@ -163,9 +163,9 @@ EOF3d(void *process)
       eigenvectors[varID].resize(n_eig);
       eigenvalues[varID].resize(nts);
 
-      for (int i = 0; i < n; i++)
+      for (size_t i = 0; i < n; i++)
         {
-          if (i < n_eig) eigenvectors[varID][i].resize(temp_size, missval);
+          if (i < (size_t)n_eig) eigenvectors[varID][i].resize(temp_size, missval);
 
           eigenvalues[varID][i].resize(1, missval);
         }
@@ -318,7 +318,7 @@ EOF3d(void *process)
 
       if (Options::cdoVerbose) cdoPrint("Processed SVD decomposition for var %d from %d x %d matrix", varID, n, n);
 
-      for (int eofID = 0; eofID < n; eofID++) eigenvalues[varID][eofID][0] = eigv[eofID];
+      for (size_t eofID = 0; eofID < n; eofID++) eigenvalues[varID][eofID][0] = eigv[eofID];
 
       if (Options::Timer) timer_stop(timer_eig);
 
@@ -408,7 +408,7 @@ EOF3d(void *process)
   int64_t vdate = 10101;
   int vtime = 0;
   auto juldate = julianDateEncode(calendar, vdate, vtime);
-  for (tsID = 0; tsID < n; tsID++)
+  for (tsID = 0; tsID < (int)n; tsID++)
     {
       juldate = julianDateAddSeconds(60, juldate);
       julianDateDecode(calendar, juldate, vdate, vtime);

src/statistic.cc

@@ -187,10 +187,10 @@ pythagoras(double a, double b)
 #define MAX_ITER 1000
 
 static void
-eigen_solution_of_triangular_matrix(Varray<double> &d, Varray<double> &e, size_t n, Varray2D<double> &a, const char *prompt)
+eigen_solution_of_triangular_matrix(Varray<double> &d, Varray<double> &e, long n, Varray2D<double> &a, const char *prompt)
 {
   constexpr double eps = 1.e-6;
-  size_t i, k, l, m, iter;
+  long i, k, l, m, iter;
   double b, c, f, g, p, r, s;
 
   for (i = 1; i < n; i++) e[i - 1] = e[i];