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_points()
'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/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.
