If you want to do volume/surface coupling, you should not use emptymesh, but extract instead. See this example FreeFem-sources/helmholtz-coupled-2d-PETSc-complex.edp at master · FreeFem/FreeFem-sources · GitHub, and in particular the restriction from volume to surface matrix, RVtoS and the ASurf distributed Mat.
1 Like