Optimising iterative solution of small system

[gdalle]

@alemelis

However this is not what I see when profiling the simulation as most of the time is spent in NR! when calling \ in

Profiling only runs the code once so it’s possible that the samples were noisy. Have you tried profiling a repeating loop instead?