hi On 21/09/17 15:15, Jed Brown wrote:

Unfortunately PetscSFSetGraph/PetscSFGetGraph need custom bindings for Fortran and they have not been written. Shouldn't be difficult, but will take a little work/testing. I would probably make the array of PetscSFNode be a PetscInt array of length 2n (or of shape (2,n)) -- that's an equivalent memory representation to what we're using now.

That sounds like the sort of thing I need. 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. - Adrian

Adrian Croucher <[email protected]> writes:

...

hi

For the dual-porosity stuff I'm working on, I'm creating a modified DMPlex. Looking at e.g. TS ex11.c (in the SplitFaces() routine), it appears I am going to have to set up a PetscSF for the new DM (from what I can gather, that is needed for managing parallel communication of partition ghost cell values).

The way I've set up the modified DMPlex, the partition ghost cells are in the same locations as they were in the original single-porosity DM. So I'm thinking I basically just have to copy the SF straight over from the original DM.

The routine DMGetPointSF() seems to work OK in Fortran, but I'm having a bit of trouble with PetscSFGetGraph().

When I try to declare an array of type PetscSFNode for the iremote parameter, it doesn't seem to know about PetscSFNode, even though I have a 'use petsc' in my code, which normally gives access to everything.

Also I'm not sure of exactly how to set up the ilocal and iremote array parameters for this function in Fortran. How should they be declared- as pointer arrays?

- 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