[Numpy-discussion] linalg.qr versus scilab qr
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