Compute the base matrix of a continuous Gaussian field. Computed as a block circulant matrix on a torus where x and y is the x and y centroids (must be equally spaced)
blockcircbase(x, y, sigma, phi, model, additionalparameters, inverse = FALSE)x centroids, an equally spaced vector
y centroids, an equally spaced vector
spatial variance parameter
spatial decay parameter
covariance model, see ?CovarianceFct
additional parameters for chosen covariance model. See ?CovarianceFct
logical. Whether to return the base matrix of the inverse covariance matrix (ie the base matrix for the precision matrix), default is FALSE
the base matrix of a block circulant matrix representing a stationary covariance function on a toral grid.