Dear Peter, Thank you for the reply!
The thing is that the input data is a FE_DGQ DoF vector, and I would rather not project it into FE_Q elements. The original plan was not to create a separate dummy DoFHandler, but your idea seems like a solution for my issue. Still, I'm a bit confused about why FE_DGQ doesn't contain face DoFs. Is the reason that in this case DoF indices would be different when looking on either sides of a face? Regards, Alexander On Sun, 2022-11-20 at 08:43 -0800, Peter Munch wrote: > 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 dealii+unsubscr...@googlegroups.com. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/54a20ebc-b6e2-479d-8eed-d236e1ee90cdn%40googlegroups.com > . -- 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/a3badbdabac248c7524fb6fd353138a8ad525e0c.camel%40gmail.com.
