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.