Mixed quadrature fails when number of points for different degrees coincide

After the fix in martinal/topic-fix-mixed-quadrature a new issue pops up:

from dolfin import *
parameters['form_compiler']['representation'] = 'uflacs'

mesh = UnitSquareMesh(1, 1)
V = FunctionSpace(mesh, 'CG', 3)
f = interpolate(Expression("x[0]*x[0]*x[0]"), V)

M = [f*dx(degree=i) for i in range(4)]
print assemble(M[3] - M[2])
print assemble(M[1] - M[0])

Comments (3)