<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, Oct 28, 2016 at 8:05 AM, leejearl <span dir="ltr"><<a href="mailto:leejearl@126.com" target="_blank">leejearl@126.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
Hi, Matt:<br>
Thank you for your help. It is the answer which I expect. <br>
But, I have an another equestion. Can I obtain the coordinates
of the ghost cells:<br>
</div></blockquote><div><br></div><div>I am not sure I understand. By default, Plex stores coordinates only for vertices. Also, ghost cells</div><div>do not have a full complement of vertices because they are only used to define a state on the other</div><div>side of a face. Does that make sense?</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div> Thanks,<br>
<br>
leejearl<br>
<br>
<div class="m_7510602955081850294moz-cite-prefix">On 2016年10月28日 19:56, Matthew Knepley
wrote:<br>
</div>
<blockquote type="cite">
<div dir="ltr">
<div class="gmail_extra">
<div class="gmail_quote">On Fri, Oct 28, 2016 at 6:32 AM,
leejearl <span dir="ltr"><<a href="mailto:leejearl@126.com" target="_blank">leejearl@126.com</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi, everyone:<br>
<br>
I have a distributed DMPlex representing a 2D finite
volume mesh. I distributed it using the function
"DMPlexDistribute()".<br>
<br>
Then, I constructed the ghost cells using
"DMPlexConstructGhostCells()".<br>
<br>
The question is that I want to determine which cell is
ghosted. I obtained all the cells using
"DMPlexGetHeightStratum()", and<br>
<br>
how can I kown which cell is ghosted?<br>
</blockquote>
<div><br>
</div>
<div>1) The ghost cells are marked as "hybrid", so you can
call</div>
<div><br>
</div>
<div> ierr = DMPlexGetHybridBounds(dm, &cEndInterior,
NULL, NULL, NULL);CHKERRQ(ierr);</div>
<div>
<div> ierr = DMPlexGetHeightStratum(dm, 0, NULL,
&cEnd);CHKERRQ(ierr);</div>
</div>
<div><br>
</div>
<div>and then the ghost cells are [cEndInterior, cEnd).</div>
<div><br>
</div>
<div>2) There is also a label "ghost" which marks the ghost
faces, and ghost cells which are non-local</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, all the helps are appreciated.<span class="m_7510602955081850294gmail-HOEnZb"><font color="#888888"><br>
<br>
<br>
leejearl<br>
<br>
<br>
<br>
</font></span></blockquote>
</div>
<br>
<br clear="all"><span class="HOEnZb"><font color="#888888">
<div><br>
</div>
-- <br>
<div class="m_7510602955081850294gmail_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>
</font></span></div>
</div>
</blockquote>
<br>
<pre class="m_7510602955081850294moz-signature" cols="72"></pre>
</div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="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>