Prevent division by zero in jaccard computation

Closes #17.
broadinstitute · Jan 11, 2023 · 67688d9 · 67688d9
1 parent 50fbed1
commit 67688d9
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/strainge/comparison.py b/src/strainge/comparison.py
@@ -36,6 +36,10 @@ def jaccard(kmers1, kmers2):
     """Computes jaccard similarity. Returns numerator and denominator
     separately."""
     intersection = kmerizer.count_common(kmers1, kmers2)
+
+    if kmers1.size == 0 and kmers2.size == 0:
+        return 0
+
     return intersection / (kmers1.size + kmers2.size - intersection)