Can once trick DifferentialEquations.jl to LU factor M + dt * K once only prior to the time stepping loop? Thx.
It already reuses the jacobian (and factorization) when possible.
When you solved, if you check sol.stats, how many times did it factorize?
Many thanks for referring to sol.stats !
I am looking further into this further.