uflacs failing in 1D

Consider the following code:

from dolfin import *

parameters["form_compiler"]["representation"] = "uflacs"

mesh = IntervalMesh(400, 0, 1)
V = FunctionSpace(mesh, "CG", 1)
u = interpolate(Constant(1), V)
v = interpolate(Constant(2), V)
form = inner(u - v, u - v)*dx

dform = assemble(derivative(form, u)) # crashes

It fails to compile, with

/tmp/tmpZqVvI6/ffc_form_0ee9a0b845857a9f82a6163e1aafc9fd36a8a174.cpp: In member function ‘virtual void ffc_form_0ee9a0b845857a9f82a6163e1aafc9fd36a8a174_cell_integral_main_otherwise::tabulate_tensor(double*, const double* const*, const double*, int) const’:
/tmp/tmpZqVvI6/ffc_form_0ee9a0b845857a9f82a6163e1aafc9fd36a8a174.cpp:82:46: error: ‘FE1_C0_D1_Q2’ was not declared in this scope
     const double J_c0 = coordinate_dofs[0] * FE1_C0_D1_Q2[0][0][0] + coordinate_dofs[1] * FE1_C0_D1_Q2[0][0][1];

Comments (3)