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
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
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&
>>
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,
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