gibbs builtin

gibbs(
    coords,
    bounds,
    variogram,
    *,
    iterations=200,
    burn_in=50,
    seed=0,
)

One Gaussian draw at coords honoring bounds ((n, 2) lower/upper).