Abstract Consideration is given to the problem of obtaining the best finite element method, in the sense of the technology that delivers a given level of accuracy with a minimum of unknowns. More specifically, this paper addresses the question: Beginning with an initial FEM approximation, how can the element sizes and spectral orders be re-distributed or changed throughout the mesh to produce the best improvement of the solution with the least number of additional unknowns?