Simulation of stationary Gaussian time series using the Davies-Harte algorithm