From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 7522 invoked by alias); 21 Aug 2008 19:40:14 -0000 Received: (qmail 7435 invoked by uid 22791); 21 Aug 2008 19:40:13 -0000 X-Spam-Check-By: sourceware.org Received: from mailout04.t-online.de (HELO mailout04.t-online.de) (194.25.134.18) by sourceware.org (qpsmtpd/0.31) with ESMTP; Thu, 21 Aug 2008 19:39:21 +0000 Received: from fwd32.aul.t-online.de by mailout04.sul.t-online.de with smtp id 1KWG06-0003nW-00; Thu, 21 Aug 2008 21:39:18 +0200 Received: from [192.168.2.20] (TJtGnuZeQhZ8zvsegWHRDFChxqiFjP2UzIEVxv1TwLSECN9mhOU48Ut97ew3Hs5gtF@[84.154.74.159]) by fwd32.aul.t-online.de with esmtp id 1KWFzn-1nLhrs0; Thu, 21 Aug 2008 21:38:59 +0200 Message-ID: <48ADC452.30707@lrz.de> Date: Thu, 21 Aug 2008 19:40:00 -0000 From: bader-reinhold@t-online.de (Reinhold Bader) Reply-To: Reinhold.Bader@lrz.de User-Agent: Thunderbird 2.0.0.14 (X11/20080421) MIME-Version: 1.0 To: gsl-discuss@sourceware.org CC: Brian Gough Subject: Documentation (2) Content-Type: text/plain; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: 7bit X-ID: TJtGnuZeQhZ8zvsegWHRDFChxqiFjP2UzIEVxv1TwLSECN9mhOU48Ut97ew3Hs5gtF X-TOI-MSGID: 705c1308-9aa9-48dd-b0ae-cc557fb0fadf Mailing-List: contact gsl-discuss-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: gsl-discuss-owner@sourceware.org X-SW-Source: 2008-q3/txt/msg00020.txt.bz2 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