dokken
September 11, 2023, 3:19pm
10
emiel_brinkhuis1:
Traceback (most recent call last):
File "/home3/s3234266/nsx/examples/scripts/readmsh.py", line 4, in <module>
mesh, cell_markers, facet_markers = gmshio.read_from_msh("../meshes/patient9cgs.msh", MPI.COMM_WORLD, gdim=3)
File "/home3/s3234266/spack/var/spack/environments/fenicsx/.spack-env/view/lib/python3.10/site-packages/dolfinx/io/gmshio.py", line 311, in read_from_msh
msh = model_to_mesh(gmsh.model, comm, rank, gdim=gdim, partitioner=partitioner)
File "/home3/s3234266/spack/var/spack/environments/fenicsx/.spack-env/view/lib/python3.10/site-packages/dolfinx/io/gmshio.py", line 221, in model_to_mesh
cell_id = cell_information[perm_sort[-1]]["id"]
IndexError: index -1 is out of bounds for axis 0 with size 0
So this error message means that you are not using phyiscal entity grouping in Gmsh, Which is adviced to avoid Getting duplicate nodes and cells in your mesh.