Hello, I am trying to implement and HDG like method from [1], but without the elimination process described in step 51. Thereby creating a block matrix.
To implement this I use as finite element a system of FE_DGQ + FE_FaceQ. On the face I need to compute integrals between the functions from the trace space (i.e., from FE_FaceQ) and the normal times gradient of the functions inside the element (i.e., from FE_DGQ). Therefore, I pass update_gradients to FE_FaceValues, but this runs into the exception that the gradients are not implemented for FE_FaceQ. Note that the lack of gradients of FE_FaceQ is no problem, as its gradients are not need and (I expect) not well defined, but I do need the gradients of FE_DGQ. Is there a way around this without eliminating the degrees of freedom for the elements (like in step 51)? Regards, Lars Corbijn [1] https://epubs.siam.org/doi/abs/10.1137/090775464 -- 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/5f675d13-d851-47f5-8beb-1ce30a125804%40googlegroups.com.
