def METHOD_NAME(self):
    x = np.array([1.0])
    y = np.array([2.0])
    phi_E = 1.0
    gamma = 2.0
    q = 0.9
    phi_G = 1.0
    s_scale = 0.001
    e1, e2 = param_util.phi_q2_ellipticity(phi_G, q)
    f_xx, f_xy, f_yx, f_yy = self.SPEMD.hessian(x, y, phi_E, gamma, e1, e2, s_scale)
    if fastell4py_bool:
        f_xx_nie, f_xy_nie, f_yx_nie, f_yy_nie = self.NIE.hessian(
            x, y, phi_E, e1, e2, s_scale
        )
        npt.assert_almost_equal(f_xx, f_xx_nie, decimal=4)
        npt.assert_almost_equal(f_yy, f_yy_nie, decimal=4)
        npt.assert_almost_equal(f_xy, f_xy_nie, decimal=4)
        npt.assert_almost_equal(f_yx, f_yx_nie, decimal=4)
    else:
        npt.assert_almost_equal(f_xx, 0, decimal=7)
        npt.assert_almost_equal(f_yy, 0, decimal=7)
        npt.assert_almost_equal(f_xy, 0, decimal=7)
    npt.assert_almost_equal(f_xy, f_yx, decimal=8)