def METHOD_NAME(fn):
    num_f = getattr(num, fn)
    np_f = getattr(np, fn)
    par = (SIZES, range(NDIMS), DTYPES)
    for size, ndims, dtype in product(*par):
        shape = ndims * [size]
        xf = num_f(shape, dtype=dtype)
        yf = np_f(shape, dtype=dtype)
        assert np.array_equal(xf, yf)
        assert xf.dtype == yf.dtype