Abstract
Multigrid methods are frequently used when solving systems of linear equations, applied either as standalone solvers or as preconditioners for iterative methods. Within each cycle, the approximation is computed using smoothing on fine levels and solving on the coarsest level.
With growth of the size of the problems that are being solved, the size of the problems on the coarsest grid is also growing and their solution can become a computational bottleneck. In practice the problems on the coarsest-grid are often solved approximately, for example by Krylov subspace methods or direct methods based on low rank approximation; see, e.g., [M. Huber, Massively parallel and fault-tolerant multigrid solvers on peta-scale systems, Ph.D. Thesis, Technical University of Munich, Germany, 2019], [Buttari et al., Numerical Linear Algebra with Applications (2021)]. The accuracy of the coarsest-grid solver is typically determined experimentally in order to balance the cost of the solves and the total number of multigrid cycles required for convergence.
In this talk, we present an approach to analyzing the effect of approximate coarsest-grid solves in the multigrid V-cycle method for symmetric positive definite problems. We discuss several stopping criteria derived based on the analysis and suggest a strategy for utilizing them in practice. The results are illustrated through numerical experiments.