Since FFC is about to move to UFL, I guess this is the best place to ask
certain questions now ;).
Will UFL support the absolute value of sums, e.g.
vector = VectorElement("Lagrange", "tetrahedron", 1)
scalar = Element("Lagrange", "tetrahedron", 1)
v = TestFunction(scalar)
u = Function(vector)
n = Function(vector)
L = v*abs(inner(u, n))*ds
This would be useful for upwinding fluxes.
Garth
_______________________________________________
UFL-dev mailing list
UFL-dev@xxxxxxxxxx
http://fenics.org/mailman/listinfo/ufl-dev