Thank you for the github link. In the end I fixed the issue by checking 
which value of the dof_index vector was not equal to 
numbers::invalid_unsigned_int. I added this line :

*unsigned int component_i = (dof_index[0] != numbers::invalid_unsigned_int) ? *
*                                    
fe_RTE.system_to_component_index(dof_index[0]).first:*
*                                    
fe_RTE.system_to_component_index(dof_index[1]).first;*

Thank you Corbin as well. I will try that next time I am in need of it !

Best,

Sylvain M.
El viernes, 23 de julio de 2021 a la(s) 20:58:25 UTC+2, Wolfgang Bangerth 
escribió:

> On 7/23/21 9:41 AM, Sylvain Mathonnière wrote:
> > 
> > As for the 4294967295 occurence, is it because I am dealing with an 
> interface 
> > which is a boundary as well ? Isn't this case handled by the 
> *boundary_worker* 
> > and not the *face_worker* (using step-12 notation) ?
>
> No, it's because you are using DG elements. I tried to clarify this here:
> https://github.com/dealii/dealii/pull/12595
>
> Best
> W.
>
> -- 
> ------------------------------------------------------------------------
> Wolfgang Bangerth email: [email protected]
> 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].
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/83b84ffd-3a84-43c6-b918-bb9cc35705bfn%40googlegroups.com.

Reply via email to