Fespaces in vectorized formulation

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?

It work and x y are the coordinate of the quadrature point use to compute the integral.

1 Like