--- a/src/ImFilter.jl Wed Dec 25 17:42:32 2019 +0200 +++ b/src/ImFilter.jl Fri Dec 27 20:59:44 2019 +0200 @@ -43,12 +43,12 @@ tmp = 0.0 it₁ = inside(i, a₁, b₁, 1, n) it₂ = inside(j, a₂, b₂, 1, m) - for p=it₁ + @inbounds for p=it₁ @simd for q=it₂ - tmp += @inbounds(kp[p-o₁, q-o₂]*b[i+p,j+q]) + tmp += kp[p-o₁, q-o₂]*b[i+p,j+q] end end - res[i, j] = tmp + @inbounds res[i, j] = tmp end end