<div dir="ltr"><div dir="ltr">On Sat, Mar 27, 2021 at 9:27 AM Nicolas Barral <<a href="mailto:nicolas.barral@math.u-bordeaux.fr">nicolas.barral@math.u-bordeaux.fr</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi all,<br>
<br>
First, I'm not sure I understand what the overlap parameter in <br>
DMPlexDistributeOverlap does. I tried the following: generate a small <br>
mesh on 1 rank with DMPlexCreateBoxMesh, then distribute it with <br>
DMPlexDistribute. At this point I have two nice partitions, with shared <br>
vertices and no overlapping cells. Then I call DMPlexDistributeOverlap <br>
with the overlap parameter set to 0 or 1, and get the same resulting <br>
plex in both cases. Why is that ?<br></blockquote><div><br></div><div>The overlap parameter says how many cell adjacencies to go out. You should not get the same</div><div>mesh out. We have lots of examples that use this. If you send your small example, I can probably</div><div>tell you what is happening.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Second, I'm wondering what would be a good way to handle two overlaps <br>
and associated local vectors. In my adaptation code, the remeshing <br>
library requires a non-overlapping mesh, while the refinement criterion <br>
computation is based on hessian computations, which require a layer of <br>
overlap. What I can do is clone the dm before distributing the overlap, <br>
then manage two independent plex objects with their own local sections <br>
etc. and copy/trim local vectors manually. Is there a more automatic way <br>
to do this ?<br></blockquote><div><br></div><div>DMClone() is a shallow copy, so that will not work. You would maintain two different Plexes, overlapping</div><div>and non-overlapping, with their own sections and vecs. Are you sure you need to keep around the non-overlapping one?</div><div>Maybe if I understood what operations you want to work, I could say something more definitive.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Thanks<br>
<br>
-- <br>
Nicolas<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>