bork.py

from firedrake import *

mesh = UnitTriangleMesh()

V = FunctionSpace(mesh, 'RT', 1)

u = TrialFunction(V)
v = TestFunction(V)
n = FacetNormal(mesh)
mu = Constant(1.0)

muM = as_tensor([[mu, 0], [0, mu]])

a = mu*(-inner(outer(v,n),grad(u)) -inner(outer(u,n),grad(v)))*ds
b = (-inner(outer(v, n), dot(grad(u), muM)) - inner(outer(u, n), dot(grad(v), muM)))*ds


print "b"
print assemble(b).M.values
print "a"
print assemble(a).M.values