Ask Your Question

LU decomposition for a matrix

asked 2015-05-04 15:59:21 +0200

rafarob32 gravatar image

updated 2023-01-09 23:59:37 +0200

tmonteil gravatar image

There is an inconsistency in the use of the LU decomposition of a matrix. If matrix A is declared in the RDF ring then P, L, U = A.LU () provides PA = LU but if A is declared in the RR ring (or QQ, ZZ,...) then A = PLU. Therefore in a case P is P.transpose () on the other.

edit retag flag offensive close merge delete

1 Answer

Sort by ยป oldest newest most voted

answered 2015-05-04 17:05:12 +0200

tmonteil gravatar image

updated 2015-05-04 17:24:09 +0200

While this is correctly documented (hence not a bug), i agree that there is some inconsistency here. Thanks for reporting, it is now trac ticket 18365.

edit flag offensive delete link more

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools



Asked: 2015-05-04 15:59:21 +0200

Seen: 2,712 times

Last updated: May 05 '15