<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Jul 9, 2015 at 7:42 AM, Alejandro D Otero <span dir="ltr"><<a href="mailto:aotero@fi.uba.ar" target="_blank">aotero@fi.uba.ar</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div>Hi, sorry if this is an obvious question, but I cannot figure out how to recover finite element nodes coordinates once I have distributed a mesh stored as a dmplex. I am using petsc4py as interface to petsc rutines.<br><br></div><div>I first created a dmplex using:<br>dm.createFromCellList()<br><br></div><div>In a sequential run I got the coordinates with:<br>Coords = dm.getCoordinates()<br></div><div><br></div>which gave a sequential vector with the coordinates of the mesh nodes.<br><br></div>When I distribute the mesh with:<br></div>dm.distribute()<br><br></div>each mpi process has it own dm but the indexing of the vector resulting from getCoordinates() or getCoordinatesLocal() seems not consistent with the local numbering of the cells and nodes.<br></div></div></div></blockquote><div><br></div><div>When the mesh is distributed, the vertices are renumbered. Thus the coordinates you get out are</div><div>for reordered local vertices, but they are consistent with the local topology (cells still contain the</div><div>right vertices) and the overlap mapping (SF still connects the shared vertices).</div><div><br></div><div>What do you need it to do?</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 dir="ltr"><div><div></div>Which is the correct way of doing this in PETSc philosophy? <br><br></div>Thanks in advance, <br><div><div><div>Alejandro<br></div></div></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">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></div>