as per title, see https://docs.rs/diffsol/latest/diffsol/nonlinear_solver/trait.NonLinearSolver.html#tymethod.solve_in_place
This would allow the solvers to directly do the solve in each column of the diff matrix, which would mean that the old_state
does not need to be stored, a significant storage reduction for the bdf and sdirk solvers