The example on the sunode readme page uses the sundials solvers through Aesara in a PyMC model.
It doesn’t use pm.ode.DifferentialEquation if that’s your question.
Depending on your model, sunode or jax.experimental.odeint will be the fastest, but for a fair comparison you must run them with the same tolerance settings!!