Add comment

Author: szhan

python/tests/beagle_numba.py

@@ -562,6 +562,7 @@ def interpolate_allele_probs(
         for a in alleles:
             is_a_in_ref_h = ref_h[i, :] == a
             if np.sum(is_a_in_ref_h) == 0:
+                # This avoids division by zero when getting a threshold adaptively below.
                 continue
             if use_threshold:
                 # TODO: Check whether this is implemented correctly. Not used by default.