This package provides an implementation of direct Gaussian simulation (a.k.a. LU simulation) as described in Alabert 1987. In this method, the full covariance matrix is built to include all locations of the simulation domain, and samples from the multivariate Gaussian are drawn via LU factorization.

The method, which is widely implemented in many packages for Gaussian processes (e.g. GaussianProcesses.jl, GaussianRandomFields.jl), is appropriate for relatively small simulation domains (e.g. 100x100 grids) where it is feasible to factorize the full covariance. For larger domains (e.g. 3D grids), other methods are available such as sequential Gaussian simulation, spectral methods, and FFT moving averages.


Get the latest stable release with Julia's package manager:

] add DirectGaussianSimulation


This package is part of the GeoStats.jl framework.

For a simple example of usage, please check this notebook.

Asking for help

If you have any questions, please open an issue.