help-gsl
[Top][All Lists]
Advanced

[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;
 



reply via email to

[Prev in Thread] Current Thread [Next in Thread]