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