From c21dfe50df1c1399190a7c5c100c5d0b14ab7ef8 Mon Sep 17 00:00:00 2001
From: Pjotr Prins
Date: Fri, 28 Sep 2018 15:22:11 +0000
Subject: Fixes regression #179

---
 src/mvlmm.cpp | 10 +++++++++-
 1 file changed, 9 insertions(+), 1 deletion(-)

(limited to 'src/mvlmm.cpp')

diff --git a/src/mvlmm.cpp b/src/mvlmm.cpp
index a3c7c85..4b8db05 100644
--- a/src/mvlmm.cpp
+++ b/src/mvlmm.cpp
@@ -256,7 +256,15 @@ double EigenProc(const gsl_matrix *V_g, const gsl_matrix *V_e, gsl_vector *D_l,
   gsl_blas_dgemm(CblasNoTrans, CblasNoTrans, 1.0, V_e_hi, VgVehi, 0.0, Lambda);
 
   // Eigen decomposition of Lambda.
-  EigenDecomp_Zeroed(Lambda, U_l, D_l, 0);
+  EigenDecomp(Lambda, U_l, D_l, 0);
+
+  for (size_t i = 0; i < d_size; i++) {
+    d = gsl_vector_get(D_l, i);
+    if (d < 0) {
+      gsl_vector_set(D_l, i, 0);
+    }
+  }
+
 
   // Calculate UltVeh and UltVehi.
   gsl_blas_dgemm(CblasTrans, CblasNoTrans, 1.0, U_l, V_e_h, 0.0, UltVeh);
-- 
cgit v1.2.3