The fastest way to check if a sparse matrix is singular or not?

26 views (last 30 days)
Dear All,
For a large sparse matrix A, how can I check in a fastest way if matrix A is singluar or not? The matrix has a dimension of 6000 x 6000. I know lu is fast, but I need to repeat this calculations many times for different matrix A.
Thanks a lot in advance.
Beson

Accepted Answer

Matt J
Matt J on 12 Apr 2019
condest(A)
  2 Comments
Benson Gou
Benson Gou on 13 Apr 2019
@Matt. Thanks for your reply. But I think there must be a threshold to judge if martrix A is singular or not. How to decide the threshold? Thanks.
Matt J
Matt J on 13 Apr 2019
Edited: Matt J on 13 Apr 2019
The condition number calculated by condest measures how much error magnification you get when you try to invert the matrix when solving a linear equation
The "threshold" that you are pursuing would be determined from the amount of error magnification that your specific application can tolerate (note that it is not possible to have zero error magnification). If the error magnification is greater than what you can tolerate, you would consider the matrix "singular".

Sign in to comment.

More Answers (0)

Categories

Find more on Sparse Matrices in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!