Re: [petsc-dev] tiny bug in DMPlexDistribute

2017-10-13 Thread Mark Adams
Thanks Vaclav, I've been seeing a valgrind warning in PlexDestroy for a long time. Mark On Fri, Oct 13, 2017 at 10:33 AM, Matthew Knepley wrote: > On Fri, Oct 13, 2017 at 10:26 AM, Vaclav Hapla > wrote: > >> Hello >> >> In DMPlexDistribute, when it is run on 1 process and sf != NULL, the >> out

Re: [petsc-dev] tiny bug in DMPlexDistribute

2017-10-13 Thread Matthew Knepley
On Fri, Oct 13, 2017 at 10:26 AM, Vaclav Hapla wrote: > Hello > > In DMPlexDistribute, when it is run on 1 process and sf != NULL, the > output *sf is undefined. This in turn leads to crash of PetscSFDestroy, > e.g. in src/dm/impls/plex/examples/tutorials/ex5.c:45. > > I created the corresponding

[petsc-dev] tiny bug in DMPlexDistribute

2017-10-13 Thread Vaclav Hapla
Hello In DMPlexDistribute, when it is run on 1 process and sf != NULL, the output *sf is undefined. This in turn leads to crash of PetscSFDestroy, e.g. in src/dm/impls/plex/examples/tutorials/ex5.c:45. I created the corresponding pull request a minute ago. It is my first PR - please, should I

Re: [petsc-dev] PetscSF in Fortran

2017-10-13 Thread Matthew Knepley
On Fri, Oct 13, 2017 at 12:25 AM, Adrian Croucher wrote: > hi Jed > > I have had a go at writing a Fortran binding for the PetscSFGetGraph() > function, along the lines of what you suggested: > > On 21/09/17 15:15, Jed Brown wrote: > Unfortunately PetscSFSetGraph/PetscSFGetGraph need custom bindi