#62 QR without the Q matrix

[vasek]

For some tasks (triangularization) the Q matrix is not really needed. Attached is an overloaded version of IT++ QR that does not compute Q.

h-file:
bool qr(const mat &A, mat &R);

.cpp see attached.

Comments wellcome.
Vasek

[vasek]

[Adam Piątyszek]

Logged In: YES
user_id=1004597
Originator: NO

Hi Vasek,

Could you prepare a proper patch against the qr.{cpp,h} files in IT++? Please do not forget about the complex version, some Doxygen documentation and extensions to the test program (qr_test).

Thanks in advance!
/Adam

[vasek]

[vasek]

Logged In: YES
user_id=1169711
Originator: YES

Hi Adam,
The attached patch is against rev.1508. Most of it is just a modification of the original routine "with Q". Let me know if I forgot something.
I did not do any modification of the case with pivoting. I am not using this version and I do not know if it would be meaningful.
Best
Vasek
File Added: qr-noQ.dif

[Adam Piątyszek]

Logged In: YES
user_id=1004597
Originator: NO

Thanks for the patch. Your contribution with minor clean-up changes and a small fix in the test program is now in the IT++ SVN repository (trunk, rev. 1509).

/Adam