← Back to team overview

dolfin team mailing list archive

compile form with normal

 

Hello.

I realize that there has been a big change lately with FFC, but will it be possible in the near future to compile the following form with FFC:

---------------------------

element = FiniteElement("Lagrange", "triangle", 1)

v = TestFunction(element)
u = TrialFunction(element)
f = Function(element)

n = FacetNormal("triangle")

a = v*u*dx + dot(grad(v), grad(u))*dx + \
(D(v, 0) * (-n[1]) + D(v, 1) * n[0])*(D(u, 0) * (-n[1]) + D(u, 1) * n[0])*ds

L = v*f*dx

----------------------------

Was it possible to do this with the old FFC?

- Shawn


References