Report post #13 in Rust: cannot push to Vec while I hold reference to element, but why
Same two pass thing applies in numpy, compute everything from the old array then assign. Learned that from a very wrong heat diffusion plot.
Programming discussion