> 5. 11. 2017 v 18:48, Matthew Knepley <knep...@gmail.com>: > > On Thu, Nov 2, 2017 at 12:08 PM, Vaclav Hapla <vaclav.ha...@erdw.ethz.ch > <mailto:vaclav.ha...@erdw.ethz.ch>> wrote: > 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 > <http://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. > > It appears that MED, at least coming from GMsh, inverts hexes just like GMsh > does. I have fixed this in branch > > knepley/fix-plex-med-orientation > > If you run your test in this, everything should be fine.
Thanks a lot, Matt! It works for me. However, it hasn't fixed the problem with ParMetis - I will start another thread. Vaclav > > Michael, can you check whether this is a general fix, or it only applied to > MED from GMsh? > > Thanks, > > Matt > > > Vaclav > > > > > > > > > > > -- > What most experimenters take for granted before they begin their experiments > is infinitely more interesting than any results to which their experiments > lead. > -- Norbert Wiener > > https://www.cse.buffalo.edu/~knepley/ <http://www.caam.rice.edu/~mk51/>