Make sure tau_inv is initialized outside loop - gets rid of compiler warning

1 files changed, 2 insertions, 1 deletions

diff --git a/src/vc.cpp b/src/vc.cpp
index 416a3e4..19590f7 100644
--- a/src/vc.cpp
+++ b/src/vc.cpp
@@ -1935,7 +1935,7 @@ void VC::CalcVCacl(const gsl_matrix *K, const gsl_matrix *W,
   size_t n1 = K->size1, n2 = K->size2;
   size_t n_vc = n2 / n1;
 
-  double d, y2_sum, tau_inv;
+  double d, y2_sum;
 
   // New matrices/vectors.
   gsl_matrix *K_scale = gsl_matrix_alloc(n1, n2);
@@ -2049,6 +2049,7 @@ void VC::CalcVCacl(const gsl_matrix *K, const gsl_matrix *W,
   // Iterate to solve tau and h's.
   size_t it = 0;
   double s = 1;
+  double tau_inv = y2_sum / (double)n1 - d / ((double)n1 * ((double)n1 - 1));
   while (abs(s) > 1e-3 && it < 100) {
 
     // Update tau_inv.