Hey there, I have a model that spends most time in functions like this function Gux!(dudx::Matrix{Numtype},u::Matrix{Numtype}) dudx[1:end-1,:] = u[2:end,:]-u[1:end-1,:] dudx[end,:] = u[1,:]-u[end,:]` end that will be evaluated over and over again with changing input matrices u. Prealloc…

Efficient finite difference operators

dpsanders July 17, 2018, 9:19pm 11

The explicit loop is much more readable IMO.

Strange performance of a loop