changed binning rng to be dependent on the seed from preseq

Author: terencewtli

src/load_data_for_complexity.cpp

@@ -751,14 +751,14 @@ SplitMappedRead(const bool VERBOSE,
 size_t
 load_coverage_counts_MR(const bool VERBOSE,
                         const string input_file_name,
+			const unsigned long int seed,
                         const size_t bin_size,
                         const size_t max_width,
                         vector<double> &coverage_hist) {
 
   srand(time(0) + getpid());
   //Runif runif(rand());
-  std::random_device rd;
-  std::mt19937 generator(rd());
+  std::mt19937 generator(seed);
 
   std::ifstream in(input_file_name.c_str());
   if (!in)
@@ -813,14 +813,14 @@ load_coverage_counts_MR(const bool VERBOSE,
 
 size_t
 load_coverage_counts_GR(const string input_file_name,
+			const unsigned long int seed,
                         const size_t bin_size,
                         const size_t max_width,
                         vector<double> &coverage_hist) {
 
   srand(time(0) + getpid());
   //Runif runif(rand());
-  std::random_device rd;
-  std::mt19937 generator(rd());
+  std::mt19937 generator(seed);
 
   std::ifstream in(input_file_name.c_str());
   if (!in)

src/load_data_for_complexity.hpp

@@ -26,13 +26,15 @@
 size_t
 load_coverage_counts_MR(const bool VERBOSE,
                         const std::string input_file_name,
+			const unsigned long int seed,
                         const size_t bin_size,
                         const size_t max_width,
                         std::vector<double> &coverage_hist);
 
 
 size_t
 load_coverage_counts_GR(const std::string input_file_name,
+			const unsigned long int seed,
                         const size_t bin_size,
                         const size_t max_width,
                         std::vector<double> &coverage_hist);

src/preseq.cpp

@@ -883,13 +883,13 @@ gc_extrap(const int argc, const char **argv) {
     if (NO_SEQUENCE) {
       if (VERBOSE)
         cerr << "BED FORMAT" << endl;
-      n_reads = load_coverage_counts_GR(input_file_name, bin_size,
+      n_reads = load_coverage_counts_GR(input_file_name, seed, bin_size,
                                         max_width, coverage_hist);
     }
     else {
       if (VERBOSE)
         cerr << "MAPPED READ FORMAT" << endl;
-      n_reads = load_coverage_counts_MR(VERBOSE, input_file_name, bin_size,
+      n_reads = load_coverage_counts_MR(VERBOSE, input_file_name, seed, bin_size,
                                         max_width, coverage_hist);
     }