1mesh Th=square(10,10);
2fespace Vh(Th,P1);
3// bug due  to optimisation process .. HARD ....  FH
4Vh u = y ? 1./y : x;