<div dir="ltr"><div dir="ltr">On Wed, Sep 29, 2021 at 5:18 PM Eric Chamberland <<a href="mailto:Eric.Chamberland@giref.ulaval.ca">Eric.Chamberland@giref.ulaval.ca</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,<br>
<br>
I come back with _almost_ the original question:<br>
<br>
I would like to add an integer information (*our* original element <br>
number, not petsc one) on each element of the DMPlex I create with <br>
DMPlexBuildFromCellListParallel.<br>
<br>
I would like this interger to be distribruted by or the same way <br>
DMPlexDistribute distribute the mesh.<br>
<br>
Is it possible to do this?<br></blockquote><div><br></div><div>I think we already have support for what you want. If you call</div><div><br></div><div>  <a href="https://petsc.org/main/docs/manualpages/DM/DMSetUseNatural.html">https://petsc.org/main/docs/manualpages/DM/DMSetUseNatural.html</a></div><div><br></div><div>before DMPlexDistribute(), it will compute a PetscSF encoding the global to natural map. You</div><div>can get it with</div><div><br></div><div>  <a href="https://petsc.org/main/docs/manualpages/DMPLEX/DMPlexGetGlobalToNaturalSF.html">https://petsc.org/main/docs/manualpages/DMPLEX/DMPlexGetGlobalToNaturalSF.html</a></div><div><br></div><div>and use it with</div><div><br></div><div>  <a href="https://petsc.org/main/docs/manualpages/DMPLEX/DMPlexGlobalToNaturalBegin.html">https://petsc.org/main/docs/manualpages/DMPLEX/DMPlexGlobalToNaturalBegin.html</a></div><div><br></div><div>Is this sufficient?</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>
Eric<br>
<br>
On 2021-07-14 1:18 p.m., Eric Chamberland wrote:<br>
> Hi,<br>
><br>
> I want to use DMPlexDistribute from PETSc for computing overlapping <br>
> and play with the different partitioners supported.<br>
><br>
> However, after calling DMPlexDistribute, I noticed the elements are <br>
> renumbered and then the original number is lost.<br>
><br>
> What would be the best way to keep track of the element renumbering?<br>
><br>
> a) Adding an optional parameter to let the user retrieve a vector or <br>
> "IS" giving the old number?<br>
><br>
> b) Adding a DMLabel (seems a wrong good solution)<br>
><br>
> c) Other idea?<br>
><br>
> Of course, I don't want to loose performances with the need of this <br>
> "mapping"...<br>
><br>
> Thanks,<br>
><br>
> Eric<br>
><br>
-- <br>
Eric Chamberland, ing., M. Ing<br>
Professionnel de recherche<br>
GIREF/Université Laval<br>
(418) 656-2131 poste 41 22 42<br>
<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>