dolfin team mailing list archive

["Garth N. Wells" <gnw20@xxxxxxxxx>]

We have an issue with MixedFunctionSpace which makes incorrect usage
very easy. For example, if I have 

    V = VectorFunctionSpace(mesh, "CG", 2)
    Q = FunctionSpace(mesh, "CG", 1)
    W = V * Q

The dof maps for Q and W.sub(1) are not the same. I've seen a number of
users using V or Q in setting boundary conditions rather than W.sub(0)
and W.sub(1).

When using 

    W = MixedFunctionSpace([ V, Q])

could we create W and then inside the MixedFunctionSpace constructor set
 
   V = W.sub(0)  
   Q = W.sub(1)

?

Garth