def METHOD_NAME(self):
    # test whether inverse ellipticity distrotions can be fully described by reduced shear distortions
    # elliptical gaussian
    x, y = util.make_grid(numPix=10, deltapix=0.1)
    gauss = GaussianEllipse()
    shear = ShearReduced()
    flux_round = gauss.function(x, y, amp=1, sigma=0.1, e1=0, e2=0)
    (
        e1,
        e2,
    ) = (
        0,
        0.1,
    )
    f_x, f_y = shear.derivatives(x, y, gamma1=-e1, gamma2=-e2)
    flux_distorted = gauss.function(
        x - f_x, y - f_y, amp=1, sigma=0.1, e1=e1, e2=e2
    )
    npt.assert_almost_equal(flux_round, flux_distorted, decimal=9)
    (
        e1,
        e2,
    ) = (
        0.2,
        0,
    )
    f_x, f_y = shear.derivatives(x, y, gamma1=-e1, gamma2=-e2)
    flux_distorted = gauss.function(
        x - f_x, y - f_y, amp=1, sigma=0.1, e1=e1, e2=e2
    )
    npt.assert_almost_equal(flux_round, flux_distorted, decimal=9)
    (
        e1,
        e2,
    ) = (
        -0.2,
        0.1,
    )
    f_x, f_y = shear.derivatives(x, y, gamma1=-e1, gamma2=-e2)
    flux_distorted = gauss.function(
        x - f_x, y - f_y, amp=1, sigma=0.1, e1=e1, e2=e2
    )
    npt.assert_almost_equal(flux_round, flux_distorted, decimal=9)