That works great, thanks!<br><br>-Chris<br><br><div class="gmail_quote">On Wed, Jul 31, 2013 at 4:59 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">Chris Eldred <<a href="mailto:chris.eldred@gmail.com">chris.eldred@gmail.com</a>> writes:<br>
> if (rank==0) then<br>
> call DMPlexCreateFromCellList(PETSC_COMM_WORLD, 2, 2, 4, 3,  1 , cell_list<br>
> , 2 , vertex_coords , dm,  ierr)<br>
<br>
</div>The problem here is that after this line, dm is garbage on rank 1.  You<br>
want to start by creating the DM on COMM_WORLD, but have all the cells<br>
on rank 0.  Try the attached.<br>
<br>
</blockquote></div><br><br clear="all"><br>-- <br>Chris Eldred<br>DOE Computational Science Graduate Fellow<br>Graduate Student, Atmospheric Science, Colorado State University<br>B.S. Applied Computational Physics, Carnegie Mellon University, 2009<br>
<a href="mailto:chris.eldred@gmail.com" target="_blank">chris.eldred@gmail.com</a> / <a href="mailto:celdred@atmos.colostate.edu" target="_blank">celdred@atmos.colostate.edu</a>