Linearizations for interpolation bases – a comparison I
share
One strategy to solve a nonlinear eigenvalue problem T(λ)x=0 is to solve a polynomial eigenvalue problem (PEP) P(λ)x=0 that approximates the original problem through interpolation. Then, this PEP is usually solved by linearization. Most of the literature about linearizations assumes that P(λ) is expressed in the monomial basis but, because of the polynomial approximation techniques, in this context, P(λ) is expressed in a non-monomial basis. The bases used with most frequency are the Chebyshev basis, the Newton basis and the Lagrange basis. In this paper we construct a family of linearizations of P(λ) that is easy to construct from the matrix coefficients of P(λ) when this polynomial is expressed in any of those three bases. We also provide recovery formulas of eigenvectors (when P(λ) is regular) and recovery formulas of minimal bases and minimal indices (when P(λ) is singular). Our ultimate goal is to compare the numerical behavior of these linearizations, within the same family (to select the best one) and with the linearizations of other families based on the location of the eigenvalues with respect to the interpolation nodes.
READ FULL TEXT