<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Dec 19, 2017 at 11:40 AM, Yann JOBIC <span dir="ltr"><<a href="mailto:yann.jobic@univ-amu.fr" target="_blank">yann.jobic@univ-amu.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br>
<br>
We want to extract the cell connectivity from a DMPlex. We have no problem for a sequential run.<br></blockquote><div><br></div><div>Do you want it on disk? If so, you can just DMView() for HDF5. That outputs the connectivity in a global numbering.</div><div>I can show you the calls I use inside if you want. I usually put</div><div><br></div><div>  DMViewFromOptions(dm, NULL, "-dm_view")</div><div><br></div><div>Then</div><div><br></div><div>  -dm_view hdf5:mesh.h5</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">
However for parallel ones, we need to get the node numbering in the global ordering, as when we distribute the mesh, we only have local nodes, and thus local numbering.<br>
<br>
It seems that we should use DMGetLocalToGlobalMapping (we are using Fortran with Petsc 3.8p3). However, we get the running error :<br>
<br>
[0]PETSC ERROR: No support for this operation for this object type<br>
[0]PETSC ERROR: DM can not create LocalToGlobalMapping<br>
<br>
Is it the right way to do it ?<br>
<br>
Many thanks,<br>
<br>
Regards,<br>
<br>
Yann<br>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><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.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div>
</div></div>