[Dune] Intersection iterator and global labels of vertices

Dedner, Andreas A.S.Dedner at warwick.ac.uk
Fri May 3 18:26:33 CEST 2013


Hi.
Definitely a dune question:

Probably the easiest way is to use the codim 1 subentity which (in the case of alberta) corresponds to the intersection, i.e., e=intersection.numberInSelf().
Then you can get the index of your corners for that edge by using the subIndex method.
Best
Andreas


________________________________
From: dune-bounces+a.s.dedner=warwick.ac.uk at dune-project.org [dune-bounces+a.s.dedner=warwick.ac.uk at dune-project.org] on behalf of Sacconi, Andrea [a.sacconi11 at imperial.ac.uk]
Sent: 03 May 2013 17:01
To: dune at dune-project.org
Subject: [Dune] Intersection iterator and global labels of vertices

Hi everyone,

I have a quick question about the Intersection Iterator. Before starting, I'm using AlbertGrid.

Suppose you have a triangle ABC, and you need to perform some operations related to its boundary edges, i.e., the three segments AB, AC, BC.
When you call the intersection iterator, you get three iterators, one for each edge. So far, so good.
But what is the correspondence? I mean, iterator -> geometry().corner(0) gives you the local first corner of the codim-1 entity. To which global vertex does it correspond? A, B or C?

I'm not even sure if this question is related to DUNE in general, or maybe should be addressed to DUNE-FEM.
Thanks for your help!

Andrea

__________________________________________________________

Andrea Sacconi
PhD student, Applied Mathematics
AMMP Section, Department of Mathematics, Imperial College London,
London SW7 2AZ, UK
a.sacconi11 at imperial.ac.uk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20130503/ebdc1d8a/attachment.htm>


More information about the Dune mailing list