A black graph with white text. Y-axis is Average time; X-axis is number of required integration steps. There are several dots showing the performance of CyRK's cysolve_ivp, pysolve_ivp, and nbsolve_ivp functions compared to SciPy's solve_ivp. In all cases CyRK is significantly faster than SciPy. The plot shows the cysolve_ivp method to be the best in all circumstances.
I wanted to share a side project I have been working on: CyRK github.com/jrenaud90/CyRK a #Python package that mimics much of the functionality of #SciPy solve_ivp but can perform #ODE integrations 100's x faster. If you are looking for a numerical integrator for your python project, check it out!