[petsc-users] Partitioning does not work

[Lawrence Mitchell]

> On 15 Apr 2017, at 22:31, Barry Smith <bsmith at 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.shibli at 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