<div dir="ltr">Thank you very much Jed, now I have a clearer explanation and another point of view of the situation. This is like take the red or blue pill for me...<br></div><div class="gmail_extra"><br><div class="gmail_quote">2017-06-28 21:06 GMT+02:00 Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">There is no single correct way to do this so you can do whatever makes<br>
the most sense for your application.  That ranges from calling ParMETIS<br>
directly and creating a numbering using any scheme you like to using<br>
PETSc functions for everything.  Note that assembled linear algebra<br>
(matrices and vectors) cares about the dof connectivity graph which is<br>
different from the element graph.  Some people (usually those using<br>
lowest order methods) don't pay attention to elements at all (even<br>
integrating elements redundantly instead of communicating again) while<br>
others focus entirely on elements and make arbitrary decisions about<br>
vertex ownership.  For every choice, you can find someone that did it<br>
one way and someone else that did it the other way -- both will be<br>
convinced that their choice was the only correct choice.<br>
<div class="HOEnZb"><div class="h5"><br>
Guido Giuntoli <<a href="mailto:giuntoli1991@gmail.com">giuntoli1991@gmail.com</a>> writes:<br>
<br>
> Hi,<br>
><br>
> I am using ParMETIS to get a partition of my mesh (a group of elements for<br>
> each process) and now I want to renumber the nodes contiguous for each<br>
> process according to their rank. Can I take advantage of the IS functions<br>
> (or others like AO) directly from this point or is crucial to do the<br>
> partition with ParMETIS using PETSc functions and apply<br>
> ISPartitionToNumbering like in the manual's example ?<br>
><br>
> Another question, which is the best way of determine the ownership of the<br>
> nodes ? for example : node 101 belongs to an element in process 0 and to an<br>
> element in process 1, how do I  take the decision ? In the past I did a<br>
> partition with METIS and he returned an array of nodes belonging like the<br>
> elements so that was easy but now my only idea is to get the new<br>
> renumbering and check if the new renumbering fall inside my range.<br>
><br>
> Thank you, Guido.<br>
</div></div></blockquote></div><br></div>