Optimizing Large-Scale ODE Simulations