Interpolation matrix with non matching meshes

Linear Algebra
21

msh_0.geometry.dofmap.links(i) returns vertex ids of the ith cell, e.g. it returns 3 integers for mesh with triangle elements. To use that with the following

cell_geom_nodes = msh_0.geometry.dofmap.links(i)
cell_geom_coords = msh_0.geometry.x[cell_geom_nodes]

In v0.8.0, you can do

cell_geom_nodes = msh_0.geometry.dofmap[i]
cell_geom_coords = msh_0.geometry.x[cell_geom_nodes]
1 Like
22

I also want to add that there’s a unit test for interpolation that includes nonmatching meshes, e.g. at dolfinx/python/test/unit/fem/test_interpolation.py at main · FEniCS/dolfinx · GitHub

1 Like
23

I am sorry, but at the moment I am still working with version 0.6.0 as well. I have built an entire set of tools for vibroacoustics with it, and my daily work relies on it. I will port all the scripts to the most recent version, but I cannot promise when I will do that.

1 Like
24

You can also have a look at the following post, where strategy and functionality are used to create a interpolation matrix in parallel.