Thank you Dr. Wolfgang! I'll try to implement a patch if I have time for it.
Regards, Mohamad On Mon, Apr 8, 2024 at 9:47 PM Wolfgang Bangerth <bange...@colostate.edu> wrote: > On 4/6/24 16:24, Mohamad Ghadban wrote: > > * > > * > > But when I changed the first FESystem into: > > *FESystem[FE_Q(1)^2-**FE_Q(1)^1**]* > > it worked fine because the number of base elements in both FESystems is > now > > the same. My question is then, is this assertion necessary? > > Julian's suggestion works, but if you have two variables that are defined > everywhere and one that is defined in only part of the domain, then what > you > are suggesting (working with only three components) makes sense to me. > > You've already found how to make that work without getting the assertion. > As > to the question of whether the assertion is necessary: Perhaps not, but > the > function in question (FESystem::compare_for_domination()) might be > substantially more complicated to implement. The way it is implemented is > to > assume that the two elements you are comparing have the same number of > base > elements and multiplicities, and in that case determining which element > dominates is a question of determining which base element dominates. If > you > wanted to remove the assertion (patches as always welcome!) you have to > break > things into smaller components and compare those for domination. Here, if > you have > FESystem[FE_Q(1)^3] > and > FESystem[FE_Q(1)^2-FE_Nothing^1] > that would not be too difficult. I bet that could be done in 25 extra > lines of > code (again: patches are very much welcome!) but the difficult case is you > have > FESystem[FE_Q(1)^3] > and > FESystem[FESystem[FE_Q(1),FE_Q(1)]-FE_Nothing^1] > for example. > > Best > Wolfgang > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: bange...@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 a topic in the > Google Groups "deal.II User Group" group. > To unsubscribe from this topic, visit > https://groups.google.com/d/topic/dealii/csAJHjYuK88/unsubscribe. > To unsubscribe from this group and all its topics, send an email to > dealii+unsubscr...@googlegroups.com. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/32538aef-9787-40ea-af0e-f230177e0c45%40colostate.edu > . > -- 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/CAP2_OwGsNMt83Gq6on9nifmyVfjMpct%3DiiuzvGq722Kc78tuLg%40mail.gmail.com.
