If k-mer length not given explicitly, calculate from pigeonhole lemma

seqan · May 9, 2023 · 5e1c59c · 5e1c59c
1 parent 3f75081
commit 5e1c59c
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/include/stellar/options/index_options.hpp b/include/stellar/options/index_options.hpp
@@ -8,7 +8,7 @@ namespace stellar
 
 struct IndexOptions
 {
-    size_t qGram{7}; // length of the q-grams
+    size_t qGram{std::numeric_limits<size_t>::max()}; // length of the q-grams
     double qgramAbundanceCut{1};
 };
 

diff --git a/include/stellar/stellar_types.hpp b/include/stellar/stellar_types.hpp
@@ -72,7 +72,7 @@ struct StellarStatistics
         assert(options.epsilon.numerator() == 0 || StellarOptions::absoluteErrors(options.epsilon, n1) == e1);
 
         kMerLength = options.qGram;
-        kMerComputed = options.qGram == (unsigned)-1;
+        kMerComputed = options.qGram == (size_t)-1;
 
         if (kMerComputed)
             kMerLength = std::min(std::max(1u, smin), 32u);