We present an eigen-based high-order expansion basis for the spectral element approach with structured elements. The new basis exhibits a numerical efficiency significantly superior , in terms of the conditioning of coefficient matrices and the number of iterations to convergence for the conjugate gradient solver, to the commonly-used Jacobi polynomial-based expansion basis. This basis results ...