aboutsummaryrefslogtreecommitdiff
path: root/src/varcov.cpp
diff options
context:
space:
mode:
authorPjotr Prins2017-10-05 09:18:59 +0000
committerPjotr Prins2017-10-05 09:19:09 +0000
commit530488454ae6a35b098b81b45581f88f02d6de0e (patch)
tree7ee31c78a678ea4fa364267a9de8564faf4b159e /src/varcov.cpp
parente90c2d25cdb77d41d6587188db9e61cb988f2c78 (diff)
downloadpangemma-530488454ae6a35b098b81b45581f88f02d6de0e.tar.gz
Adding debug statements on entering functions and added test for issue 58
Diffstat (limited to 'src/varcov.cpp')
-rw-r--r--src/varcov.cpp2
1 files changed, 2 insertions, 0 deletions
diff --git a/src/varcov.cpp b/src/varcov.cpp
index 0f87ba8..39c3523 100644
--- a/src/varcov.cpp
+++ b/src/varcov.cpp
@@ -247,6 +247,7 @@ void Calc_Cor(vector<vector<double>> &X_mat, vector<double> &cov_vec) {
// genotype file and calculate the covariance matrix for neighboring
// SNPs output values at 10000-SNP-interval.
void VARCOV::AnalyzeBimbam() {
+ debug_msg("entering");
igzstream infile(file_geno.c_str(), igzstream::in);
if (!infile) {
cout << "error reading genotype file:" << file_geno << endl;
@@ -345,6 +346,7 @@ void VARCOV::AnalyzeBimbam() {
}
void VARCOV::AnalyzePlink() {
+ debug_msg("entering");
string file_bed = file_bfile + ".bed";
ifstream infile(file_bed.c_str(), ios::binary);
if (!infile) {