[firedrake] Mesh coordinates for higher degree of FE

[Luca Mechelli]

*******************
This email from originates from outside Imperial. Do not click on links and attachments unless you recognise the sender. 
If you trust the sender, add them to your safe senders list https://spam.ic.ac.uk/SpamConsole/Senders.aspx to disable email stamping for this address.
*******************
Dear all,

I am wondering if it is possible to access the list of all points used by piecewise quadratic (or higher order) elements on a regular mesh.
For example, let's say I use the simple example of the Poisson equation at https://www.firedrakeproject.org/demos/poisson.py.html

If I print mesh.coordinates.dat.data_ro I can see the verteces of the mesh, but not the middle points used to build the quadratic elements, then if I print u.dat.data_ro the array has not the same length of mesh.coordinates.dat.data_ro because the middle points nodes are missing from the second array. Is it possible to obtain the full list of coordinates used? If yes, how?

Thank you in advance,
best regards,
Luca Mechelli


--
Dr. Luca Mechelli
Wissenschaftlicher Mitarbeiter
Arbeitsgruppe Numerik
Fachbereich Mathematik und Statistik

Universität Konstanz
Fach D 198
D-78464 Konstanz
Telefon: +49-(0) 7531-88-3646