From 71553f5e5626e1d791b5be24c84ea6b17ae81cc7 Mon Sep 17 00:00:00 2001
From: Pjotr Prins
Date: Sun, 22 Aug 2021 11:48:43 +0200
Subject: Disable round-off correction

---
 src/lmm.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

(limited to 'src/lmm.cpp')

diff --git a/src/lmm.cpp b/src/lmm.cpp
index 2055c45..cc2bd33 100644
--- a/src/lmm.cpp
+++ b/src/lmm.cpp
@@ -846,7 +846,7 @@ double LogRL_f(double l, void *params) {
   }
   index_ww = GetabIndex(n_cvt + 2, n_cvt + 2, n_cvt);
   double P_yy = gsl_matrix_safe_get(Pab, nc_total, index_ww);
-  if (P_yy < 0.00000001) P_yy = 0.00000001;
+  // if (P_yy > 0.0 && P_yy < 0.00000001) P_yy = 0.00000001; // control potential round-off
 
   double c = 0.5 * df * (safe_log(df) - safe_log(2 * M_PI) - 1.0);
   f = c - 0.5 * logdet_h - 0.5 * logdet_hiw - 0.5 * df * safe_log(P_yy);
-- 
cgit v1.2.3