ValueError: matrices are not aligned

The error

ValueError: matrices are not aligned

is produced wehn using the vertex quadrature rule.

MWE:

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}
              )

Comments (8)