From e760e53fac224f252db47dfe7c37f76d6bebb4dc Mon Sep 17 00:00:00 2001
From: Pjotr Prins
Date: Sun, 26 Aug 2018 09:24:26 +0000
Subject: Make sure tau_inv is initialized outside loop - gets rid of compiler
 warning

---
 src/vc.cpp | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

(limited to 'src')

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.
-- 
cgit v1.2.3