Hello everbody, I am quite new to field and I face the problem that the ALE.move command produces unexpected results in case of free meshes. The 1st problems looks fine:
from dolfin import *
from mshr import *
mesh = UnitSquareMesh(4, 4)
bmesh = RectangleMesh(Point(1, 1), Point(3, 3), 4, 4)
bmesh = BoundaryMesh(bmesh, 'exterior')
u = ALE.move(mesh, bmesh)
which means, that the initial mesh is moved in order to match the new boundaries. But for the 2nd case:
from dolfin import *
from mshr import *
mesh = UnitSquareMesh(4, 4)
geom = Rectangle(Point(1, 1), Point(3, 3))
bmesh = generate_mesh(geom,4)
bmesh = BoundaryMesh(bmesh, 'exterior')
u = ALE.move(mesh, bmesh)
the resulting mesh looks somehow different but that case might sometimes the more practically relevant one.
I would be very glad if you have any suggestions.
Best Regards
