[SciPy-dev] NumPy_for_Matlab_Addicts

[Nils Wagner]

Hi,

The QR decomposition is already available in scipy.

qr(a, overwrite_a=0, lwork=None)
     QR decomposition of an M x N matrix a.

     Description:

       Find a unitary matrix, q, and an upper-trapezoidal 
matrix r
       such that q * r = a

     Inputs:

       a -- the matrix
       overwrite_a=0 -- if non-zero then discard the 
contents of a,
                      i.e. a is used as a work array if 
possible.

       lwork=None -- >= shape(a)[1]. If None (or -1) 
compute optimal
                     work array size.

     Outputs:

       q, r -- matrices such that q * r = a

Nils