This forum has been archived. All content is frozen. Please use KDE Discuss instead.

UmfPack Support: factorization fails if matrix is singular

Tags: None
(comma "," separated)
joaoruileal
Registered Member
Posts
18
Karma
0
OS
Hello,

UmfPackLU factorization fails when the matrix is singular, which might not be a problem when all it is needed are the LU matrices.
If this factorization is not used for the solution of a system then everything is OK.

Perhaps one could have access to the last status code returned by umfpack_numeric to check if the status is:
UMFPACK_WARNING_singular_matrix: Numeric factorization was successful, but the matrix is singular. (from the manual)

Any chance of adding this to UmfPackLU?

Best regards,
João Leal
User avatar
ggael
Moderator
Posts
3447
Karma
19
OS
Why not. Feel free to propose a patch doing so.
joaoruileal
Registered Member
Posts
18
Karma
0
OS
I implemented the changes in order to have access to the return status of the factorization.
I've also added access to the control vector used by umfpack in case anyone would like to customize some options.

Here is the code:
http://pastebin.com/dFivurz1
joaoruileal
Registered Member
Posts
18
Karma
0
OS
Is this the appropriate place to propose a patch ???
User avatar
ggael
Moderator
Posts
3447
Karma
19
OS
It is better to put it on the bug tracker (eigen.tuxfamily.org/bz/): it will become more visible to other developers and there is no risk we forgot about it.


Bookmarks



Who is online

Registered users: Bing [Bot], Google [Bot], Yahoo [Bot]