how to access the interpolation nodes of a finite element grid for order p >1

Question

Dear All

I need to implement a boundary condition for the sum of the finite element function
and another globally defined function.
For this I need to know the interpolation node coordinates, not the vertices only.

regards

Moritz

Answer 1

Take a look at the GenericDofMap class. It has a function for getting the spatial coordinates of dofs for a chosen cell.