 Analogously, splitting a hexahedron into eight subelements and varying their polynomial degrees by at most two yields 3^8 = 6,561 refinement candidates.Clearly, standard FEM error estimates providing one constant number per element are not enough to guide automatic hp-adaptivity.In standard FEM one only works with shape functions associated with grid
vertices (the so-called "vertex functions").