Recursive Krylov‐based multigrid cycles

Abstract

AbstractWe consider multigrid (MG) cycles based on the recursive use of a two‐grid method, in which the coarse‐grid system is solved by µ⩾1 steps of a Krylov subspace iterative method. The approach is further extended by allowing such inner iterations only at the levels of given multiplicity, whereas V‐cycle formulation is used at all other levels. For symmetric positive definite systems and symmetric MG schemes, we consider a flexible (or generalized) conjugate gradient method as Krylov subspace solver for both inner and outer iterations. Then, based on some algebraic (block matrix) properties of the V‐cycle MG viewed as a preconditioner, we show that the method can have optimal convergence properties if µ is chosen to be sufficiently large. We also formulate conditions that guarantee both, optimal complexity and convergence, bounded independently of the number of levels. Our analysis shows that the method is, at least, as effective as the standard W‐cycle, whereas numerical results illustrate that it can be much faster than the latter, and actually more robust than predicted by the theory. Copyright © 2007 John Wiley & Sons, Ltd.