Mercurial > repos > public > sbplib_julia

comparison diffOp.jl @ 83:b795ec7f9ca0 cell_based_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 24 Jan 2019 17:46:57 +0100
parents 45dece5e4928 fbf7398f8154
children 8d505e9bc715
comparison
equal deleted inserted replaced
82:45dece5e4928 83:b795ec7f9ca0
69 li = LinearIndices(L.grid.numberOfPointsPerDim) 69 li = LinearIndices(L.grid.numberOfPointsPerDim)
70 ci = CartesianIndices(L.grid.numberOfPointsPerDim) 70 ci = CartesianIndices(L.grid.numberOfPointsPerDim)
71 I = ci[i] 71 I = ci[i]
72 72
73 # 2nd x-derivative 73 # 2nd x-derivative
74 vx = uview(v, uview(li,:,I[2])) 74 @inbounds vx = uview(v, uview(li,:,I[2]))
75 uᵢ = apply(L.op, h[1], vx , I[1]) 75 uᵢ = apply(L.op, h[1], vx , I[1])
76 # 2nd y-derivative 76 # 2nd y-derivative
77 vy = uview(v, uview(li,I[1],:)) 77 @inbounds vy = uview(v, uview(li,I[1],:))
78 uᵢ += apply(L.op, h[2], vy, I[2]) 78 uᵢ += apply(L.op, h[2], vy, I[2])
79 79
80 return uᵢ 80 return uᵢ
81 end 81 end