[Dune] Iterator over Ghosts Nodes comes to the same node multiple times

Rebecca Neumann rebecca.neumann at iwr.uni-heidelberg.de
Thu Oct 13 21:05:35 CEST 2011


Dear Patrick,
That's right, in ALUGrid each border face has his "own" ghosts. 
Therefore you have several ghosts with the same global id (and so the 
same coordinates) but different local ids.

We had this topic last year, you can find it here:
http://lists.dune-project.org/pipermail/dune/2010-October/007407.html

Best,
Rebecca

On 13.10.2011 20:49, Patrick Leidenberger wrote:
> Hi all Dune-Developers,
>
> I am a little bit confused by the behavior of
>
> typedef typename GridType::template Codim<3>::
>   template Partition<Dune::Ghost_Partition>::LeafIterator
>     IteratorGhost;
>
> iterator. It comes to the same node multiple times! This is not the 
> behavior I expected.
> So my question is: Is this a bug or a feature?
>
> My system:
> Current dune from SVN, last co 1h ago.
> ALUgrid 1.23, with simplices in 3d
> Standart Ubuntu 11.04, 10.10 with standart gcc, openmpi & co.
>
> That you can reproduce my results, I modified the test-alugrid.cc from 
> dune-grid/dune/grid/test and attached the file. Run it e.g. with
> mpirun -np 3 ./test-alugrid 3dsimp
> It iterates over InteriorBoarder_Partition and GhostPartition. For 
> both cases it collects all coordinates of the nodes, put them in a 
> list (separate ones) and sort the lists. You get e.g.
> ...
> rank=1, Ghost,   -0.848528, -1.84853, 6.875
> rank=1, Ghost,   -0.494975, -1.49497, 6.875
> rank=1, Ghost,   -0.494975, -1.49497, 6.875
> rank=1, Ghost,   -0.494975, -1.49497, 6.875
> rank=1, Ghost,   -0.494975, -1.49497, 6.875
> rank=1, Ghost,   -5.14352e-16, 0.2, 9.25
> rank=1, Ghost,   -3.00038e-16, -0.3, 9.25
> rank=1, Ghost,   -3.00038e-16, -0.3, 9.25
> rank=1, Ghost,   -3.00038e-16, -0.3, 9.25
> rank=1, Ghost,   -3.00038e-16, -0.3, 9.25
> rank=1, Ghost,   -3.00038e-16, -0.3, 9.25
> rank=1, Ghost,   -8.57253e-17, -0.8, 8.75
> rank=1, Ghost,   -8.57253e-17, -0.8, 8.75
> ...
>
> Thanks a lot for some tips, best regards
>
> Patrick
>
>
> _______________________________________________
> Dune mailing list
> Dune at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20111013/e5860080/attachment.htm>


More information about the Dune mailing list