Mercurial > repos > public > sbplib_julia

comparison diffOp.jl @ 60:7ab0043aca53 cell_based_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Only do inner stencil
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 17 Jan 2019 09:41:25 +0100
parents 178a203f3e6d
children 17b971a0b852
comparison
equal deleted inserted replaced
57:178a203f3e6d 60:7ab0043aca53
32 error("not implemented") 32 error("not implemented")
33 end 33 end
34 34
35 # DiffOp must have a grid!!! 35 # DiffOp must have a grid!!!
36 function apply!(D::DiffOp, u::AbstractVector, v::AbstractVector) 36 function apply!(D::DiffOp, u::AbstractVector, v::AbstractVector)
37 for i ∈ 1:Grid.numberOfPoints(D.grid) 37 li = LinearIndices(D.grid.numberOfPointsPerDim)
38
39 Nx, Ny = D.grid.numberOfPointsPerDim
40
41 for i ∈ view(li, 5:Nx-4, 5:Ny-4)
38 u[i] = apply(D, v, i) 42 u[i] = apply(D, v, i)
39 end 43 end
40 44
41 return nothing 45 return nothing
42 end 46 end