Computing ODE sensitivities with sparse Jacobians in solver

You’ll need to switch to SparsepakFactorization(), which it should do automatically?