<div dir="ltr">Hi Matthew and Lawrence,<div><br></div><div>Thank you for your responses. The second approach involving DMPlexDistributeOverlap is viable. However there are two issues: 1) I have lots of points to interpolate at and 2) The points are not always in the next nearest neighbor process but can extend beyond, so essentially creating an overlap that extends beyond the nearest neighbor process may not be good idea?</div><div><br></div><div>Then the earlier approach which I thought was to extract and store the mesh locally and create local communicator to communicate the fields needed for interpolation. This is a bit of an effort to code, but may be more efficient. This step is also a critical part of my algorithm, so I am trying to code it efficiently. How can a region of dmplex be extracted and stored locally?</div><div><br></div><div>Inputs and ideas are always welcome.</div><div><br></div><div>Sincerely,</div><div>SG</div></div>