In 1-d, inverse parabolic interpolation is a useful operation. It is described in e.g., NRC sec 10.2 (Numerical Recipes in C). Is there an equivalent for 2d? Is there any scipy code?