boost/multiprecision
2
0
Fork 0
mirror of https://github.com/boostorg/multiprecision.git synced 2026-01-27 07:02:10 +00:00
Code Issues Projects Releases Wiki Activity

These error values for + - * / √ were tested 2e7 times and are passing.

Browse Source
This commit is contained in:
Janek Kozicki
2021-08-04 21:56:19 +02:00
parent 6e50915a99
commit 659af61946
1 changed files with 13 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

20
test/test_cpp_double_float_arithmetic.cpp
View File

@@ -29,6 +29,7 @@
#include <boost/multiprecision/cpp_double_float.hpp>
#include <boost/multiprecision/cpp_dec_float.hpp>
#include <boost/random/uniform_real_distribution.hpp>
#include <boost/multiprecision/traits/max_digits10.hpp>
#include <boost/core/demangle.hpp>
#if defined(__clang__)
@@ -51,14 +52,19 @@ namespace local
{
using float_type = FloatingPointConstituentType;
static constexpr int digits = 2 * (std::numeric_limits<float_type>::digits - 1);
static constexpr int digits10 = (int) (float(digits - 1) * 0.301F);
static constexpr int digits = 2 * std::numeric_limits<float_type>::digits;
static constexpr int digits10 = boost::multiprecision::detail::calc_digits10<digits>::value;
static constexpr int max_digits10 = boost::multiprecision::detail::calc_max_digits10<digits>::value;
static unsigned seed_prescaler;
using double_float_type = boost::multiprecision::number<boost::multiprecision::backends::cpp_double_float<float_type>, boost::multiprecision::et_off>;
using control_float_type = boost::multiprecision::number<boost::multiprecision::cpp_dec_float<(2 * std::numeric_limits<double_float_type>::digits10) + 1>, boost::multiprecision::et_off>;
static_assert( digits == std::numeric_limits<double_float_type>::digits , "" );
static_assert( digits10 == std::numeric_limits<double_float_type>::digits10 , "" );
static_assert( max_digits10 == std::numeric_limits<double_float_type>::max_digits10 , "" );
template<const std::size_t DigitsToGet = digits10>
static void get_random_fixed_string(std::string& str, const bool is_unsigned = false)
{
@@ -163,7 +169,7 @@ namespace local
{
bool result_is_ok = true;
const control_float_type MaxError = ldexp(control_float_type(1), -std::numeric_limits<double_float_type>::digits + 0);
const control_float_type MaxError = ldexp(control_float_type(1), 2 - std::numeric_limits<double_float_type>::digits);
for(std::uint32_t i = 0U; ((i < count) && result_is_ok); ++i)
{
@@ -196,7 +202,7 @@ namespace local
{
bool result_is_ok = true;
const control_float_type MaxError = ldexp(control_float_type(1), -std::numeric_limits<double_float_type>::digits + 0);
const control_float_type MaxError = ldexp(control_float_type(1), 2 - std::numeric_limits<double_float_type>::digits);
for(std::uint32_t i = 0U; ((i < count) && result_is_ok); ++i)
{
@@ -229,7 +235,7 @@ namespace local
{
bool result_is_ok = true;
const control_float_type MaxError = ldexp(control_float_type(1), -std::numeric_limits<double_float_type>::digits + 2);
const control_float_type MaxError = ldexp(control_float_type(1), 4 - std::numeric_limits<double_float_type>::digits);
for(std::uint32_t i = 0U; ((i < count) && result_is_ok); ++i)
{
@@ -262,7 +268,7 @@ namespace local
{
bool result_is_ok = true;
const control_float_type MaxError = ldexp(control_float_type(1), -std::numeric_limits<double_float_type>::digits + 2);
const control_float_type MaxError = ldexp(control_float_type(1), 3 - std::numeric_limits<double_float_type>::digits);
for(std::uint32_t i = 0U;((i < count) && result_is_ok); ++i)
{
@@ -295,7 +301,7 @@ namespace local
{
bool result_is_ok = true;
const control_float_type MaxError = ldexp(control_float_type(1), -std::numeric_limits<double_float_type>::digits + 1);
const control_float_type MaxError = ldexp(control_float_type(1), 3 - std::numeric_limits<double_float_type>::digits);
for(std::uint32_t i = 0U; ((i < count) && result_is_ok); ++i)
{