That works great, thanks! -Chris On Wed, Jul 31, 2013 at 4:59 PM, Jed Brown <[email protected]> wrote:

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.

-- Chris Eldred DOE Computational Science Graduate Fellow Graduate Student, Atmospheric Science, Colorado State University B.S. Applied Computational Physics, Carnegie Mellon University, 2009 [email protected] / [email protected]