Hi Wolfgang, Thanks for the response. I think FE_FaceQ is exactly what I am looking for, at least in the case of representing the normal trace of RT_0. But I have trouble using the FEFieldFunction with the FE_FaceQ element. I tried the following:
*Functions::FEFieldFunction<dim,DoFHandler<dim>,Vector<double>> fe_interface_data(dof_handler_face,solution_vector);* * MappingQGeneric<dim> mapping_generic(1);* * DoFTools::map_dofs_to_support_points(mapping_generic,dof_handler_face,support_points); fe_interface_data.value_list(support_points,values);* And got all values to be NaN. Am I doing something wrong? Sam. On Wed, Oct 16, 2019 at 2:41 PM Wolfgang Bangerth <bange...@colostate.edu> wrote: > On 10/16/19 12:15 PM, Sam wrote: > > I would like to define a FE defined by lagrange interpolants on the > faces of > > the cell. The objective is to setup a FE corresponding to the normal > trace of > > RT(degree=k). > > For example, for RT(degree=0), the dofs are given by face integrals on > faces > > and I need a DGQ finite element which has support points exactly at the > center > > of each face. > > In general how can I construct a DGQ finite element which is continuous > only > > across faces of the cell. > > Before we think about implementing a new element, have you checked whether > the > FE_FaceQ and FEFaceP classes do what you need? > > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: bange...@colostate.edu > 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 dealii+unsubscr...@googlegroups.com. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/6bc97617-70c0-395e-1f36-2fc48b720dbf%40colostate.edu > . > -- 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/CAP8kMBEYZTg3KzhhhetgTTbz%3D_GoVC-7%3D4CWTR_QL6EcAAqeQw%40mail.gmail.com.
