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]
