[petsc-dev] DMPlexCreateMedFromFile broken

[Vaclav Hapla]

> 5. 11. 2017 v 18:48, Matthew Knepley <knepley at gmail.com>:
> 
> On Thu, Nov 2, 2017 at 12:08 PM, Vaclav Hapla <vaclav.hapla at erdw.ethz.ch <mailto:vaclav.hapla at 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/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20171106/fcfde14a/attachment.html>