Direct Gaussian simulation solver for the GeoStats.jl framework
Author JuliaEarth
0 Stars
Updated Last
3 Months Ago
Started In
March 2018


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.

Used By Packages

No packages found.