Solution value with respective vertex coordinate

[DaSchn...@gmail.com]

Dear all,

i have a question as already mentioned in the headline.
My work is based on the step- 44 tutorial and I am trying to get the solution of a vertex with the respective coordinate.
I wrote a quite circumstantial code to get every vertex coordinate at the boundary (only once) by using cell->face->vertex. For the solution I used an ElementIterator with the relevant DOFs.
But the coordinates are stored by using a cell iterator and so the solution doesn't fit necessarily to the coordinate.

The FE System has only generalized_support_points. I tried also to make use of map_dofs_to_support points with the generalized support points, but without any success.

Maybe someone has another or a better way to get rid of this.
Any help is appreciated.
Best regards and many thanks in advance,
David

[Jean-Paul Pelteret]

Dear David,

This is quite a frequently asked question (well, some permutation of it), and so I refer you to the list of answers for this search query on the forum instead of repeating the discussion here. This reminds me though that I should make a FAQ entry for it.

If those forum posts don't satisfy your question, then you can also take a look at this function in one of the code-gallery contributions.

I hope that this help!

Regards,

Jean-Paul

[daschn...@googlemail.com]

Dear Jean-Paul,

Thank you very much for the fast reply, this was very helpful. I just missed the keyword vertex_dof_index.

Regards!