<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, May 19, 2015 at 10:14 AM, Mark Lohry <span dir="ltr"><<a href="mailto:mlohry@princeton.edu" target="_blank">mlohry@princeton.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">I'm interested in trying DMPlex for a parallel 3D unstructured DG solver.  Is it appropriate for this, and is it considerably more work than for the simple DMDACreate3D examples? The manual example seems fairly clear on how to set up a mesh in serial, but does this play nicely with METIS for partitioning, and DMDA for parallel arrays/communication? I'm typically using CGNS files, so is DMPlexCreateCGNS working in parallel and appropriate for this?<br></blockquote><div><br></div><div>Plex can certainly handle 3D unstructured meshes, read from CGNS, and can partition them and repartition them.</div><div><br></div><div>You can easily get started by just calling DMPlexDistribute() on a serial mesh.</div><div><br></div><div>To use it to handle fields, you construct a PetscSection that models your data layout for the DG field. This just means</div><div>tell it how many dofs sit on each cell, face, edge, and vertex. Then it can handle all the parallel communication for you.</div><div>Note that you need to specify what kind of adjacency you have for it to automatically construct the comm pattern:</div><div><br></div><div>  <a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DM/DMPlexSetAdjacencyUseClosure.html">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DM/DMPlexSetAdjacencyUseClosure.html</a></div><div>  <a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DM/DMPlexSetAdjacencyUseCone.html#DMPlexSetAdjacencyUseCone">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DM/DMPlexSetAdjacencyUseCone.html#DMPlexSetAdjacencyUseCone</a></div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
Any examples/guidance related to this would be appreciated.</blockquote><div><br></div><div>I do cG FEM in SNES ex12 and ex62.</div><div><br class="">  Thanks,</div><div><br></div><div>     Matt </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><span class=""><font color="#888888"><br>
-Mark<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>