Store function with meshio

[dokken]

I would suggest saving your data to xdmf, as explained here: https://jsdokken.com/converted_files/tutorial_pygmsh.html#second

I’ve also made a similar version coming from cell-sets (pygmsh)

Subdomain definition based on input file

How about: import meshio import numpy as np mesh = meshio.read("mesh.inp") cells = mesh.get_cells_type("triangle") cell_sets = mesh.cell_sets_dict cell_data = np.zeros(len(cells)) for marker, set in cell_sets.items(): for type, entities in set.items(): if type == "triangle": # Create marker with int from integer from name of input cell_data[entities] = int(marker[-1]) out_mesh = meshio.Mesh(points=mesh.points, cells={ "triangle": cells…

and

Transitioning from mesh.xml to mesh.xdmf, from dolfin-convert to meshio

I had a play with pygmsh over the weekend, and what seems to be done in the later versions is to save the physical markers to cell sets. I have made the following function for this: def create_entity_mesh(mesh, cell_type, prune_z=False, remove_unused_points=False): """ Given a meshio mesh, extract mesh and physical markers for a given entity. We assume that all unused points are at the end of the mesh.points (this happens when we use physical markers with…