[Fwd: QRDecomposition]

[Ron Boisvert <boisvert@nist.gov>]

Sender: Sione <sionep@xtra.co.nz>
Subject: QRDecomposition



Hi All,

Does anyone know of how to modify the  QRDecomposition  class so   to be 
able to do what  Matlab  "qr"  command does,  such as the following 
description in Matlab.

[Q,R,E] = QR(A) produces unitary Q, upper triangular R and a
permutation matrix E so that A*E = Q*R. The column permutation E is
chosen so that ABS(DIAG(R)) is decreasing.

I want to modify the  QRDecomposition  class so that it also calculates  
"E" , the permutation matrix.

Any help would be appreciated.

Cheers,
Sione.