Thanks for the reply. I tested the change and it did not work, as get_generalized_face_support_points() is also empty. The cause seems that FESystem doesn't initialize the generalized_face_support_points. To test that I added some code to initialize it when initialize_unit_face_support_points() 'fails', to work similar but to choose between unit_face_support_points and generalized_face_support_points.
Then my test with step-20 seems to not crash anymore, whether the normal flux enforcement works I did not test. Best, Lars On Mon, Feb 20, 2017 at 6:15 PM, Wolfgang Bangerth <bange...@colostate.edu> wrote: > > After a bit of digging, I think the following happens. The basis for the >> problem is a system of RaviarThomas + DGQ. The RaviarThomas element has no >> support points, but has dofs_per_face. Therefore the same applies to the >> FESystem, which therefore triggers the above mentioned error. >> > > Interesting. The RT element has *generalized* support points, though. What > happens if you change the line in vector_tools.templates.h that calls > unit_support_points = fe_collection[i].get_unit_face_support_points(); > to > unit_support_points = fe_collection[i].get_generaliz > ed_face_support_points(); > ? > > 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/fo > rum/dealii?hl=en > --- You received this message because you are subscribed to a topic in the > Google Groups "deal.II User Group" group. > To unsubscribe from this topic, visit https://groups.google.com/d/to > pic/dealii/jFjGBuEhdXQ/unsubscribe. > To unsubscribe from this group and all its topics, send an email to > dealii+unsubscr...@googlegroups.com. > For more options, visit https://groups.google.com/d/optout. > -- 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.
