> On 15 Apr 2017, at 22:31, Barry Smith <bsm...@mcs.anl.gov> wrote: > > > The DMPlex experts and point you in the correct direction. > > >> On Apr 15, 2017, at 4:28 PM, Orxan Shibliyev <orxan.shi...@gmail.com> wrote: >> >> About partitioning of DMPlex, manual says that, >> >> In exactly the same way as MatPartitioning and MatOrdering, we encode the >> results of a partition or order in an IS. However, the graph we are dealing >> with now is not the adjacency graph of the problem Jacobian, but the mesh >> itself. >> >> I could not find a function related to partitioning to feed in DM object. I >> know how to partition with adjacency matrix but not with DM. Would you refer >> me to an example in which DMPlex is partitioned or to a proper function?
You want DMPlexDistribute. See, for example snes/examples/tutorials/ex12.c, which solves Poisson with simplicial elements. Lawrence
