Auto-diff through transfer function interpolation and integration

[fjebaker]

This means reworking the caches to use dual caches and fetching the right thing depending on whether either one of the interpolation arguments or the grid itself is of a dual number type.