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; _______________________________________________ Help-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/help-gsl
