DifferentialEquations.jl: LU factor ( M + dt * K ) once prior to time-stepping loop

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.