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.

1 files changed, 1 insertions, 1 deletions

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);
       }