[petsc-dev] Segfault in DMPlexDistribute

[Jed Brown]

Chris Eldred <chris.eldred at gmail.com> writes:
> if (rank==0) then
> call DMPlexCreateFromCellList(PETSC_COMM_WORLD, 2, 2, 4, 3,  1 , cell_list
> , 2 , vertex_coords , dm,  ierr)

The problem here is that after this line, dm is garbage on rank 1.  You
want to start by creating the DM on COMM_WORLD, but have all the cells
on rank 0.  Try the attached.

-------------- next part --------------
A non-text attachment was scrubbed...
Name: t.F90
Type: text/x-fortran
Size: 1101 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130801/8ddbe7c2/attachment.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130801/8ddbe7c2/attachment.sig>