[Dune] face loop

Marco Cisternino marco.cisternino at optimad.it
Tue Jun 5 17:27:43 CEST 2012


Hi,
I'm looking at the tutorial codes, specifically, at evolve.hh. I see 
that a loop over the intersections is nested in a loop over the elements.
I would like to loop over all the faces (old intersections) to compute 
the fluxes but I would need to know which elements the face belong to.
How could I do?
I defined the face iterator on the leaves like this:
typedef typename GridView::template Codim<1>::template 
Partition<Dune::All_Partition>::Iterator FLeafIterator;
Is there a way to ask for the elements of codim=0 which own the face 
pointed by FLeafIterator??
Thanks for any help.

Cheers,
Marco




More information about the Dune mailing list