From ab434bbd38a61941bf32afcf91b7cc397462a53b Mon Sep 17 00:00:00 2001
From: Pjotr Prins
Date: Sun, 30 Nov 2025 08:53:10 +0100
Subject: Write out K

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

(limited to 'src/gemma.cpp')

diff --git a/src/gemma.cpp b/src/gemma.cpp
index 6b42763..1cced77 100644
--- a/src/gemma.cpp
+++ b/src/gemma.cpp
@@ -1915,7 +1915,16 @@ void GEMMA::BatchRun(PARAM &cPar) {
     cout << "Calculating Relatedness Matrix ... " << endl;
 
     if (cPar.is_mdb) {
-      cPar.MdbCalcKin();
+      time_start = clock();
+      auto K = cPar.MdbCalcKin();
+      cPar.time_G = (clock() - time_start) / (double(CLOCKS_PER_SEC) * 60.0);
+      if (cPar.a_mode == M_KIN) {
+        cPar.WriteMatrix(K, "cXX");
+      } else { // M_KIN2
+        cPar.WriteMatrix(K, "sXX");
+      }
+
+      gsl_matrix_safe_free(K);
       return;
     }
 
-- 
cgit 1.4.1