On 08/04/2017 07:07 AM, [email protected] wrote:
Yes, that makes sense to me. But now if we have neighboring faces with
differrent refinement levels (e.g. the face element on one side is only half
as long as the one on the other side) the situation gets a little more
complicated, right? Since then you cannot just copy the coefficients 1 to 1,
or do I get that wrong?
Yes. But then make_hanging_node_constraints() knows what to do with this. The
hp_*_identities() functions aren't used for this.
So what you suggest is that one could actually go on and implement
hp_*_identities for FE_FaceP, in the way you describe? Do you know of some
other element that already implements this, not using interpolation, which one
could use as reference?
I don't know for sure if that element implements the hp constraints, but the
modal RT element has a similar structure in that the degrees of freedom are
located on the faces and are defined via modes or integrals, rather than
interpolation. (There is also an interpolatory RT element; both of them are
available in deal.II.)
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].
For more options, visit https://groups.google.com/d/optout.
