Chris Eldred <[email protected]> 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).