<div dir="ltr"><div dir="ltr">On Wed, Dec 23, 2020 at 11:29 AM Hui LIU <<a href="mailto:h.liu@whu.edu.cn">h.liu@whu.edu.cn</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"><div lang="ZH-CN"><div class="gmail-m_-577384413770344926WordSection1"><p class="MsoNormal"><span lang="EN-US">Dear the maintainer of the PETSc: <u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-US">I am implementing the CPU parallel finite element calculation of the unstructured grids by using the great PETSc platform. <u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">The unstructured grids are imported from other open access software, for example, Gmsh. <u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">Herein we call the node numbering for the mesh generated by the Gmsh as the original ordering, before distributing the mesh to each processor. <u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">After distributing this unstructured mesh to each process, is it possible to find the mapping between the original ordering (before distributing) and the PETSc/global ordering (after distributing)? <u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">If yes, how to extract or find it? <u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"></span></p></div></div></blockquote><div><br></div><div>It is possible using <a href="https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DMPLEX/DMPlexNaturalToGlobalBegin.html">https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DMPLEX/DMPlexNaturalToGlobalBegin.html</a></div><div>However, I do not recommend you do this. It is not scalable and makes your code fragile.</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"><div lang="ZH-CN"><div class="gmail-m_-577384413770344926WordSection1"><p class="MsoNormal"><span lang="EN-US">If no, could you give me some suggestions for better applying the displacement constraints and load boundary conditions?</span></p></div></div></blockquote><div><br></div><div>Sure. You should mark your boundaries with a marker that indicates the kind of condition to apply. This is easy in GMsh, and those</div><div>markers will be translated to DMLabel objects when you import the mesh. These markers will be preserved when the mesh is distributed,</div><div>or redistributed for load balance.</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"><div lang="ZH-CN"><div class="gmail-m_-577384413770344926WordSection1"><p class="MsoNormal"><span lang="EN-US"> <u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">Another question is how to distinguish the nodes subject to displacement boundary constraints and the nodes subject to the external force loads.</span></p></div></div></blockquote><div><br></div><div>Use a different marker.</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"><div lang="ZH-CN"><div class="gmail-m_-577384413770344926WordSection1"><p class="MsoNormal"><span lang="EN-US">Thank you very much. Looking forward to your reply.<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-US">Best regards,<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-US">Hui Liu<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">Wuhan University<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p></div></div></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>