Fast and exact synthesis of stationary multivariate Gaussian time series using circulant embedding