Merge branch 'master' of github.com:marbl/MashMap

marbl · Nov 30, 2023 · cc49b9f · cc49b9f
2 parents 1a405d2 + 9afc465
commit cc49b9f
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/src/map/include/computeMap.hpp b/src/map/include/computeMap.hpp
@@ -1638,6 +1638,12 @@ namespace skch
                   ) / it->n_merged; // this would scale directly by the number of mappings in the chain
 
 
+              //Mean kmer complexity of all mappings in the chain
+              it->kmerComplexity = ( std::accumulate(
+                                      it, it_end, 0.0,
+                                      [](double x, MappingResult &e){ return x + e.kmerComplexity; })
+                  ) / it->n_merged; // this would scale directly by the number of mappings in the chain
+
               //Discard other mappings of this chain
               std::for_each( std::next(it), it_end, [&](MappingResult &e){ e.discard = 1; });