Orthogonal polynomials of the [formula omitted]-linear generalized minimal residual method

Abstract

The speed of convergence of the R-linear GMRES method is bounded in terms of a polynomial approximation problem on a finite subset of the spectrum. This result resembles the classical GMRES convergence estimate except that the matrix involved is assumed to be condiagonalizable. The bounds obtained are applicable to the CSYM method, in which case they are sharp. Then a new three term recurrence for generating a family of orthogonal polynomials is shown to exist, yielding a natural link with complex symmetric Jacobi matrices. This shows that a mathematical framework analogous to the one appearing with the Hermitian Lanczos method exists in the complex symmetric case. The probability of being condiagonalizable is estimated with random matrices.