Hi Alexander, You are using FE_DGQ, which per definition only has DoFs in the cell. Switch to quadratic FE_Q and then your idea with the dummy DoFHandler to globally uniquely identify faces should work! Not sure why you create a system!?
Peter On Friday, 18 November 2022 at 18:30:38 UTC+1 Alexander Kiselyov wrote: > Dear deal.II users, > > I'm currently developing a way to uniquely identify a boundary face on > a mesh. My idea has been to gather global DoF indices on a face and > hash them together. However, `get_dof_indices` method of a face > accessor doesn't return any useful information, and > `FESystem::max_dofs_per_face()`, surprisingly, returns 0. See the > attachment for a minimal example. The problem persists when running on > a single process. The version of deal.II used is 9.4.0. > > What am I doing wrong here? Is there a better way to globally identify > a face? > > Best regards, > Alexander Kiselyov > -- 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 [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/54a20ebc-b6e2-479d-8eed-d236e1ee90cdn%40googlegroups.com.
