[petsc-users] [DMCreate from a Gmsh with mixed mesh crash for DMPlex]

[Mike Michell]

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.

Below functions are used to mesh distribution.
call DMCreate(PETSC_COMM_WORLD, dm_g, ierr);CHKERRA(ierr)
call DMSetType(dm_g, DMPLEX, ierr);CHKERRA(ierr)
call DMSetFromOptions(dm_g, ierr);CHKERRA(ierr)

Below error messages got from run.
[0]PETSC ERROR: #1 DMPlexGetRawFaces_Internal() at
/home/Mike/Library/petsc_partition/src/dm/impls/plex/plexinterpolate.c:312
[0]PETSC ERROR: #2 DMPlexInterpolateFaces_Internal() at
/home/Mike/Library/petsc_partition/src/dm/impls/plex/plexinterpolate.c:350
[0]PETSC ERROR: #3 DMPlexInterpolate() at
/home/Mike/Library/petsc_partition/src/dm/impls/plex/plexinterpolate.c:1327
[0]PETSC ERROR: #4 DMPlexCreateGmsh() at
/home/Mike/Library/petsc_partition/src/dm/impls/plex/plexgmsh.c:1634
[0]PETSC ERROR: #5 DMPlexCreateGmshFromFile() at
/home/Mike/Library/petsc_partition/src/dm/impls/plex/plexgmsh.c:1418
[0]PETSC ERROR: #6 DMPlexCreateFromFile() at
/home/Mike/Library/petsc_partition/src/dm/impls/plex/plexcreate.c:4721
[0]PETSC ERROR: #7 DMPlexCreateFromOptions_Internal() at
/home/Mike/Library/petsc_partition/src/dm/impls/plex/plexcreate.c:3212
[0]PETSC ERROR: #8 DMSetFromOptions_Plex() at
/home/Mike/Library/petsc_partition/src/dm/impls/plex/plexcreate.c:3433
[0]PETSC ERROR: #9 DMSetFromOptions() at
/home/Mike/Library/petsc_partition/src/dm/interface/dm.c:887
[0]PETSC ERROR: #10 User provided function() at User file:0
Abort(63) on node 0 (rank 0 in comm 16): application called
MPI_Abort(MPI_COMM_SELF, 63) - process 0

Can I get any comments on that?

Thanks,
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20220713/71c5a565/attachment.html>