Error with boundary integral and VectorFunctionSpace

[marazzaf]

Hello,

The following code gets me a PETSc Error 77 whereas the second does not. I do not get why.

from firedrake import *

N = 5
mesh = UnitSquareMesh(N, N)

V = VectorFunctionSpace(mesh, 'CG', 1)
u = TrialFunction(V)
v = TestFunction(V)
a = inner(u, v) * ds(1)
A = assemble(a)

from firedrake import *

N = 5
mesh = UnitSquareMesh(N, N)

V = FunctionSpace(mesh, 'CG', 1)
u = TrialFunction(V)
v = TestFunction(V)
a = inner(u, v) * ds(1)
A = assemble(a)

Any idea?

[connorjward]

This seems to be a bug. I've created an issue to track it.