<div dir="ltr">Tim, Iulian,<div><br></div><div>Thanks for the pointers. I think going from the highest-dimension down and removing shared entities from the list makes sense. I have a partial implementation that already does this and will check whether this algorithm yields what I need. </div>

<div><br></div><div>I was originally thinking that Skinner:: find_geometric_skin would give me the boundary vertices but doing it by the process of elimination might be the right way since there is no reliance on GEOM_DIMENSION tag being explicitly defined.</div>

<div><br></div><div>Shall update soon.</div><div>Vijay</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Aug 14, 2013 at 9:50 PM, Tim Tautges <span dir="ltr"><<a href="mailto:tautges@mcs.anl.gov" target="_blank">tautges@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">Note, I've added this to the MOAB faq, <a href="http://trac.mcs.anl.gov/projects/ITAPS/wiki/ParallelGeomSkin" target="_blank">http://trac.mcs.anl.gov/<u></u>projects/ITAPS/wiki/<u></u>ParallelGeomSkin</a>.<span class="HOEnZb"><font color="#888888"><br>


<br>
- tim</font></span><div class="im HOEnZb"><br>
<br>
On 08/14/2013 09:31 PM, Tim Tautges wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
This is tricky.  As you say, if you get the local skin, that'll include interior entities shared with/ghosted from other<br>
procs.  If you get the skin and remove the shared entities, that'll remove the vertices that are on the geometric skin<br>
but also shared with other procs.  The right way to do this is:<br>
<br>
- get the highest-dimensional entities on the skin (1d if in 2d, 2d if in 3d)<br>
- remove any shared entities from that range (because if it's shared, then by definition there will be d+1 - dimensional<br>
entities on the other proc)<br>
- if you're in 3d and you also want skin edges, get the edges adj to the entities in the range (remember the UNION arg<br>
on get_adjacencies), put in tmp range, then do the same for vertices, also put in tmp range<br>
- merge the tmp range with the original range<br>
<br>
Like I said, tricky, but welcome to parallel mesh.<br>
</blockquote>
<br></div><div class="HOEnZb"><div class="h5">
-- <br>
==============================<u></u>==============================<u></u>====<br>
"You will keep in perfect peace him whose mind is<br>
  steadfast, because he trusts in you."               Isaiah 26:3<br>
<br>
             Tim Tautges            Argonne National Laboratory<br>
         (<a href="mailto:tautges@mcs.anl.gov" target="_blank">tautges@mcs.anl.gov</a>)      (telecommuting from UW-Madison)<br>
 phone (gvoice): <a href="tel:%28608%29%20354-1459" value="+16083541459" target="_blank">(608) 354-1459</a>      1500 Engineering Dr.<br>
            fax: <a href="tel:%28608%29%20263-4499" value="+16082634499" target="_blank">(608) 263-4499</a>      Madison, WI 53706<br>
<br>
</div></div></blockquote></div><br></div>