Squeezing performance out of ODE solves with interpolations

To reduce the number of f calls, have you specified the Jacobian sparsities?