skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Parallel accelerated cyclic reduction preconditioner for three-dimensional elliptic PDEs with variable coefficients

Journal Article · · Journal of Computational and Applied Mathematics

We present a robust and scalable preconditioner for the solution of large-scale linear systems that arise from the discretization of elliptic PDEs amenable to rank compression. The preconditioner is based on hierarchical low-rank approximations and the cyclic reduction method. The setup and application phases of the preconditioner achieve log-linear complexity in memory footprint and number of operations, and numerical experiments exhibit good weak and strong scalability at large processor counts in a distributed memory environment. Numerical experiments with linear systems that feature symmetry and nonsymmetry, definiteness and indefiniteness, constant and variable coefficients demonstrate the preconditioner applicability and robustness. Furthermore, it is possible to control the number of iterations via the accuracy threshold of the hierarchical matrix approximations and their arithmetic operations, and the tuning of the admissibility condition parameter. Together, these parameters allow for optimization of the memory requirements and performance of the preconditioner.

Research Organization:
Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)
Sponsoring Organization:
USDOE Office of Science (SC)
Grant/Contract Number:
AC02-05CH11231
OSTI ID:
1580345
Journal Information:
Journal of Computational and Applied Mathematics, Vol. 344, Issue C; ISSN 0377-0427
Publisher:
ElsevierCopyright Statement
Country of Publication:
United States
Language:
English
Citation Metrics:
Cited by: 4 works
Citation information provided by
Web of Science

References (36)

A Fast Direct Solution of Poisson's Equation Using Fourier Analysis journal January 1965
Parallel -Matrix Arithmetics on Shared Memory Systems journal December 2004
$${{\fancyscript{H}}} $$ H -LU factorization on many-core systems journal June 2013
On the Numerical Rank of the Off-Diagonal Blocks of Schur Complements of Discretized Elliptic PDEs journal January 2010
The Multifrontal Solution of Indefinite Sparse Symmetric Linear journal September 1983
A bibliography on semiseparable matrices* journal December 2005
A Fast $ULV$ Decomposition Solver for Hierarchically Semiseparable Representations journal January 2006
Superfast Multifrontal Method for Large Structured Linear Systems of Equations journal January 2010
Fast algorithms for hierarchically semiseparable matrices journal November 2010
Randomized Sparse Direct Solvers journal January 2013
An Efficient Multicore Implementation of a Novel HSS-Structured Multifrontal Solver Using Randomized Sampling journal January 2016
A Parallel Geometric Multifrontal Solver Using Hierarchically Semiseparable Structure journal May 2016
An $$\mathcal O (N \log N)$$ O ( N log N )   Fast Direct Solver for Partial Hierarchically Semi-Separable Matrices: With Application to Radial Basis Function Interpolation journal April 2013
A fast, memory efficient and robust sparse preconditioner based on a multifrontal approach with applications to finite‐element matrices journal February 2016
Improving Multifrontal Methods by Means of Block Low-Rank Representations journal January 2015
Hierarchical Cholesky decomposition of sparse matrices arising from curl–curl-equation journal January 2007
Parallel black box $$\mathcal {H}$$ -LU preconditioning for elliptic boundary value problems journal April 2008
Domain decomposition based $${\mathcal H}$$ -LU preconditioning journal March 2009
Robust Approximate Cholesky Factorization of Rank-Structured Symmetric Positive Definite Matrices journal January 2010
A Sparse Matrix Arithmetic Based on $\Cal H$ -Matrices. Part I: Introduction to ${\Cal H}$ -Matrices journal April 1999
On Direct Methods for Solving Poisson’s Equations journal December 1970
Parallel Distributed FFT-Based Solvers for 3-D Poisson Problems in Meso-Scale Atmospheric Simulations journal February 2001
Sweeping preconditioner for the Helmholtz equation: Hierarchical matrix representation journal February 2011
A Parallel Sweeping Preconditioner for Heterogeneous 3D Helmholtz Equations journal January 2013
Preconditioning by incomplete block cyclic reduction journal May 1984
On the Approximate Cyclic Reduction Preconditioner journal January 1999
A Cyclic Reduction Approach to the Numerical Solution of Boundary Value ODEs journal January 1997
A Parallel Algorithm for Solving Tridiagonal Linear Systems on Distributed-Memory Multiprocessors journal September 1994
Predicting Performance and Selecting Modes of Parallelism: A Case Study Using Cyclic Reduction on Three Parallel Machines journal November 1993
A Parallel and Vector Variant of the Cyclic Reduction Algorithm journal July 1988
Cyclic Reduction Tridiagonal Solvers on GPUs Applied to Mixed-Precision Multigrid journal January 2011
Methods of conjugate gradients for solving linear systems journal December 1952
A generic grid interface for parallel and adaptive scientific computing. Part I: abstract framework journal April 2008
SuperLU_DIST: A scalable distributed-memory sparse direct solver for unsymmetric linear systems journal June 2003
PetIGA: A framework for high-performance isogeometric analysis journal August 2016
Hierarchical Matrices book January 2015

Cited By (2)

Distributed-memory lattice H-matrix factorization journal January 2019
Hierarchical algorithms on hierarchical architectures
  • Keyes, D. E.; Ltaief, H.; Turkiyyah, G.
  • Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, Vol. 378, Issue 2166 https://doi.org/10.1098/rsta.2019.0055
journal January 2020

Similar Records

Massively parallel solvers for elliptic partial differential equations in numerical weather and climate prediction
Journal Article · Thu Jan 16 00:00:00 EST 2014 · Quarterly Journal of the Royal Meteorological Society · OSTI ID:1580345

An Algebraic Multilevel Preconditioner with Low-Rank Corrections for Sparse Symmetric Matrices
Journal Article · Thu Mar 03 00:00:00 EST 2016 · SIAM Journal on Matrix Analysis and Applications · OSTI ID:1580345

A fast, memory efficient and robust sparse preconditioner based on a multifrontal approach with applications to finite‐element matrices
Journal Article · Tue Feb 02 00:00:00 EST 2016 · International Journal for Numerical Methods in Engineering · OSTI ID:1580345