[MOAB-dev] Boundary nodes and faces

[Tim Tautges]

Note, I've added this to the MOAB faq, http://trac.mcs.anl.gov/projects/ITAPS/wiki/ParallelGeomSkin.

- tim

On 08/14/2013 09:31 PM, Tim Tautges wrote:
> This is tricky.  As you say, if you get the local skin, that'll include interior entities shared with/ghosted from other
> procs.  If you get the skin and remove the shared entities, that'll remove the vertices that are on the geometric skin
> but also shared with other procs.  The right way to do this is:
>
> - get the highest-dimensional entities on the skin (1d if in 2d, 2d if in 3d)
> - remove any shared entities from that range (because if it's shared, then by definition there will be d+1 - dimensional
> entities on the other proc)
> - 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
> on get_adjacencies), put in tmp range, then do the same for vertices, also put in tmp range
> - merge the tmp range with the original range
>
> Like I said, tricky, but welcome to parallel mesh.

-- 
================================================================
"You will keep in perfect peace him whose mind is
   steadfast, because he trusts in you."               Isaiah 26:3

              Tim Tautges            Argonne National Laboratory
          (tautges at mcs.anl.gov)      (telecommuting from UW-Madison)
  phone (gvoice): (608) 354-1459      1500 Engineering Dr.
             fax: (608) 263-4499      Madison, WI 53706