mirror of
https://github.com/boostorg/math.git
synced 2025-05-11 21:33:52 +00:00
09e1350f6f
* Update tanh_sinh handling of boundaries. So that we don't accidentality end up at an end point, even when arithmetic is inexact, and the FP type has no denomrms. Fixes: https://github.com/boostorg/math/issues/893 * Add licence/copyright to new file. * Remove tabs.
33 lines
910 B
C++
33 lines
910 B
C++
// Copyright John Maddock, 2022
|
|
// 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)
|
|
|
|
#define BOOST_TEST_MODULE issue893
|
|
|
|
#include <iostream>
|
|
#include <sstream>
|
|
#include <boost/test/included/unit_test.hpp>
|
|
#include <boost/multiprecision/cpp_bin_float.hpp>
|
|
#include <boost/math/quadrature/tanh_sinh.hpp>
|
|
|
|
using boost::math::quadrature::tanh_sinh;
|
|
|
|
#include <iostream>
|
|
|
|
BOOST_AUTO_TEST_CASE(issue893) {
|
|
typedef boost::multiprecision::cpp_bin_float_100 real;
|
|
|
|
auto fun = [](real x) -> real {
|
|
return 1.0;
|
|
};
|
|
|
|
tanh_sinh<real> integrator;
|
|
const real a = 0.0;
|
|
const real b = -0.9999995515592481132478776023609116290187750667053638330158486516399489191171270344610533516275817076;
|
|
real y = integrator.integrate(fun, -b, a);
|
|
|
|
BOOST_CHECK(y < 1);
|
|
}
|