Localized bases of eigensubspaces and operator compression

Given a complex local operator, such as the generator of a Markov chain on a large network, a differential operator, or a large sparse matrix that comes from the discretization of a differential operator, we would like to find its best finite dimensional approximation with a given dimension. The answer to this question is often given simply by the projection of the original operator to its eigensubspace of the given dimension that corresponds to the smallest or largest eigenvalues, depending on the setting. The representation of such subspaces, however, is far from being unique and our interest is to find the most localized bases for these subspaces. The reduced operator using these bases would have sparsity features similar to that of the original operator. We will discuss different ways of obtaining localized bases, and we will give an explicit characterization of the decay rate of these basis functions. We will also discuss efficient numerical algorithms for finding such basis functions and the reduced (or compressed) operator.

[1]  Y. Saad,et al.  Numerical Methods for Large Eigenvalue Problems , 2011 .

[2]  E Weinan,et al.  Linear-scaling subspace-iteration algorithm with optimally localized nonorthogonal wave functions for Kohn-Sham density functional theory , 2009 .

[3]  E Weinan,et al.  Orbital minimization with localization , 2008 .

[4]  Tiejun Li,et al.  Image segmentation by using the localized subspace iteration algorithm , 2008 .

[5]  E. Weinan,et al.  A sub-linear scaling algorithm for computing the electronic structure of materials , 2007 .

[6]  Othmar Koch,et al.  Dynamical Low-Rank Approximation , 2007, SIAM J. Matrix Anal. Appl..

[7]  Petros Drineas,et al.  FAST MONTE CARLO ALGORITHMS FOR MATRICES II: COMPUTING A LOW-RANK APPROXIMATION TO A MATRIX∗ , 2004 .

[8]  Robert M. Gray,et al.  Toeplitz and Circulant Matrices: A Review , 2005, Found. Trends Commun. Inf. Theory.

[9]  Doron Levy,et al.  On Wavelet-Based Numerical Homogenization , 2005, Multiscale Model. Simul..

[10]  Jinchao Xu,et al.  On an energy minimizing basis for algebraic multigrid methods , 2004 .

[11]  R. Martin,et al.  Electronic Structure: Basic Theory and Practical Methods , 2004 .

[12]  D. Sorensen.  Numerical methods for large eigenvalue problems , 2002, Acta Numerica.

[13]  Tony F. Chan,et al.  An Energy-minimizing Interpolation for Robust Multigrid Methods , 1999, SIAM J. Sci. Comput..

[14]  S. Goedecker.  Linear scaling electronic structure methods , 1999 .

[15]  D. Muoio,et al.  Effect of sex steroids on beta-endorphin levels at rest and during submaximal treadmill exercise in anovulatory and ovulatory runners. , 1999, Fertility and sterility.

[16]  Chi-Wang Shu,et al.  On the Gibbs Phenomenon and Its Resolution , 1997, SIAM Rev..

[17]  N. Marzari,et al.  Maximally localized generalized Wannier functions for composite energy bands , 1997, cond-mat/9707145.

[18]  Thomas Y. Hou,et al.  A Multiscale Finite Element Method for Elliptic Problems in Composite Materials and Porous Media , 1997 .

[19]  V. Zhikov,et al.  Homogenization of Differential Operators and Integral Functionals , 1994 .

[20]  R. Coifman,et al.  Fast wavelet transforms and numerical algorithms I , 1991 .

[21]  C. T. Mullis,et al.  Inverse eigenvalue problem for real symmetric Toeplitz matrices , 1988, ICASSP-88., International Conference on Acoustics, Speech, and Signal Processing.

[22]  I. Babuska,et al.  Generalized Finite Element Methods: Their Performance and Their Relation to Mixed Methods , 1983 .

[23]  G. Golub.  Matrix computations , 1983 .

[24]  A. Bensoussan,et al.  Asymptotic analysis for periodic structures , 1979 .

[25]  P. Heywood.  Trigonometric Series , 1968, Nature.

[26]  John G. Kemeny,et al.  Finite Markov chains , 1960 .

[27]  Walter Kohn,et al.  Analytic Properties of Bloch Waves and Wannier Functions , 1959 .

[28]  G. Wannier.  The Structure of Electronic Excitation Levels in Insulating Crystals , 1937 .