Dear Praveen, I think FEInterfaceValues is not implemented for FESystem. In my case, it actually just returned zeros without throwing an assertion. It is however more than two years ago I experimented with it.
An alternative way to implement a face_worker without FEInterfaceValues can be found in Step 12b, i.e. on https://github.com/dealii/dealii/issues/15839. Notice, it is probably best to focus on the face_worker only, because the tutorial seem to contain outdated code, see https://github.com/dealii/dealii/issues/15839 Regards, Nils Praveen C <prav...@math.tifrbng.res.in> writes: > Hello > > In a continuous Galerkin method, I want to assemble terms involving jump in > solution gradient and shape function gradient. > > For this, FEInterfacesValues seems like a good choice. But the function > > get_jump_in_function_gradients > > https://www.dealii.org/current/doxygen/deal.II/classFEInterfaceValues.html#aef23f7eee71cc8cbf22ca2cfde9dcdd9 > > seems to work only for scalar FE, while I have an FESystem. > > Is this not implemented for vector FE, and what other options are there to > assemble such terms ? > > best > praveen -- Nils Schween PhD Student Phone: +49 6221 516 557 Mail: nils.schw...@mpi-hd.mpg.de PGP-Key: 4DD3DCC0532EE96DB0C1F8B5368DBFA14CB81849 Max Planck Institute for Nuclear Physics Astrophysical Plasma Theory (APT) Saupfercheckweg 1, D-69117 Heidelberg https://www.mpi-hd.mpg.de/mpi/en/research/scientific-divisions-and-groups/independent-research-groups/apt -- 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/87bk5xx5zh.fsf%40mpi-hd.mpg.de.
smime.p7s
Description: S/MIME cryptographic signature
