Hello, The cell numbering between gmsh elementTags and doflinx original_cell_indexis not consistent for the following with this mesh.msh : import dolfinx.io.gmshio import numpy as np from mpi4py import MPI from dolfinx import fem from dolfinx.plot import vtk_mesh import pyvista import gmsh from mpi…

Ok, I understood how to proceed. Here are the updates: #[...] phys_grps = gmsh.model.getPhysicalGroups() numbering = [] for dim, tag in phys_grps: if dim == 2: entities = gmsh.model.getEntitiesForPhysicalGroup(dim, tag) for entity in entities: (entity_types, entity_…

Thanks for your reply. I already spotted this, but it is not exactly the same issue: here it is the numbering of cells instead of vertices (which [vertices] is actually ok in my case). In my case, it works for simple geometries with few physical entities – by using the mapping mesh_data.mesh.topol…

My question is: how to feed dolfinx data array from gmsh array following the elementTags numbering, i.e. tri_tags[0] <=> f_arr[:]: gmsh.view.addModelData(view_tag, 0, '', "ElementData", tri_tags[0], f_arr.reshape(-1, 1)) gmsh.view.addModelData

[image] pa-beaufort: here it is the numbering of cells instead I specifically mention cell indices: [image] Renumbered .msh imported data giving wrong subdomains Secondly, we do not associate the original_cell_index with the index it had in GMSH, as in the msh format, there might be cells…

ok, then I simply don’t understand the workflow you proposed

topologies_2 is the nodal definition of the cells. What’s the next step? I essentially need a numpy.array(dtype=int) (dolfinx2gmsh) mapping doflinx dofs to gmsh ones: gmsh.view.addModelData(view_tag, 0, '', "ElementData", tri_tags[0], f_arr.reshape(-1, 1)[dolfinx2gmsh]) equivalent to gmsh.view.…

I tried topologies_2 = [] entities = gmsh.model.getEntities() entities_sorted = np.sort([e[1] for e in entities]) for obj in entities_sorted: (entity_types, entity_tags, entity_topologies) = gmsh.model.mesh.getElements(2, tag=obj) assert len(entity_types) == 1 properties = gmsh.model.me…

[image] Renumbered .msh imported data giving wrong subdomains entities_sorted = np.sort([e[1] for e in entities]) for obj in entities_sorted: (entity_types, entity_tags, entity_topologies) = gmsh.model.mesh.getElements(2, tag=obj) assert len(entity_types) == 1 properties = gmsh.model.…

Cell numbering: inconsistent gmsh <=> dolfinx numbering

mesh

dokken January 22, 2025, 2:43pm 2

See: Renumbered .msh imported data giving wrong subdomains - #4 by dokken

Topic		Replies	Views
Renumbered .msh imported data giving wrong subdomains dolfinx mesh	3	182	October 9, 2024
Unable to import extruded geometry mesh mesh dolfinx	5	159	September 30, 2025
Tranfinite gmsh mesh and read_meshtags lead to unexpected crash Errors io , dolfinx , gmsh	14	833	March 20, 2024
Gmshio.model_to_mesh - error mesh	3	960	August 19, 2023
Regarding equivalent representation of dolfin.Mesh module expressions in dolfinx.gmshio dolfinx mesh , boundary-conditions	2	248	December 19, 2023

Cell numbering: inconsistent gmsh <=> dolfinx numbering

Related topics