Dear all,
I use a change of coordinate to solve my equations, and I then save the solutions in the initial domain.
The transform is \hat{x} = L(x-R)/(L-R).
In 1D, my code to save the data in the initial (not transformed) mesh is:
----------------------------
# mesh on which I solve the equations (\hat{x})
mesh_transformed = mesh.coordinates.vector().get_local()
# back to the initial mesh (applying the inverse transform) (x)
mesh.coordinates.vector().set_local(mesh_transformed*(L-R)/L + R)
# save data on the initial mesh
save_h << h
save_psi << psi
# Back to the transformed mesh to solve the WF (\hat{x})
mesh.coordinates.vector().set_local(mesh_transformed)
----------------------------
This works well in 1D. However, I would like to extend my model in 2D, where the transform is still
only in x. Therefore, I would like to apply the previous manipulation to the x coordinate only, and not on the 2D mesh. The previous code changes both x and y.
How can I change only x ?
Thanks,
Floriane