Mercurial > repos > public > sbplib_julia

comparison diffOp.jl @ 61:8a7a537f54e5 default_perf_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:44:06 +0100
parents c62ea0112d4d
children
comparison
equal deleted inserted replaced
59:60d575e4a2d8 61:8a7a537f54e5
64 n_x, n_y = L.grid.numberOfPointsPerDim 64 n_x, n_y = L.grid.numberOfPointsPerDim
65 65
66 66
67 # For each x 67 # For each x
68 temp = zeros(eltype(u), n_y) 68 temp = zeros(eltype(u), n_y)
69 for i ∈ 1:n_x 69 for i ∈ 5:n_x-4
70 70
71 v_i = view(v, li[i,:]) 71 v_i = view(v, li[i,:])
72 apply!(L.op, temp, v_i, h[2]) 72 apply!(L.op, temp, v_i, h[2])
73 73
74 u[li[i,:]] += temp 74 u[li[i,:]] += temp
75 end 75 end
76 76
77 # For each y 77 # For each y
78 temp = zeros(eltype(u), n_x) 78 temp = zeros(eltype(u), n_x)
79 for i ∈ 1:n_y 79 for i ∈ 5:n_y-4
80 v_i = view(v, li[:,i]) 80 v_i = view(v, li[:,i])
81 apply!(L.op, temp, v_i, h[1]) 81 apply!(L.op, temp, v_i, h[1])
82 82
83 u[li[:,i]] += temp 83 u[li[:,i]] += temp
84 end 84 end