Instead of constructing a sparse matrix in the inner loop it would be more 
efficient to write an in place stencil kernel function to perform the 
equivalent operation.

Reply via email to