<div dir="ltr"><div class="gmail_quote"><div dir="ltr">On Thu, Jan 17, 2019 at 3:34 AM leejearl 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:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi all Petscer,<br>
        I have ask helps for some questions. Thanks for the replies<br>
from developer. I have known more about the petsc. I have also sought<br>
helps in the mail lists, and I find that there are many subjects<br>
focused on such a problem. Some subjects are listed as follow<br>
<br>
<br>
<a href="https://lists.mcs.anl.gov/pipermail/petsc-users/2019-January/037425.html" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/pipermail/petsc-users/2019-January/037425.html</a><br>
,<br>
<br>
<br>
<a href="https://lists.mcs.anl.gov/mailman/htdig/petsc-users/2015-January/024068.html" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/htdig/petsc-users/2015-January/024068.html</a><br>
,<br>
<br>
<br>
<a href="https://lists.mcs.anl.gov/mailman/htdig/petsc-users/2008-November/003633.html" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/htdig/petsc-users/2008-November/003633.html</a><br>
,<br>
<br>
<br>
<a href="https://lists.mcs.anl.gov/pipermail/petsc-users/2019-January/037436.html" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/pipermail/petsc-users/2019-January/037436.html</a><br>
<br>
......<br>
<br>
The problem can be summarized as<br>
1. Getting the value of vector from other processor<br>
2. How to set the value of ghost cell, the value of the ghost cell<br>
might be as same as a very cell in other processor. <br>
<br>
I think the above problem is very popular. I suffer such a problem when<br>
I want to treat the periodic boundary for the FVM method.</blockquote><div><br></div><div>Ah, I think I am now seeing the problem. In DMPlex, we do not implement</div><div>periodicity of the mesh by putting in extra communication. We make a</div><div>periodic mesh directly. For example, in 1D what you describe would mesh</div><div>a line segment, and then communicate the values from one end to the other.</div><div>Instead, in Plex we just mesh a circle.</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"> After the dm<br>
object is distributed, the donor cell of bound might be in other<br>
processor. Since the donor cell must be matched as a pair correctly,<br>
using the routines DMGlobaltoLocal and DMLocaltoGlobal can not get the<br>
expected results. <br>
<br>
In fact, such a problem is not very difficult for a mpi program. But we<br>
are coding with the petsc, we always think that whether we can<br>
implement more easily.<br>
<br>
Whether our developer can create a demo for us to following. I think<br>
such a demo is very useful for user.<br>
<br>
Thanks<br>
<br>
leejeal<br>
<br>
<br>
<br>
<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>