diff options
author | Pjotr Prins | 2020-10-01 09:57:34 +0100 |
---|---|---|
committer | Pjotr Prins | 2020-10-01 09:57:34 +0100 |
commit | 40a9d521f4ecc9d4762a5f8c1953ec5f321755f0 (patch) | |
tree | 2215dc71fdf994ea30c5ec00be3e91f65ae58228 /src/gemma_io.cpp | |
parent | 85272f8d045e9ecb0a72ad5b2738eb135708bb01 (diff) | |
download | pangemma-40a9d521f4ecc9d4762a5f8c1953ec5f321755f0.tar.gz |
These variables are introduced in DEBUG mode
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 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)); |