Pymc inference for a high dimensional ODE system

Ahh you have an odeint there, what library is that?

There’s sunode for PyMC models that doesn’t need any as_op for sure