Added the performance script I used.

Author: pelson

stratify/_vinterp.pyx

@@ -392,7 +392,6 @@ cdef class _PythonExtrapKernel(ExtrapKernel):
         with gil:
             self.extrap_kernel(direction, z_src, fz_src, level, fz_target)
 
-
 cdef class _TestableDirectionExtrapKernel(ExtrapKernel):
     @cython.boundscheck(False)
     @cython.wraparound(False)

stratify/tests/performance.py

@@ -0,0 +1,25 @@
+"""
+Functions that may be used to measure performance of a component.
+
+"""
+import numpy as np
+import stratify
+
+
+def src_data(shape=(400, 500, 100)):
+    z = np.tile(np.linspace(0, 100, shape[-1]),
+                np.prod(shape[:2])).reshape(shape)
+    fz = np.arange(np.prod(shape)).reshape(shape)
+    return z, fz 
+
+
+def interp_and_extrap(shape,
+                      interp=stratify.INTERPOLATE_LINEAR,
+                      extrap=stratify.EXTRAPOLATE_NEAREST):
+    z, fz = src_data(shape)
+    stratify.interpolate(np.linspace(-20, 120, 50), z, fz,
+                         interpolation=interp, extrapolation=extrap)
+
+
+if __name__ == '__main__':
+    interp_and_extrap(shape=(500, 600, 100))