Dear all, As indicated in the documentation <http://dealii.org/developer/doxygen/deal.II/classFE__DGP.html>, the finite elements FE_DGP and FE_FaceP are not based on nodal interpolation (as eg. FE_DGQ and FE_FaceQ), but on projection. Apparently, this prevents these elements from being used in combination with a hp::DoFHandler, since there seems to be no easy way (or even no way?) to implement the functions hp_*_identities for FE_FaceP (which is needed to support hp) when no nodal interpolation (ie. no unit_support_points) is available.
I would really like to use the FE_DGP/FE_FaceP combination in a hp FEM computation, as this requires fewer DoFs than the FE_DGQ/FE_FaceQ couple. Does anyone have an idea how to fix this problem? I imagine that either: - one could implement hp_*_identities, somehow, based on the existing code with projection, or - one could implement a different version of FE_FaceP with nodal interpolation but i don't really know whether these approaches could work and where to start. Best, Samuel -- 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. For more options, visit https://groups.google.com/d/optout.