PSelInv -- A Distributed Memory Parallel Algorithm for Selected Inversion : the Symmetric Case

Abstract

We describe an efficient parallel implementation of the selected inversion algorithm for distributed memory computer systems, which we call \texttt{PSelInv}. The \texttt{PSelInv} method computes selected elements of a general sparse matrix $A$ that can be decomposed as $A=LU$, where $L$ is lower triangular and $U$ is upper triangular. The implementation described in this paper focuses on the case of sparse symmetric matrices. It contains an interface that is compatible with the distributed memory parallel sparse direct factorization \texttt{SuperLU\_DIST}. However, the underlying data structure and design of \texttt{PSelInv} allows it to be easily combined with other factorization routines such as \texttt{PARDISO}. We discuss general parallelization strategies such as data and task distribution schemes. In particular, we describe how to exploit the concurrency exposed by the elimination tree associated with the $LU$ factorization of $A$. We demonstrate the efficiency and accuracy of \texttt{PSelInv} by presenting a number of numerical experiments. In particular, we show that \texttt{PSelInv} can run efficiently on more than $4,000$ cores for a modestly sized matrix. We also demonstrate how \texttt{PSelInv} can be used to accelerate large-scale electronic structure calculations.

References in corpus (9)

• The SIESTA method for ab initio order-N materials simulation
• Efficient index handling of multidimensional periodic boundary conditions
• Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework I: Total energy calculation
• Accelerating Atomic Orbital-based Electronic Structure Calculation via Pole Expansion and Selected Inversion
• Pole-based approximation of Fermi-Dirac function
• Solving Fluctuation-Enhanced Poisson-Boltzmann Equations
• Extension and optimization of the FIND algorithm: computing Green's and less-than Green's functions (with technical appendix)
• A Nested Dissection Approach to Modeling Transport in Nanodevices: Algorithms and Applications
• Enhancing the scalability and load balancing of the parallel selected inversion algorithm via tree-based asynchronous communication

Cited by in corpus (3)

• ELSI: A Unified Software Interface for Kohn-Sham Electronic Structure Solvers
• ELSI -- An Open Infrastructure for Electronic Structure Solvers
• Enhancing the scalability and load balancing of the parallel selected inversion algorithm via tree-based asynchronous communication