Re: [petsc-users] DMPlexDistributeField

2019-07-10 Thread Matthew Knepley via petsc-users
Crap! The wrong thing got pushed. Thanks, Matt On Wed, Jul 10, 2019 at 9:49 PM Adrian Croucher wrote: > Probably don't want the two extra PetscSFView() calls (lines 1687, 1688) > though- presumably they were just for temporary debugging? > > - Adrian > On 11/07/19 2:18 PM, Adrian Crouche

Re: [petsc-users] DMPlexDistributeField

2019-07-10 Thread Adrian Croucher via petsc-users
Probably don't want the two extra PetscSFView() calls (lines 1687, 1688) though- presumably they were just for temporary debugging? - Adrian On 11/07/19 2:18 PM, Adrian Croucher wrote: On 10/07/19 3:40 PM, Matthew Knepley wrote: Sorry this took so long. I put in a PR for this: https://bit

Re: [petsc-users] DMPlexDistributeField

2019-07-10 Thread Adrian Croucher via petsc-users
On 10/07/19 3:40 PM, Matthew Knepley wrote: Sorry this took so long. I put in a PR for this: https://bitbucket.org/petsc/petsc/pull-requests/1858/knepley-fix-plex-distribute-overlap/diff I think it fixes your problem. Thanks very much Matt, that does fix it. Good stuff! - Adrian -- Dr Adr

Re: [petsc-users] DMPlexDistributeField

2019-07-07 Thread Adrian Croucher via petsc-users
hi Matt, On 29/06/19 12:07 AM, Matthew Knepley wrote: Okay, I ran through this. I am attaching my C version which was easier for me to play with in the debugger, but it should not be hard to put the few lines into yours. The problem comes from using overlap in the redistribution. You can see

Re: [petsc-users] DMPlexDistributeField

2019-06-27 Thread Adrian Croucher via petsc-users
On 28/06/19 10:09 AM, Zhang, Junchao wrote: Check how the graph is created and then whether the parameters to PetscSFSetGraph() are correct. Yes, unfortunately I don't have a good enough understanding of how DMPlexDistribute() works to see what the problem is. - Adrian -- Dr Adrian Crou

Re: [petsc-users] DMPlexDistributeField

2019-06-27 Thread Zhang, Junchao via petsc-users
On Thu, Jun 27, 2019 at 4:50 PM Adrian Croucher mailto:a.crouc...@auckland.ac.nz>> wrote: hi On 28/06/19 3:14 AM, Zhang, Junchao wrote: > You can dump relevant SFs to make sure their graph is correct. Yes, I'm doing that, and the graphs don't look correct. Check how the graph is created and t

Re: [petsc-users] DMPlexDistributeField

2019-06-27 Thread Adrian Croucher via petsc-users
hi On 28/06/19 3:14 AM, Zhang, Junchao wrote: You can dump relevant SFs to make sure their graph is correct. Yes, I'm doing that, and the graphs don't look correct. - Adrian -- Dr Adrian Croucher Senior Research Fellow Department of Engineering Science University of Auckland, New Zealand em

Re: [petsc-users] DMPlexDistributeField

2019-06-27 Thread Zhang, Junchao via petsc-users
On Wed, Jun 26, 2019 at 11:12 PM Adrian Croucher mailto:a.crouc...@auckland.ac.nz>> wrote: hi On 27/06/19 4:07 PM, Zhang, Junchao wrote: Adrian, I am working on SF but know nothing about DMPlexDistributeField. Do you think SF creation or communication is wrong? If yes, I'd like to know the

Re: [petsc-users] DMPlexDistributeField

2019-06-26 Thread Adrian Croucher via petsc-users
hi On 27/06/19 4:07 PM, Zhang, Junchao wrote:  Adrian, I am working on SF but know nothing about DMPlexDistributeField. Do you think SF creation or communication is wrong? If yes, I'd like to know the detail.  I have a branch jczhang/sf-more-opts, which adds some optimizations to SF.  It pr

Re: [petsc-users] DMPlexDistributeField

2019-06-26 Thread Zhang, Junchao via petsc-users
On Mon, Jun 24, 2019 at 6:23 PM Adrian Croucher via petsc-users mailto:petsc-users@mcs.anl.gov>> wrote: hi Thanks Matt for the explanation about this. I have been trying a test which does the following: 1) read in DMPlex from file 2) distribute it, with overlap = 1, using DMPlexDistribute()

Re: [petsc-users] DMPlexDistributeField

2019-06-26 Thread Adrian Croucher via petsc-users
hi On 25/06/19 11:22 AM, Adrian Croucher wrote: I know you originally recommended using overlap = 0 for the initial distribution and only adding overlap for the redistribution. But then Stefano indicated that it should work with overlap now. And it would simplify my code if I could use overl

Re: [petsc-users] DMPlexDistributeField

2019-06-24 Thread Adrian Croucher via petsc-users
hi Thanks Matt for the explanation about this. I have been trying a test which does the following: 1) read in DMPlex from file 2) distribute it, with overlap = 1, using DMPlexDistribute() 3) create FVM cell and face geometry vectors using DMPlexComputeGeometryFVM() 4) re-distribute, again

Re: [petsc-users] DMPlexDistributeField

2019-06-20 Thread Adrian Croucher via petsc-users
I have been trying to get this FVM geometry data re-distribution to work using DMPlexDistributeField(). It seems to be working OK for the cell geometry data (cell volumes and centroids). But it is making a mess of the face geometry data (face normals and centroids). Should I even expect DMPl

[petsc-users] DMPlexDistributeField

2019-06-05 Thread Adrian Croucher via petsc-users
hi I have some questions about using the DMPlexDistributeField() function. I have finite volume mesh geometry data stored in two local vectors created using DMPlexComputeGeometryFVM(), and I need to redistribute these after calling DMPlexDistribute() to redistribute my mesh. (I need the geome