diff options
author | Pjotr Prins | 2017-10-13 12:28:26 +0000 |
---|---|---|
committer | Pjotr Prins | 2017-10-13 15:27:24 +0000 |
commit | a610dd723a233aed1abe31aa32e3137b23b5f983 (patch) | |
tree | 66172cf10d0d85f5a61219c1b7b074cbadb45575 /src/lmm.cpp | |
parent | 25ad9756ebecfdb2f01b6f87c35bd731e3a3186d (diff) | |
download | pangemma-a610dd723a233aed1abe31aa32e3137b23b5f983.tar.gz |
OpenBlas: preparing for dgemm use
Diffstat (limited to 'src/lmm.cpp')
-rw-r--r-- | src/lmm.cpp | 5 |
1 files changed, 4 insertions, 1 deletions
diff --git a/src/lmm.cpp b/src/lmm.cpp index a49c8c5..71aa184 100644 --- a/src/lmm.cpp +++ b/src/lmm.cpp @@ -39,6 +39,7 @@ #include "gsl/gsl_vector.h" #include "eigenlib.h" + #include "gzstream.h" #include "io.h" #include "lapack.h" @@ -1269,6 +1270,7 @@ void LMM::AnalyzeGene(const gsl_matrix *U, const gsl_vector *eval, return; } + void LMM::Analyze(std::function< SnpNameValues(size_t) >& fetch_snp, const gsl_matrix *U, const gsl_vector *eval, const gsl_matrix *UtW, const gsl_vector *Uty, @@ -1358,7 +1360,7 @@ void LMM::Analyze(std::function< SnpNameValues(size_t) >& fetch_snp, }; const auto num_snps = indicator_snp.size(); - const size_t progress_step = (num_snps/20>d_pace ? num_snps/20 : d_pace); + const size_t progress_step = (num_snps/50>d_pace ? num_snps/50 : d_pace); for (size_t t = 0; t < num_snps; ++t) { if (t % progress_step == 0 || t == (num_snps - 1)) { @@ -1429,6 +1431,7 @@ void LMM::Analyze(std::function< SnpNameValues(size_t) >& fetch_snp, batch_compute(msize); } batch_compute(c % msize); + ProgressBar("Reading SNPs", num_snps - 1, num_snps - 1); // cout << "Counted SNPs " << c << " sumStat " << sumStat.size() << endl; cout << endl; |