aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorPjotr Prins2022-08-05 13:17:06 +0200
committerPjotr Prins2022-08-05 13:17:06 +0200
commit9ea4215bd1c2dbd616049538180a67d5fb65766c (patch)
tree9fde4bd824e66e84beb6e91faa23644bcc6e736a
parent1c02316f9cc44b092d72cd361964fcaac04be357 (diff)
downloadpangemma-9ea4215bd1c2dbd616049538180a67d5fb65766c.tar.gz
Bump version and minor output
-rw-r--r--src/gemma_io.cpp11
-rw-r--r--src/version.h6
2 files changed, 11 insertions, 6 deletions
diff --git a/src/gemma_io.cpp b/src/gemma_io.cpp
index 45d1eb4..962e0c0 100644
--- a/src/gemma_io.cpp
+++ b/src/gemma_io.cpp
@@ -150,7 +150,7 @@ std::istream &safeGetline(std::istream &is, std::string &t) {
// Read SNP file. A single column of SNP names.
bool ReadFile_snps(const string file_snps, set<string> &setSnps) {
- debug_msg("entered");
+ debug_msg("enter ReadFile_snps");
setSnps.clear();
igzstream infile(file_snps.c_str(), igzstream::in);
@@ -329,6 +329,9 @@ bool ReadFile_anno(const string &file_anno, map<string, string> &mapRS2chr,
mapRS2bp[rs] = b_pos;
mapRS2cM[rs] = cM;
}
+ // for (auto& [key, value] : mapRS2bp) {
+ // cerr << key << endl;
+ //}
infile.close();
infile.clear();
@@ -693,6 +696,9 @@ bool ReadFile_geno(const string &file_geno, const set<string> &setSnps,
file_pos = 0;
auto count_warnings = 0;
auto infilen = file_geno.c_str();
+ // for (auto& [key, value] : mapRS2bp) {
+ // cerr << key << endl;
+ // }
while (!safe_get_line(infile, line).eof()) {
ch_ptr = strtok_safe2((char *)line.c_str(), " ,\t",infilen);
rs = ch_ptr;
@@ -714,8 +720,7 @@ bool ReadFile_geno(const string &file_geno, const set<string> &setSnps,
if (mapRS2bp.count(rs) == 0) {
if (is_debug_mode() && count_warnings++ < 10) {
- std::string msg = "Can't figure out position for ";
- msg += rs;
+ std::string msg = "Can't figure out position for <" + rs + ">";
debug_msg(msg);
if (count_warnings == 10)
debug_msg("Skipping similar warnings");
diff --git a/src/version.h b/src/version.h
index 53ab379..25ccfcd 100644
--- a/src/version.h
+++ b/src/version.h
@@ -1,5 +1,5 @@
// version.h generated by GEMMA scripts/gen_version_info.sh
-#define GEMMA_VERSION "0.98.5"
-#define GEMMA_DATE "2021-08-25"
-#define GEMMA_YEAR "2021"
+#define GEMMA_VERSION "0.98.6"
+#define GEMMA_DATE "2022-08-05"
+#define GEMMA_YEAR "2022"
#define GEMMA_PROFILE "/gnu/store/8rvid272yb53bgascf5c468z0jhsyflj-profile"