Hi Maurice!

On Monday, November 16, 2020 at 8:09:14 AM UTC-7 maurice....@googlemail.com 
wrote:

> Looking at the doc of `collect_periodic_faces` (This function will collect 
> periodic face pairs on the coarsest mesh level of the given mesh (a 
> Triangulation 
> <https://www.dealii.org/9.0.0/doxygen/deal.II/classTriangulation.html> or 
> DoFHandler 
> <https://www.dealii.org/9.0.0/doxygen/deal.II/classDoFHandler.html>) and 
> add them to the vector matched_pairs leaving the original contents 
> intact.)  the resulting vector will only contain parent cells, on which the 
> call `is_artifiical()` is not possible. 
>

It appears to me that `GridTools::get_active_child_cells()` is the function 
you are looking for (see here 
<https://www.dealii.org/current/doxygen/deal.II/namespaceGridTools.html#ab5c182c1901bb0d09046271809f918c6>).
 
In your case, you would get all active children of the coarsest mesh cells 
if applied on the results of `GridTools::collect_periodic_faces()`. You can 
iterate over those and identify which one of them are located at the 
boundary.

However, I am not sure whether `GridTools::get_active_child_cells()` works 
with `parallel::distributed::Triangulations`. If it does though, the 
results may contain locally owned, ghost, and/or artificial cells. As the 
results will be all active cells, you are able to check for these 
attributes, i.e., check which ones are locally relevant.

Hope this helps!

Marc

-- 
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- 
You received this message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/ea516d6f-1eeb-411d-9ca6-c7bf483efe7fn%40googlegroups.com.

Reply via email to