[Numpy-discussion] linalg.qr versus scilab qr

[R. David]

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