aboutsummaryrefslogtreecommitdiff
path: root/src/gemma_io.cpp
diff options
context:
space:
mode:
authorPjotr Prins2020-10-01 09:57:34 +0100
committerPjotr Prins2020-11-29 08:50:44 +0000
commitb298c466171006dfb1196c74d7017d0dc8e45492 (patch)
tree5c532828ffe437536ec1b8ba7c916efe2c17ceb5 /src/gemma_io.cpp
parent1a431e1319e2ff3855f8617f446a33a3f931bad6 (diff)
downloadpangemma-b298c466171006dfb1196c74d7017d0dc8e45492.tar.gz
These variables are introduced in DEBUG mode
Diffstat (limited to 'src/gemma_io.cpp')
-rw-r--r--src/gemma_io.cpp2
1 files changed, 2 insertions, 0 deletions
diff --git a/src/gemma_io.cpp b/src/gemma_io.cpp
index 569c79b..81182f8 100644
--- a/src/gemma_io.cpp
+++ b/src/gemma_io.cpp
@@ -1557,6 +1557,7 @@ bool BimbamKin(const string file_geno, const set<string> ksnps,
// FIXME: the following is not so slow but appears to generate an
// identical matrix
+ /*
for (size_t i = 0; i < ni_total; ++i) {
for (size_t j = 0; j < i; ++j) {
double d = gsl_matrix_get(matrix_kin, j, i);
@@ -1564,6 +1565,7 @@ bool BimbamKin(const string file_geno, const set<string> ksnps,
}
}
write(matrix_kin,"K rotated");
+ */
// GSL is faster - and there are even faster methods
// enforce_gsl(gsl_matrix_transpose(matrix_kin));