def METHOD_NAME(self):
    def f(a, b):
        ac = at_complex(a[0], a[1])
        return abs((ac * b) ** 2).sum()
    rng = np.random.default_rng(9333)
    aval = np.asarray(rng.standard_normal((2, 5)))
    bval = rng.standard_normal(5)
    try:
        utt.verify_grad(f, [aval, bval])
    except GradientError as e:
        print(e.num_grad.gf)
        print(e.analytic_grad)
        raise