Mpirun performs mesh generation

[dokken]

The example above moves the internal nodes based on the data in a dolfinx.fem.Function, which I assume is what you will have as input to your own problem. It doesn’t have any assumptions that would keep the exterior boundary fixed, that is just due to the specifics of the prescribed functions expr and expr2

@dokken

def expr(x):
        move_1 = x[1] * (0.5 - x[1]) * (x[1] < 0.5)
        move_2 = (x[1] - 1) * (x[1] - 0.5) * (x[1] > 0.5)
        return (np.zeros(x.shape[1]), move_1 + move_2)

    u = dolfinx.fem.Function(V)
    u.interpolate(expr)
    perturb_mesh(u)

    # Create the dof -> vertex map and vertex->node map
    with dolfinx.io.XDMFFile(mesh.comm, "moved_mesh.xdmf", "w") as xdmf:
        xdmf.write_mesh(mesh)

    # Reset mesh
    u.x.array[:] *= -1
    perturb_mesh(u)

    def expr2(x, eps=1e-8):
        xi = np.floor(x[0] * N)
        yi = np.floor(x[1] * N)
        move_x = (
            0.2
            / N
            * np.cos(10 * np.pi * x[0])
            * (x[1] > eps)
            * (x[1] < 1 - eps)
            * (-1) ** yi
        )
        move_y = (
            0.2
            / N
            * np.cos(10 * np.pi * x[1])
            * (x[0] > eps)
            * (x[0] < 1 - eps)
            * (-1) ** xi
        )
        return (move_y, move_x)

If u is a solution from a elasticity probloem, you can simply pass that to perturb_mesh(u)