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

Alternatively, it occurred to me I could probably just use DMCreateDefaultSF() to create the point SF on the new DM- would there be any disadvantage in doing it that way?

However, at the moment I can't, because it appears the Fortran interface for DMCreateDefaultSF() is missing. I would have thought the interface for that should be simple enough, as it just takes a DM and a couple of PetscSections.

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? - 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