Hello,
the routine gsl_linalg_QRPT_update appears to be also incorrectly documented;
in this case I think the comments in the source are also partially incorrect.
My own testing implies that the update formula is
Q'R' = Q ( R + w v^T ) P
and not Q'R' = Q R + w v^T, as the 1.10 manual presently says, or
Q'R' = Q ( R + w v^T P), as mentioned in the qrpt source file.
Regards
--
Dr. Reinhold Bader
Leibniz-Rechenzentrum, Abt. Hochleistungssysteme | Tel. +49 89 35831 8825
Boltzmannstr. 1, 85748 Garching | Fax +49 89 35831 9700
