I would like to define a VectorH1 space with different Dirichlet BC conditions for different components, currently, I am only doing the following, and then I have to define grad/div manually, which is cumbersome…
V1 = H1(mesh, dirichlet="top")
V2 = H1(mesh, dirichelt="bottom")
V = FESpace([V1, V2])
Is there a way to use VectorH1 to define the above space?