Hi,
Unfortunately no-one has implemented FEInterfaceValues for hp yet (I think it
would be great if someone did). The only workaround I know is to NOT use
FEInterfaceValues and instead:
1. Use 2 hp::FEFaceValues objects, one initialized with cell 0 and one
initialized with cell 1.
2. Split the jump term in 4 terms: [u][v] = u_0*v_0 - u0*v_1 - u_1*v_0 + u_1*v_1
3. Assemble each term in a separate local matrix (A_00, A_01, A_10, A_11), where
A_01 - Rows correspond to local dofs on cell 0, columns to local dofs on
cell 1.
4. Use global_dof_indices of cell 0 and cell 1 to map each local matrix to the
global matrix.
Best,
Simon
On 09/08/2022 00:58, Marco Feder wrote:
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/
<http://www.dealii.org/>
For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en
<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]
<mailto:[email protected]>.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/2d71b802-132c-4f8a-a20d-03542939036an%40googlegroups.com
<https://groups.google.com/d/msgid/dealii/2d71b802-132c-4f8a-a20d-03542939036an%40googlegroups.com?utm_medium=email&utm_source=footer>.
--
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/23e81ebd-cf4e-60ad-ea66-ad825581b533%40gmail.com.
