Assembly of rectangular matrices broken

Issue #352 resolved
Jan Blechta created an issue

Assembly of rectangular matrices is broken as pointed out in Q4289. MWE:

from dolfin import *
mesh = UnitSquareMesh(6, 6)
V = VectorFunctionSpace(mesh, "CG", 2)
Q = FunctionSpace(mesh, "CG", 1)
v = TestFunction(V)
p = TrialFunction(Q)
assemble(inner(grad(p), v)*dx)

raises

RuntimeError: *** Error: All dofmaps must have same block size (for now)

Comments (8)

  1. Log in to comment