diff --git a/src/oblique_stencil.jl b/src/oblique_stencil.jl index 076b41a..e4a10c4 100644 --- a/src/oblique_stencil.jl +++ b/src/oblique_stencil.jl @@ -10,12 +10,12 @@ 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, singularity_cutout=4, direction=(0, 0, 1), gridsize=(64, 64, 64)) +function stencil(; st=27, res=(1.0, 1.0, 1.0), singularity_cutout=4, 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] + coord = [((1:gridsize[i]) .- middle[i]) for i in 1:3] - d = [dipole(x, y, z, singularity_cutout * res0, direction) for x in coord[1], y in coord[2], z in coord[3]] + d = [dipole(x, y, z, singularity_cutout, direction) for x in coord[1], y in coord[2], z in coord[3]] d_mask = isfinite.(d) # stencil mask