You don’t need a special class, it should be very straightforward for P1
geometry meshes.
def move(mesh, u):
x = mesh.geometry.x
gdim = mesh.geometry.dim
u_x = u.compute_point_values()
x[:,:gdim] += u_x
You don’t need a special class, it should be very straightforward for P1
geometry meshes.
def move(mesh, u):
x = mesh.geometry.x
gdim = mesh.geometry.dim
u_x = u.compute_point_values()
x[:,:gdim] += u_x