Re: [petsc-users] Understanding DMPlexDistribute overlap

2017-04-06 Thread Francesco Caimmi
Hi Michael, thank you for verifying that the test works. I actually was on `maint`, as per https://www.mcs.anl.gov/petsc/download/index.html. Now that I switched to `master ` everything works as expected, both my test example and the DMPlex test ex12, which I can now build without issue. I

Re: [petsc-users] Understanding DMPlexDistribute overlap

2017-04-06 Thread Michael Lange
Hi Francesco, Ok, I can confirm that your test runs fine for me with the latest master branch. I'm attaching the log for two processes up to overlap 7, where the entire mesh is effectively replicated on each partition. The command I ran was: for OL in 1 2 3 4 5 6 7; do mpiexec -np 2 python

Re: [petsc-users] Understanding DMPlexDistribute overlap

2017-04-06 Thread Francesco Caimmi
Dear Matt, thanks for your reply. On mercoledì 5 aprile 2017 21:09:43 CEST Matthew Knepley wrote: > On Wed, Apr 5, 2017 at 6:03 AM, Francesco Caimmi >[...] > > the program fails with the error message captured in the attached file > > error.log. Changing the number

Re: [petsc-users] Understanding DMPlexDistribute overlap

2017-04-05 Thread Matthew Knepley
On Wed, Apr 5, 2017 at 6:03 AM, Francesco Caimmi wrote: > Hi Michael, > > thanks for the prompt reply! > > While I am happy I mostly got it right, this means I have some kind of > problem > I cannot solve on my own. :( > > I have this very simple 2D mesh I am

Re: [petsc-users] Understanding DMPlexDistribute overlap

2017-04-05 Thread Francesco Caimmi
Hi Michael, thanks for the prompt reply! While I am happy I mostly got it right, this means I have some kind of problem I cannot solve on my own. :( I have this very simple 2D mesh I am experimenting with: a rectangle with 64 vertexes and 45 cells (attached in exodus format as cantilever.e);

Re: [petsc-users] Understanding DMPlexDistribute overlap

2017-04-05 Thread Michael Lange
Hi Francesco, Your description is almost correct: the overlap defines the topological depth of shared entities as counted in "neighboring cells", where a cell counts as a neighbor of an owned cell according to the defined adjacency style. So for overlap=0 only facets, edges and vertices may

[petsc-users] Understanding DMPlexDistribute overlap

2017-04-05 Thread Francesco Caimmi
Dear all, I was playing with DMPlex objects and I was trying to exactly figure out what the `overlap` parameter in DMPlexDistribute does. From the tutorial "Flexible, Scalable Mesh and Data Management using PETSc DMPlex" (slide 10) and from the work by Knepley et al. "Unstructured Overlapping