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
108
109
110
|
/* SPDX-License-Identifier: LGPL-2.1-or-later */
#include <float.h>
#include "math-util.h"
#include "tests.h"
TEST(iszero_safe) {
/* zeros */
assert_se(iszero_safe(0.0));
assert_se(iszero_safe(-0.0));
assert_se(iszero_safe(0e0));
assert_se(iszero_safe(-0e0));
assert_se(iszero_safe(0e+0));
assert_se(iszero_safe(0e-0));
assert_se(iszero_safe(-0e-0));
assert_se(iszero_safe(-0e000));
assert_se(iszero_safe(0e000));
/* non-zero normal values */
assert_se(!iszero_safe(42.0));
assert_se(!iszero_safe(M_PI));
assert_se(!iszero_safe(DBL_MAX));
assert_se(!iszero_safe(-DBL_MAX));
assert_se(!iszero_safe(DBL_MIN));
assert_se(!iszero_safe(-DBL_MIN));
assert_se(!iszero_safe(1 / DBL_MAX));
/* subnormal values */
assert_se(!iszero_safe(DBL_MIN / 2));
assert_se(!iszero_safe(-DBL_MIN / 42));
assert_se(!iszero_safe(1 / DBL_MAX / 2));
/* too small values which cannot be in subnormal form */
assert_se( iszero_safe(DBL_MIN / DBL_MAX));
assert_se( iszero_safe(DBL_MIN / -DBL_MAX));
assert_se( iszero_safe(-DBL_MIN / DBL_MAX));
assert_se( iszero_safe(-DBL_MIN / -DBL_MAX));
/* NaN or infinity */
assert_se(!iszero_safe(NAN));
assert_se(!iszero_safe(INFINITY));
assert_se(!iszero_safe(-INFINITY));
assert_se(!iszero_safe(1 / NAN));
/* inverse of infinity */
assert_se( iszero_safe(1 / INFINITY));
assert_se( iszero_safe(1 / -INFINITY));
assert_se( iszero_safe(-1 / INFINITY));
assert_se( iszero_safe(-1 / -INFINITY));
assert_se( iszero_safe(42 / -INFINITY));
assert_se( iszero_safe(-42 / -INFINITY));
assert_se( iszero_safe(DBL_MIN / INFINITY));
assert_se( iszero_safe(DBL_MIN / -INFINITY));
assert_se( iszero_safe(DBL_MAX / INFINITY / 2));
assert_se( iszero_safe(DBL_MAX / -INFINITY * DBL_MAX));
/* infinity / infinity is NaN */
assert_se(!iszero_safe(INFINITY / INFINITY));
assert_se(!iszero_safe(INFINITY * 2 / INFINITY));
assert_se(!iszero_safe(INFINITY / DBL_MAX / INFINITY));
}
TEST(fp_equal) {
/* normal values */
assert_se( fp_equal(0.0, -0e0));
assert_se( fp_equal(3.0, 3));
assert_se(!fp_equal(3.000001, 3));
assert_se( fp_equal(M_PI, M_PI));
assert_se(!fp_equal(M_PI, -M_PI));
assert_se( fp_equal(DBL_MAX, DBL_MAX));
assert_se(!fp_equal(DBL_MAX, -DBL_MAX));
assert_se(!fp_equal(-DBL_MAX, DBL_MAX));
assert_se( fp_equal(-DBL_MAX, -DBL_MAX));
assert_se( fp_equal(DBL_MIN, DBL_MIN));
assert_se(!fp_equal(DBL_MIN, -DBL_MIN));
assert_se(!fp_equal(-DBL_MIN, DBL_MIN));
assert_se( fp_equal(-DBL_MIN, -DBL_MIN));
/* subnormal values */
assert_se( fp_equal(DBL_MIN / 10, DBL_MIN / 10));
assert_se(!fp_equal(DBL_MIN / 10, -DBL_MIN / 10));
assert_se(!fp_equal(-DBL_MIN / 10, DBL_MIN / 10));
assert_se( fp_equal(-DBL_MIN / 10, -DBL_MIN / 10));
assert_se(!fp_equal(DBL_MIN / 10, DBL_MIN / 15));
assert_se(!fp_equal(DBL_MIN / 10, DBL_MIN / 15));
/* subnormal difference */
assert_se(!fp_equal(DBL_MIN / 10, DBL_MIN + DBL_MIN / 10));
assert_se( fp_equal(3.0, 3.0 + DBL_MIN / 2)); /* 3.0 + DBL_MIN / 2 is truncated to 3.0 */
/* too small values */
assert_se( fp_equal(DBL_MIN / DBL_MAX, -DBL_MIN / DBL_MAX));
/* NaN or infinity */
assert_se(!fp_equal(NAN, NAN));
assert_se(!fp_equal(NAN, 0));
assert_se(!fp_equal(NAN, INFINITY));
assert_se(!fp_equal(INFINITY, INFINITY));
assert_se(!fp_equal(INFINITY, -INFINITY));
assert_se(!fp_equal(-INFINITY, INFINITY));
assert_se(!fp_equal(-INFINITY, -INFINITY));
/* inverse of infinity */
assert_se( fp_equal(0, 1 / INFINITY));
assert_se( fp_equal(42 / INFINITY, 1 / -INFINITY));
assert_se(!fp_equal(42 / INFINITY, INFINITY / INFINITY));
}
DEFINE_TEST_MAIN(LOG_DEBUG);
|
