Abstract:On current computer architectures, GMRES' performance can be limited by its communication cost to generate orthonormal basis vectors of the Krylov subspace. To address this performance bottleneck, its $s$-step variant orthogonalizes a block of $s$ basis vectors at a time, potentially reducing the communication cost by a factor of $s$. Unfortunately, for a large step size $s$, the solver can generate extremely ill-conditioned basis vectors, and to maintain stability in practice, a conservatively small step size is used, which limits the performance of the $s$-step solver. To enhance the performance using a small step size, in this paper, we introduce a two-stage block orthogonalization scheme. Similar to the original scheme, the first stage of the proposed method operates on a block of $s$ basis vectors at a time, but its objective is to maintain the well-conditioning of the generated basis vectors with a lower cost. The orthogonalization of the basis vectors is delayed until the second stage when enough basis vectors are generated to obtain higher performance. Our analysis shows the stability of the proposed two-stage scheme. The performance is improved because while the same amount of computation as the original scheme is required, most of the communication is done at the second stage of the proposed scheme, reducing the overall communication requirements. Our performance results with up to 192 NVIDIA V100 GPUs on the Summit supercomputer demonstrate that when solving a 2D Laplace problem, the two-stage approach can reduce the orthogonalization time and the total time-to-solution by the respective factors of up to $2.6\times$ and $1.6\times$ over the original $s$-step GMRES, which had already obtained the respective speedups of $2.1\times$ and $1.8\times$ over the standard GMRES. Similar speedups were obtained for 3D problems and for matrices from the SuiteSparse Matrix Collection.

Near-optimal convergence of the full orthogonalization method

Global FOM and GMRES algorithms for a class of complex matrix equations.

Two-Stage Block Orthogonalization to Improve Performance of $s$-step GMRES

Restarted weighted full orthogonalization method for shifted linear systems

Global Convergence of Triangularized Orthogonalization-free Method

Optimal Solutions of Well-Posed Linear Systems via Low-Precision Right-Preconditioned GMRES with Forward and Backward Stabilization

A Kaczmarz-Inspired Method for Orthogonalization

Orthogonal Directions Constrained Gradient Method: from non-linear equality constraints to Stiefel manifold

On IGMRES: an Incomplete Generalized Minimal Residual Method for Large Unsymmetric Linear Systems

On IOM(q): the Incomplete Orthogonalization Method for Large Unsymmetric Linear Systems

Guaranteed $ε$-optimal solutions with the linear optimizer ART3+O

Improving the Performance of the GMRES Method using Mixed-Precision Techniques

Certified Multi-Fidelity Zeroth-Order Optimization

Parallelizable Algorithms for Optimization Problems with Orthogonality Constraints

Optimal Convergence Rates for the Orthogonal Greedy Algorithm

A Nonlinear GMRES Optimization Algorithm for Canonical Tensor Decomposition

Convergence of Generalized AOR Iterative Method for Linear Systems with Strictly Diagonally Dominant Matrices

On the backward stability of s-step GMRES

Analyzing Vector Orthogonalization Algorithms

Gl-QFOM and Gl-QGMRES: two efficient algorithms for quaternion linear systems with multiple right-hand sides