Re: [petsc-users] Partitioning does not work

2017-04-18 Thread Lawrence Mitchell
> On 15 Apr 2017, at 22:31, Barry Smith wrote: > > > The DMPlex experts and point you in the correct direction. > > >> On Apr 15, 2017, at 4:28 PM, Orxan Shibliyev wrote: >> >> About partitioning of DMPlex, manual says that, >> >> In exactly

Re: [petsc-users] Partitioning does not work

2017-04-15 Thread Barry Smith
The DMPlex experts and point you in the correct direction. > On Apr 15, 2017, at 4:28 PM, Orxan Shibliyev wrote: > > About partitioning of DMPlex, manual says that, > > In exactly the same way as MatPartitioning and MatOrdering, we encode the > results of a

Re: [petsc-users] Partitioning does not work

2017-04-15 Thread Barry Smith
> On Apr 15, 2017, at 2:48 PM, Orxan Shibliyev wrote: > > Yes I just realized that before your reply. I would like to ask a few more > questions. > • MatCreateMPIAdj() takes "number of local rows" as argument. Well > that's what normally I don't know prior to

Re: [petsc-users] Partitioning does not work

2017-04-15 Thread Barry Smith
Put a MatView() on the mesh matrix. The mesh information in your modified version is not a list of cells of triangles or quads and hence the conversion to dual doesn't find anything. Each row of the mesh matrix (for trianglular mesh) needs three entries indicating the vertices of the

Re: [petsc-users] Partitioning does not work

2017-04-15 Thread Barry Smith
> On Apr 15, 2017, at 12:13 PM, Orxan Shibliyev wrote: > > I modified ex11.c: How did you modify it, what exactly did you change? > Tests MatMeshToDual() in order to partition the unstructured grid provided in > Petsc documentation, page 71. The resulting code is

[petsc-users] Partitioning does not work

2017-04-15 Thread Orxan Shibliyev
I modified ex11.c: Tests MatMeshToDual() in order to partition the unstructured grid provided in Petsc documentation , page 71. The resulting code is as