[petsc-dev] Will DMPlexDistribute only ever distribute Cells?

[Jacob Faibussowitsch]

Hello All,

TL;DR: Is it possible now, or is it a planned feature, for plex to distribute over anything but points with height = 0? 

If I understand this correctly when plex currently partitions a graph, points with height 0 are only owned by a single process, but all other points can be co-owned by multiple procs.  For example for a 2D plex with 8 vertices, 12 edges, and 6 cells over 2 procs these are the global-local IS’s for all points on processes (negative values indicate ownership by another proc) the final IS corresponding to cells will always have positive values as each proc is the sole owner of its cells.

VERTICES
IS Object: 2 MPI processes
  type: general
[0] Number of indices in set 7
[0] 0 -2
[0] 1 0
[0] 2 -3
[0] 3 -4
[0] 4 -5
[0] 5 -6
[0] 6 -8
[1] Number of indices in set 7
[1] 0 1
[1] 1 2
[1] 2 3
[1] 3 4
[1] 4 5
[1] 5 6
[1] 6 7

EDGES
IS Object: 2 MPI processes
  type: general
[0] Number of indices in set 9
[0] 0 0
[0] 1 1
[0] 2 -4
[0] 3 -5
[0] 4 -6
[0] 5 2
[0] 6 -7
[0] 7 -9
[0] 8 -11
[1] Number of indices in set 9
[1] 0 3
[1] 1 4
[1] 2 5
[1] 3 6
[1] 4 7
[1] 5 8
[1] 6 9
[1] 7 10
[1] 8 11

CELLS
IS Object: 2 MPI processes
  type: general
[0] Number of indices in set 3
[0] 0 0
[0] 1 1
[0] 2 2
[1] Number of indices in set 3
[1] 0 3
[1] 1 4
[1] 2 5 

Best regards,

Jacob Faibussowitsch
(Jacob Fai - booss - oh - vitch)
Cell: (312) 694-3391

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20200416/56940985/attachment.html>