def METHOD_NAME(self):
    array = np.array([[0.0, 1.0, 2.0], [0.0, 1.0, 2.0], [0.0, 1.0, 2.0]])
    axis_values = np.array([0.0, 1.0, 2.0])
    f, dfdx, dfdy = interp_smoothed_array_and_derivatives(
        array, axis_values, axis_values, 0, 0, indexing="xy"
    )
    self.assertArraysAlmostEqual(
        [f([0.0, 1.0])[0], dfdx([0.0, 1.0])[0], dfdy([0.0, 1.0])[0]],
        [array[1][0], 1.0, 0.0],
    )