From 4d1a44cb90affbe9c73a27767bae4a4b69c72402 Mon Sep 17 00:00:00 2001
From: Pjotr Prins
Date: Sat, 25 Aug 2018 09:48:03 +0000
Subject: segfpe handling

---
 src/gemma_io.cpp | 2 ++
 1 file changed, 2 insertions(+)

(limited to 'src/gemma_io.cpp')

diff --git a/src/gemma_io.cpp b/src/gemma_io.cpp
index 818c5e8..fad3e47 100644
--- a/src/gemma_io.cpp
+++ b/src/gemma_io.cpp
@@ -1483,12 +1483,14 @@ bool BimbamKin(const string file_geno, const set<string> ksnps,
   // and transpose
   // FIXME: the following is very slow
 
+  debug_msg("begin transpose");
   for (size_t i = 0; i < ni_total; ++i) {
     for (size_t j = 0; j < i; ++j) {
       d = gsl_matrix_get(matrix_kin, j, i);
       gsl_matrix_set(matrix_kin, i, j, d);
     }
   }
+  debug_msg("end transpose");
   // GSL is faster - and there are even faster methods
   // enforce_gsl(gsl_matrix_transpose(matrix_kin));
 
-- 
cgit v1.2.3