I will try to do so somewhere this week, but it depends a bit on how much free evening time I have this week. Furthermore, I don't know what the 'fall-out' is of having to initialize the generalized support points, ass far as I see now I will have to implement and learn about the interpolation functions of FiniteElement.
For the change to VectorTools, I can do it. But I don't have an idea if the code following it does assume in any way that it has normal support points and not generalized ones. If you (or somebody else with more experience) can tell me it doesn't I will also make a patch for it. Lars On Mon, Feb 20, 2017 at 10:56 PM, Wolfgang Bangerth <bange...@colostate.edu> wrote: > On 02/20/2017 07:22 PM, Lars Corbijn van Willenswaard wrote: > >> 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_p >> oints() >> 'fails', to work similar but to choose between unit_face_support_points >> and >> generalized_face_support_points. >> > > Ah, interesting. > > Would you want to submit a patch to fix the initialization of generalized > support points in FESystem? That seems like an obvious oversight that > FESystem doesn't currently do it. > > Let us know if you'd like help in figuring out how to submit a patch! > > > Then my test with step-20 seems to not crash anymore, whether the normal >> flux >> enforcement works I did not test. >> > > That would then be a second bug that it doesn't use the generalized > support points. We'd take a (separate) patch for that as well. > > > 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.
