From fe713dc5a75951bfd114d0be272180c1a081f95a Mon Sep 17 00:00:00 2001 From: Elias Carvalho Date: Tue, 5 Dec 2023 09:33:04 -0300 Subject: [PATCH] Update to Meshes.jl v0.37 --- Project.toml | 2 +- src/field/gaussian.jl | 4 ++-- src/field/gaussian/fft.jl | 5 +++-- 3 files changed, 6 insertions(+), 5 deletions(-) diff --git a/Project.toml b/Project.toml index b5591ce..ae0c63c 100644 --- a/Project.toml +++ b/Project.toml @@ -43,7 +43,7 @@ GeoStatsModels = "0.2" GeoTables = "1.9" ImageQuilting = "0.23" LinearAlgebra = "1.9" -Meshes = "0.36" +Meshes = "0.37" ProgressMeter = "1.9" Random = "1.9" Statistics = "1.9" diff --git a/src/field/gaussian.jl b/src/field/gaussian.jl index e0c1b11..dbebec7 100644 --- a/src/field/gaussian.jl +++ b/src/field/gaussian.jl @@ -25,8 +25,8 @@ include("gaussian/seq.jl") function defaultmethod(::GaussianProcess, setup::RandSetup) dom = setup.domain - udom, _ = unview(dom) - if udom isa Grid + pdom = parent(dom) + if pdom isa Grid FFTMethod() elseif nelements(dom) < 10000 LUMethod() diff --git a/src/field/gaussian/fft.jl b/src/field/gaussian/fft.jl index 62212d4..328aa72 100644 --- a/src/field/gaussian/fft.jl +++ b/src/field/gaussian/fft.jl @@ -55,7 +55,7 @@ function randprep(::AbstractRNG, process::GaussianProcess, method::FFTMethod, se dom = setup.domain data = setup.geotable - grid, _ = unview(dom) + grid = parent(dom) dims = size(grid) center = CartesianIndex(dims .รท 2) cindex = LinearIndices(dims)[center] @@ -104,7 +104,8 @@ function randsingle(rng::AbstractRNG, process::GaussianProcess, method::FFTMetho # retrieve domain info dom = setup.domain data = setup.geotable - grid, inds = unview(dom) + grid = parent(dom) + inds = parentindices(dom) dims = size(grid) # retrive variogram model and mean