I wonder whether this shouldn't be rolled together with entropy, with a symmetric keyword for example
https://github.com/scipy/scipy/blob/c254b4ff5def78c2a5e2d9decbd35882a2f0f650/scipy/stats/_distn_infrastructure.py#L2292

Otherwise we should get the same broadcasting behavior, i.e. makes sum wrt axis=0