Hello,
I am trying to find the easiest way to get the coordinates of dofs for a piecewise quadratic function u1 over a 3D mesh:
mesh3 Th3 = tetg(ThSurf, regionlist=domain);
fespace Vh2(Th3,P2, periodic=[[XminVal, y, z], [XmaxVal, y, z], [ZminVal, x, y], [ZmaxVal, x, y]]);
Vh2 u1;
I know that having access to Th3, Vh2 and u1[][i] is enough to get the coordinates of the i-th dof, but how do I do it?
Thank you.