mesh.mf_map[(a,b)] will contain a mesh function with the intersection of mesh with mesh (id=a) and mesh (id=b) and dx_map will be the Measure() used to integrate