def METHOD_NAME(m, geo):
    mesh_name, dim, e_meas, e_s_1, e_s_2, e_s_bdy=geo
    m2d= feelpp.load(m, mesh_name, 0.1)
    Xh = feelpp.functionSpace(mesh=m2d)
    if feelpp.Environment.isMasterRank():
        print("Xh basisname: ", Xh.basisName())
        print("Xh nDof: ", Xh.nDof())
        print("Xh nLocalDof: ", Xh.nLocalDof())
        print("Xh nLocalDofWithGhost: ", Xh.nLocalDofWithGhost())
        print("Xh nLocalDofWithoutGhost: ", Xh.nLocalDofWithoutGhost())
    m3 = Xh.mesh()
    assert m3 == m2d
    b = feelpp.backend(worldcomm=feelpp.Environment.worldCommPtr())
    v = b.newVector(dm=Xh.mapPtr())
    # test constant
    v.setConstant(1.0)
    # create u which use the same memory storage as v
    u = Xh.element(v)
    w = (2*u-u)-u
    l2_w = feelpp.normL2(range=feelpp.elements(m2d), expr=w)
    print("norm(w)", l2_w)
    assert abs(l2_w) < 1e-12
    uu = Xh.element()
    uu.on(range=feelpp.elements(m2d), expr=feelpp.expr("x:x"))
    vv = uu.to_petsc()
    u.on(range=feelpp.elements(m2d), expr=feelpp.expr("x:x"))# now v contains also x
    w=uu-u
    l2_w = feelpp.normL2(range=feelpp.elements(m2d), expr=w)
    print("norm(uu-u)={}, uu.l2={}, u.l2={}".format(l2_w, uu.l2Norm(), u.l2Norm()))
    assert abs(l2_w) < 1e-12
    vv-=v
    l2_w = w.l2Norm()
    print("norm(vv_petsc-v_petsc)={}, vv.l2={}, v.l2={}".format(l2_w, vv.l2Norm(), v.l2Norm()))
    assert abs(l2_w) < 1e-12