You have a quantity named nTonEdge
assciated to each edge, that contains the number of triangles touching the edge. Its value is thus 2 for internal edge, and 1 for boundary edge.
You can just then mutiply your expression in intalledges
by (nTonEdge-1)
, this will give zero for the boundary terms, and 1 for internal terms. For example
+intalledges(Th)(mean(u)*jump(v)*(nTonEdge-1))
3 Likes