The FiniteElementDomination logic in the codim=0 case would indeed make up a cheap a priori check in this context.
In case a FE_Nothing has been configured to dominate, the solution should be continuous on the interface if I understood correctly, i.e., zero on the face. I will write a few tests to see if this is actually automatically the case in user applications. If so, this check for domination will help other users to avoid this pitfall :) Marc On Monday, December 28, 2020 at 4:13:40 PM UTC-7 Wolfgang Bangerth wrote: > > > The problem here is that the solution is not continuous across the face > of a > > FE_Q and a FE_Nothing element. If a FE_Nothing is turned into a FE_Q > element, > > the solution is suddenly expected to be continuous, and we have no rule > in > > deal.II yet how to continue in the situation. In my opinion, we should > throw > > an assertion in this case. > > Denis actually thought of that already a while back. It wasn't well > documented, but see here now: > https://github.com/dealii/dealii/pull/11430 > > Does that also address what you wanted to do in your patch? > > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: bang...@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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/5c35ba6e-6deb-42a2-8ef1-aa45bba097dan%40googlegroups.com.
