Jit compiler error when creating Quadrature FunctionSpace for Quadrilateral Grids

Consider the following short python code

import dolfin as dl

nx = 1
ny = 1
qdegree = 1
mesh = dl.UnitSquareMesh.create(nx, ny, dl.CellType.Type_quadrilateral)
element = dl.FiniteElement("Quadrature", mesh.ufl_cell(), qdegree, quad_scheme="default")
Vh = dl.FunctionSpace(mesh, element)

which leads to the following error when trying to create the FunctionSpace Vh:

** Reason:  ffc.jit failed with message:
Traceback (most recent call last):
  File "/usr/lib/python3/dist-packages/dolfin/compilemodules/jit.py", line 142, in jit
    result = ffc.jit(ufl_object, parameters=p)
  File "/usr/local/lib/python3.5/dist-packages/ffc/jitcompiler.py", line 218, in jit
    module = jit_build(ufl_object, module_name, parameters)
  File "/usr/local/lib/python3.5/dist-packages/ffc/jitcompiler.py", line 134, in jit_build
    generate=jit_generate)
  File "/usr/local/lib/python3.5/dist-packages/dijitso/jit.py", line 167, in jit
    header, source, dependencies = generate(jitable, name, signature, params["generator"])
  File "/usr/local/lib/python3.5/dist-packages/ffc/jitcompiler.py", line 67, in jit_generate
    prefix=module_name, parameters=parameters, jit=True)
  File "/usr/local/lib/python3.5/dist-packages/ffc/compiler.py", line 150, in compile_element
    prefix, parameters, jit)
  File "/usr/local/lib/python3.5/dist-packages/ffc/compiler.py", line 190, in compile_ufl_objects
    ir = compute_ir(analysis, prefix, parameters, jit)
  File "/usr/local/lib/python3.5/dist-packages/ffc/representation.py", line 171, in compute_ir
    for e in elements]
  File "/usr/local/lib/python3.5/dist-packages/ffc/representation.py", line 171, in <listcomp>
    for e in elements]
  File "/usr/local/lib/python3.5/dist-packages/ffc/representation.py", line 226, in _compute_element_ir
    ir["evaluate_basis"] = _evaluate_basis(ufl_element, fiat_element, parameters["epsilon"])
  File "/usr/local/lib/python3.5/dist-packages/ffc/representation.py", line 706, in _evaluate_basis
    e.get_coeffs()
  File "/usr/local/lib/python3.5/dist-packages/FIAT/tensor_product.py", line 434, in get_coeffs
    raise self.element.get_coeffs()
AttributeError: 'QuadratureElement' object has no attribute 'get_coeffs'
.
*** Where:   This error was encountered inside jit.py.
*** Process: 0
***
*** DOLFIN version: 2017.2.0
*** Git changeset:  0baf73825079a581e43ab1705370043040aa213d
*** -------------------------------------------------------------------------

Aborted
$

The same code works fine with triangular grids (dl.CellType.Type_triangle).

Thanks in advance!

Umberto

Comments (1)