I would like to create a vector valued finite element space. I know you can do VectorH1, but what if the number of components is not equal to the spatial dimension?
For example, suppose I want an H1 vector-valued function (defined on a 3-D mesh), with 5 components. How should I do it? I guess I could do this:
V = H1(mesh,order=k, dirichlet="whatever")
X = FESpace([V,V,V,V,V])
Would that work?
Also, is there a tensor-valued finite element space option? Say for symmetric matrices? That are trace-free?