87 lines
1.8 KiB
C++
87 lines
1.8 KiB
C++
/*
|
|
* 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, 135);
|
|
}
|
|
|
|
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();
|
|
}
|