On 9/22/21 2:37 AM, Sebastian Kinnewig wrote:
When I implement the computation of the embedding matrix , as done in
FE_Nedelec, in FE_NedelecSZ the program compiles and runs without any
warnings or errors (even in debug mode).
But from the results it appears that there still is a problem with the
hanging nodes.
Alexander knows this element much better than I do, so I have little to
add. The current way we deal with constraints lives by the assumption
that on every cell, the shape functions are mapped from a reference cell
in the same way. I believe that that is not actually the case for the SZ
element because, if I recall correctly, the sign (vector direction) of
shape functions depends on the node numbers. A check for you might be to
see whether at hanging nodes the *magnitude* of the solution is correct
but the *sign* is wrong.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
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 [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/e1cc5e10-0269-72cb-7352-1076c94c100f%40colostate.edu.
