Malkovsky

import numpy as np
from scipy.stats import binom
NEG_INF = -np.inf
LOG2 = np.log(2.0)
def logsubexp(log_a: float, log_b: float) -> float:
    """
    Compute log(exp(log_a) - exp(log_b)) stably, assuming log_a >= log_b.
    """
    if np.isneginf(log_b):
        return log_a