Abstract:We present three methods for distributed memory parallel inverse factorization of block-sparse Hermitian positive definite matrices. The three methods are a recursive variant of the AINV inverse Cholesky algorithm, iterative refinement, and localized inverse factorization, respectively. All three methods are implemented using the Chunks and Tasks programming model, building on the distributed sparse quad-tree matrix representation and parallel matrix-matrix multiplication in the publicly available Chunks and Tasks Matrix Library (CHTML). Although the algorithms are generally applicable, this work was mainly motivated by the need for efficient and scalable inverse factorization of the basis set overlap matrix in large scale electronic structure calculations. We perform various computational tests on overlap matrices for quasi-linear Glutamic Acid-Alanine molecules and three-dimensional water clusters discretized using the standard Gaussian basis set STO-3G with up to more than 10 million basis functions. We show that for such matrices the computational cost increases only linearly with system size for all the three methods. We show both theoretically and in numerical experiments that the methods based on iterative refinement and localized inverse factorization outperform previous parallel implementations in weak scaling tests where the system size is increased in direct proportion to the number of processes. We show also that compared to the method based on pure iterative refinement the localized inverse factorization requires much less communication.

Impact of two factors on several domain decomposition based parallel incomplete factorizations for the meso-scale simulation of concrete

An Efficient Domain Decomposition Method in Multi-physics Simulation of Package Systems

Parallelization of incomplete factorization preconditioning of block tridiagonal linear systems with 1-D domain decomposition

Domain decomposition approach for parallel improvement of tetrahedral meshes.

Parallel Data Partitioning Strategy In Solving Large Scale Electromagnetic Scattering Problems

Parallel incomplete factorization preconditioning of block tridiagonal linear systems with 2-D domain decomposition

A Parallelization Technique Based on Factor Combination and Graph Partitioning for General Incomplete LU Factorization

An Adaptive Lu Factorization Algorithm For Parallel Circuit Simulation

Block iterative methods for linear algebraic equation and domain decomposition method for incompressible viscous flow

Scalable Hierarchical Parallel Algorithm for the Solution of Super Large-Scale Sparse Linear Equations

Efficient Noninteractive Outsourcing of Large-Scale QR and LU Factorizations

Parallelization and scalability analysis of inverse factorization using the Chunks and Tasks programming model

Block Size Selection of Parallel LU and QR on PVP-based and RISC-based Supercomputers

Iteration-merged Solver for Linear Banded Matrices Based on Domain Decomposition and MPI

Parallel Factorizations in Numerical Analysis

Selection of the Ideal Number of Sub-Domains for Domain Decomposition in MPP

Coupling parallel adaptive mesh refinement with a nonoverlapping domain decomposition solver

Computing Large-Scale Matrix and Tensor Decomposition with Structured Factors: A Unified Nonconvex Optimization Perspective

Accelerating Large-Scale Sparse LU Factorization for RF Circuit Simulation.

Implementation of the Domain Decomposition Method-Based Parallel Finite Element Method in Geotechnical Engineering