Computation of the multivariate Gaussian rate-distortion-perception function