Diffrax vs solve_ivp: what actually got faster
A cosmologist moved his likelihood from SciPy to Diffrax and got roughly seven times faster solves plus exact gradients. I go through what the numbers really show, the three silent porting traps, and the cases where I would keep solve_ivp.