summaryrefslogtreecommitdiffstats
path: root/src/boost/libs/multiprecision/test/bug12039.cpp
diff options
context:
space:
mode:
authorDaniel Baumann <daniel.baumann@progress-linux.org>2024-04-21 11:54:28 +0000
committerDaniel Baumann <daniel.baumann@progress-linux.org>2024-04-21 11:54:28 +0000
commite6918187568dbd01842d8d1d2c808ce16a894239 (patch)
tree64f88b554b444a49f656b6c656111a145cbbaa28 /src/boost/libs/multiprecision/test/bug12039.cpp
parentInitial commit. (diff)
downloadceph-e6918187568dbd01842d8d1d2c808ce16a894239.tar.xz
ceph-e6918187568dbd01842d8d1d2c808ce16a894239.zip
Adding upstream version 18.2.2.upstream/18.2.2
Signed-off-by: Daniel Baumann <daniel.baumann@progress-linux.org>
Diffstat (limited to 'src/boost/libs/multiprecision/test/bug12039.cpp')
-rw-r--r--src/boost/libs/multiprecision/test/bug12039.cpp40
1 files changed, 40 insertions, 0 deletions
diff --git a/src/boost/libs/multiprecision/test/bug12039.cpp b/src/boost/libs/multiprecision/test/bug12039.cpp
new file mode 100644
index 000000000..f6e601045
--- /dev/null
+++ b/src/boost/libs/multiprecision/test/bug12039.cpp
@@ -0,0 +1,40 @@
+///////////////////////////////////////////////////////////////////////////////
+// Copyright 2016 John Maddock. Distributed under 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 <boost/multiprecision/cpp_bin_float.hpp>
+
+int main()
+{
+ typedef boost::multiprecision::number<boost::multiprecision::backends::cpp_bin_float<256> > ext_float_t;
+ typedef boost::multiprecision::number<boost::multiprecision::backends::cpp_bin_float<2046> > long_ext_float_t;
+
+ ext_float_t x = 5e15;
+ x += 0.5;
+ ext_float_t x1 = x + 255.0 / (1 << 20); // + 2^-12 - eps
+ ext_float_t x2 = x + 257.0 / (1 << 20); // + 2^-12 + eps
+ double d1 = x1.convert_to<double>();
+ double d2 = x2.convert_to<double>();
+
+ std::cout << std::setprecision(18) << d1 << std::endl;
+ std::cout << std::setprecision(18) << d2 << std::endl;
+
+ x = 1e7 + 0.5;
+ x1 = x + ldexp(255.0, -38); // + 2^-30 - eps
+ x2 = x + ldexp(257.0, -38); // + 2^-30 + eps
+ float f1 = x1.convert_to<float>();
+ float f2 = x2.convert_to<float>();
+
+ std::cout << std::setprecision(9) << f1 << std::endl;
+ std::cout << std::setprecision(9) << f2 << std::endl;
+
+ long_ext_float_t lf(1);
+ lf += std::numeric_limits<long_ext_float_t>::epsilon();
+ lf += std::numeric_limits<float>::epsilon() / 2;
+ BOOST_MP_ASSERT(lf != 1);
+ float f3 = lf.convert_to<float>();
+ std::cout << std::setprecision(9) << f3 << std::endl;
+
+ return (d1 == d2) && (f1 == f2) && (f3 != 1) ? 0 : 1;
+}