UFL input for the Poisson bilinear formΒΆ

The bilinear form for a stiffness matrix (Poisson):

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

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

a = dot(grad(v), grad(u))*dx