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