[Dune] ALUGrid: double ghostnodes and communication

[Markus Blatt]

Hi Robert,

On Tue, Oct 19, 2010 at 07:05:23PM +0200, Robert Kloefkorn wrote:
> 
> The local ids should be the same, as well as the global ids.

Well, the local ids are not the same. Rebecca tested it for one
problem where she outputs the Local+GlobalIDs using
(alu_cube.localIdSet().id(*it) und alu_cube.globalIdSet().id(*it),
respectively. it is a leaf iterator here.

On processor 0 the output reads:

Center 0.5 1 0 ghost LocalId 900000031 GlobalId ([7,-1,-1,-1],1,3)
Center 0 1 0 ghost LocalId 900000032 GlobalId ([6,-1,-1,-1],1,3)
Center 0 1 0.5 ghost LocalId 900000033 GlobalId ([15,-1,-1,-1],1,3)
Center 0.5 1 0.5 ghost LocalId 900000034 GlobalId ([16,-1,-1,-1],1,3)
Center 1 1 0 ghost LocalId 900000035 GlobalId ([8,-1,-1,-1],1,3)
Center 0.5 1 0 ghost LocalId 900000036 GlobalId ([7,-1,-1,-1],1,3)
Center 0.5 1 0.5 ghost LocalId 900000037 GlobalId ([16,-1,-1,-1],1,3)

Apparently GlobalId 16 is there twice, but each time with a different LocalId.

Unfortunately, having them twice with the same local id renders my parallel
code unusable.

Cheers,

Markus

-- 
Markus Blatt, Interdisciplinary Center for Scientific Computing,
University Heidelberg, Im Neuenheimer Feld 368, D-69120 Heidelberg

'Mathematics knows no races or geographic boundaries; for
mathematics, the cultural world is one country' - David Hilbert