Mercurial > repos > public > sbplib_julia

comparison diffOp.jl @ 34:bb841977d198

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Move stencil operator application to its own function
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 10 Jan 2019 15:49:44 +0100
parents e8d7137b3f07
children ef060ab3b035
comparison
equal deleted inserted replaced
33:e8d7137b3f07 34:bb841977d198
33 N = closureSize(L.op) 33 N = closureSize(L.op)
34 M = length(v) 34 M = length(v)
35 35
36 h = scaling(L.grid) 36 h = scaling(L.grid)
37 37
38 for i ∈ 1:N 38 apply!(L.op, u, v, grid.spacings(L.grid)[1], 1, L.grid.numberOfPointsPerDim, stride=1)
39 u[i] = apply(L.op.closureStencils[i], v, i)/h^2
40 end
41
42 for i ∈ N+1:M-N
43 u[i] = apply(L.op.innerStencil, i)/h^2
44 end
45
46 for i ∈ M:-1:M-N+1
47 u[i] = apply(flip(L.op.closureStencils[M-i+1]), v, i)/h^2
48 end
49
50 return nothing 39 return nothing
51 end 40 end