<div dir="ltr"><div class="gmail_quote"><div dir="ltr">On Wed, Dec 12, 2018 at 9:01 AM Artur Palha Da Silva Clérigo - LR via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">



<div style="word-wrap:break-word;line-break:after-white-space">
Dear All,
<div><br>
</div>
<div>I having been trying to add periodic domain functionality to my code. I am able to generate a periodic mesh with gmsh and load it using dmplex. The problem is that I am unable to have a numbering of the degrees of freedom that reflects the periodicity
 (they are numbered as if there was no periodicity).</div></div></blockquote><div><br></div><div>Why would there be a difference? They should be numbered in the same way.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word;line-break:after-white-space">
<div>I read this thread: <a href="https://lists.mcs.anl.gov/mailman/htdig/petsc-users/2018-October/036539.html" target="_blank">https://lists.mcs.anl.gov/mailman/htdig/petsc-users/2018-October/036539.html</a></div>
<div><br>
</div>
<div>There it is mentioned that DMLocalizeCoordinates is required after loading the mesh.</div></div></blockquote><div><br></div><div>Yes, this makes cell-wise coordinates, which have a jump at the periodic boundary.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word;line-break:after-white-space">
<div>My questions are the following:</div>
<div><br>
</div>
<div>1. Is it correct that I need to run DMLocalizeCoordinates on my DM?</div></div></blockquote><div><br></div><div>Yes.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word;line-break:after-white-space">
<div>2. I am currently using petsc4py but I am unable to find DMLocalizeCoordinates. Is this functionality missing or is there an alternative to it?\</div></div></blockquote><div><br></div><div>Its possible no wrapper has been written for this yet. Lisandro, is it missing? (I am at a meeting, or I would check).</div><div><br></div><div>  Thanks,</div><div><br></div><div>    Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word;line-break:after-white-space">
<div>Thank you for your time.</div>
<div><br>
</div>
<div>Kind regards,</div>
<div><br>
</div>
<div>-artur palha</div>
</div>

</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="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>