<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, May 19, 2015 at 10:32 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Great thanks, I'll look more into this.<br>
<br>
Question on the dofs that wasn't clear from the manual - are the dofs on cells, faces, edges, and vertices *all* actual degrees of freedom / solution unknowns in the system to solve, or does one say the dofs on the cells represent actual solution unknowns, and face dofs are just things to be communicated for inter-cell flux computation? e.g. say i was doing the simplest cellwise-constant finite volume scheme for a scalar; do i have 1 dof per cell, and zero dofs anywhere else?<br></blockquote><div><br></div><div>There are many ways you could imagine structuring the data. The thing about PetscSection is that is easily allows</div><div>you to pick any one of those.</div><div><br></div><div>In a FV example for PETSc, TS ex11, we choose to store only cell data. Then in loops over the faces, you</div><div><br></div><div>  1) Call DMPlexGetSupport(dm, face, &cells) to get the two cells</div><div><br></div><div>  2) Call DMPlexGetPointGlobalRead(dm, cell[i], solArray, &val) to get the cell data</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
You can easily get started by just calling DMPlexDistribute() on a serial mesh.<br>
<br>
</blockquote>
<br>
perfect, i'll start with that.<br>
<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
<br>
I do cG FEM in SNES ex12 and ex62.<br>
<br>
</blockquote>
<br>
Great thanks, didn't see these linked from the DMPlex doc pages.<br>
</blockquote></div><br>They are in the dev docs I believe.</div><div class="gmail_extra"><br></div><div class="gmail_extra">  Thanks,</div><div class="gmail_extra"><br></div><div class="gmail_extra">    Matt<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>