eigs with sparse matrices

6 次查看(过去 30 天)
Liz
Liz 2011-7-13
Hi all,
I'm working with very sparse positive-definite symmetric matrices, and I would like to solve for the eigenvectors/values (I use eigs of course) . Unfortunately, I cannot solve for the eigenvectors for many of the matrices I construct. I get an ARPACK error from the routine dneupd saying that no eigenvalues could be found to sufficient accuracy.
I was wondering two things:
1. is there a more reliable way (i.e. doesn't exit with errors) to estimate eigenvalues/eigenvectors with very sparse matrices?
2. is there a faster way to tell if you can solve for a matrix's eigenvalues? It takes a long time for eigs to decide that it cannot find any eigenvalues.
I would appreciate any pointers,tips, and/or advice. Thanks in advance.

回答(2 个)

the cyclist
the cyclist 2011-7-13
It would be really helpful if you could provide an example that gives the error.
If the error is related to the condition number (singular matrix), then you may be able to use the cond() function.

Andrew Knyazev
Andrew Knyazev 2011-8-1

类别

Help CenterFile Exchange 中查找有关 Linear Algebra 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by