On 9/23/19 11:12 PM, Phạm Ngọc Kiên wrote: > The function fe.has_support_points() returns false. This means that the > FE_Raviart_Thomas does not have support points. > I think generalized_support_points are used for this type of finite element.
You are right -- the FE_RT is defined in terms of integrals on the faces, not point values. > For FE_RaviartThomas > <https://www.dealii.org/current/doxygen/deal.II/classFE__RaviartThomas.html#adc8e64d0cdb38af6e29dc9ac0279f121>(0), > > are the dofs defined at the centers of the faces of a cell? No, in fact they're not defined through interpolation at points at all, but the degree of freedom is the integral along an entire face. But at least for RT(0), that's only a constant multiple of the point value at the center of the cell. > Is the order of the dofs the same as that of the faces. In that case, I think > I can loop over the faces to evaluate the shape function at a point on each > surface for the purpose of finding the direction. Yes, same order as the order of the faces of the cell. 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/67370078-f7bb-64b9-53f8-9f31d73db260%40colostate.edu.
