Quoting Fredrik Bengzon <fredrik.bengzon@xxxxxxxxxxx>:
Hello
I'm wondering if there is some simple way (i.e., other than component
wise) of writing the term
integral ( [ outer(n, v) ] : < grad u > ) dS
If I do the following:
element = FiniteElement("BDM", "triangle", 1)
v = TrialFunction(element)
u = TestFunction(element)
n = FacetNormal("triangle")
print "jump", jump(outer(n,v))
print "avg", avg(grad(u))
a = dot(jump(outer(n,v)), avg(grad(u)))*dS
print "a", a
at least the output (the print statements) looks OK as far as I can tell, but
have a look. The above form will however produce an error in simplify.py. Maybe
you can have a look and see what is going wrong?
