[Numpy-discussion] linalg.qr versus scilab qr

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


Hello,

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 ?

	Romaric


More information about the Numpy-discussion mailing list