def METHOD_NAME(self):
    """Test the `w_f_approx()` method with values computed using
    `scipy.special.wofz()`.
    :return:
    :rtype:
    """
    x = np.logspace(-3.0, 3.0, 100)
    y = np.logspace(-3.0, 3.0, 100)
    X, Y = np.meshgrid(x, y)
    w_f_app = self.gaussian_kappa_ellipse.w_f_approx(X + 1j * Y)
    w_f_scipy = wofz(X + 1j * Y)
    npt.assert_allclose(w_f_app.real, w_f_scipy.real, rtol=4e-5, atol=0)
    npt.assert_allclose(w_f_app.imag, w_f_scipy.imag, rtol=4e-5, atol=0)
    # check `derivatives()` method with and without `scipy.special.wofz()`
    x = 1.0
    y = 1.0
    e1, e2 = 5e-5, 0
    sigma = 1.0
    amp = 2.0
    # with `scipy.special.wofz()`
    gauss_scipy = GaussianEllipseKappa(use_scipy_wofz=True)
    f_x_sp, f_y_sp = gauss_scipy.derivatives(x, y, amp, sigma, e1, e2)
    # with `GaussEllipseKappa.w_f_approx()`
    gauss_approx = GaussianEllipseKappa(use_scipy_wofz=False)
    f_x_ap, f_y_ap = gauss_approx.derivatives(x, y, amp, sigma, e1, e2)
    npt.assert_almost_equal(f_x_sp, f_x_ap, decimal=4)
    npt.assert_almost_equal(f_y_sp, f_y_ap, decimal=4)