[Numpy-discussion] linalg.qr versus scilab qr

R. David david at icps.u-strasbg.fr
Mon Jan 8 04:09:59 EST 2007


In scilab, the QR decomposition can produce  "sorted" Q and R matrices, where
diagonal elements of R are decreasing. Thus, a permutation vector E is returned,
so that X*E = Q*R.

Is there a way to compute such a permutation QR in numpy ? It seems that
numpy.linalg.qr does not sort the diagonal R elements.

Do I have to wrap around lapack's zgeqpf or zgeqp3 ?


More information about the NumPy-Discussion mailing list