|
Locally Optimal Block Preconditioned Conjugate Gradient (LOBPCG) is a matrix-free method for finding the largest (or smallest) eigenvalues and the corresponding eigenvectors of a symmetric positive definite generalized eigenvalue problem : for a given pair of complex Hermitian or real symmetric matrices, where the matrix is also assumed positive-definite. ==Algorithm== The method performs an iterative maximization (or minimization) of the generalized Rayleigh quotient : which results in finding largest (or smallest) eigenpairs of The direction of the steepest ascent, which is the gradient, of the generalized Rayleigh quotient is positively proportional to the vector : called the eigenvector residual. If a preconditioner is available, it is applied to the residual giving vector : called the preconditioned residual. Without preconditioning, we set and so . An iterative method : or, in short, : : : is known as preconditioned steepest ascent (or descent), where the scalar is called the step size. The optimal step size can be determined by maximizing the Rayleigh quotient, i.e., : (use in case of minimizing). The maximization/minimization of the Rayleigh quotient in a 3-dimensional subspace can be performed numerically by the Rayleigh–Ritz method. As the iterations converge, the vectors and become nearly linearly dependent, making the Rayleigh–Ritz method numerically unstable in the presence of round-off errors. It is possible to substitute the vector with an explicitly computed difference making the Rayleigh–Ritz method more stable; see. This is a single-vector version of the LOBPCG method. It is one of possible generalization of the preconditioned conjugate gradient linear solvers to the case of symmetric eigenvalue problems. Even in the trivial case and the resulting approximation with will be different from that obtained by the Lanczos algorithm, although both approximations will belong to the same Krylov subspace. Iterating several approximate eigenvectors together in a block in a similar locally optimal fashion, gives the full block version of the LOBPCG. It allows robust computation of eigenvectors corresponding to nearly-multiple eigenvalues. 抄文引用元・出典: フリー百科事典『 ウィキペディア(Wikipedia)』 ■ウィキペディアで「LOBPCG」の詳細全文を読む スポンサード リンク
|