How can the Cholesky decomposition step in eigs() be avoided without passing a matrix to eigs that is a Cholesky decomposition?

Question

Allo 2012-1-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/26812-how-can-the-cholesky-decomposition-step-in-eigs-be-avoided-without-passing-a-matrix-to-eigs-that-i

回答： Andrew Knyazev 2015-5-15

Hello,

I have been looking at the following set of notes:

http://people.inf.ethz.ch/arbenz/ewp/Lnotes/chapter10.pdf

and specifically this quote in those notes:

"If SIGMA is a real or complex scalar including 0, EIGS finds the eigenvalues closest to SIGMA. For scalar SIGMA, and when SIGMA = ’SM’, B need only be symmetric (or Hermitian) positive semi-definite since it is not Cholesky factored as in the other cases."

I have a Hermitian positive-semidefinite matrix A, of which I want to find the 3 smallest eigenvalues. The Cholesky-decomposition is too memory intensive for the matrices I am working with. Please, is there a way to use eigs() without having to perform the Cholesky decomposition either in eigs() or outside of it?

Thank you very much.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Walter Roberson 2012-1-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/26812-how-can-the-cholesky-decomposition-step-in-eigs-be-avoided-without-passing-a-matrix-to-eigs-that-i#answer_34909

在 MATLAB Online 中打开

Try

eigs(YourArray, 3, 'SM')

However, note that this requires that you be seeking the 3 eigenvalues with smallest absolute magnitude. If you need to find the smallest magnitude (e.g., -11.49 being smaller than -1.149) then you will not be able to use this option.

2 个评论
显示无隐藏无

Allo 2012-1-23

Hello and thank you for responding.

I tried this

D=eigs(M,6,'SM');

Because M=S*S, where S is my original Hermitian but not positive semi-definite matrix, because I wanted to submit a Hermitian positive semi-definite matrix to the eigs() function for the sake of exploiting the information I found above. It doesn't seem to be yielding the decrease in memory consumption that I expected. I thought the Cholesky decomposition of eigs() was supposed to require far more memory than the rest of eigs(), so I'm a bit confused.

Allo 2012-1-23