I would like to put a change into the calls to Intel MKL to speed up simulation. Specifically, this is to reuse the symbolic factorization from the first Newton iteration for the remaining iterations.

The advantages is a significant speed increase, since the symbolic factorization is no longer done for each nonlinear iteration. The regression tests still seem to run to completion.

However, there are slight differences in the results, as the symbolic factorization is based only at the simulation matrix from the first iteration. Note that the per iteration matrix values are still considered for the numerical factorization.

Should this change should be the default, or should it be enabled by the user on a per simulation basis?