Mapping between meshes in dolfinx

[dokken]

If the meshes are meshed separately you need to map the function or derived quantity from one mesh to the other, see for instance

Interpolate on boundary

Consider the minimal working example: from mpi4py import MPI import dolfinx import numpy as np L = 1 N = 5 def interface(x): return np.isclose(x[0], L) mesh_1 = dolfinx.mesh.create_box( comm=MPI.COMM_WORLD, points=[(0, 0, 0), (L, L, L)], n=[N, N, N] ) V = dolfinx.fem.functionspace(mesh_1, ("Lagrange", 2)) u = dolfinx.fem.Function(V) u.interpolate(lambda x: x[0] + np.cos(x[1]) + np.sin(np.pi * x[2])) Nx, Ny, Nz = 12, 13, 15 mesh_2 = dolfinx.mesh.create_box( comm=MPI.COMM_WORLD…

which illustrates this.