ADIOS2 (2.8.3 and 2.9) BP4 files written with VTXWriter contain the mesh partition boundaries when viewed in ParaView (with opacity<1). Can this be avoided?
MWE: run with mpirun -n 2
import dolfinx
from dolfinx import fem, io
from mpi4py import MPI
mesh = dolfinx.mesh.create_unit_cube(MPI.COMM_WORLD, 16, 16, 16)
V = fem.FunctionSpace(mesh, ("CG", 1))
u = fem.Function(V, name="u")
with io.VTXWriter(MPI.COMM_WORLD, "test.bp4", [u], engine="BP4") as f:
f.write(0.0)
The boundary has to exist when you use ADIOS2 VTXWriter as the mesh is written in partitions, rather than as a global array with global numbering. This is how the ADIOS people want to have their mesh writer, ref VTK support for unstructured grids with global arrays. · Issue #3688 · ornladios/ADIOS2 · GitHub
where I propose the usage of global arrays, which would hide these partition boundaries (and add some features that would be useful for checkpointing).
Too bad, doesn’t look like they’re very open to that proposal.
I guess we’ll stick with XDMF and just hope it’s not going to be phased out very soon.
