Mercurial > repos > public > sbplib_julia

changeset 82:45dece5e4928 cell_based_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Use unsafe views
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 24 Jan 2019 17:42:51 +0100
parents 4640839b1616
children b795ec7f9ca0
files diffOp.jl
diffstat 1 files changed, 4 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/diffOp.jl	Fri Jan 18 14:26:04 2019 +0100
+++ b/diffOp.jl	Thu Jan 24 17:42:51 2019 +0100
@@ -60,6 +60,8 @@
     return uᵢ
 end
 
+using UnsafeArrays
+
 # u = L*v
 function apply(L::Laplace{2}, v::AbstractVector, i::Int)
     h = Grid.spacings(L.grid)
@@ -69,10 +71,10 @@
     I = ci[i]
 
     # 2nd x-derivative
-    vx = view(v, view(li,:,I[2]))
+    vx = uview(v, uview(li,:,I[2]))
     uᵢ  = apply(L.op, h[1], vx , I[1])
     # 2nd y-derivative
-    vy = view(v, view(li,I[1],:))
+    vy = uview(v, uview(li,I[1],:))
     uᵢ += apply(L.op, h[2], vy, I[2])
 
     return uᵢ