Mercurial > repos > public > sbplib_julia

changeset 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 60d575e4a2d8
children 5b09c16d669d
files diffOp.jl sbpD2.jl
diffstat 2 files changed, 2 insertions(+), 10 deletions(-) [+]
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])
 
--- a/sbpD2.jl	Wed Jan 16 12:40:20 2019 +0100
+++ b/sbpD2.jl	Thu Jan 17 09:44:06 2019 +0100
@@ -4,20 +4,12 @@
     N = length(v)
     cSize = closureSize(op)
 
-    for i ∈ range(1; length=cSize)
-        u[i] = apply(op.closureStencils[i], v, i)/h^2
-    end
-
     innerStart = 1 + cSize
     innerEnd = N - cSize
     for i ∈ range(innerStart, stop=innerEnd)
         u[i] = apply(op.innerStencil, v, i)/h^2
     end
 
-    for i ∈ range(innerEnd+1, length=cSize)
-        u[i] = Int(op.parity)*apply(flip(op.closureStencils[N-i+1]), v, i)/h^2
-    end
-
     return nothing
 end