[petsc-dev] Segfault in DMPlexDistribute

[Jed Brown]

Chris Eldred <chris.eldred at gmail.com> writes:
> Here is my relevant code:
>
> cell_list = (/ 0 , 1 , 2  , 1 , 3 , 2 /)
> vertex_coords = 0.0
>
> call DMPlexCreateFromCellList(PETSC_COMM_WORLD, 2, 2, 4, 3,  1 , cell_list
> , 2 , vertex_coords , dm,  ierr)
> CHKERRQ(ierr)

Which process executes this line?  DMPlexDistribute expects rank 0 to
hold the mesh to be distributed.  If this doesn't fix the problem, can
you send a test code so we can debug and add it to our test suite?

> call DMPlexDistribute(dm,'chaco',0,dm_distrib,ierr)
> CHKERRQ(ierr)
>
> I am running the code with mpirun -np 2 ./test (it also segfaults with
> higher values for np).
>
> Any ideas? The code works fine with mpirun -np 1 (presumably since it is
> not actually attempting to distribute the mesh).

-------------- 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/030f4ef7/attachment.sig>