<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Sat, Nov 16, 2013 at 7:22 AM, Michael Lange <span dir="ltr"><<a href="mailto:michael.lange@imperial.ac.uk" target="_blank">michael.lange@imperial.ac.uk</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Hi,<br>
<br>
I notice that, when creating the point SF for the parallel partition in DMPlexDistribute, cells are assigned to procs according to the original partition but vertices aren't. Was this done by design or is this a bug?<br>
</blockquote><div><br></div><div>If this were true, there would be no communication for the P1 test of SNES ex12. Here is running it with</div><div>-interpolate 1 and -dm_view ::ascii_info_detail</div><div><br></div><div>
<div>PetscSF Object: 2 MPI processes</div><div>  type: basic</div><div>    sort=rank-order</div><div>  [0] Number of roots=19, leaves=5, remote ranks=1</div><div>  [0] 4 <- (1,6)</div><div>  [0] 5 <- (1,8)</div><div>
  [0] 7 <- (1,9)</div><div>  [0] 10 <- (1,13)</div><div>  [0] 11 <- (1,17)</div><div>  [1] Number of roots=19, leaves=0, remote ranks=0</div><div>  [0] Roots referenced by my leaves, by rank</div><div>  [0] 1: 5 edges</div>
<div>  [0]    4 <- 6</div><div>  [0]    5 <- 8</div><div>  [0]    7 <- 9</div><div>  [0]    10 <- 13</div><div>  [0]    11 <- 17</div><div>  [1] Roots referenced by my leaves, by rank</div></div><div><br></div>
<div>So there are 3 vertices and 2 edges in the point SF.</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">

In case it is a bug, I have attached a patch that fixes this by using the closure of the original partition instead.<br>
<br>
Thanks and kind regards<span class=""><font color="#888888"><br>
Michael<br>
</font></span></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>