dolfinx icon indicating copy to clipboard operation
dolfinx copied to clipboard

Interpolation on a submesh

Open jpdean opened this issue 2 years ago • 0 comments

Interpolation on a mesh where a process has no cells currently causes a segfault. Returning an empty array of interpolation coordinates on processes with no cells (see https://github.com/FEniCS/dolfinx/tree/jpdean/submesh_function_fix) fixes this. However, if a process happens to own a dof that is not in a cell on that process, it's value will not be set on that process during intepolation. A work around is to send ghost values back to their owners, but this seems like a hack. This issue should only affect point evaluation dofs on the mesh vertices.

jpdean avatar Apr 13 '22 09:04 jpdean