mirror of
https://github.com/boostorg/math.git
synced 2026-01-23 05:42:15 +00:00
fae96bf542
* Remove tests we don't need right now. !!!REVERT THIS COMMIT BEFORE MERGING!!! * Add s390x testing to drone. * Correct drone file. * Correct drone file (again) * Prevent complete cancellation in bessel_jy logic. * Correct testing for 128-bit floats. * Make some more tests 128-bit long double safe. * Make more tests 128-bit float safe. * Fix some more 128-bit testing issues. * More 128-bit float fixes. * Make more tests 128-bit float safe. * Fix up remaining tests for 128-bit floats. * Yet more 128-bit float test case fixes. * Fix up more tests for 128-bit floats and non-intel platforms. * Fix up more tests to be 128-bit long double safe. * More test case adjustments. * More 128-bit float error rate adjustments. * Fixes for autodiff tests * Two more test fixes. * Fix up daubechies_scaling_test.cpp and reinstate full CI. Co-authored-by: Matt Borland <matt@mattborland.com>
89 lines
1.9 KiB
C++
89 lines
1.9 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>
|
|
#include <cmath>
|
|
|
|
using boost::math::tools::cohen_acceleration;
|
|
using boost::multiprecision::cpp_bin_float_100;
|
|
using boost::math::constants::pi;
|
|
using std::log;
|
|
|
|
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, (std::numeric_limits<Real>::digits > 100 ? 350 : 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();
|
|
}
|