[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[lmi-commits] [lmi] master 8b486fc 4/5: Show TOMS 748 result for a probl
From: |
Greg Chicares |
Subject: |
[lmi-commits] [lmi] master 8b486fc 4/5: Show TOMS 748 result for a problematic test |
Date: |
Mon, 27 Sep 2021 17:46:55 -0400 (EDT) |
branch: master
commit 8b486fcbda1b7bfcae6787bcc7d41f19fea3fa6b
Author: Gregory W. Chicares <gchicares@sbcglobal.net>
Commit: Gregory W. Chicares <gchicares@sbcglobal.net>
Show TOMS 748 result for a problematic test
Document and show how to reproduce the value calculated today using the
TOMS algorithm as translated by 'f2c':
0.022317679157464680173 (x86_64-pc-linux-gnu)
That matches the value published by the TOMS 748 authors in 'testout',
to the number of decimal places specified in that file:
0.022317679157465 published by TOMS 748 authors
---
zero_test.cpp | 8 ++++++++
1 file changed, 8 insertions(+)
diff --git a/zero_test.cpp b/zero_test.cpp
index 1819e7e..1e0b742 100644
--- a/zero_test.cpp
+++ b/zero_test.cpp
@@ -2355,6 +2355,14 @@ void test_alefeld_examples(int alefeld_count, double tol)
++g;
show_details();
+#if 0
+# if defined TEST_ALEFELD
+ // Result is 0.022317679157464680173 for x86_64-pc-linux-gnu:
+ r = toms748_root(f13n26, bound0, bound1, bias_none, 1000);
+ std::cout << "#13: " << r.root << std::endl;
+# endif // defined TEST_ALEFELD
+#endif // 0
+
// Table I #14 = FORTRAN #27
e = 14; f = 27;
auto f14n27 = [&n](double x)