def METHOD_NAME(self):
    x = np.array([(0,0), (-0.5,-0.5), (-0.5,0.5), (0.5, 0.5), (0.25, 0.3)],
                 dtype=np.double)
    y = (np.arange(x.shape[0], dtype=np.double)[:,None]
         + np.array([0,1])[None,:])
    for method in ('nearest', 'linear', 'cubic'):
        for rescale in (True, False):
            msg = repr((method, rescale))
            yi = griddata((x[:,0], x[:,1]), y, (x[:,0], x[:,1]), method=method,
                          rescale=rescale)
            assert_allclose(y, yi, atol=1e-14, err_msg=msg)