import quincy beta 17.1.0

[ceph.git] / ceph / src / boost / libs / math / test / cohen_acceleration_test.cpp

diff --git a/ceph/src/boost/libs/math/test/cohen_acceleration_test.cpp b/ceph/src/boost/libs/math/test/cohen_acceleration_test.cpp
new file mode 100644 (file)
index 0000000..291cb85
--- /dev/null
+++ b/ceph/src/boost/libs/math/test/cohen_acceleration_test.cpp
@@ -0,0 +1,86 @@
+/*
+ * Copyright Nick Thompson, 2020
+ * Use, modification and distribution are subject to the
+ * Boost Software License, Version 1.0. (See accompanying file
+ * LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
+ */
+
+#include "math_unit_test.hpp"
+#include <boost/math/tools/cohen_acceleration.hpp>
+#include <boost/math/constants/constants.hpp>
+#ifdef BOOST_HAS_FLOAT128
+#include <boost/multiprecision/float128.hpp>
+using boost::multiprecision::float128;
+#endif
+#include <boost/multiprecision/cpp_bin_float.hpp>
+
+using boost::math::tools::cohen_acceleration;
+using boost::multiprecision::cpp_bin_float_100;
+using boost::math::constants::pi;
+
+template<typename Real>
+class G {
+public:
+    G(){
+        k_ = 0;
+    }
+    
+    Real operator()() {
+        k_ += 1;
+        return 1/(k_*k_);
+    }
+
+private:
+    Real k_;
+};
+
+template<typename Real>
+void test_pisq_div12()
+{
+    auto g = G<Real>();
+    Real x = cohen_acceleration(g);
+    CHECK_ULP_CLOSE(pi<Real>()*pi<Real>()/12, x, 3);
+}
+
+template<typename Real>
+class Divergent {
+public:
+    Divergent(){
+        k_ = 0;
+    }
+
+    // See C3 of: https://people.mpim-bonn.mpg.de/zagier/files/exp-math-9/fulltext.pdf
+    Real operator()() {
+        using std::log;
+        k_ += 1;
+        return log(k_);
+    }
+
+private:
+    Real k_;
+};
+
+template<typename Real>
+void test_divergent()
+{
+    auto g = Divergent<Real>();
+    Real x = -cohen_acceleration(g);
+    CHECK_ULP_CLOSE(log(pi<Real>()/2)/2, x, 80);
+}
+
+int main()
+{
+    test_pisq_div12<float>();
+    test_pisq_div12<double>();
+    test_pisq_div12<long double>();
+
+    test_divergent<float>();
+    test_divergent<double>();
+    test_divergent<long double>();
+
+    #ifdef BOOST_HAS_FLOAT128
+    test_pisq_div12<float128>();
+    test_divergent<float128>();
+    #endif
+    return boost::math::test::report_errors();
+}