I have been trying to get a hang on the JuAFEM package and running the Heat Equation example from the documentation I have noticed that the ordering of the resulting vector appears to not correspond to the ordering of nodes I find in the mesh.
This makes sense to me, as it seems to be fairly standard to reorder the nodes using a reverse cuthill-mckee. However, I would like to be able to
Get the dof node coordinates. For 1st order elements these correspond to the mesh node coordinates, but for higher order elements they won’t.
For the case of 1st order elements I would like the permutation mapping the elements in the solution vector u to the node ordering in the mesh.