Mercurial > repos > public > sbplib_julia

diff diffOp.jl @ 70:e4fa13137d12 cell_based_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Use double view
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 18 Jan 2019 13:51:32 +0100
parents 543b7a5ab831
children fbf7398f8154 45dece5e4928
line wrap: on
line diff
--- a/diffOp.jl	Thu Jan 17 16:04:07 2019 +0100
+++ b/diffOp.jl	Fri Jan 18 13:51:32 2019 +0100
@@ -69,9 +69,11 @@
     I = ci[i]
 
     # 2nd x-derivative
-    uᵢ  = apply(L.op, h[1], view(v, li[:,I[2]]), I[1])
+    vx = view(v, view(li,:,I[2]))
+    uᵢ  = apply(L.op, h[1], vx , I[1])
     # 2nd y-derivative
-    uᵢ += apply(L.op, h[2], view(v, li[I[1],:]), I[2])
+    vy = view(v, view(li,I[1],:))
+    uᵢ += apply(L.op, h[2], vy, I[2])
 
     return uᵢ
 end