Is there a way to findout the "coordinates" of a degree of freedom? I want to strongly enforce all degrees of freedom on a boundary be 0 in a DG dof_handler. I know this is not normal for DG. It is done to avoid an edge that otherwise contains a singularity. However, the degrees of freedom are not logically on the boundary and as far as I understand they aren't even on the face of their respective cells but the interior. I can determine the cells that contain degrees of freedom that are physically on this outer boundary but not "logically" there within Deal.II. I just need to be able to differentiate between the other degrees of freedom within the cell now so that I can set the constraints.
Thanks for any time spent answering this,
Sean Johnson