[lmi-commits] [lmi] master 8b486fc 4/5: Show TOMS 748 result for a problematic test

[Greg Chicares]

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)