def METHOD_NAME():
    A = np.random.rand(4, 4).astype(np.float32)
    B = np.random.rand(1).astype(np.float32)
    sdfg = create_reduce_sdfg("lambda a,b: a+b", (0, 1), "reduction_sum_all_axis", A, B, dace.float32)
    from dace.libraries.standard import Reduce
    Reduce.default_implementation = "FPGAPartialReduction"
    sdfg.expand_library_nodes()
    sdfg(A=A, B=B)
    assert np.allclose(B, np.sum(A, axis=(0, 1)))
    return sdfg