Mercurial > repos > public > sbplib_julia
diff sbpD2.jl @ 60:7ab0043aca53 cell_based_test
Only do inner stencil
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Thu, 17 Jan 2019 09:41:25 +0100 |
| parents | 27a8d3021a1c |
| children |
line wrap: on
line diff
--- a/sbpD2.jl Tue Jan 15 19:51:16 2019 +0100 +++ b/sbpD2.jl Thu Jan 17 09:41:25 2019 +0100 @@ -4,13 +4,7 @@ cSize = closureSize(op) N = length(v) - if i ∈ range(1; length=cSize) - uᵢ = apply(op.closureStencils[i], v, i)/h^2 - elseif i ∈ range(N - cSize+1, length=cSize) - uᵢ = Int(op.parity)*apply(flip(op.closureStencils[N-i+1]), v, i)/h^2 - else - uᵢ = apply(op.innerStencil, v, i)/h^2 - end + uᵢ = apply(op.innerStencil, v, i)/h^2 return uᵢ end