lobpcg.m
lobpcg.m, MATLAB implementation of the locally optimal block preconditioned conjugate gradient method: Toward the optimal preconditioned eigensolver: Locally optimal block preconditioned conjugate gradient method. We describe new algorithms of the locally optimal block preconditioned conjugate gradient (LOBPCG) method for symmetric eigenvalue problems, based on a local optimization of a three-term recurrence, and suggest several other new methods. To be able to compare numerically different methods in the class, with different preconditioners, we propose a common system of model tests, using random preconditioners and initial guesses. As the “ideal” control algorithm, we advocate the standard preconditioned conjugate gradient method for finding an eigenvector as an element of the null-space of the corresponding homogeneous system of linear equations under the assumption that the eigenvalue is known. We recommend that every new preconditioned eigensolver be compared with this “ideal” algorithm on our model test problems in terms of the speed of convergence, costs of every iteration, and memory requirements. We provide such comparison for our LOBPCG method. Numerical results establish that our algorithm is practically as efficient as the “ideal” algorithm when the same preconditioner is used in both methods. We also show numerically that the LOBPCG method provides approximations to first eigenpairs of about the same quality as those by the much more expensive global optimization method on the same generalized block Krylov subspace. We propose a new version of block Davidson’s method as a generalization of the LOBPCG method. Finally, direct numerical comparisons with the Jacobi-Davidson method show that our method is more robust and converges almost two times faster.
Keywords for this software
References in zbMATH (referenced in 144 articles , 1 standard article )
Showing results 1 to 20 of 144.
Sorted by year (- Carson, Erin; Lund, Kathryn; Rozložník, Miroslav; Thomas, Stephen: Block Gram-Schmidt algorithms and their stability properties (2022)
- Eisenmann, Henrik; Nakatsukasa, Yuji: Solving two-parameter eigenvalue problems using an alternating method (2022)
- Ma, Xijun; Shen, Chungen; Wang, Li; Zhang, Lei-Hong; Li, Ren-Cang: A self-consistent-field iteration for MAXBET with an application to multi-view feature extraction (2022)
- Nigam, Nilima; Pollock, Sara: A simple extrapolation method for clustered eigenvalues (2022)
- Rontsis, Nikitas; Goulart, Paul; Nakatsukasa, Yuji: Efficient semidefinite programming with approximate ADMM (2022)
- Wang, Wei; Zhang, Zhimin: Spectral element methods for eigenvalue problems based on domain decomposition (2022)
- Yin, Jianyuan; Zhang, Lei; Zhang, Pingwen: Solution landscape of the Onsager model identifies non-axisymmetric critical points (2022)
- An, Dong; Lin, Lin; Xu, Ze: Split representation of adaptively compressed polarizability operator (2021)
- Giani, Stefano; Grubišić, Luka; Heltai, Luca; Mulita, Ornela: Smoothed-adaptive perturbed inverse iteration for elliptic eigenvalue problems (2021)
- Guarín-Zapata, Nicolás; Gomez, Juan; Hadjesfandiari, Ali Reza; Dargush, Gary F.: Variational principles and finite element Bloch analysis in couple stress elastodynamics (2021)
- Jolivet, Pierre; Roman, Jose E.; Zampini, Stefano: KSPHPDDM and PCHPDDM: extending PETSc with advanced Krylov methods and robust multilevel overlapping Schwarz preconditioners (2021)
- Krumnow, Christian; Pfeffer, Max; Uschmajew, André: Computing eigenspaces with low rank constraints (2021)
- Li, Ji; Cai, Jian-Feng; Zhao, Hongkai: Scalable incremental nonconvex optimization approach for phase retrieval (2021)
- Li, Tianyi; Abdelmoula, Radhi: Gradient damage analysis of a cylinder under torsion: bifurcation and size effects (2021)
- Miao, Cun-Qiang; Wu, Wen-Ting: On relaxed filtered Krylov subspace method for non-symmetric eigenvalue problems (2021)
- Pollock, Sara; Scott, L. Ridgway: Extrapolating the Arnoldi algorithm to improve eigenvector convergence (2021)
- Williamson, Kevin; Cho, Heyrim; Sousedík, Bedřich: Application of adaptive ANOVA and reduced basis methods to the stochastic Stokes-Brinkman problem (2021)
- Ezvan, Olivier; Zeng, Xiaoshu; Ghanem, Roger; Gencturk, Bora: Multiscale modal analysis of fully-loaded spent nuclear fuel canisters (2020)
- Ferrari, Federico; Sigmund, Ole: Towards solving large-scale topology optimization problems with buckling constraints at the cost of linear analyses (2020)
- Kelley, C. T.; Bernholc, J.; Briggs, E. L.; Hamilton, Steven; Lin, Lin; Yang, Chao: Mesh independence of the generalized Davidson algorithm (2020)