Hi,
Is it correct to have different fespaces and solve a vectorized equation as:
fespace Vh(Omega, P2);
fespace Ph(Omega, P1);
Vh wr,wz,urs,uzs;
Ph hs,q;
problem Myproblem([hs,urs,uzs],[q,wr,wz])
=int2d(Omega)((y*[hs,urs,uzs])'*[q,wr,wz]))…
if yes, how do the coordinates x or y work?