Dear Wolfgang, Thank you for your response. I identified the correct dofs by trial and error which worked (9 dofs for a scalar quantity with quadratic polyorder). However, when I use *"fe_face.system_to_base_index().first.first" *to identify, which base element these dofs belong to, it does not at all fit the theory as they should all belong to base 1, i.e. the electric potential (I copied the output from my terminal)
The base group of dof 3 is: 1 The base group of dof 10 is: 1 The base group of dof 17 is: 1 The base group of dof 24 is: 1 The base group of dof 31 is: 1 The base group of dof 41 is: 2 The base group of dof 51 is: 0 The base group of dof 61 is: 2 The base group of dof 71 is: 2 This explains, why my former code did not work, however it raises the question, why the dofs do not belong to the (correct) base element. Do you have any idea? Thank you and kind regards, Markus Wolfgang Bangerth schrieb am Dienstag, 15. Februar 2022 um 04:32:09 UTC+1: > On 2/11/22 12:12, 'Markus Mehnert' via deal.II User Group wrote: > > *** Caution: EXTERNAL Sender *** > > > > Dear Wolfgang, > > > > Thank you for your quick response. The dofs per face are defined as / > > / > > / > > / > > /const unsigned int dofs_per_face = fe_cell.dofs_per_face;/ > > / > > / > > where fe_cell is the/FESystem/ that consists of > > > > fe_cell(FE_Q<dim> (poly_order), dim, //Displacement > > FE_Q<dim> (poly_order), 1, //Electric Potential > > FE_Q<dim> (poly_order+1), dim, //Vectorial Electronic order > > parameter > > FE_DGPMonomial<dim>(poly_order - 1), 1, // Dilatation > > FE_DGPMonomial<dim>(poly_order - 1), 1), // Pressure > > > > I checked the number of dofs and it seems to be correct. > > Then it's time to count how many constraints you expect (on a piece of > paper), > how many you actually enter into the constraint object, and how many the > constraints object actually holds :-) > > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: [email protected] > www: http://www.math.colostate.edu/~bangerth/ > > -- 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/106f344a-45f5-4b19-81bf-8e92463dbb24n%40googlegroups.com.
