[SciPy-User] permutation matrix in LU decomposition

Jason Grout jason-sage@creativetrax....
Mon Nov 26 11:21:37 CST 2012


I just noticed that the permutation matrix P returned in scipy's LU 
decomposition is the inverse (i.e., transpose) of the permutation matrix 
returned from Matlab's LU function.  In scipy, A = PLU, but in Matlab, 
PA=LU [1].  It's just a convention, of course, but perhaps a note about 
this could be added to the relevant table row of 
http://www.scipy.org/NumPy_for_Matlab_Users?

Thanks,

Jason

[1] http://www.mathworks.com/help/matlab/ref/lu.html


More information about the SciPy-User mailing list