[Dune] vote: what should the codim<3> iterator do in the curvilinear case

Aleksejs Fomins aleksejs.fomins at lspr.ch
Thu Jan 15 13:33:43 CET 2015


-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

Dear Dune,

I just realized a slightly puzzling thing:

In curvilinear geometry, an element has a set of interpolatory
vertices. For a tetrahedron, there are always 4 corners, plus a bunch
of internal vertices based on interpolation order of the element.

Currently, I have implemented it in the grid such that all
interpolatory vertices are codim<3> entities. This means that
currently, if one iterates over codim<3> entities, the interpolatory
vertices will also participate.

Having said that, I personally think it does not make sense. If
somebody wishes to implement a nodal-based FEM code, they probably
only care about the corners. I do not see a reason why should a user
ever care about the non-corner interpolatory vertices, unless he
directly accesses the Geometry class.

My proposal:
1) All vertices are codim<3> entities
2) standard codim<3> iterator will only iterate over corners
3) an additional InterpolatoryIterator will be provided which will be
able to iterate over all codim<3> entities
4) DataHandle communication interface will only explicitly be allowed
to communicate over corners.

What do people think about this approach

Best,
Aleksejs
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.12 (GNU/Linux)

iQIcBAEBAgAGBQJUt7OmAAoJEDkNM7UOEwMZzycP/3wt0/eAuFO1flvfP8uxKEvJ
FSig076Zt68DYiziveHOe540AJ5Si4X7psk08r9Vdu7VyDSn+xZZ8rzoAwMexLvV
etj0F0y2zDMyfjrgyzGxUss/Z82BFmLMKKxnCT8LRpjCCznEInfqu9SpFLj2dUXu
bSw0fYRib85vJsRtvYdAUcGIM9SEEJVnl/5tcga7+vUR1C/eN2zKdkCG+UkqdKMV
3VxYyz8Ayv9khY75MTjHd+BVX8dpNKwYDqtLJzMFiz4Rcm1hhjISORszhqSOTDuT
+rwJMcUk9JhVrJ5IUNBa3ePBz2Aj/lrMruag2nanIUTh6SNaQ64QIN4fRWFLMPh5
fo+IrSjJoFFR1VsBCetBCuxTWekCzoPDwX3BGQqDFB/u9bq4eE0+SI+Bah+YgBO8
v8NlUn83+ypWFk2ejvP0opSiI9Kum3iqPbG2df9rHKZNxG97F6OQS91sPXR6295e
OX653s9UhRE4GQ+oox9FZOevB1HeZmXBx97PW45NyB205fqPN2dQKtpsJixoLNOn
zoL2tZIlfcRKHf3B9qag+sIVsT2fOj6GoP93FJede7/BhNBHE49RYE2vrGzicZdw
aZyFIU66BQl8HqUwhZsCu74f0R6T+N5cIg2Y6YoY4m87Vk8Vw+50JQRlNrjoN8xs
B41PIBR6eqgfZT5SWaAX
=ng3b
-----END PGP SIGNATURE-----



More information about the Dune mailing list