Plot mesh parallel, show partitioning

[dokken]

I am not quite sure what you are trying to achieve here.
You have ignored the crucial:

# Map to global dof indices
global_dofs = V.dofmap.index_map.local_to_global(topology[topology_dofs].copy())
# Overwrite topology
topology[topology_dofs] = global_dofs

How to collect the global mesh without writing a file

Here is an example code of how to do it: # Copyright (C) 2023 Jørgen S. Dokken # # SPDX-License-Identifier: MIT # # Code to plot mesh with pyvista in parallel by gather data on a root process # # Last changed: 2023-10-23 import dolfinx from mpi4py import MPI import pyvista import numpy as np pyvista.start_xvfb(0.2) mesh = dolfinx.mesh.create_unit_square(MPI.COMM_WORLD, 10, 10) V = dolfinx.fem.functionspace(mesh, ("Lagrange", 2)) u = dolfinx.fem.Function(V) u.interpolate(lambda x: x[0] + 2*…

This is crucial as one maps the local topology indices to a single global numbering.
For more about this see for instance section 6.2 of:

DOLFINx: The next generation FEniCS problem solving environment

which tries to explain global ordering.