I can get around this error by first creating a bubble-enriched scalar element, then combining several copies of it into a MixedElement
(which is equivalent to a VectorElement
):
from dolfin import *
mesh = UnitSquareMesh(200,200)
#P1 = VectorElement("P",triangle,1)
#B1 = VectorElement("B",triangle,3)
#V = VectorFunctionSpace(mesh,P1+B1)
P1_scalar = FiniteElement("P",triangle,1)
B1_scalar = FiniteElement("B",triangle,3)
VE_scalar = P1_scalar + B1_scalarb
VE = MixedElement(2*[VE_scalar])
V = FunctionSpace(mesh,VE)