[petsc-users] DMPlex overlap

Nicolas Barral nicolas.barral at math.u-bordeaux.fr
Wed Mar 31 10:59:30 CDT 2021


On 31/03/2021 17:47, Matthew Knepley wrote:
> On Sat, Mar 27, 2021 at 9:27 AM Nicolas Barral 
> <nicolas.barral at math.u-bordeaux.fr 
> <mailto:nicolas.barral at math.u-bordeaux.fr>> wrote:
> 
>     Hi all,
> 
> 
> I promise to read this today. While I am doing that, I have a branch
> 
> https://gitlab.com/petsc/petsc/-/commits/knepley/feature-plex-vertex-mapping
> 
> which I think does the face labeling you want, but I have not tested it yet.
> 
Thanks Matt, I'll have a look later today :) I also have working code 
for what I need - maybe itś  time to share my branch even if it's still 
in progress.

Thanks

-- 
Nicolas

>    Thanks,
> 
>      Matt
> 
>     First, I'm not sure I understand what the overlap parameter in
>     DMPlexDistributeOverlap does. I tried the following: generate a small
>     mesh on 1 rank with DMPlexCreateBoxMesh, then distribute it with
>     DMPlexDistribute. At this point I have two nice partitions, with shared
>     vertices and no overlapping cells. Then I call DMPlexDistributeOverlap
>     with the overlap parameter set to 0 or 1, and get the same resulting
>     plex in both cases. Why is that ?
> 
>     Second, I'm wondering what would be a good way to handle two overlaps
>     and associated local vectors. In my adaptation code, the remeshing
>     library requires a non-overlapping mesh, while the refinement criterion
>     computation is based on hessian computations, which require a layer of
>     overlap. What I can do is clone the dm before distributing the overlap,
>     then manage two independent plex objects with their own local sections
>     etc. and copy/trim local vectors manually. Is there a more automatic
>     way
>     to do this ?
> 
>     Thanks
> 
>     -- 
>     Nicolas
> 
> 
> 
> -- 
> 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/>


More information about the petsc-users mailing list