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)