Submesh from parallel mesh

[dokken]

Not really sure what you are trying to do, please clarify.

In the previous post I showed how to transfer meshtags to the sub-mesh. Then you can output that for visualization/debugging.

If you want to use pyvista to plot in parallel, you should look at:

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*…