FEniCS / dolfinx

Next generation FEniCS problem solving environment
https://fenicsproject.org
GNU Lesser General Public License v3.0
777 stars 181 forks source link

Read mixed topology mesh from file #3231

Open chrisrichardson opened 5 months ago

chrisrichardson commented 5 months ago

Add functionality to read a mixed topology mesh (e.g. triangle + quad, hex+prism+pyramid+tet) from file.

garth-wells commented 5 months ago
  • file format? VTK, XDMF, ADIOS2

I'd suggest VTKHDF, with maybe VTK first because it's easy to read.

jorgensd commented 5 months ago

If we go for VTKHDF, please consider: https://discourse.vtk.org/t/vtkhdf-point-ownership/13913/4 as it is to me unclear how to actually make use of the partitioning offsets in VTKHDF. https://gist.github.com/jorgensd/6b3db748009ac1017929d7e7f3e83dc2 shows how the data-structure is written for a single cell type, which is fairly easy reverse engineer to reading and/or mixed topology.

chrisrichardson commented 5 months ago

Thanks @jorgensd - I will take a look.