Compute autocovariance or autocorrelation function via Wiener-Khinchin theorem using Fast Fourier Transform
ac_fft(x, demean = TRUE)
matrix with time (iteration number) along the rows and variables along the columns.
whether to subtract from each column its mean.
A matrix of the same size as x with autocovariances at all lags from 1 to the number of rows.