Mercurial > repos > public > sbplib_julia

diff 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
line wrap: on
line diff
--- a/diffOp.jl	Wed Jan 16 12:40:20 2019 +0100
+++ b/diffOp.jl	Thu Jan 17 09:44:06 2019 +0100
@@ -66,7 +66,7 @@
 
     # For each x
     temp = zeros(eltype(u), n_y)
-    for i ∈ 1:n_x
+    for i ∈ 5:n_x-4
 
         v_i = view(v, li[i,:])
         apply!(L.op, temp, v_i, h[2])
@@ -76,7 +76,7 @@
 
     # For each y
     temp = zeros(eltype(u), n_x)
-    for i ∈ 1:n_y
+    for i ∈ 5:n_y-4
         v_i = view(v, li[:,i])
         apply!(L.op, temp, v_i, h[1])