Hi,
Thansk for your reply. So, could you please help to further elaborate this hierarchical basis function assignment flow?
How to find what the el.ndofs (where i >=4) 's related vertices are. ( i know the first 4 are just the vertices of the tetrahedral cell (in 3D case) and the other 6 ndofs (e.g, case of order=2) are computed from the first 4). and also, are these (4+6) ndofs mapped in serially in GridFunction(fes) ?
Futhermore, suppose i have a solution (dimension Nsol), and i want to load/map it into fes, how do i do it correctly?
Best,
Alex