[petsc-users] DMPlex, multilevel join to find face from vertices

[Matthew Knepley]

On Tue, Jul 16, 2019 at 11:33 AM Jason Monschke via petsc-users <
petsc-users at mcs.anl.gov> wrote:

> Hello,
>
> I have a DMPlex 3D mesh containing cells, faces, edges, and vertices. What
> is the best way to get the face which contains a set of vertices in this
> mesh? I can find the edge which contains a set of vertices using
> DMPlexGetJoin. However, the documentation explicitly says that it does not
> work for multilevel joins.
>
> My current idea is to:
> 1. Find the supporting edges of each vertex using DMPlexGetSupport and add
> them to a list of edges
> 2. Remove all edges that only appear once in the list
> 3. Use DMPlexGetJoin to find the face that is the join of these edges
>
> Does this sound like it will work? Is there a better way to do this?
>

https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DMPLEX/DMPlexGetFullJoin.html

I could not think of a better name.

  Thanks,

    Matt


> Cheers,
> Jason
>


-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener

https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190716/2359738a/attachment.html>