Dofmap[icell] and links(icell) return different local connectivities in parallel

[dokken]

The geometry and topology dofmap does not necessarily use the same ghost ordering.

You need to use the function dolfinx.mesh.entities_to_geometry to map vertices to geoemtry nodes. See for instance

Locate_entities returning incorrect points with MPI

Since entities can be edges, facets or cells, it would return (num entities, num nodes per entity). Since your case is a vertex, there is only one node (1-1 map). I would call: idx = entities_to_geometry(mesh, 0, entity_idx).flatten() Rather than doing it after