--- 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