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

We present a Krylov subspace-type projection method for a quadratic matrix polynomial lambda(2)I - lambdaA - 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 reformulation 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.