Re: [petsc-users] DMPlex: Ghost points after DMRefine

2015-12-01 Thread Morten Nobel-Jørgensen
I found a solution to my problem by using the global section instead. I still don’t quite understand what my problem ISLocalToGlobalMapping was. > Yes, that is the solution I use. Thanks – good to hear that I’m on the right track :) > I think I ignore nonlocal indices in the l2g mapping rather

Re: [petsc-users] DMPlex: Ghost points after DMRefine

2015-12-01 Thread Matthew Knepley
MappingGetIndices(ltogm, &array); > return array[point]<0; > } > > > From: Morten Nobel-Jørgensen > Date: Monday 30 November 2015 at 17:24 > To: Matthew Knepley > > Cc: "petsc-users@mcs.anl.gov" > Subject: Re: [petsc-users] DMPlex: Ghost poin

Re: [petsc-users] DMPlex: Ghost points after DMRefine

2015-12-01 Thread Morten Nobel-Jørgensen
rom, &cellsTo); for (PetscInt i=cellsFrom;imailto:knep...@gmail.com>> Date: Monday 30 November 2015 at 14:08 To: Morten Nobel-Jørgensen mailto:m...@mek.dtu.dk>> Cc: "petsc-users@mcs.anl.gov<mailto:petsc-users@mcs.anl.gov>" mailto:petsc-users@mcs.anl.gov&

Re: [petsc-users] DMPlex: Ghost points after DMRefine

2015-11-30 Thread Morten Nobel-Jørgensen
>> Date: Monday 30 November 2015 at 14:08 To: Morten Nobel-Jørgensen mailto:m...@mek.dtu.dk>> Cc: "petsc-users@mcs.anl.gov<mailto:petsc-users@mcs.anl.gov>" mailto:petsc-users@mcs.anl.gov>> Subject: Re: [petsc-users] DMPlex: Ghost points after DMRefine On Mon,

Re: [petsc-users] DMPlex: Ghost points after DMRefine

2015-11-30 Thread Matthew Knepley
On Mon, Nov 30, 2015 at 7:01 AM, Morten Nobel-Jørgensen wrote: > I have a very simple unstructured mesh composed of two triangles (four > vertices) with one shared edge using a DMPlex: > > /|\ > / | \ > \ | / > \|/ > > After distributing this mesh to two processes, each process owns a > triangl