[petsc-users] How to use DMPlex_Distribute(_Overlap) without redistribution of the mesh

[Matthew Knepley]

On Tue, Oct 1, 2019 at 7:39 PM Fabian.Jakub via petsc-users <
petsc-users at mcs.anl.gov> wrote:

> Dear Petsc team,
> dear Matt,
>
> I have a 3D, distributed DMPlex with 0 overlap in FVM adjacency.
> Now, I would like to add overlap but not change the local portions.
> I was wondering if you could suggest how to do that best/easiest.
>
> My problem is that just using DMPlexDistribute changes the partitioning,
> i.e. cells change owners.
>
> I guess I could trace the steps you took in
> DMPlexDistributeOverlap
>

DistributeOverlap() should not change the partition. Are you sure it does?
I think you just call that.

  Thanks,

    Matt


> and the following stuff you have in
> DMPlexDistribute() { if overlap > 0 }
>
> Another thought was to use a custom partitioner so as not to change the
> local layouts. Is this a good idea?
>
> In case that you suggest to go for first steps (the duplication of the
> DMPlexDistributeOverlap stuff) could you add a fortran stub for it?
>
> Many thanks,
>
> Fabian
>


-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener

https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20191002/1d36642f/attachment.html>