QuadratureElement crashed when using a P2 element

[Garth Wells]

The code

    element = FiniteElement("Lagrange", triangle, 2)
    QE = FiniteElement("Quadrature", triangle, 3)

    v = TestFunction(element)
    u = TrialFunction(element)
    u0 = Function(element)
    C = Function(QE)
    a = v.dx(i)*C*u.dx(i)*dx + v.dx(i)*2*u0*u*u0.dx(i)*dx

crashes with the message

    --- CHECK 1 ---
    Here I think there are four points (on the cell)
    self.__num_quad_points = 4
    ---------------

    Compiler stage 2: Computing form representation(s)
    --------------------------------------------------
      Computing quadrature representation.

    --- CHECK 3 ---
    This is where the check fails
    len(points) = 16
    self.__num_quad_points = 4
    ---------------

    Points must be equal to coordinates of quadrature points

    *** FFC: Points must be equal to coordinates of quadrature points
    *** FFC: To get more information about this error, rerun FFC with --debug.

I doesn't seem to work for any QuadratureFunction order.