Skip to content

Commit

Permalink
set default stencil grid to 64^3;
Browse files Browse the repository at this point in the history
reduces time from 15s to 1s, might have accuracy implications
  • Loading branch information
korbinian90 committed Dec 23, 2023
1 parent 6251318 commit 0eb4326
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion src/oblique_stencil.jl
Original file line number Diff line number Diff line change
Expand Up @@ -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]
Expand Down

0 comments on commit 0eb4326

Please sign in to comment.