Failure to assemble surface integrals of a FacetElement

[pbrubeck]

I have this code:

from firedrake import *

mesh = ExtrudedMesh(UnitSquareMesh(1, 1, quadrilateral=True), 1)
element = FiniteElement("Lagrange", cell=mesh.ufl_cell(), degree=3)
element = FacetElement(element)
V = FunctionSpace(mesh, element)
u = Function(V)

assemble(u*ds_v)

which produces the following error: NotImplementedError: No rule for factorising expressions of this kind.

The code just works if I change ds_v to ds_t or ds_b

[pbrubeck]

The bug is also present in 2D:

from firedrake import *
mesh = UnitSquareMesh(1, 1, quadrilateral=True)
element = FiniteElement("Lagrange", cell=mesh.ufl_cell(), degree=3)
element = FacetElement(element)
V = FunctionSpace(mesh, element)
u = Function(V)
assemble(u*ds)

[ksagiyam]

I think https://github.com/firedrakeproject/tsfc/pull/282 will fix this.

[pbrubeck]

I think firedrakeproject/tsfc#282 will fix this.

@ksagiyam I just tested on your branch, and it seems to work, thanks!