|[eigen] [Sparse][LU] scale factors|
[ Thread Index |
| More lists.tuxfamily.org/eigen Archives
I'm using SuperLU and UmfPack backends for Sparse matrix LU
Those backends expose the P and Q permutations, but not the scaling
(I believe that SuperLU backend gets them but does not give access,
while UmfPack backend does not read the
scaling factors , for UmfPack, one would also need acces to the
Are there some strong reasons not to provide those access ?
PS: Is there a "nice" way to convert the permutations vectors to sparse
matrices such that P* and *Q do the permutation ?
I'm well aware that this would not be the proper way to perform those
permutations, but I need thoses sparse matrices anyway
, and my current code reading the values in the vectors and setting the
1. by hand in the sparse matrix feels ugly.
PS2: I saw that Scipy patched SuperLU , and while I'm against
embedding libraries, those patches seem serious ☹