def METHOD_NAME(dtype, outtype, op):
    shape = (1, 2, 3)
    out_shape = (6,)
    arr_np = mk_0to1_array(np, shape)
    arr_num = mk_0to1_array(num, shape)
    out_np = np.empty(out_shape, dtype=outtype)
    out_num = num.empty(out_shape, dtype=outtype)
    # NumPy ndarray doesn't have nancumprod, use module level API instead
    getattr(np, op)(arr_np, dtype=dtype, out=out_np)
    getattr(arr_num, op)(dtype=dtype, out=out_num)
    assert np.allclose(out_np, out_num)