Sandia National Laboratories

Polynomial preconditioning for linear solvers is well-known but not frequently used in current scientific applications.  Furthermore, polynomial preconditioning has long been touted as well-suited for parallel computing; does this claim still hold in our new world of GPU-dominated machines?  We give details of the GMRES polynomial preconditioner and discuss its simple implementation, its properties such as eigenvalue remapping, and choices such as the starting vector and added roots.  We compare polynomial preconditioned GMRES to related methods such as FGMRES and full GMRES without restarting. We conclude with initial evaluations of the polynomial preconditioner for parallel and GPU computing, further discussing how polynomial preconditioning can become useful to real-word applications.

A link for this talk will be sent to our mailing list a day or two in advance.  If you are not on the list and wish to be sent a link, please contact @email.

We show that Sorensen's (1992) implicitly restarted Arnoldi method

(IRAM) (including its block extension) is non-stationary simultaneous

iteration in disguise. By using the geometric convergence theory for

non-stationary simultaneous iteration due to Watkins and Elsner (1991)

we prove that an implicitly restarted Arnoldi method can achieve a

super-linear rate of convergence to the dominant invariant subspace of

a matrix. We conclude with some numerical results the demonstrate the

efficiency of IRAM.