[SciPy-User] permutation matrix in LU decomposition

Jason Grout jason-sage at creativetrax.com
Mon Nov 26 12:21:37 EST 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