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
111
112
113
114
115
116
117
118
119
120
121
122
123
124
|
/* { dg-do run } */
/* { dg-skip-if "PR68356 no math-errno on darwin" { "*-*-darwin*" } } */
/* { dg-add-options ieee } */
/* { dg-require-effective-target fenv_exceptions } */
/* { dg-skip-if "fenv" { powerpc-ibm-aix* } } */
#include <fenv.h>
#include <math.h>
#include <errno.h>
extern void abort (void) __attribute__ ((noreturn));
#define LARGE_NEG_MAYBE_ERANGE 0x01
#define LARGE_NEG_ERANGE 0x02
#define LARGE_POS_ERANGE 0x04
#define LARGE_NEG_EDOM 0x08
#define LARGE_POS_EDOM 0x10
#define LARGE_ERANGE (LARGE_NEG_ERANGE | LARGE_POS_ERANGE)
#define LARGE_EDOM (LARGE_NEG_EDOM | LARGE_POS_EDOM)
#define POWER_ERANGE (LARGE_NEG_MAYBE_ERANGE | LARGE_POS_ERANGE)
#define TEST(CALL, FLAGS) (CALL, tester (FLAGS))
volatile double d;
volatile int i;
static void (*tester) (int);
void
check_quiet_nan (int flags __attribute__ ((unused)))
{
if (fetestexcept (FE_ALL_EXCEPT))
abort ();
if (errno)
abort ();
}
void
check_large_neg (int flags)
{
if (flags & LARGE_NEG_MAYBE_ERANGE)
return;
int expected_errno = (flags & LARGE_NEG_ERANGE ? ERANGE
: flags & LARGE_NEG_EDOM ? EDOM
: 0);
if (expected_errno != errno)
abort ();
errno = 0;
}
void
check_large_pos (int flags)
{
int expected_errno = (flags & LARGE_POS_ERANGE ? ERANGE
: flags & LARGE_POS_EDOM ? EDOM
: 0);
if (expected_errno != errno)
abort ();
errno = 0;
}
void
test (void)
{
TEST (acos (d), LARGE_EDOM);
TEST (asin (d), LARGE_EDOM);
TEST (acosh (d), LARGE_NEG_EDOM);
TEST (atanh (d), LARGE_EDOM);
TEST (cosh (d), LARGE_ERANGE);
TEST (sinh (d), LARGE_ERANGE);
TEST (log (d), LARGE_NEG_EDOM);
#if defined (__sun__) && defined (__unix__)
/* Disabled due to a bug in Solaris libm. */
if (0)
#endif
TEST (log2 (d), LARGE_NEG_EDOM);
TEST (log10 (d), LARGE_NEG_EDOM);
#if defined(__GLIBC__) && (__GLIBC__ < 2 || (__GLIBC__ == 2 && __GLIBC_MINOR__ < 22))
/* Disabled due to glibc PR 6792, fixed in glibc 2.22. */
if (0)
#endif
TEST (log1p (d), LARGE_NEG_EDOM);
TEST (exp (d), POWER_ERANGE);
#if (defined (__sun__) || defined(__hppa__)) && defined (__unix__)
/* Disabled due to a bug in Solaris libm. HP PA-RISC libm doesn't support
ERANGE for exp2. */
if (0)
#endif
{
TEST (exp2 (d), POWER_ERANGE);
#if defined(__GLIBC__) && (__GLIBC__ < 2 || (__GLIBC__ == 2 && __GLIBC_MINOR__ < 11))
/* Disabled due to glibc PR 6788, fixed in glibc 2.11. */
if (0)
#endif
TEST (expm1 (d), POWER_ERANGE);
}
TEST (sqrt (d), LARGE_NEG_EDOM);
TEST (pow (100.0, d), POWER_ERANGE);
TEST (pow (i, d), POWER_ERANGE);
}
int
main (void)
{
errno = 0;
i = 100;
d = __builtin_nan ("");
tester = check_quiet_nan;
feclearexcept (FE_ALL_EXCEPT);
test ();
d = -1.0e80;
tester = check_large_neg;
errno = 0;
test ();
d = 1.0e80;
tester = check_large_pos;
errno = 0;
test ();
return 0;
}
|
