def METHOD_NAME(self):
    from lenstronomy.LightModel.Profiles.nie import NIE as NIE_Light
    nie_light = NIE_Light()
    kwargs = {"e1": 0.3, "e2": -0.05, "s_scale": 0.5}
    x, y = util.make_grid(numPix=10, deltapix=0.1)
    f_xx, f_xy, f_yx, f_yy = self.nie.hessian(x, y, theta_E=1, **kwargs)
    kappa = 1 / 2.0 * (f_xx + f_yy)
    flux = nie_light.function(x, y, amp=1, **kwargs)
    npt.assert_almost_equal(kappa / np.sum(kappa), flux / np.sum(flux), decimal=5)