From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 15426 invoked by alias); 16 Nov 2002 19:05:50 -0000 Mailing-List: contact gsl-discuss-help@sources.redhat.com; run by ezmlm Precedence: bulk List-Subscribe: List-Archive: List-Post: List-Help: , Sender: gsl-discuss-owner@sources.redhat.com Received: (qmail 15411 invoked from network); 16 Nov 2002 19:05:48 -0000 Received: from unknown (HELO blueyonder.co.uk) (195.188.53.121) by sources.redhat.com with SMTP; 16 Nov 2002 19:05:48 -0000 Received: from pcow035o.blueyonder.co.uk ([127.0.0.1]) by blueyonder.co.uk with Microsoft SMTPSVC(5.5.1877.757.75); Sat, 16 Nov 2002 19:05:47 +0000 Received: from localhost (unverified [62.30.73.59]) by pcow035o.blueyonder.co.uk (Content Technologies SMTPRS 4.2.9) with ESMTP id ; Sat, 16 Nov 2002 19:05:47 +0000 Received: by network-theory.co.uk via sendmail from stdin id (Debian Smail3.2.0.102) for P.J.P.Roche@damtp.cam.ac.uk; Sat, 16 Nov 2002 18:06:12 +0000 (GMT) From: Brian Gough MIME-Version: 1.0 Content-Type: text/plain; charset="us-ascii" Content-Transfer-Encoding: 7bit Message-ID: <15830.35092.202901.707582@debian> Date: Sat, 16 Nov 2002 11:05:00 -0000 To: Peter Roche Cc: gsl-discuss@sources.redhat.com, bug-gsl@gnu.org Subject: Re: LU Decomposition error In-Reply-To: References: X-SW-Source: 2002-q4/txt/msg00153.txt.bz2 Peter Roche writes: > I hope this isn't a repeat of an already dealt with problem, but I > couldn't see anything in the archives. I have encountered some problems > with the gsl_linalg_LU_decomp rountine. The routine ran without exiting > with an error message, but I found that if I multiply the decomposed > Lower and Upper triangular matrices I do not get the original matrix I > started with before decomposing. I have included a section of my code > and also a section of the original matrix.... Please send a complete example program which reproduces the problem. Thanks.