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.
1 Like
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.