On 08/18/2017 06:19 AM, Eldar Khattatov wrote:
It tries to compare FE_RaviartThomas. This element indeed does not implement
compare_for_face_domination() function,
Correct. Though you could think about implementing it.
so I tried using FE_RaviartThomasNodal
instead, as it provides an implementation for this function. Unfortunately,
from what I understood from the implementation in RT_Nodal, it works only if
another element in the comparison is also of the same type RT_Nodal, otherwise
it has same Assert(false, ...) before return statement,
Yes, indeed. But it shouldn't be very difficult to extend it for the
particular case you are considering: that the other element is an FE_Nothing.
The FE_Nothing describes a field that is constant zero, and so the appropriate
interface condition is that the FE_RT_Nodal should also be zero at that
interface, or be left unconstrained, depending on how you want the FE_Nothing
to be interpreted. That should be easy to implement.
Take a look at how the FE_Q element does this, in one of the classes that sit
below class FiniteElement and class FE_Q in the class inheritance diagram. I
think you'll see pretty quickly how you want to implement things for your case.
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.
