Dear all, I'm using an hp::FECollection object to describe different FE spaces in my grid. In particular, I'm in the following scenario:
x- - - - - -x | | | 0 | | | x - - F - -x | | | 1 | | | x- - - - - -x On cell number 0 I have a FESystem with FE_Q(1) and FE_Q(1), which corresponds to fe_collection[0]. On cell number 1, I have again a FESystem with FE_Q(1) and FE_Nothing, corresponding to fe_collection[1]. I need to compute the jump [d_n u_h]_F, being F the face shared by the two cells. Unfortunately, I can't find a proper way to initialize FEInterfaceValues so that it understands that I need the jump of the first components only. Indeed, if I use fe_collection[0] (or fe_collection[1]) as first argument in the constructor, I have a runtime error in reinit(cell) saying (as expected): * The FiniteElement you provided to FEValues and the FiniteElement that belongs to the DoFHandler that provided the cell iterator do not match.* since we have different spaces on cell0 and cell1. Is there a possible workaround to this? Best, Marco -- 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]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/2d71b802-132c-4f8a-a20d-03542939036an%40googlegroups.com.
