Updating an array inside a for loop- comparison with MATLAB

[ChrisRackauckas]

@Stephen_Vavasis

Explicit loops reduce the readability of the code, but this can be partly hidden using macro packages like devectorize or einsum.

In my finite element code, I use explicit loops and scalar ops for all dense matrix-vector operations on smaller matrices and vectors. I use matrix/vector ops for operations on big matrices like the global stiffness matrix. I’m reasonably happy with the performance.

You should checkout InplaceOps.jl for making this easily in-place.