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