@josef-pkt wrt to the broadcast behaviour: Unlike entropy, which makes sense for a single sample, this is only meaningful in a pairwise fashion. Thus, I think the interface like the ones in ``scipy.spatial.distance`` is more appropriate. If others disagree, then I can change the PR, but as a user, it'd probably take me a while to figure out what the interface is.