Stochastic variance-reduced Gaussian variational inference on the Bures-Wasserstein manifold