diff options
Diffstat (limited to 'src/lmm.cpp')
| -rw-r--r-- | src/lmm.cpp | 2 |
1 files changed, 1 insertions, 1 deletions
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); |