Debugged for a case with small number of cpgs

smithlabcode · Dec 19, 2023 · b38eecd · b38eecd
1 parent a5c618d
commit b38eecd
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/src/amrfinder/amrfinder.cpp b/src/amrfinder/amrfinder.cpp
@@ -330,7 +330,9 @@ process_chrom(const bool verbose, const uint32_t n_threads,
 
   const auto n_blocks = n_threads*blocks_per_thread;
 
-  const uint64_t lim = n_cpgs - window_size + 1;
+  const uint64_t lim = (n_cpgs >= window_size) ? 
+                       n_cpgs - window_size + 1 : 0;
+                                                    ;
   const auto blocks = get_block_bounds(static_cast<uint64_t>(0),
                                        lim, lim/n_blocks);