diff options
| author | Pjotr Prins | 2017-10-05 09:18:59 +0000 |
|---|---|---|
| committer | Pjotr Prins | 2017-10-05 09:19:09 +0000 |
| commit | 530488454ae6a35b098b81b45581f88f02d6de0e (patch) | |
| tree | 7ee31c78a678ea4fa364267a9de8564faf4b159e /src/lmm.cpp | |
| parent | e90c2d25cdb77d41d6587188db9e61cb988f2c78 (diff) | |
| download | pangemma-530488454ae6a35b098b81b45581f88f02d6de0e.tar.gz | |
Adding debug statements on entering functions and added test for issue 58
Diffstat (limited to 'src/lmm.cpp')
| -rw-r--r-- | src/lmm.cpp | 5 |
1 files changed, 5 insertions, 0 deletions
diff --git a/src/lmm.cpp b/src/lmm.cpp index eb76265..37f2f5b 100644 --- a/src/lmm.cpp +++ b/src/lmm.cpp @@ -1319,6 +1319,7 @@ void LMM::AnalyzeBimbam(const gsl_matrix *U, const gsl_vector *eval, const gsl_matrix *UtW, const gsl_vector *Uty, const gsl_matrix *W, const gsl_vector *y, const set<string> gwasnps) { + debug_msg("entering"); clock_t time_start = clock(); // LOCO support @@ -1482,6 +1483,7 @@ void LMM::AnalyzeBimbam(const gsl_matrix *U, const gsl_vector *eval, void LMM::AnalyzePlink(const gsl_matrix *U, const gsl_vector *eval, const gsl_matrix *UtW, const gsl_vector *Uty, const gsl_matrix *W, const gsl_vector *y) { + debug_msg("entering"); string file_bed = file_bfile + ".bed"; ifstream infile(file_bed.c_str(), ios::binary); if (!infile) { @@ -1676,6 +1678,7 @@ void LMM::AnalyzePlink(const gsl_matrix *U, const gsl_vector *eval, void LMM::Analyzebgen(const gsl_matrix *U, const gsl_vector *eval, const gsl_matrix *UtW, const gsl_vector *Uty, const gsl_matrix *W, const gsl_vector *y) { + debug_msg("entering"); string file_bgen = file_oxford + ".bgen"; ifstream infile(file_bgen.c_str(), ios::binary); if (!infile) { @@ -2282,6 +2285,7 @@ void LMM::AnalyzeBimbamGXE(const gsl_matrix *U, const gsl_vector *eval, const gsl_matrix *UtW, const gsl_vector *Uty, const gsl_matrix *W, const gsl_vector *y, const gsl_vector *env) { + debug_msg("entering"); igzstream infile(file_geno.c_str(), igzstream::in); if (!infile) { cout << "error reading genotype file:" << file_geno << endl; @@ -2430,6 +2434,7 @@ void LMM::AnalyzePlinkGXE(const gsl_matrix *U, const gsl_vector *eval, const gsl_matrix *UtW, const gsl_vector *Uty, const gsl_matrix *W, const gsl_vector *y, const gsl_vector *env) { + debug_msg("entering"); string file_bed = file_bfile + ".bed"; ifstream infile(file_bed.c_str(), ios::binary); if (!infile) { |