Hi everyone. I’m trying to read a mesh generated by Ansys in FEniCSx. The mesh is in .msh. I change the format with meshio convert fluent.msh venturi_1600.xdmf but when I try to open import dolfinx #################################################################################################…

I have tried that code with the following change: meshio.write("venturi_1600.xdmf, create_mesh, (mesh, "quadrilateral", True)) and I get the next error KeyError: 'quadrilateral' I find the name of elements here: https://fenicsproject.org/olddocs/dolfinx/dev/python/generated/dolfinx.cpp.mesh.ht…

[image] Karlosjusus: meshio.write("venturi_1600.xdmf, create_mesh, (mesh, "quadrilateral", True)) This does not make sense, it should be: meshio.write("venturi_1600.xdmf, create_mesh(mesh, "quad", True)) if mesh is a meshio mesh. See https://github.com/nschloe/meshio/blob/0138cc8692b806b44…

Thank you @dokken , but now I try def create_mesh(mesh, cell_type, prune_z=False): cells = mesh.get_cells_type(cell_type) cell_data = mesh.get_cell_data("gmsh:physical", cell_type) points = mesh.points[:,:2] if prune_z else mesh.points out_mesh = meshio.Mesh(points=points, cells={c…

I think you need to remove the cell_data part and try like - def create_mesh(mesh, cell_type, prune_z=False): cells = mesh.get_cells_type(cell_type) points = mesh.points[:,:2] if prune_z else mesh.points out_mesh = meshio.Mesh(points=points, cells={cell_type: cells}) return out_mesh…

sorry for the delay. When I run it, I get this error. Warnig: py.warning:/urs/local/lib/python3.10/dist-packages/h5py/__init__.py:36 UserWarning: h5py is running against HDF5 1.10.7 when it was bulit against 1.12.2, this may cause problems and the kernel dies… I tried to install these version bu…

See The kernel appears to have died. It will restart automatically - #13 by dokken

I followed the steps, uninstalled h5py and ran in terminalI followed the steps, uninstalled h5py and ran in terminal python3 -m pip install --no-binary=h5py h5py meshio and I get Defaulting to user installation because normal site-packages is not writeable Collecting h5py Using cached h5py-3.9.…

Could you check what cython version you have on your system, as they just released a new version 2 days ago (it might be that you need to pin it to an earlier version of cython). [image] Releases · cython/cython The most widely used Python to C compiler. Contribute to cython/cython devel…

Problem with read a external mesh

General

dokken July 6, 2023, 2:55pm 2

See: Can dolfinx read a GMSH file containing 1D mesh of an ellipse? - #2 by dokken

Topic		Replies	Views
Importing mesh gives error dolfinx	12	3372	December 9, 2021
Read a mesh by meshio mesh	12	3187	April 5, 2023
The kernel appears to have died. It will restart automatically General	25	2688	November 24, 2022
Problem when reading .h5 mesh in dolfin General	3	529	July 4, 2023
Fluent.msh to XDMF; error mesh	12	3733	January 22, 2020

Problem with read a external mesh

Related topics