Here’s a situation where we replaced an integrate_1d call with the ode integrator: Difference in behavior between integrate_1d and integrate
Not great but it worked. Try the integrate_1d quadrature in your fancier case before you toss it though. It’s pretty problem dependent.