Hello,
I have a very basic question about facet integrals. For example, in 2D
do you use a gauss quadrature in 1D or in 2D?
For example, with the following 2 lines
wh = Function(DG0)
a = assemble(wh('+')*dS)
what I compute?
When I use the metadata
dS1 = dS(metadata={'quadrature_degree': 2})
and then I compute
a = assemble(wh('+')*dS1)
how I can print the quadrature points and the weights of the quadrature rule?
Thank you in advance,
Fotini
