When performing a refinement, the old set of vertices is a subset of the set of vertices after refinement.
Is there any way of identifying the newly created vertices?
Thank you very much!
The new vertices are appended at the end of the vertex array. Store the old number of vertices.
Or you can create an H1(order=1) fe-space, and ask for the prolongation operator, it knows the space dimension at every level:
fes.Prolongation().LevelDofs(level)
Joachim
wow! that’s kind of nice!
thank you very much!