From 0eb4326a922152640c7ba4618427c4ef77a7d975 Mon Sep 17 00:00:00 2001 From: Korbinian Eckstein Date: Sun, 24 Dec 2023 09:22:33 +1000 Subject: [PATCH] set default stencil grid to 64^3; reduces time from 15s to 1s, might have accuracy implications --- src/oblique_stencil.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/oblique_stencil.jl b/src/oblique_stencil.jl index 86faa9d..9d4fc5b 100644 --- a/src/oblique_stencil.jl +++ b/src/oblique_stencil.jl @@ -10,7 +10,7 @@ function dipole(x, y, z, r_treshold, direction=(0, 0, 1)) return kappa end -function stencil(; st=27, res=(1.0, 1.0, 1.0), res0=1.0, direction=(0, 0, 1), gridsize=(128, 128, 128)) +function stencil(; st=27, res=(1.0, 1.0, 1.0), res0=1.0, direction=(0, 0, 1), gridsize=(64, 64, 64)) middle = floor.(Int, gridsize ./ 2) .+ 1 coord = [((1:gridsize[i]) .- middle[i]) * res0 for i in 1:3]