[Numpy-discussion] Confused with qr decomposition function

Virgil Stokes vs at it.uu.se
Mon Nov 19 18:17:53 EST 2012

I am using the latest versions of numpy (from 
numpy-1.7.0b2-win32-superpack-python2.7.exe) and scipy (from 
scipy-0.11.0-win32-superpack-python2.7.exe ) on a windows 7 (32-bit) 

I have used

import numpy as np
q,r = np.linalg.qr(A)

and compared the results obtained  from MATLAB (R2010B)

[q,r] = qr(A)

The q,r returned from numpy are both the negative of theq,r returned 
from MATLAB for the same matrix A. I believe that the q,r returned from 
MATLAB are correct. Why am I getting their negative from numpy?

Note, I have tried this on several different matrices ---  numpy always 
gives the negative of MATLAB's q,r values.
[I mistakenly have already sent a similar email to the scipy list --- 
please excuse this mistake.]

More information about the NumPy-Discussion mailing list