SUMMARYRecently the hp version of the finite element method, in which adaptivity occurs in both the size, h, of the elements and in the order, p, of the approximating piecewise polynomials, has received increasing attention. It is desirable to combine this optimal order discretization method with an optimal order algebraic solution method like multigrid. An intriguing notion is to use the values of p as the levels of a multilevel method. In this paper we present such a method, known as hp-multigrid, for high-order finite elements and hp-adaptive grids. We present a survey of the development of p-multigrid and hp-multigrid, define an hp-multigrid algorithm based on the p-hierarchical basis for the p levels and h-hierarchical basis for an h-multigrid solution of the p = 1 'coarse grid' equations, and present numerical convergence results using hp-adaptive grids. The numerical results suggest that the method has a convergence rate of