FFC generates code which fails to compile for Quadrature elements of degree greater than one

The assembly of an inner product works, but not the derivative.

Fails with uflacs and tsfc. Other backends not tested.

from dolfin import *
mesh = UnitSquareMesh(4,4)

for ii in range(3):
    print('')
    print(ii)
    elmt = FiniteElement('Quadrature', mesh.cell_name(), degree=ii, quad_scheme='default')
    V = FunctionSpace(mesh, elmt)
    u = Function(V)
    ee = inner(u, u) * dx(degree=elmt.degree())
    print(assemble(ee)) # works
    print(assemble(derivative(ee, u, TestFunction(V)))) # fails for ii == 2

Is this use case intended at all?

Comments (0)