A Krylov subspace method for quadratic matrix polynomials with application to constrained least squares problems

Abstract. We present a Krylov subspace–type projection method for a quadratic matrix poly-nomial λ2I − λA − B that works directly with A and B without going through any linearization. We discuss a special case when one matrix is a low rank perturbation of the other matrix. We also apply the method to solve quadratically constrained linear least squares problem through a refor-mulation of Gander, Golub, and von Matt as a quadratic eigenvalue problem, and we demonstrate the effectiveness of this approach. Numerical examples are given to illustrate the efficiency of the algorithms