<div dir="ltr">Hi, DMCreate() is used to load/distribute grid built from gmsh, and the function crashes when a mixed mesh of tetra and pyramids in 3D. It looks PETSc can handle the cell types of Tetra/Hexa/Prism/Pyramids/Polygon/Polyhedra in 3D. Thus I am unsure why it crashes when pyramids included in gmsh file. <br><br>Below functions are used to mesh distribution.<br>call DMCreate(PETSC_COMM_WORLD, dm_g, ierr);CHKERRA(ierr)<br>call DMSetType(dm_g, DMPLEX, ierr);CHKERRA(ierr)<br>call DMSetFromOptions(dm_g, ierr);CHKERRA(ierr)<br><br>Below error messages got from run.<br>[0]PETSC ERROR: #1 DMPlexGetRawFaces_Internal() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexinterpolate.c:312<br>[0]PETSC ERROR: #2 DMPlexInterpolateFaces_Internal() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexinterpolate.c:350<br>[0]PETSC ERROR: #3 DMPlexInterpolate() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexinterpolate.c:1327<br>[0]PETSC ERROR: #4 DMPlexCreateGmsh() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexgmsh.c:1634<br>[0]PETSC ERROR: #5 DMPlexCreateGmshFromFile() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexgmsh.c:1418<br>[0]PETSC ERROR: #6 DMPlexCreateFromFile() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexcreate.c:4721<br>[0]PETSC ERROR: #7 DMPlexCreateFromOptions_Internal() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexcreate.c:3212<br>[0]PETSC ERROR: #8 DMSetFromOptions_Plex() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexcreate.c:3433<br>[0]PETSC ERROR: #9 DMSetFromOptions() at /home/Mike/Library/petsc_partition/src/dm/interface/dm.c:887<br>[0]PETSC ERROR: #10 User provided function() at User file:0<br>Abort(63) on node 0 (rank 0 in comm 16): application called MPI_Abort(MPI_COMM_SELF, 63) - process 0<br><br>Can I get any comments on that? <br><br>Thanks,<br><br><br></div>