def METHOD_NAME(self):
    x = np.linspace(0, 5, 10)
    y = np.linspace(0, 5, 10)
    amp = 1.0 * 2 * np.pi
    center_x = 0.0
    center_y = 0.0
    sigma = 1.0
    f_xx, f_xy, f_yx, f_yy = self.gaussian_kappa.hessian(
        x, y, amp, sigma, center_x, center_y
    )
    kappa = 1.0 / 2 * (f_xx + f_yy)
    amp_3d = self.gaussian_kappa._amp2d_to_3d(amp, sigma, sigma)
    density_2d = self.gaussian_kappa.density_2d(
        x, y, amp_3d, sigma, center_x, center_y
    )
    npt.assert_almost_equal(kappa[1], density_2d[1], decimal=5)
    npt.assert_almost_equal(kappa[2], density_2d[2], decimal=5)