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?