Hello, dear participants. I solve the problem of elasticity in 2D using varf. I get the solution as a vector [u1,u2]. Is it possible to get directly the arrays u1 and u2 separately using some command in FreeFem++?

Look at the manual or search on the forum, e.g., A small puzzle about the length of vector valued fespace

The simple way is just to copy your solution u1 u2 in an other fe function w1 and w2

Where

fespace Wh(Th,P2);

Wh w2=u1,w2=u2;

Now w1 is the vector with only u1 composante.