def METHOD_NAME(m, geo):
    tic()
    mesh_name, dim, e_meas, e_s_1, e_s_2, e_s_bdy=geo
    m2d= feelpp.load(m, mesh_name, 0.1)
    toc("Mesh")
    tic()
    Xh = feelpp.functionSpace(mesh=m2d)
    toc("functionSpace")
    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
    tic()
    u = Xh.element()
    toc("FunctionSpace::Element")
    tic()
    u.on(range=feelpp.elements(m2d), expr=feelpp.expr("x:x"))
    toc("FunctionSpace::Element::on")
    tic()
    w = (2*u-u)-u
    l2_w = feelpp.normL2(range=feelpp.elements(m2d), expr=w)
    toc("FunctionSpace::Element::normL2")
    print("norm(w)", l2_w)
    assert abs(l2_w) < 1e-12
    assert u.functionSpace() == Xh
    assert u.size() == Xh.nDof()