[petsc-dev] DMPlexCreateMedFromFile broken

[Vaclav Hapla]

It seems that DMPlexCreateMedFromFile leaves out some mesh elements. I found it out when investigating why ParMetis redistribution crashes.

I attach the datafile $PETSC_DIR/share/petsc/datafiles/meshes/blockcylinder-50.exo converted to GMSH and MED format.
The conversion EXO to GMSH was done by meshio (github.com/nschloe/meshio), and GMSH to MED by GMSH itself.

I did:

cd $PETSC_DIR/src/dm/impls/plex/examples/tutorials
make ex5
FILE=blockcylinder-50.exo && ./ex5 -filename $FILE -dm_view hdf5:$FILE.h5 && $PETSC_DIR/bin/petsc_gen_xdmf.py $FILE.h5
FILE=blockcylinder-50.msh && ./ex5 -filename $FILE -dm_view hdf5:$FILE.h5 && $PETSC_DIR/bin/petsc_gen_xdmf.py $FILE.h5
FILE=blockcylinder-50.med && ./ex5 -filename $FILE -dm_view hdf5:$FILE.h5 && $PETSC_DIR/bin/petsc_gen_xdmf.py $FILE.h5

While the output from blockcylinder-50.exo and blockcylinder-50.msh looks OK, that from blockcylinder-50.med is corrupted.

I also attach my HDF5/XDMF outputs and screenshots from ParaView.

Vaclav

-------------- next part --------------
A non-text attachment was scrubbed...
Name: inputs.zip
Type: application/zip
Size: 6336 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20171102/1aad1998/attachment-0003.zip>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: outputs.zip
Type: application/zip
Size: 11801 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20171102/1aad1998/attachment-0004.zip>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: screenshots.zip
Type: application/zip
Size: 49663 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20171102/1aad1998/attachment-0005.zip>
-------------- next part --------------