Large Sparse Nonsymmetric Linear Systems Research Articles

Adaptive procedure for estimating parameters for the nonsymmetric Tchebychev iteration

An iteration based upon the Tchebychev polynomials in the complex plane can be used to solve large sparse nonsymmetric linear systems whose eigenvalues lie in the right half plane. The iteration depends upon two parameters which can be chosen from knowledge of the convex hull of the spectrum of the linear operator. This paper deals with a procedure based upon the power method for dynamically estimating the convex hull of the spectrum. The stability of the procedure is discussed in terms of the field of values of the operator. Results show the adaptive procedure to be an effective method of determining parameters. The Tchebychev iteration compares favorably with several competing iterative methods.