All of what Steven said, and then the next big performance thing would be to switch to DataInterpolations (preferably something with a continuous derivative like a cubic spline). DataInterpolations has some caching features so that repeated “nearby” interpolations should be about an order of magnitude faster, which is a special case that is common when used with ODE integration that the library is prepared for.