<div dir="ltr">On Thu, Mar 7, 2013 at 3:56 PM, Garnet Vaz <span dir="ltr"><<a href="mailto:garnet.vaz@gmail.com" target="_blank">garnet.vaz@gmail.com</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
For an small distributed unstructured mesh  with overlap 1<br>
which is *not* interpolated yet, if I do a DMPlexGetCellNumbering()<br>
I obtain an IS as follows:<br>
<br>
[0] Number of indices in set 6<br>
[0] 0 -6<br>
[0] 1 0<br>
[0] 2 -8<br>
[0] 3 1<br>
[0] 4 2<br>
[0] 5 3<br>
[1] Number of indices in set 6<br>
[1] 0 4<br>
[1] 1 5<br>
[1] 2 6<br>
[1] 3 -1<br>
[1] 4 7<br>
[1] 5 -4.<br>
<br>
Running the same example but by DMPlexInterpolate() and then<br>
DMPlexGetCellNumbering() provides instead the output:<br>
<br>
[0] Number of indices in set 6<br>
[0] 0 0<br>
[0] 1 1<br>
[0] 2 2<br>
[0] 3 3<br>
[0] 4 4<br>
[0] 5 5<br>
[1] Number of indices in set 6<br>
[1] 0 6<br>
[1] 1 7<br>
[1] 2 8<br>
[1] 3 9<br>
[1] 4 10<br>
[1] 5 11<br></blockquote><div><br></div><div style>There is something wrong with this result. It has to do with the order in which you created</div><div style>things. Right now, I do not have DMInterpolate() update the communication structures, since</div>
<div style>it would involve coordination between processes. Thus you should call DMInterpolate()</div><div style>before DMDistribute().</div><div style><br></div><div style>  Thanks,</div><div style><br></div><div style>     Matt</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
The former result I find very useful since it provides the number of<br>
ghost cells on each processor (the negative IS values) and moreover<br>
the negative IS values also indicate the global numbering of the<br>
ghost cells. Eg: On the first processor the cell 0 is a ghost cell with<br>
its global numbering 5.<br>
<br>
The behavior of DMPlexGetCellNumbering() after the interpolation does<br>
not bother me but I would like to know in the former case if the output<br>
is the right output i.e, it is not merely a side effect which might get<br>
changed later on.<br>
<br>
Thank You.<br>
<br>
--<br>
Regards,<br>
Garnet<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>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>