Cubic Interpolation With Derivatives In Numpy

A recent immigrant to Python and scientific computing with Python. This is a question to avoid any duplication of code that might already exist. I have a field that is sampled as a

Solution 1:

The scipy.interpolate.interp2d.__call__ method has options dx and dy to evaluate higher derivatives at point (at least since version 0.14.0).