Mercurial > repos > public > sbplib_julia

diff diffOp.jl @ 34:bb841977d198

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Move stencil operator application to its own function
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 10 Jan 2019 15:49:44 +0100
parents e8d7137b3f07
children ef060ab3b035
line wrap: on
line diff
--- a/diffOp.jl	Thu Jan 10 14:46:15 2019 +0100
+++ b/diffOp.jl	Thu Jan 10 15:49:44 2019 +0100
@@ -35,17 +35,6 @@
 
     h = scaling(L.grid)
 
-    for i ∈ 1:N
-        u[i] = apply(L.op.closureStencils[i], v, i)/h^2
-    end
-
-    for i ∈ N+1:M-N
-        u[i] = apply(L.op.innerStencil, i)/h^2
-    end
-
-    for i ∈ M:-1:M-N+1
-        u[i] = apply(flip(L.op.closureStencils[M-i+1]), v, i)/h^2
-    end
-
+    apply!(L.op, u, v, grid.spacings(L.grid)[1], 1, L.grid.numberOfPointsPerDim, stride=1)
     return nothing
 end