Hello,
how can I change the order of integration for internal facets integrals? The following does not work
dS2 = dS({"integration_order":1}) print assemble(a_uh('-')*v_h('-')*dS2).array()
thanks a lot, Fotini
Hi, try instead
dS1 = dS(metadata={'quadrature_degree': 1})