Incomplete selected inversion for linear-scaling electronic structure calculations

01/17/2020
by   Simon Etter, et al.
0

Pole Expansion and Selected Inversion (PEXSI) is an efficient scheme for evaluating selected entries of functions of large sparse matrices as required e.g. in electronic structure algorithms. We show that the triangular factorizations computed by the PEXSI scheme exhibit a localization property similar to that of matrix functions, and we present a modified PEXSI algorithm which exploits this observation to achieve linear scaling. To the best of our knowledge, the resulting incomplete PEXSI (iPEXSI) algorithm is the first linear-scaling algorithm which scales provably better than cubically even in the absence of localization, and we hope that this will help to further lower the critical system size where linear-scaling algorithms begin to outperform the diagonalization algorithm.

READ FULL TEXT
research
04/09/2016

A Left-Looking Selected Inversion Algorithm and Task Parallelism on Shared Memory Systems

Given a sparse matrix A, the selected inversion algorithm is an efficien...
research
08/14/2017

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

This paper generalizes the parallel selected inversion algorithm called ...
research
01/08/2021

On the Recognition of Strong-Robinsonian Incomplete Matrices

A matrix is incomplete when some of its entries are missing. A Robinson ...
research
11/25/2020

Quantum algorithms for matrix scaling and matrix balancing

Matrix scaling and matrix balancing are two basic linear-algebraic probl...
research
04/16/2021

Enabling Electronic Structure-Based Ab-Initio Molecular Dynamics Simulations with Hundreds of Millions of Atoms

We push the boundaries of electronic structure-based ab-initio molecular...
research
12/15/2014

Finding a sparse vector in a subspace: Linear sparsity using alternating directions

Is it possible to find the sparsest vector (direction) in a generic subs...

Please sign up or login with your details

Forgot password? Click here to reset