Hello, I have a discontinuous field U
W = functionspace(mesh, ("DG", 0, (mesh.geometry.dim, )))
U = Function(W)
I wish to integrate over a curve consisting of mesh edges.
This curve bounds a region, which are both tagged, using gmsh physical tags.
I have to decide whether to integrate U("+")
, or U("-")
using the following measure
dS = Measure("dS", domain=mesh, subdomain_data=facet_tags, sudomain_id=curve2integrate)
Is "+"
pointing inward or outward the region bounded by the curve?