<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">On Thu, May 9, 2019 at 2:07 PM Praveen C via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear all<br>
<br>
I am trying to understand partitioning using DMPlexDistribute.<br>
<br>
I see an option “overlap”. Does this determine how many ghost cells are added ? Are these determined based on face neighboring cells or vertex neighboring cells ? <br></blockquote><div><br></div><div>Its user settable: <a href="https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DM/DMSetAdjacency.html">https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DM/DMSetAdjacency.html</a></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
In<br>
<br>
<a href="https://www.mcs.anl.gov/petsc/petsc-current/src/ts/examples/tutorials/ex11.c.html" rel="noreferrer" target="_blank">https://www.mcs.anl.gov/petsc/petsc-current/src/ts/examples/tutorials/ex11.c.html</a><br>
<br>
I see a call to DMPlexConstructGhostCells after DMPlexDistribute. Can you explain the purpose of this and how it is related to overlap ?<br></blockquote><div><br></div><div>This is a different thing. It is a way to enforce boundary conditions in finite volume methods. You make a fake cell on the other side of every boundary face so that you can put a state there which creates the right face flux.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
After DMPlexDistribute, how can I identify locally owned cells and ghost cells ?<br></blockquote><div><br></div><div>1) Ghost cells are special, not shared with other processors, and marked with a label</div><div><br></div><div>2) Overlap cells are regular cells, and shared with other processes. You can determine if you own the cell by</div><div>    checking whether it is present in the point PetscSF: <a href="https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DM/DMGetPointSF.html">https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DM/DMGetPointSF.html</a></div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Thanks<br>
praveen<br>
<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div></div></div>