From mboxrd@z Thu Jan  1 00:00:00 1970
Return-Path: <gsl-discuss-return-2653-listarch-gsl-discuss=sourceware.cygnus.com@sources.redhat.com>
Received: (qmail 15441 invoked by alias); 5 Sep 2002 15:37:53 -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 15432 invoked from network); 5 Sep 2002 15:37:52 -0000
Received: from unknown (HELO escudo.nhh.no) (158.37.96.28)
  by sources.redhat.com with SMTP; 5 Sep 2002 15:37:52 -0000
Received: from nhh.no (colon.nhh.no [158.37.98.22])
	by escudo.nhh.no (8.12.2/8.12.2) with SMTP id g85FbpFM077211
	for <gsl-discuss@sources.redhat.com>; Thu, 5 Sep 2002 17:37:51 +0200 (CEST)
Received: from d226-24-170.home.cgocable.net ([24.226.24.170])
        (SquirrelMail authenticated user sameos)
        by webmail.nhh.no with HTTP;
        Thu, 5 Sep 2002 11:37:51 -0400 (EDT)
Message-ID: <50246.24.226.24.170.1031240271.squirrel@webmail.nhh.no>
Date: Thu, 05 Sep 2002 08:37:00 -0000
Subject: matrix decompositions fails on Max os X 10.2
From: "=?iso-8859-1?Q?Erik_=D8._S=F8rensen?=" <erik.sorensen@nhh.no>
To: <gsl-discuss@sources.redhat.com>
MIME-Version: 1.0
Content-Type: text/plain; charset=iso-8859-1
Content-Transfer-Encoding: 8bit
X-SW-Source: 2002-q3/txt/msg00188.txt.bz2

Matrix decompositions seem to fail:

PASS:   bidiag_decomp m(5,3)
(  9,  7)[0,0]:    1.00000095367431574                        1
(  9,  7)[0,1]:   0.500000476837157981                      0.5
(  9,  7)[0,3]:   0.250000238418579046                     0.25
(  9,  7)[2,1]:   0.250000238418578935                     0.25
(  9,  7)[3,4]:    0.12500011920928944                    0.125
(  9,  7)[4,3]:   0.125000119209289468                    0.125
(  9,  7)[5,2]:   0.125000119209289468                    0.125
(  9,  7)[6,1]:    0.12500011920928944                    0.125
FAIL:   bidiag_decomp m(9,7)
(  2,  2)[0,1]:   0.500000476837157759                      0.5
(  2,  2)[1,0]:   0.500000476837157759                      0.5
FAIL:   bidiag_decomp hilbert(2)
(  3,  3)[0,0]:    1.00000095367431574                        1

and it goes on, some of the decompositions passing, many of them failing,
seems to be mainly in the 7th  to  9th significant digits. Would be very
grateful for any help with what causes this -- I have never had any
problems with gsl on the linux laptop which has suffered massive hardware
failure.

Thanks,
Erik

-- 
Erik Øiolf Sørensen  <mailto:erik.sorensen@nhh.no>
phd student <http://www.geocities.com/erik_oiolf/>
Department of Economics, Norwegian School of Economics
Visiting Queen's University, Kingston, Ontario.
phone, office: (+1) 613-533-6000, ext. 74030.