Thank you for the answer, it was indeed very easy to modify the RT_Nodal according to your comments. If this patch could be of use to the library/community I can commit my changes to the repo.
Best, Eldar On Friday, August 18, 2017 at 6:23:31 PM UTC-4, Wolfgang Bangerth wrote: > > 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] > <javascript:> > 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.
