def METHOD_NAME():
    m = UnitSquareMesh(3, 3)
    x = SpatialCoordinate(m)
    V = FunctionSpace(m, 'RT', 1)
    f1 = Function(V)
    f2 = Function(V)
    project(as_vector((x[0], x[1])), f1)
    domain = "{[i]: 0 <= i < f1.dofs}"
    instructions = """
    for i
        f2[i, 0] = f1[i, 0]
    end
    """
    par_loop((domain, instructions), dS, {'f1': (f1, READ), 'f2': (f2, WRITE)})
    par_loop((domain, instructions), ds, {'f1': (f1, READ), 'f2': (f2, WRITE)})
    assert errornorm(f1, f2, degree_rise=0) < 1e-10