[Dune] Reference tetrahedron, again

[Oliver Sander]

>
> Well, if nobody else has a preference, I suggest we adopt mine ;-).

Does your numbering have the required consistency?  That is, is


my_element->subIndex( refElement(i, 1, j, dim) )

always equal to

my_element->entity<1>(i)->subIndex<dim>(j)

?

If so, I have no objections to your numbering.

>
> (Having the same orientation on all faces makes sense if you have a complete 
> set of entities (i.e. you actually store faces, egdes etc.). In this case, 
> you can implement a generic routine to calculate the twist of a face relative 
> to an element's face. This is of course not needed in Dune, but I think it is 
> a nice feature nevertheless.)
>
> Adi
>
> Peter Bastian wrote:
>
>> Robert Kloefkorn schrieb:
>> 
>>> Peter Bastian wrote:
>>> 
>>> 
>>>> Hi,
>>>> 
>>>> what is the class ReferenceTopologySet for anyway? I thought, on the
>>>> last dune meeting we decided that the reference elements are those that
>>>> are in the class ReferenceElement. This class contains a superset of the
>>>> information in ReferenceTopologySet and is implemented for more cases.
>>>> 
>>> Wir haben beschlossen, dass diese es sein sollen. Da zum Zeitpunkt des
>>> Beschlusses aber
>>> noch nicht klar war, wie die Referenzsimplices aussehen, müssen wir
>>> diese Diskussion eben jetzt nachholen. Ich denke bei der Nummerierung
>>> der Knoten und Flächen sind wir uns einig, dass es die sein soll, die
>>> auch im web steht. Bei der Nummerierung der Knoten der Flächen wäre es
>>> allerdings schön, wenn diese so wäre, wie Adrian es vorgeschalgen hat,
>>> da dann z.B. alle Normalen in eine Richtung (nach innen) zeigen. Bei der
>>> von Heidelberg vorgeschlagenen Variante ist dies nicht der Fall.
>>> Deshalb wäre es schön, mal die Vorteile der von Heidelberg
>>> vorgeschlagenen Variante zu hören, damit wir abwägen können, welche
>>> Varainte ins Interface aufgenommen sind.
>>> 
>>> 
>>> 
>> I do not have any preference for any enumeration. The thing with the
>> normals is nice but
>> I dont see where this is needed because you get the normals from the
>> intersection iterator and don't compute them for yourself. What I would
>> prefer is a numbering that can be defined for all simplices in n
>> dimensions (although we do not have a grid manager that can do that).
>> Might be interesting for the future.
>> 
>> -- Peter
>> ------------------------------------------------------------------
>> Peter Bastian, IWR,Uni Heidelberg, INF 348,R 020, 69120 Heidelberg
>> email: Peter.Bastian at iwr.uni-heidelberg.de   Tel: +49 6221 54 4984
>> WWW: http://www.iwr.uni-heidelberg.de/~Peter.Bastian Fax: ... 8860
>> 
>> 
>> _______________________________________________
>> Dune mailing list
>> Dune at hal.iwr.uni-heidelberg.de
>> http://hal.iwr.uni-heidelberg.de/cgi-bin/mailman/listinfo/dune
>> 
>
> _______________________________________________
> Dune mailing list
> Dune at hal.iwr.uni-heidelberg.de
> http://hal.iwr.uni-heidelberg.de/cgi-bin/mailman/listinfo/dune
>