A method for simulating a stationary Gaussian process on a fine rectangular grid in [0, 1] 𝑑 ⊂ ℝ 𝑑 is described. It is assumed that the process is stationary with respect to translations of R d, but ...