From 2b85e9dc5111a8f51170c745c449697d2c77df1e Mon Sep 17 00:00:00 2001
From: Artyom Bologov
Date: Fri, 16 Aug 2024 01:58:59 +0400
Subject: gemma.cpp: Do eigendecomposition the same way as with LMM.

This way, one can export the eigen* data that matches the one used
while doing LMM. Otherwise the exported eigenvectors are sliiiiightly
off and break the computations depending on them in Mgamma.
---
 src/gemma.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

(limited to 'src')

diff --git a/src/gemma.cpp b/src/gemma.cpp
index 2525075..4dfd7fd 100644
--- a/src/gemma.cpp
+++ b/src/gemma.cpp
@@ -2617,7 +2617,7 @@ void GEMMA::BatchRun(PARAM &cPar) {
       time_start = clock();
 
       if (cPar.a_mode == M_EIGEN) {
-        cPar.trace_G = EigenDecomp_Zeroed(G, U, eval, 1);
+        cPar.trace_G = EigenDecomp_Zeroed(G, U, eval, 0);
       } else {
         cPar.trace_G = EigenDecomp_Zeroed(G, U, eval, 0);
       }
-- 
cgit v1.2.3