Relaxation strategies for nested Krylov methods

  • Authors:

  • Jasper van den Eshof;Gerard L. G. Sleijpen;Martin B. van Gijzen

  • Affiliations:

  • Department of Mathematics, University of Düsseldorf Universitätsstr, 1, D-40224, Düsseldorf, Germany;Department of Mathematics, Utrecht University, P.O. Box 80.010, NL-3508 TA Utrecht, The Netherlands;CERFACS, 42 Avenue Gaspard Coriolis, 31057 Toulouse, Cedex 01, France

  • Venue:
  • Journal of Computational and Applied Mathematics
  • Year:
  • 2005

Quantified Score

Hi-index 7.29

Visualization

Abstract

This paper studies computational aspects of Krylov methods for solving linear systems where the matrix-vector products dominate the cost of the solution process because they have to be computed via an expensive approximation procedure. In recent years, so-called relaxation strategies for tuning the precision of the matrix-vector multiplications in Krylov methods have proved to be effective for a range of problems. In this paper, we will argue that the gain obtained from such strategies is often limited. Another important strategy for reducing the work in the matrix-vector products is preconditioning the Krylov method by another iterative Krylov method. Flexible Krylov methods are Krylov methods designed for this situation. We combine these two approaches for reducing the work in the matrix-vector products. Specifically, we present strategies for choosing the precision of the matrix-vector products in several flexible Krylov methods as well as for choosing the accuracy of the variable preconditioner such that the overall method is as efficient as possible. We will illustrate this computational scheme with a Schur-complement system that arises in the modeling of global ocean circulation.