[petsc-users] DMPlexDistributeField

[Adrian Croucher]

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 overlap for the initial distribution 
> (because if dual porosity cells are not being used, then there is no 
> second redistribution).
>
Actually I have also tried this test with overlap = 0 for the initial 
distribution, and overlap = 1 for the redistribution, and that doesn't 
seem to give correct results for the face geometry either. So that would 
suggest that, if there is a problem, it probably doesn't lie in the 
stuff that removes overlap cells from the original DM.

I am running PETSc master 3.11.2. I see 3.11.3 is just out, not sure if 
there are any changes there which might affect this?

- Adrian

> Is this a bug or is there something I'm doing wrong?
>
> - Adrian
>
-- 
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.croucher at auckland.ac.nz
tel: +64 (0)9 923 4611