Alternating Anderson-Richardson method: an efficient alternative to preconditioned Krylov methods for large, sparse linear systems

From MaRDI portal
Publication:6043327

DOI10.1016/J.CPC.2018.07.007arXiv1606.08740WikidataQ129511157 ScholiaQ129511157MaRDI QIDQ6043327FDOQ6043327


Authors: Phanish Suryanarayana, Phanisri P. Pratapa, John E. Pask Edit this on Wikidata


Publication date: 5 May 2023

Published in: Computer Physics Communications (Search for Journal in Brave)

Abstract: We present the Alternating Anderson-Richardson (AAR) method: an efficient and scalable alternative to preconditioned Krylov solvers for the solution of large, sparse linear systems on high performance computing platforms. Specifically, we generalize the recently proposed Alternating Anderson-Jacobi (AAJ) method (Pratapa et al., J. Comput. Phys. (2016), 306, 43--54) to include preconditioning, discuss efficient parallel implementation, and provide serial MATLAB and parallel C/C++ implementations. In serial applications to nonsymmetric systems, we find that AAR is comparably robust to GMRES, using the same preconditioning, while often outperforming it in time to solution; and find AAR to be more robust than Bi-CGSTAB for the problems considered. In parallel applications to the Helmholtz and Poisson equations, we find that AAR shows superior strong and weak scaling to GMRES, Bi-CGSTAB, and Conjugate Gradient (CG) methods, using the same preconditioning, with consistently shorter times to solution at larger processor counts. Finally, in massively parallel applications to the Poisson equation, on up to 110,592 processors, we find that AAR shows superior strong and weak scaling to CG, with shorter minimum time to solution. We thus find that AAR offers a robust and efficient alternative to current state-of-the-art solvers, with increasing advantages as the number of processors grows.


Full work available at URL: https://arxiv.org/abs/1606.08740







Cites Work


Cited In (5)





This page was built for publication: Alternating Anderson-Richardson method: an efficient alternative to preconditioned Krylov methods for large, sparse linear systems

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6043327)