origin of 2D-meshes appear at [0,0] instead of at the chosen coordinate

[VolkerRaab]

When I execute the following script:

mesh Th = square(100,100, [10*(x-1/2), 20*(y-1/2)]);
fespace F(Th, P1);
F sinsin = sin(x)*sin(y);
plot(sinsin);

I get a tiling with a knot in the bottom left corner instead of in the center of the plot. The mesh seems to respect the different extensions along x and y but not the shifted origin. Is this by design?

[prj-]

Your mesh definition is wrong, you are performing integer divisions (whereas you want floating-point divisions).

mesh Th = square(100,100, [10.0*(x-1.0/2.0), 20.0*(y-1.0/2.0)]);

[VolkerRaab]

Dear Pierre, thanks a lot for the quick reply 👍 and sorry for that somewhat stupid error 🥴.