I noticed that the default ODE tuning parameters for `pmx_solve_bdf`

are 1e-6, 1e-6, 1e+6; while the defaults for `integrate_ode_bdf`

is 1e-10, 1e-10, 1e+8. Is the difference deliberate?

+1 for having the algebraic solver tuning parameters for the steady states.