diff options
author | Pjotr Prins | 2021-08-16 10:56:02 +0200 |
---|---|---|
committer | Pjotr Prins | 2021-08-16 10:56:02 +0200 |
commit | 0dbe0670fe437ddadd713f8abc404a4be215ac7f (patch) | |
tree | 75d425d21560fa8580b2c4ec0331feca9e3a5351 /src | |
parent | 42de3667c90964206c3764be1867527d8eeb8cfc (diff) | |
download | pangemma-0dbe0670fe437ddadd713f8abc404a4be215ac7f.tar.gz |
README.md
README
Diffstat (limited to 'src')
-rw-r--r-- | src/mvlmm.cpp | 6 |
1 files changed, 6 insertions, 0 deletions
diff --git a/src/mvlmm.cpp b/src/mvlmm.cpp index 4b8db05..51ed901 100644 --- a/src/mvlmm.cpp +++ b/src/mvlmm.cpp @@ -3666,6 +3666,8 @@ void MVLMM::AnalyzePlink(const gsl_matrix *U, const gsl_vector *eval, t_last++; } for (vector<SNPINFO>::size_type t = 0; t < snpInfo.size(); ++t) { + cout << snpInfo[t].rs_number << endl; + if (t % d_pace == 0 || t == snpInfo.size() - 1) { ProgressBar("Reading SNPs", t, snpInfo.size() - 1); } @@ -3752,6 +3754,10 @@ void MVLMM::AnalyzePlink(const gsl_matrix *U, const gsl_vector *eval, gsl_matrix_set_zero(Xlarge); for (size_t i = 0; i < l; i++) { + double t = (double(clock()-time_start)/(double(CLOCKS_PER_SEC))); + printf("%.2f ",t); + if (i>10) exit(1); + cout << i << "," << l << endl; gsl_vector_view UtXlarge_col = gsl_matrix_column(UtXlarge, i); gsl_vector_memcpy(&X_row.vector, &UtXlarge_col.vector); |