We develop an algorithm to simulate a Gaussian stochastic process that is non-δ-correlated in both space and time coordinates. The colored noise obeys a linear reaction-diffusion Langevin equation with Gaussian white noise. This equation is exactly simulated in a discrete Fourier space.

Peer Reviewed