[Dune] Intersection iterator problem

Christian Engwer christi at uni-hd.de
Tue Aug 4 15:26:11 CEST 2009


Dear Aleksejs,

please try to be more specific. Your description leaves me clueless
because it lacks a lot of information.

> I have got a problem with intersection iterator.
> Shortly the program runs as follows:
> 1) Load the grid containing 2 adjacent tetrahedrons from a file (this 
> part is certainly correct)
> 2) Use leaf iterator to loop over all elements of the grid
>     2.1) Use intersection iterator to loop through all faces of each 
> elements

Please note, that the intersection iterators do not interate over
intersections, but over intersections.

>        2.1.1) For each face it prints out its local number, global 
> number, and whether it is a neighbor.
> I have previously declared corresponding mappers to obtain the global 
> numbers.
> 
> The problem is that the neighboring elements are not the ones with the 
> same global number.
> 
> Could you please suggest what might be the problem.

Just some suggestions on what might be helpful information, if you
want an answer:

* which grid are you using?
* how does your grid look like (in detail)?
* which result exactly did you expect, which result did you obtain?
* perhaps everything with a stripped down snipped from your code.

>From your mails I guess that you are just starting to use Dune and
also that you didn't write much code yourself yet.

At PSI you have many colleagues which might help you in getting the
basic right and which can help you learning how to search for a
bug. We are trying to help everybody with DUNE problems, but we also
don't have super natural powers to do thought reading, so you have to
give a proper problem description or you are on your own. Please ask
someone at PSI to help you. You two might have a look at your code,
most often it helps already to try to properly state the problem and
suddenly you find the bug yourself. If you two can't find the problem,
he/she can help you to collect all necessary information to state the
question on the list.

Cheers
Christian





More information about the Dune mailing list