Mercurial > repos > public > sbplib_julia

changeset 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 d485da6e3a77
children 18d0d794d3bb fbf7398f8154
files diffOp.jl
diffstat 1 files changed, 4 insertions(+), 2 deletions(-) [+]
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