src/common.hpp: adding the report_histogram to the common sources

Author: andrewdavidsmith

src/common.cpp

@@ -334,7 +334,7 @@ write_predicted_complexity_curve(const string &outfile, const double c_level,
                                  const vector<double> &yield_upper_ci_lognorm) {
   std::ofstream of;
   if (!outfile.empty())
-    of.open(outfile.c_str());
+    of.open(outfile);
   std::ostream out(outfile.empty() ? std::cout.rdbuf() : of.rdbuf());
 
   out << "TOTAL_READS\tEXPECTED_DISTINCT\t"

src/common.hpp

@@ -23,6 +23,9 @@
 
 #include <cstddef>  // std::size_t
 #include <cstdint>  // std::uint64_t
+#include <fstream>
+#include <iostream>
+#include <ostream>
 #include <random>
 #include <stdexcept>
 #include <string>
@@ -113,4 +116,18 @@ multinomial(std::mt19937 &gen, const std::vector<double> &mult_probs,
     throw std::runtime_error("multinomial sampling failed");
 }
 
+template <typename H>
+void
+report_histogram(const std::string &outfile, const H &h) {
+  std::ofstream of;
+  if (!outfile.empty())
+    of.open(outfile);
+  std::ostream o(outfile.empty() ? std::cerr.rdbuf() : of.rdbuf());
+  o << "OBSERVED COUNTS (" << std::size(h) << ")" << std::endl;
+  for (auto i = 0u; i < std::size(h); ++i)
+    if (h[i] > 0)
+      o << i << '\t' << static_cast<std::uint32_t>(h[i]) << std::endl;
+  o << std::endl;
+}
+
 #endif  // SRC_COMMON_HPP_