Can someone help me how to write the second term of the bilinear form??
Also, given that theta= 0 on boundary of the domain.
The problem is if w is not constant a boundary edge the set \partial E_{-}^w is a full edge.
other wise you can add term like :
([w1,w2]'*N < 0)
Dear professor please check if the following FreeFem++ code is correct code of the second term of the bilinear form according to your suggesion:
Blockquote
+intalledges(Th)((abs(N.xmean(v1old)+N.ymean(v2old))jump(v1)jump(vv1))/nTonEdge(nTonEdge-1))
+int1d(Th)(abs(N.xv1old+N.yv2old)v1vv1)
+intalledges(Th)((abs(N.xmean(v1old)+N.ymean(v2old))jump(v2)jump(vv2))/nTonEdge(nTonEdge-1))
+int1d(Th)(abs(N.xv1old+N.yv2old)v2vv2)
Blockquote
Here, v is known ,v=(v1old,v2old), z=(v1,v2) is unkown, theta=(vv1,vv2) is test function.