1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
|
// Copyright (C) 2020-2023 Free Software Foundation, Inc.
//
// This file is part of the GNU ISO C++ Library. This library is free
// software; you can redistribute it and/or modify it under the
// terms of the GNU General Public License as published by the
// Free Software Foundation; either version 3, or (at your option)
// any later version.
//
// This library is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License along
// with this library; see the file COPYING3. If not see
// <http://www.gnu.org/licenses/>.
// only: float|double|ldouble * * *
// expensive: * [1-9] * *
#include "bits/main.h"
template <typename V>
void
test()
{
vir::test::setFuzzyness<float>(0);
vir::test::setFuzzyness<double>(0);
using T = typename V::value_type;
constexpr T inf = std::__infinity_v<T>;
constexpr T nan = std::__quiet_NaN_v<T>;
constexpr T denorm_min = std::__denorm_min_v<T>;
constexpr T norm_min = std::__norm_min_v<T>;
constexpr T max = std::__finite_max_v<T>;
#if defined __LONG_DOUBLE_IBM128__
// On POWER with IBM128 long double, 1+eps and 2-eps is not a constant
// expression. Until this is fixed, just use const instead of constexpr.
// (error: '(1.0e+0l + 4.94065645841246544176568792868221e-324l)' is not a
// constant expression)
const T after_one = 1 + std::__epsilon_v<T>;
const T before_one = (2 - std::__epsilon_v<T>) / 2;
#else
constexpr T after_one = 1 + std::__epsilon_v<T>;
constexpr T before_one = (2 - std::__epsilon_v<T>) / 2;
#endif
const std::initializer_list<T>
input_values = {+0.,
0.5,
-0.5,
before_one,
-before_one,
after_one,
-after_one,
1.5,
-1.5,
2 * before_one,
-2 * before_one,
2 * after_one,
-2 * after_one,
2.5,
-2.5,
0x1.fffffffffffffp52,
-0x1.fffffffffffffp52,
0x1.ffffffffffffep52,
-0x1.ffffffffffffep52,
0x1.ffffffffffffdp52,
-0x1.ffffffffffffdp52,
0x1.fffffep21,
-0x1.fffffep21,
0x1.fffffcp21,
-0x1.fffffcp21,
0x1.fffffep22,
-0x1.fffffep22,
0x1.fffffcp22,
-0x1.fffffcp22,
0x1.fffffep23,
-0x1.fffffep23,
0x1.fffffcp23,
-0x1.fffffcp23,
0x1.8p23,
-0x1.8p23,
inf,
-inf,
-0.,
nan,
denorm_min,
norm_min / 3,
norm_min,
max};
test_values<V>(input_values, {10000}, MAKE_TESTER(erf), MAKE_TESTER(erfc),
MAKE_TESTER(tgamma), MAKE_TESTER(lgamma), MAKE_TESTER(ceil),
MAKE_TESTER(floor), MAKE_TESTER(trunc), MAKE_TESTER(round),
MAKE_TESTER(lround), MAKE_TESTER(llround),
MAKE_TESTER(nearbyint), MAKE_TESTER(rint), MAKE_TESTER(lrint),
MAKE_TESTER(llrint), MAKE_TESTER(ilogb));
// sqrt(x) on x87 is precise in 80 bits, but the subsequent rounding can be
// wrong (up to 1 ULP)
#if __FLT_EVAL_METHOD__ == 1
vir::test::setFuzzyness<float>(1);
vir::test::setFuzzyness<double>(0);
#elif __FLT_EVAL_METHOD__ == 2
vir::test::setFuzzyness<float>(1);
vir::test::setFuzzyness<double>(1);
#endif
test_values<V>(input_values, {10000}, MAKE_TESTER(sqrt));
}
|
