This is the mail archive of the
gsl-discuss@sourceware.org
mailing list for the GSL project.
Documentation (2)
- From: bader-reinhold at t-online dot de (Reinhold Bader)
- To: gsl-discuss at sourceware dot org
- Cc: Brian Gough <bjg at network-theory dot co dot uk>
- Date: Thu, 21 Aug 2008 21:38:58 +0200
- Subject: Documentation (2)
- Reply-to: Reinhold dot Bader at lrz dot de
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