On 25/09/17 13:37, Adrian Croucher wrote:

Actually it looks like DMGetDefaultSF() does already work in Fortran though, and could be used to do the same thing.

Do you think that would be a reasonable way to do it?

Hmm, I just tried it and it looks like the default SF produced by this function doesn't have any overlap. So that won't work for me. Guess I might need a Fortran interface for PetscSFSetGraph/PetscSFGetGraph after all. If you could put one it that would be great. Another simpler possibility that would work for my case is if there was a DMCopySF() function for copying the point SF from one DM to another. Up to you if you think that would be of enough general usefulness to warrant inclusion. - Adrian -- Dr Adrian Croucher Senior Research Fellow Department of Engineering Science University of Auckland, New Zealand email: [email protected] tel: +64 (0)9 923 4611