ValueError: matrices are not aligned

[Nico Schlömer]

When trying to assemble a matrix with the vertex-quadrature rule,

form_compiler_parameters={'quadrature_rule': 'vertex', 'quadrature_degree': 1}

the error

ValueError: matrices are not aligned

is thrown.

Minimal example that reproduces the behavior:

from dolfin import *

n = 30

mesh = UnitIntervalMesh(n)

V = FunctionSpace(mesh, 'CG', 1)

u = TrialFunction(V)
v = TestFunction(V)

a1 = u*v*dx

A1 = assemble(a1,
              form_compiler_parameters={'quadrature_rule': 'vertex',
                                        'quadrature_degree': 1}
              )

[Prof Garth Wells]

The error comes from FIAT. Please test further to test if report should be made under FIAT or FFC.