def METHOD_NAME(self, samples):
    """
    Calculates :meth:`per mode approximate KL divergence <kl_divergence>`
    then sums these.
    See :meth:`pints.toy.ToyLogPDF.distance()`.
    """
    kl = self.kl_divergence(samples)
    return np.sum(kl)