faster diagonalgaussian likelihoods using einsum