--- a/diffOp.jl	Thu Jan 24 14:33:49 2019 +0100
+++ b/diffOp.jl	Thu Jan 24 14:58:22 2019 +0100
@@ -48,7 +48,7 @@
 end
 
 # u = L*v
-function apply!(L::Laplace{2}, u::AbstractVector, v::AbstractVector)
+@inline function apply!(L::Laplace{2}, u::AbstractVector, v::AbstractVector)
     u .= 0*u
     h = Grid.spacings(L.grid)
 
@@ -60,19 +60,19 @@
     temp = zeros(eltype(u), n_y)
     for i ∈ 1:n_x
 
-        v_i = view(v, li[i,:])
-        apply!(L.op, temp, v_i, h[2])
+        @inbounds v_i = view(v, li[i,:])
+        @inbounds apply!(L.op, temp, v_i, h[2])
 
-        u[li[i,:]] += temp
+        @inbounds u[li[i,:]] += temp
     end
 
     # For each y
     temp = zeros(eltype(u), n_x)
     for i ∈ 1:n_y
-        v_i = view(v, li[:,i])
-        apply!(L.op, temp, v_i, h[1])
+        @inbounds v_i = view(v, li[:,i])
+        @inbounds apply!(L.op, temp, v_i, h[1])
 
-        u[li[:,i]] += temp
+        @inbounds u[li[:,i]] += temp
     end
 
     u .= L.a*u