Show simple item record

dc.contributor.authorChavez Chavez, Gustavo Ivan
dc.contributor.authorTurkiyyah, George
dc.contributor.authorZampini, Stefano
dc.contributor.authorLtaief, Hatem
dc.contributor.authorKeyes, David E.
dc.date.accessioned2017-12-21T13:57:03Z
dc.date.available2017-12-21T13:57:03Z
dc.date.issued2017-12-15
dc.identifier.citationChávez G, Turkiyyah G, Zampini S, Ltaief H, Keyes D (2017) Accelerated Cyclic Reduction: A Distributed-Memory Fast Solver for Structured Linear Systems. Parallel Computing. Available: http://dx.doi.org/10.1016/j.parco.2017.12.001.
dc.identifier.issn0167-8191
dc.identifier.doi10.1016/j.parco.2017.12.001
dc.identifier.urihttp://hdl.handle.net/10754/626403
dc.description.abstractWe present Accelerated Cyclic Reduction (ACR), a distributed-memory fast solver for rank-compressible block tridiagonal linear systems arising from the discretization of elliptic operators, developed here for three dimensions. Algorithmic synergies between Cyclic Reduction and hierarchical matrix arithmetic operations result in a solver that has O(kNlogN(logN+k2)) arithmetic complexity and O(k Nlog N) memory footprint, where N is the number of degrees of freedom and k is the rank of a block in the hierarchical approximation, and which exhibits substantial concurrency. We provide a baseline for performance and applicability by comparing with the multifrontal method with and without hierarchical semi-separable matrices, with algebraic multigrid and with the classic cyclic reduction method. Over a set of large-scale elliptic systems with features of nonsymmetry and indefiniteness, the robustness of the direct solvers extends beyond that of the multigrid solver, and relative to the multifrontal approach ACR has lower or comparable execution time and size of the factors, with substantially lower numerical ranks. ACR exhibits good strong and weak scaling in a distributed context and, as with any direct solver, is advantageous for problems that require the solution of multiple right-hand sides. Numerical experiments show that the rank k patterns are of O(1) for the Poisson equation and of O(n) for the indefinite Helmholtz equation. The solver is ideal in situations where low-accuracy solutions are sufficient, or otherwise as a preconditioner within an iterative method.
dc.description.sponsorshipWe thank the anonymous reviewers for their detailed comments and suggestions for this manuscript. The authors would also like to thank Ronald Kriemann from the Max-Planck-Institute for Mathematics in the Sciences for development and continuous support of HLibPro, Alexander Litvinenko from the King Abdullah University of Science and Technology (KAUST) for the enlightening discussions and advice, and Pieter Ghysels from the Lawrence Berkeley National Laboratory for his recommendations on the use of STRUMPACK. Support from the KAUST Supercomputing Laboratory and access to Shaheen is gratefully acknowledged. The work of all authors was supported by the Extreme Computing Research Center at KAUST.
dc.publisherElsevier BV
dc.relation.urlhttp://www.sciencedirect.com/science/article/pii/S0167819117302041
dc.rightsNOTICE: this is the author’s version of a work that was accepted for publication in Parallel Computing. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. A definitive version was subsequently published in Parallel Computing, [, , (2017-12-15)] DOI: 10.1016/j.parco.2017.12.001 . © 2017. This manuscript version is made available under the CC-BY-NC-ND 4.0 license http://creativecommons.org/licenses/by-nc-nd/4.0/
dc.subjectCyclic reduction
dc.subjectHierarchical matrices
dc.subjectFast direct solvers
dc.subjectElliptic equations
dc.titleAccelerated Cyclic Reduction: A Distributed-Memory Fast Solver for Structured Linear Systems
dc.typeArticle
dc.contributor.departmentApplied Mathematics and Computational Science Program
dc.contributor.departmentComputer Science Program
dc.contributor.departmentComputer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division
dc.contributor.departmentExtreme Computing Research Center
dc.identifier.journalParallel Computing
dc.eprint.versionPost-print
dc.contributor.institutionDepartment of Computer Science, American University of Beirut (AUB), Beirut, Lebanon
dc.identifier.arxivid1701.00182
kaust.personChavez Chavez, Gustavo Ivan
kaust.personZampini, Stefano
kaust.personLtaief, Hatem
kaust.personKeyes, David E.
refterms.dateFOA2019-12-15T00:00:00Z
dc.date.published-online2017-12-15
dc.date.published-print2018-05
dc.date.posted2017-01-01


Files in this item

Thumbnail
Name:
1-s2.0-S0167819117302041-main.pdf
Size:
1.585Mb
Format:
PDF
Description:
Accepted Manuscript

This item appears in the following Collection(s)

Show simple item record