<div dir="ltr"><div dir="ltr">On Mon, Apr 13, 2020 at 6:14 AM Pierre Jolivet <<a href="mailto:pierre.jolivet@enseeiht.fr">pierre.jolivet@enseeiht.fr</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello,<br>
I’m trying to figure out why the result of DMGetNeighbors when a DMPlex is distributed with no overlap can be nonsymmetric.<br>
</blockquote><div><br></div><div>Right now, the Plex function only calls PetscSFGetRootRanks(). It should probably call PetscSFGetLeafRanks() as well and merge the lists.</div><div>Have time to make the fix? You could commit this nice test :)</div><div><br></div><div>  Thanks,</div><div><br></div><div>    Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Is this intended? If so, is there some other easy way to get the equivalent list as in the overlapping case (I’m interested in the list of ranks sharing at least one point)?<br>
Here is a MWE.<br>
$ patch -p1 < patch.txt<br>
$ cd src/dm/impls/plex/tests/<br>
$ make ex1<br>
$ mpirun -n 2 ./ex1 -dim 2 -cell_simplex 0 -interpolate -domain_shape box -domain_box_sizes 3,3 -test_redistribute<br>
1 <— how come rank 0 has 1 neighbor<br>
0 <— but rank 1 has 0 neighbor?<br>
$ mpirun -n 2 ./ex1 -dim 2 -cell_simplex 0 -interpolate -domain_shape box -domain_box_sizes 3,3 -overlap 1<br>
1<br>
1<br>
Overlap: 1<br>
<br>
Thanks,<br>
Pierre<br>
<br>
diff --git a/src/dm/impls/plex/tests/ex1.c b/src/dm/impls/plex/tests/ex1.c<br>
index 4a5d051a66..ac2ce8eb20 100644<br>
--- a/src/dm/impls/plex/tests/ex1.c<br>
+++ b/src/dm/impls/plex/tests/ex1.c<br>
@@ -3,0 +4 @@ static char help[] = "Tests various DMPlex routines to construct, refine and dis<br>
+#include <petscsf.h><br>
@@ -420,0 +422,8 @@ PetscErrorCode CreateMesh(MPI_Comm comm, AppCtx *user, DM *dm)<br>
+      PetscInt nranks;<br>
+      const PetscInt* ranks;<br>
+      PetscSF sf;<br>
+      ierr = DMGetPointSF(distributedMesh, &sf);CHKERRQ(ierr);<br>
+      ierr = PetscSFSetUp(sf);CHKERRQ(ierr);<br>
+      ierr = DMGetNeighbors(distributedMesh, &nranks, &ranks);CHKERRQ(ierr);<br>
+      ierr = PetscSynchronizedPrintf(PETSC_COMM_WORLD,"%D\n",nranks);CHKERRQ(ierr);<br>
+      ierr = PetscSynchronizedFlush(PETSC_COMM_WORLD,PETSC_STDOUT);CHKERRQ(ierr);<br>
@@ -434,0 +444,8 @@ PetscErrorCode CreateMesh(MPI_Comm comm, AppCtx *user, DM *dm)<br>
+      PetscInt nranks;<br>
+      const PetscInt* ranks;<br>
+      PetscSF sf;<br>
+      ierr = DMGetPointSF(overlapMesh, &sf);CHKERRQ(ierr);<br>
+      ierr = PetscSFSetUp(sf);CHKERRQ(ierr);<br>
+      ierr = DMGetNeighbors(overlapMesh, &nranks, &ranks);CHKERRQ(ierr);<br>
+      ierr = PetscSynchronizedPrintf(PETSC_COMM_WORLD,"%D\n",nranks);CHKERRQ(ierr);<br>
+      ierr = PetscSynchronizedFlush(PETSC_COMM_WORLD,PETSC_STDOUT);CHKERRQ(ierr);</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><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.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>