Open
Description
Summarize the issue
One can create a mesh with "shared_vertex" as input, and it will not throw an error.
How to reproduce the bug
See below. Originally reported at
https://fenicsproject.discourse.group/t/unexpected-behaviour-of-ghostmode-shared-vertex/17342/2
Minimal Example (Python)
from dolfinx import mesh
from mpi4py import MPI
domain = mesh.create_unit_square(
MPI.COMM_WORLD, 4, 4, mesh.CellType.triangle,
ghost_mode=mesh.GhostMode.shared_vertex # or shared_facet
)
tdim = domain.topology.dim
num_local = domain.topology.index_map(tdim).size_local
num_ghost = domain.topology.index_map(tdim).num_ghostsOutput (Python)
Version
main branch
DOLFINx git commit
No response
Installation
No response
Additional information
No response
Metadata
Metadata
Assignees
Labels
No labels