Daniel, Thanks a lot for your reply, it answers my question fully!
Regards, Alexander On Tue, 2022-11-22 at 10:12 -0500, Daniel Arndt wrote: > Alexander, > > Yes, we are identifying DoFs on faces and thus FE_DGQ doesn't have > any DoFs on faces. > > Best, > Daniel > > On Tue, Nov 22, 2022 at 6:48 AM Alexander Kiselyov > <[email protected]> wrote: > > 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 [email protected]. > > > 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 [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/1550a49daa0c16f3079f7c6c99916cd34936e856.camel%40gmail.com.
