def METHOD_NAME():
    mesh = UnitIntervalMesh(10)
    d = Function(FunctionSpace(mesh, "DG", 0))
    consts = []
    for i in range(20):
        consts.append(Constant(i, domain=mesh))
    arg = {}
    if mesh.comm.rank == 0:
        arg['d'] = (d, WRITE)
        for i, c in enumerate(consts):
            arg["c%d" % i] = (c, READ)
    else:
        arg['d'] = (d, WRITE)
        for i, c in enumerate(reversed(consts)):
            arg["c%d" % (len(consts) - i - 1)] = (c, READ)
    domain = "{[i]: 0 <= i < d.dofs}"
    instructions = """
    for i
        d[i, 0] = c10[0]
    end
    """
    par_loop((domain, instructions), dx, arg)
    assert np.allclose(d.dat.data, consts[10].dat.data)