Ask Your Question

LU decomposition for a matrix

asked 2015-05-04 08:59:21 -0600

rafarob32 gravatar image

updated 2017-01-05 15:17:17 -0600

FrédéricC 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 10:05:12 -0600

tmonteil gravatar image

updated 2015-05-04 10:24:09 -0600

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

1 follower


Asked: 2015-05-04 08:59:21 -0600

Seen: 392 times

Last updated: May 04 '15