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.