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

[Jason Monschke]

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?

Cheers,
Jason
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190716/8a9bbe7f/attachment.html>