amrfinder.cpp, hmr.cpp and pmd.cpp: summary output now uses 64-bit fo…

…r total sizes and is unsigned. This should help with MongoDB integer values accumulating without wrapping and going negative only to be too large as unsigned
smithlabcode · Nov 30, 2023 · ee91acf · ee91acf
1 parent 108831b
commit ee91acf
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 3 deletions.
diff --git a/src/amrfinder/amrfinder.cpp b/src/amrfinder/amrfinder.cpp
@@ -58,7 +58,7 @@ using bamxx::bgzf_file;
 struct amr_summary {
   amr_summary(const vector<GenomicRegion> &amrs) {
     amr_count = size(amrs);
-    amr_total_size = accumulate(cbegin(amrs), cend(amrs), 0u,
+    amr_total_size = accumulate(cbegin(amrs), cend(amrs), 0ul,
                                 [](const uint64_t t, const GenomicRegion &p) {
                                   return t + p.get_width();
                                 });

diff --git a/src/analysis/hmr.cpp b/src/analysis/hmr.cpp
@@ -55,7 +55,7 @@ using bamxx::bgzf_file;
 struct hmr_summary {
   hmr_summary(const vector<GenomicRegion> &hmrs) {
     hmr_count = hmrs.size();
-    hmr_total_size = accumulate(cbegin(hmrs), cend(hmrs), 0,
+    hmr_total_size = accumulate(cbegin(hmrs), cend(hmrs), 0ul,
                                 [](const uint64_t t, const GenomicRegion &p) {
                                   return t + p.get_width(); });
     hmr_mean_size =

diff --git a/src/analysis/pmd.cpp b/src/analysis/pmd.cpp
@@ -60,7 +60,7 @@ template<typename T> using num_lim = std::numeric_limits<T>;
 struct pmd_summary {
   pmd_summary(const vector<GenomicRegion> &pmds) {
     pmd_count = pmds.size();
-    pmd_total_size = accumulate(cbegin(pmds), cend(pmds), 0,
+    pmd_total_size = accumulate(cbegin(pmds), cend(pmds), 0ul,
                                 [](const uint64_t t, const GenomicRegion &p) {
                                   return t + p.get_width(); });
     pmd_mean_size =