Dear Prof. Bangerth, Thanks a lot for your explanations. Indeed, the response was correct, is just the output a bit confusing.
Have a nice day, Longying On Wednesday, 17 August 2022 at 20:29:01 UTC+3 Wolfgang Bangerth wrote: > > > I am debugging a program using the function > > *'get_generalized_support_points()' *( where**has_support_points()=0, > > while has_generalized_support_points()=1*)*. My FE system is defined as > > *'FESystem<3> fe(FE_Nedelec<3>(0), 2);'*, therefore, each active > cell > > has 12*2 dofs. So I would also expect > 'get_generalized_support_points()' can > > return the support points with a vector of size 24 (of course the value > will > > repeat once). However, it only has 12 valid Point<3> values, the other > 12 are > > zero or some crazy number. > > > > My question is, is this reasonable, or there is sth wrong with my > understand > > on this. > > Longying & Jochen: > the latter. You assume that the array you are reading from has 12*2 > entries, > but the array really only has 12 entries and as a consequence when you > output > 24 elements, the latter half is accessing invalid memory. > > Generalized support points are the ones at which you need to know the > values > of a function to compute some kind of interpolant or projection. It is > enough > to know a function (which would have 2 components) at 12 points to > determine > the 24 coefficients of the interpolant because in your specific case, the > two > elements that describe the two components have support points at the same > location. > > 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/54cbbbe2-1858-4c1a-a19b-7a6f66680128n%40googlegroups.com.
