[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Help-gsl] Re: gsl-1.14 linalg test: FAIL: QRPT_update m(5,3)
|
From: |
Brian Gough |
|
Subject: |
Re: [Help-gsl] Re: gsl-1.14 linalg test: FAIL: QRPT_update m(5,3) |
|
Date: |
Thu, 23 Dec 2010 14:10:56 +0000 |
|
User-agent: |
Wanderlust/2.15.6 (Almost Unreal) Emacs/23.2 Mule/6.0 (HANACHIRUSATO) |
At Wed, 22 Dec 2010 15:43:48 -0500,
Leo Razoumov wrote:
>
> I think I found a fix. The problem is due to infamous "extended
> precision (80 bits)" registers and 64-bit doubles.
>
Hi. It sounds like the test tolerance is too small. Does the
following fix it?
=== modified file 'linalg/test.c'
--- linalg/test.c 2010-06-13 15:00:43 +0000
+++ linalg/test.c 2010-12-23 14:08:57 +0000
@@ -1552,7 +1552,7 @@
gsl_test(f, " QRPT_update m(3,5)");
s += f;
- f = test_QRPT_update_dim(m53, 2 * 512.0 * GSL_DBL_EPSILON);
+ f = test_QRPT_update_dim(m53, 2 * 1024.0 * GSL_DBL_EPSILON);
gsl_test(f, " QRPT_update m(5,3)");
s += f;