nvm, hard cap dist at 1 instead

marbl · Jul 30, 2018 · 512bbbe · 512bbbe
1 parent aa809ba
commit 512bbbe
Showing 1 changed file with 6 additions and 1 deletion.
diff --git a/src/mash/CommandDistance.cpp b/src/mash/CommandDistance.cpp
@@ -383,9 +383,14 @@ void compareSketches(CommandDistance::CompareOutput::PairOutput * output, const
     {
         //distance = log(double(common + 1) / (denom + 1)) / log(1. / (denom + 1));
         distance = -log(2 * jaccard / (1. + jaccard)) / kmerSize;
+
+        if ( distance > 1 )
+        {
+        	distance = 1;
+        }
     }
 
-    if ( distance > maxDistance && maxDistance < 1. )
+    if ( distance > maxDistance )
     {
         return;
     }