From mboxrd@z Thu Jan  1 00:00:00 1970
Return-Path: <gsl-discuss-return-2617-listarch-gsl-discuss=sourceware.cygnus.com@sources.redhat.com>
Received: (qmail 30229 invoked by alias); 19 Aug 2002 18:32:02 -0000
Mailing-List: contact gsl-discuss-help@sources.redhat.com; run by ezmlm
Precedence: bulk
List-Subscribe: <mailto:gsl-discuss-subscribe@sources.redhat.com>
List-Archive: <http://sources.redhat.com/ml/gsl-discuss/>
List-Post: <mailto:gsl-discuss@sources.redhat.com>
List-Help: <mailto:gsl-discuss-help@sources.redhat.com>, <http://sources.redhat.com/ml/#faqs>
Sender: gsl-discuss-owner@sources.redhat.com
Received: (qmail 30222 invoked from network); 19 Aug 2002 18:32:01 -0000
Received: from unknown (HELO blueyonder.co.uk) (195.188.53.94)
  by sources.redhat.com with SMTP; 19 Aug 2002 18:32:01 -0000
Received: from pcow057o.blueyonder.co.uk ([127.0.0.1]) by blueyonder.co.uk  with Microsoft SMTPSVC(5.5.1877.757.75);
	 Mon, 19 Aug 2002 19:32:00 +0100
Received: from localhost (unverified [62.30.73.52]) by pcow057o.blueyonder.co.uk
 (Content Technologies SMTPRS 4.2.9) with ESMTP id <T5cd068cc8aac1785b311f@pcow057o.blueyonder.co.uk>;
 Mon, 19 Aug 2002 19:32:00 +0100
Received: by network-theory.co.uk
	via sendmail from stdin
	id <m17grJe-0000eBC@localhost> (Debian Smail3.2.0.102)
	for huberph@infomaniak.ch; Mon, 19 Aug 2002 19:31:50 +0100 (BST) 
From: Brian Gough <bjg@network-theory.co.uk>
MIME-Version: 1.0
Content-Type: text/plain; charset="us-ascii"
Content-Transfer-Encoding: 7bit
Message-ID: <15713.14742.437923.195861@debian>
Date: Mon, 19 Aug 2002 11:32:00 -0000
To: Philippe Huber <huberph@infomaniak.ch>
Cc: gsl-discuss@sources.redhat.com
Subject: Re: Problem with matrices bigger than 2x2 using gsl_linalg_LU_invert
In-Reply-To: <001101c2476d$18a7c250$7e28c281@osiris>
References: <001101c2476d$18a7c250$7e28c281@osiris>
X-SW-Source: 2002-q3/txt/msg00152.txt.bz2

Philippe Huber writes:
 > First of all, thank you for having created such a great library. I have
 > a small problem with gsl_linalg_LU_invert. I am using a 4x4 matrix which
 > is not singular. I made a LU decomposition using gsl_linalg_LU_decomp.
 > Then, I send it to gsl_linalg_LU_invert and it crashes. I tried with a
 > 2x2 matrix and it works perfectly well.  I verified that the program also
 > crashes with a 3x3 matrix. Is it a bug or did I make something wrong.
 > Please apologize if I don't send my code, it is well too big.

Please send a small example which reproduces the bug, with details of 
your operating system, compiler, etc.
regards
Brian Gough