1 |
689 |
jeremybenn |
/* Copyright (C) 2006 Free Software Foundation.
|
2 |
|
|
|
3 |
|
|
Verify that built-in math function constant folding of constant
|
4 |
|
|
arguments is correctly performed by the compiler.
|
5 |
|
|
|
6 |
|
|
Origin: Kaveh R. Ghazi, October 23, 2006. */
|
7 |
|
|
|
8 |
|
|
/* { dg-do link } */
|
9 |
|
|
|
10 |
|
|
/* Define "e" with as many bits as found in builtins.c:dconste. */
|
11 |
|
|
#define M_E 2.7182818284590452353602874713526624977572470936999595749669676277241
|
12 |
|
|
|
13 |
|
|
/* All references to link_error should go away at compile-time. */
|
14 |
|
|
extern void link_error(int);
|
15 |
|
|
|
16 |
|
|
/* Return TRUE if the sign of X != sign of Y. This is important when
|
17 |
|
|
comparing signed zeros. */
|
18 |
|
|
#define CKSGN_F(X,Y) \
|
19 |
|
|
(__builtin_copysignf(1.0F,(X)) != __builtin_copysignf(1.0F,(Y)))
|
20 |
|
|
#define CKSGN(X,Y) \
|
21 |
|
|
(__builtin_copysign(1.0,(X)) != __builtin_copysign(1.0,(Y)))
|
22 |
|
|
#define CKSGN_L(X,Y) \
|
23 |
|
|
(__builtin_copysignl(1.0L,(X)) != __builtin_copysignl(1.0L,(Y)))
|
24 |
|
|
|
25 |
|
|
/* Test that FUNC(ARG) == (RES). */
|
26 |
|
|
#define TESTIT(FUNC,ARG,RES) do { \
|
27 |
|
|
if (__builtin_##FUNC##f(ARG##F) != RES##F \
|
28 |
|
|
|| CKSGN_F(__builtin_##FUNC##f(ARG##F),RES##F)) \
|
29 |
|
|
link_error(__LINE__); \
|
30 |
|
|
if (__builtin_##FUNC(ARG) != RES \
|
31 |
|
|
|| CKSGN(__builtin_##FUNC(ARG),RES)) \
|
32 |
|
|
link_error(__LINE__); \
|
33 |
|
|
if (__builtin_##FUNC##l(ARG##L) != RES##L \
|
34 |
|
|
|| CKSGN_L(__builtin_##FUNC##l(ARG##L),RES##L)) \
|
35 |
|
|
link_error(__LINE__); \
|
36 |
|
|
} while (0)
|
37 |
|
|
|
38 |
|
|
/* Range test, check that (LOW) < FUNC(ARG) < (HI). */
|
39 |
|
|
#define TESTIT_R(FUNC,ARG,LOW,HI) do { \
|
40 |
|
|
if (__builtin_##FUNC##f(ARG) <= (LOW) || __builtin_##FUNC##f(ARG) >= (HI)) \
|
41 |
|
|
link_error(__LINE__); \
|
42 |
|
|
if (__builtin_##FUNC(ARG) <= (LOW) || __builtin_##FUNC(ARG) >= (HI)) \
|
43 |
|
|
link_error(__LINE__); \
|
44 |
|
|
if (__builtin_##FUNC##l(ARG) <= (LOW) || __builtin_##FUNC##l(ARG) >= (HI)) \
|
45 |
|
|
link_error(__LINE__); \
|
46 |
|
|
} while (0)
|
47 |
|
|
|
48 |
|
|
/* Test that FUNC(ARG1, ARG2) == (RES). */
|
49 |
|
|
#define TESTIT2(FUNC,ARG1,ARG2,RES) do { \
|
50 |
|
|
if (__builtin_##FUNC##f(ARG1##F, ARG2##F) != RES##F \
|
51 |
|
|
|| CKSGN_F(__builtin_##FUNC##f(ARG1##F,ARG2##F),RES##F)) \
|
52 |
|
|
link_error(__LINE__); \
|
53 |
|
|
if (__builtin_##FUNC(ARG1, ARG2) != RES \
|
54 |
|
|
|| CKSGN(__builtin_##FUNC(ARG1,ARG2),RES)) \
|
55 |
|
|
link_error(__LINE__); \
|
56 |
|
|
if (__builtin_##FUNC##l(ARG1##L, ARG2##L) != RES##L \
|
57 |
|
|
|| CKSGN_L(__builtin_##FUNC##l(ARG1##L,ARG2##L),RES##L)) \
|
58 |
|
|
link_error(__LINE__); \
|
59 |
|
|
} while (0)
|
60 |
|
|
|
61 |
|
|
/* Range test, check that (LOW) < FUNC(ARG1,ARG2) < (HI). */
|
62 |
|
|
#define TESTIT2_R(FUNC,ARG1,ARG2,LOW,HI) do { \
|
63 |
|
|
if (__builtin_##FUNC##f(ARG1, ARG2) <= (LOW) \
|
64 |
|
|
|| __builtin_##FUNC##f(ARG1, ARG2) >= (HI)) \
|
65 |
|
|
link_error(__LINE__); \
|
66 |
|
|
if (__builtin_##FUNC(ARG1, ARG2) <= (LOW) \
|
67 |
|
|
|| __builtin_##FUNC(ARG1, ARG2) >= (HI)) \
|
68 |
|
|
link_error(__LINE__); \
|
69 |
|
|
if (__builtin_##FUNC##l(ARG1, ARG2) <= (LOW) \
|
70 |
|
|
|| __builtin_##FUNC##l(ARG1, ARG2) >= (HI)) \
|
71 |
|
|
link_error(__LINE__); \
|
72 |
|
|
} while (0)
|
73 |
|
|
|
74 |
|
|
/* Test that FUNC(ARG1, ARG2, ARG3) == (RES). */
|
75 |
|
|
#define TESTIT3(FUNC,ARG1,ARG2,ARG3,RES) do { \
|
76 |
|
|
if (__builtin_##FUNC##f(ARG1##F, ARG2##F, ARG3##F) != RES##F \
|
77 |
|
|
|| CKSGN_F(__builtin_##FUNC##f(ARG1##F,ARG2##F,ARG3##F),RES##F)) \
|
78 |
|
|
link_error(__LINE__); \
|
79 |
|
|
if (__builtin_##FUNC(ARG1, ARG2, ARG3) != RES \
|
80 |
|
|
|| CKSGN(__builtin_##FUNC(ARG1,ARG2,ARG3),RES)) \
|
81 |
|
|
link_error(__LINE__); \
|
82 |
|
|
if (__builtin_##FUNC##l(ARG1##L, ARG2##L, ARG3##L) != RES##L \
|
83 |
|
|
|| CKSGN_L(__builtin_##FUNC##l(ARG1##L,ARG2##L,ARG3##L),RES##L)) \
|
84 |
|
|
link_error(__LINE__); \
|
85 |
|
|
} while (0)
|
86 |
|
|
|
87 |
|
|
/* Test that for FUNC(ARG, &ARG_S, &ARG_C);
|
88 |
|
|
assert (ARG_S == RES_S && ARG_C == RES_C);. */
|
89 |
|
|
#define TESTIT_2P(FUNC,ARG,ARG_S,ARG_C,RES_S,RES_C) do { \
|
90 |
|
|
__builtin_##FUNC##f(ARG##F, &ARG_S##f, &ARG_C##f); \
|
91 |
|
|
if (ARG_S##f != (RES_S##F) || ARG_C##f != (RES_C##F)) \
|
92 |
|
|
link_error(__LINE__); \
|
93 |
|
|
__builtin_##FUNC(ARG, &ARG_S, &ARG_C); \
|
94 |
|
|
if (ARG_S != (RES_S) || ARG_C != (RES_C)) \
|
95 |
|
|
link_error(__LINE__); \
|
96 |
|
|
__builtin_##FUNC##l(ARG##L, &ARG_S##l, &ARG_C##l); \
|
97 |
|
|
if (ARG_S##l != (RES_S##L) || ARG_C##l != (RES_C##L)) \
|
98 |
|
|
link_error(__LINE__); \
|
99 |
|
|
} while (0)
|
100 |
|
|
|
101 |
|
|
/* Test that for FUNC(ARG, &ARG_S, &ARG_C);
|
102 |
|
|
assert (LOW_S < ARG_S < HI_S && LOW_C < ARG_C < HI_C);. */
|
103 |
|
|
#define TESTIT_2P_R(FUNC,ARG,ARG_S,ARG_C,LOW_S,HI_S,LOW_C,HI_C) do { \
|
104 |
|
|
__builtin_##FUNC##f(ARG##F, &ARG_S##f, &ARG_C##f); \
|
105 |
|
|
if (ARG_S##f <= (LOW_S##F) || ARG_S##f >= (HI_S##F) \
|
106 |
|
|
|| ARG_C##f <= (LOW_C##F) || ARG_C##f >= (HI_C##F)) \
|
107 |
|
|
link_error(__LINE__); \
|
108 |
|
|
__builtin_##FUNC(ARG, &ARG_S, &ARG_C); \
|
109 |
|
|
if (ARG_S <= (LOW_S) || ARG_S >= (HI_S) \
|
110 |
|
|
|| ARG_C <= (LOW_C) || ARG_C >= (HI_C)) \
|
111 |
|
|
link_error(__LINE__); \
|
112 |
|
|
__builtin_##FUNC##l(ARG##L, &ARG_S##l, &ARG_C##l); \
|
113 |
|
|
if (ARG_S##l <= (LOW_S##L) || ARG_S##l >= (HI_S##L) \
|
114 |
|
|
|| ARG_C##l <= (LOW_C##L) || ARG_C##l >= (HI_C##L)) \
|
115 |
|
|
link_error(__LINE__); \
|
116 |
|
|
} while (0)
|
117 |
|
|
|
118 |
|
|
int main (void)
|
119 |
|
|
{
|
120 |
|
|
#ifdef __OPTIMIZE__
|
121 |
|
|
float sf, cf, oneF = 1.0F;
|
122 |
|
|
double s, c, one = 1.0;
|
123 |
|
|
long double sl, cl, oneL = 1.0L;
|
124 |
|
|
#endif
|
125 |
|
|
|
126 |
|
|
TESTIT_R (asin, -1.0, -3.15/2.0, -3.14/2.0); /* asin(-1) == -pi/2 */
|
127 |
|
|
TESTIT (asin, 0.0, 0.0); /* asin(0) == 0 */
|
128 |
|
|
TESTIT (asin, -0.0, -0.0); /* asin(-0) == -0 */
|
129 |
|
|
TESTIT_R (asin, 1.0, 3.14/2.0, 3.15/2.0); /* asin(1) == pi/2 */
|
130 |
|
|
|
131 |
|
|
TESTIT_R (acos, -1.0, 3.14, 3.15); /* acos(-1) == pi */
|
132 |
|
|
TESTIT_R (acos, 0.0, 3.14/2.0, 3.15/2.0); /* acos(0) == pi/2 */
|
133 |
|
|
TESTIT_R (acos, -0.0, 3.14/2.0, 3.15/2.0); /* acos(-0) == pi/2 */
|
134 |
|
|
TESTIT (acos, 1.0, 0.0); /* acos(1) == 0 */
|
135 |
|
|
|
136 |
|
|
TESTIT_R (atan, -1.0, -3.15/4.0, -3.14/4.0); /* atan(-1) == -pi/4 */
|
137 |
|
|
TESTIT (atan, 0.0, 0.0); /* atan(0) == 0 */
|
138 |
|
|
TESTIT (atan, -0.0, -0.0); /* atan(-0) == -0 */
|
139 |
|
|
TESTIT_R (atan, 1.0, 3.14/4.0, 3.15/4.0); /* atan(1) == pi/4 */
|
140 |
|
|
|
141 |
|
|
TESTIT_R (asinh, -1.0, -0.89, -0.88); /* asinh(-1) == -0.881... */
|
142 |
|
|
TESTIT (asinh, 0.0, 0.0); /* asinh(0) == 0 */
|
143 |
|
|
TESTIT (asinh, -0.0, -0.0); /* asinh(-0) == -0 */
|
144 |
|
|
TESTIT_R (asinh, 1.0, 0.88, 0.89); /* asinh(1) == 0.881... */
|
145 |
|
|
|
146 |
|
|
TESTIT (acosh, 1.0, 0.0); /* acosh(1) == 0. */
|
147 |
|
|
TESTIT_R (acosh, 2.0, 1.31, 1.32); /* acosh(2) == 1.316... */
|
148 |
|
|
|
149 |
|
|
TESTIT_R (atanh, -0.5, -0.55, -0.54); /* atanh(-0.5) == -0.549... */
|
150 |
|
|
TESTIT (atanh, 0.0, 0.0); /* atanh(0) == 0 */
|
151 |
|
|
TESTIT (atanh, -0.0, -0.0); /* atanh(-0) == -0 */
|
152 |
|
|
TESTIT_R (atanh, 0.5, 0.54, 0.55); /* atanh(0.5) == 0.549... */
|
153 |
|
|
|
154 |
|
|
TESTIT_R (sin, -1.0, -0.85, -0.84); /* sin(-1) == -0.841... */
|
155 |
|
|
TESTIT (sin, 0.0, 0.0); /* sin(0) == 0 */
|
156 |
|
|
TESTIT (sin, -0.0, -0.0); /* sin(-0) == -0 */
|
157 |
|
|
TESTIT_R (sin, 1.0, 0.84, 0.85); /* sin(1) == 0.841... */
|
158 |
|
|
|
159 |
|
|
TESTIT_R (cos, -1.0, 0.54, 0.55); /* cos(-1) == 0.5403... */
|
160 |
|
|
TESTIT (cos, 0.0, 1.0); /* cos(0) == 1 */
|
161 |
|
|
TESTIT (cos, -0.0, 1.0); /* cos(-0) == 1 */
|
162 |
|
|
TESTIT_R (cos, 1.0, 0.54, 0.55); /* cos(1) == 0.5403... */
|
163 |
|
|
|
164 |
|
|
TESTIT_R (tan, -1.0, -1.56, 1.55); /* tan(-1) == -1.557... */
|
165 |
|
|
TESTIT (tan, 0.0, 0.0); /* tan(0) == 0 */
|
166 |
|
|
TESTIT (tan, -0.0, -0.0); /* tan(-0) == -0 */
|
167 |
|
|
TESTIT_R (tan, 1.0, 1.55, 1.56); /* tan(1) == 1.557... */
|
168 |
|
|
|
169 |
|
|
#ifdef __OPTIMIZE__
|
170 |
|
|
/* These tests rely on propagating the variables s, c and one, which
|
171 |
|
|
happens only when optimization is turned on. */
|
172 |
|
|
TESTIT_2P_R (sincos, -1.0, s, c, -0.85, -0.84, 0.54, 0.55); /* (s==-0.841..., c==0.5403...) */
|
173 |
|
|
TESTIT_2P (sincos, 0.0, s, c, 0.0, 1.0); /* (s==0, c==1) */
|
174 |
|
|
TESTIT_2P (sincos, -0.0, s, c, -0.0, 1.0); /* (s==-0, c==1) */
|
175 |
|
|
TESTIT_2P_R (sincos, 1.0, s, c, 0.84, 0.85, 0.54, 0.55); /* (s==0.841..., c==0.5403...) */
|
176 |
|
|
TESTIT_2P_R (sincos, one, s, c, 0.84, 0.85, 0.54, 0.55); /* (s==0.841..., c==0.5403...) */
|
177 |
|
|
TESTIT_2P_R (sincos, -one, s, c, -0.85, -0.84, 0.54, 0.55); /* (s==-0.841..., c==0.5403...) */
|
178 |
|
|
#endif
|
179 |
|
|
|
180 |
|
|
TESTIT_R (sinh, -1.0, -1.18, -1.17); /* sinh(-1) == -1.175... */
|
181 |
|
|
TESTIT (sinh, 0.0, 0.0); /* sinh(0) == 0 */
|
182 |
|
|
TESTIT (sinh, -0.0, -0.0); /* sinh(-0) == -0 */
|
183 |
|
|
TESTIT_R (sinh, 1.0, 1.17, 1.18); /* sinh(1) == 1.175... */
|
184 |
|
|
|
185 |
|
|
TESTIT_R (cosh, -1.0, 1.54, 1.55); /* cosh(-1) == 1.543... */
|
186 |
|
|
TESTIT (cosh, 0.0, 1.0); /* cosh(0) == 1 */
|
187 |
|
|
TESTIT (cosh, -0.0, 1.0); /* cosh(-0) == 1 */
|
188 |
|
|
TESTIT_R (cosh, 1.0, 1.54, 1.55); /* cosh(1) == 1.543... */
|
189 |
|
|
|
190 |
|
|
TESTIT_R (tanh, -1.0, -0.77, -0.76); /* tanh(-1) == -0.761... */
|
191 |
|
|
TESTIT (tanh, -0.0, -0.0); /* tanh(-0) == -0 */
|
192 |
|
|
TESTIT (tanh, 0.0, 0.0); /* tanh(0) == 0 */
|
193 |
|
|
TESTIT_R (tanh, 1.0, 0.76, 0.77); /* tanh(1) == 0.761... */
|
194 |
|
|
|
195 |
|
|
TESTIT_R (exp, -1.0, 0.36, 0.37); /* exp(-1) == 1/e */
|
196 |
|
|
TESTIT (exp, -0.0, 1.0); /* exp(-0) == 1 */
|
197 |
|
|
TESTIT (exp, 0.0, 1.0); /* exp(0) == 1 */
|
198 |
|
|
TESTIT_R (exp, 1.0, 2.71, 2.72); /* exp(1) == e */
|
199 |
|
|
|
200 |
|
|
TESTIT (exp2, -1.0, 0.5); /* exp2(-1) == 1/2 */
|
201 |
|
|
TESTIT (exp2, -0.0, 1.0); /* exp2(-0) == 1 */
|
202 |
|
|
TESTIT (exp2, 0.0, 1.0); /* exp2(0) == 1 */
|
203 |
|
|
TESTIT (exp2, 1.0, 2.0); /* exp2(1) == 2 */
|
204 |
|
|
|
205 |
|
|
TESTIT (exp10, -1.0, 0.1); /* exp10(-1) == 1/10 */
|
206 |
|
|
TESTIT (exp10, -0.0, 1.0); /* exp10(-0) == 1 */
|
207 |
|
|
TESTIT (exp10, 0.0, 1.0); /* exp10(0) == 1 */
|
208 |
|
|
TESTIT (exp10, 1.0, 10.0); /* exp10(1) == 10 */
|
209 |
|
|
|
210 |
|
|
TESTIT (pow10, -1.0, 0.1); /* pow10(-1) == 1/10 */
|
211 |
|
|
TESTIT (pow10, -0.0, 1.0); /* pow10(-0) == 1 */
|
212 |
|
|
TESTIT (pow10, 0.0, 1.0); /* pow10(0) == 1 */
|
213 |
|
|
TESTIT (pow10, 1.0, 10.0); /* pow10(1) == 10 */
|
214 |
|
|
|
215 |
|
|
TESTIT_R (expm1, -1.0, -0.64, -0.63); /* expm1(-1) == 1/e - 1 */
|
216 |
|
|
TESTIT (expm1, -0.0, -0.0); /* expm1(-0) == 0 */
|
217 |
|
|
TESTIT (expm1, 0.0, 0.0); /* expm1(0) == 0 */
|
218 |
|
|
TESTIT_R (expm1, 1.0, 1.71, 1.72); /* expm1(1) == e - 1 */
|
219 |
|
|
|
220 |
|
|
TESTIT (log, 1.0, 0.0); /* log(1) == 0 */
|
221 |
|
|
TESTIT_R (log, M_E, 0.99, 1.01); /* log(e) == 1.000... */
|
222 |
|
|
TESTIT_R (log, M_E*M_E, 1.99, 2.01); /* log(e*e) == 2.000... */
|
223 |
|
|
|
224 |
|
|
TESTIT (log2, 1.0, 0.0); /* log2(1) == 0 */
|
225 |
|
|
TESTIT (log2, 2.0, 1.0); /* log2(2) == 1 */
|
226 |
|
|
TESTIT (log2, 4.0, 2.0); /* log2(4) == 2 */
|
227 |
|
|
|
228 |
|
|
TESTIT (log10, 1.0, 0.0); /* log10(1) == 0 */
|
229 |
|
|
TESTIT (log10, 10.0, 1.0); /* log10(10) == 1 */
|
230 |
|
|
TESTIT (log10, 100.0, 2.0); /* log10(100) == 2 */
|
231 |
|
|
|
232 |
|
|
TESTIT (log1p, 0.0, 0.0); /* log1p(0) == 0 */
|
233 |
|
|
TESTIT (log1p, -0.0, -0.0); /* log1p(-0) == -0 */
|
234 |
|
|
TESTIT_R (log1p, M_E-1, 0.99, 1.01); /* log1p(e-1) == 1.000... */
|
235 |
|
|
TESTIT_R (log1p, M_E*M_E-1, 1.99, 2.01); /* log1p(e*e-1) == 2.000... */
|
236 |
|
|
|
237 |
|
|
TESTIT (cbrt, -0.0, -0.0); /* cbrt(-0) == -0 */
|
238 |
|
|
TESTIT (cbrt, 0.0, 0.0); /* cbrt(0) == 0 */
|
239 |
|
|
TESTIT (cbrt, 1.0, 1.0); /* cbrt(1) == 1 */
|
240 |
|
|
TESTIT (cbrt, -1.0, -1.0); /* cbrt(-1) == -1 */
|
241 |
|
|
TESTIT (cbrt, 8.0, 2.0); /* cbrt(8) == 2 */
|
242 |
|
|
TESTIT (cbrt, -8.0, -2.0); /* cbrt(-8) == -2 */
|
243 |
|
|
|
244 |
|
|
TESTIT (erf, -0.0, -0.0); /* erf(-0) == -0 */
|
245 |
|
|
TESTIT (erf, 0.0, 0.0); /* erf(0) == 0 */
|
246 |
|
|
TESTIT_R (erf, 1.0, 0.84, 0.85); /* erf(1) == 0.842... */
|
247 |
|
|
TESTIT_R (erf, -1.0, -0.85, -0.84); /* erf(-1) == -0.842... */
|
248 |
|
|
|
249 |
|
|
TESTIT (erfc, -0.0, 1.0); /* erfc(-0) == 1 */
|
250 |
|
|
TESTIT (erfc, 0.0, 1.0); /* erfc(0) == 1 */
|
251 |
|
|
TESTIT_R (erfc, 1.0, 0.15, 0.16); /* erfc(1) == 0.157... */
|
252 |
|
|
TESTIT_R (erfc, -1.0, 1.84, 1.85); /* erfc(-1) == 1.842... */
|
253 |
|
|
|
254 |
|
|
TESTIT_R (tgamma, -4.5, -0.061, -0.060); /* tgamma(-4.5) == -0.06001... */
|
255 |
|
|
TESTIT_R (tgamma, -3.5, 0.27, 0.28); /* tgamma(-3.5) == 0.27008... */
|
256 |
|
|
TESTIT_R (tgamma, -2.5, -0.95, -0.94); /* tgamma(-2.5) == -0.945... */
|
257 |
|
|
TESTIT_R (tgamma, -1.5, 2.36, 2.37); /* tgamma(-1.5) == 2.363... */
|
258 |
|
|
TESTIT_R (tgamma, -0.5, -3.55, -3.54); /* tgamma(-0.5) == -3.544... */
|
259 |
|
|
TESTIT_R (tgamma, 0.5, 1.77, 1.78); /* tgamma(0.5) == 1.772... */
|
260 |
|
|
TESTIT (tgamma, 1.0, 1.0); /* tgamma(1) == 1 */
|
261 |
|
|
TESTIT_R (tgamma, 1.5, 0.88, 0.89); /* tgamma(1.5) == 0.886... */
|
262 |
|
|
TESTIT (tgamma, 2.0, 1.0); /* tgamma(2) == 1 */
|
263 |
|
|
TESTIT_R (tgamma, 2.5, 1.32, 1.33); /* tgamma(2.5) == 1.329... */
|
264 |
|
|
TESTIT (tgamma, 3.0, 2.0); /* tgamma(3) == 2 */
|
265 |
|
|
TESTIT_R (tgamma, 3.5, 3.32, 3.33); /* tgamma(3.5) == 3.323... */
|
266 |
|
|
TESTIT (tgamma, 4.0, 6.0); /* tgamma(4) == 6 */
|
267 |
|
|
TESTIT_R (tgamma, 4.5, 11.63, 11.64); /* tgamma(4.5) == 11.631... */
|
268 |
|
|
|
269 |
|
|
TESTIT2 (pow, 3.0, 4.0, 81.0); /* pow(3,4) == 81 */
|
270 |
|
|
TESTIT2 (pow, -3.0, 5.0, -243.0); /* pow(-3,5) == -243 */
|
271 |
|
|
TESTIT2 (pow, 16.0, 0.25, 2.0); /* pow(16,1/4) == 2 */
|
272 |
|
|
TESTIT2 (pow, 4.0, -2.0, 0.0625); /* pow(4,-2) == 1/16 */
|
273 |
|
|
TESTIT2 (pow, -2.0, -3.0, -0.125); /* pow(-2,-3) == -1/8 */
|
274 |
|
|
TESTIT2_R (pow, -1.5, -3.0, -0.297, -0.296); /* pow(-1.5,-3) == -1/3.375 */
|
275 |
|
|
|
276 |
|
|
TESTIT2 (hypot, 0.0, 0.0, 0.0); /* hypot(0,0) == 0 */
|
277 |
|
|
TESTIT2 (hypot, -0.0, 0.0, 0.0); /* hypot(-0,0) == 0 */
|
278 |
|
|
TESTIT2 (hypot, 0.0, -0.0, 0.0); /* hypot(0,-0) == 0 */
|
279 |
|
|
TESTIT2 (hypot, -0.0, -0.0, 0.0); /* hypot(-0,-0) == 0 */
|
280 |
|
|
TESTIT2 (hypot, 3.0, 4.0, 5.0); /* hypot(3,4) == 5 */
|
281 |
|
|
TESTIT2 (hypot, -3.0, 4.0, 5.0); /* hypot(-3,4) == 5 */
|
282 |
|
|
TESTIT2 (hypot, 3.0, -4.0, 5.0); /* hypot(3,-4) == 5 */
|
283 |
|
|
TESTIT2 (hypot, -3.0, -4.0, 5.0); /* hypot(-3,-4) == 5 */
|
284 |
|
|
TESTIT2_R (hypot, 4.0, 5.0, 6.40, 6.41); /* hypot(4,5) == 6.403... */
|
285 |
|
|
|
286 |
|
|
TESTIT2 (atan2, 0.0, 0.0, 0.0); /* atan2(0,0) == 0 */
|
287 |
|
|
TESTIT2 (atan2, -0.0, 0.0, -0.0); /* atan2(-0,0) == -0 */
|
288 |
|
|
TESTIT2_R (atan2, 0.0, -0.0, 3.14, 3.15); /* atan2(0,-0) == pi */
|
289 |
|
|
TESTIT2_R (atan2, -0.0, -0.0, -3.15, -3.14); /* atan2(-0,-0) == -pi */
|
290 |
|
|
TESTIT2_R (atan2, 0.0, -1.0, 3.14, 3.15); /* atan2(0,-1) == pi */
|
291 |
|
|
TESTIT2_R (atan2, -0.0, -1.0, -3.15, -3.14); /* atan2(-0,-1) == -pi */
|
292 |
|
|
TESTIT2 (atan2, 0.0, 1.0, 0.0); /* atan2(0,1) == 0 */
|
293 |
|
|
TESTIT2 (atan2, -0.0, 1.0, -0.0); /* atan2(-0,1) == -0 */
|
294 |
|
|
TESTIT2_R (atan2, -1.0, 0.0, -1.58, -1.57); /* atan2(-1,0) == -pi/2 */
|
295 |
|
|
TESTIT2_R (atan2, 1.0, 0.0, 1.57, 1.58); /* atan2(1,0) == pi/2 */
|
296 |
|
|
|
297 |
|
|
TESTIT2 (fdim, 0.0, 0.0, 0.0); /* fdim(0,0) == 0 */
|
298 |
|
|
TESTIT2 (fdim, -0.0, 0.0, 0.0); /* fdim(-0,0) == 0 */
|
299 |
|
|
TESTIT2 (fdim, 0.0, -0.0, 0.0); /* fdim(0,-0) == 0 */
|
300 |
|
|
TESTIT2 (fdim, -0.0, -0.0, 0.0); /* fdim(-0,-0) == 0 */
|
301 |
|
|
TESTIT2 (fdim, 5.0, 5.0, 0.0); /* fdim(5,5) == 0 */
|
302 |
|
|
TESTIT2 (fdim, 5.0, 6.0, 0.0); /* fdim(5,6) == 0 */
|
303 |
|
|
TESTIT2 (fdim, 6.0, 5.0, 1.0); /* fdim(6,5) == 1 */
|
304 |
|
|
TESTIT2 (fdim, -5.0, -6.0, 1.0); /* fdim(-5,-6) == 1 */
|
305 |
|
|
TESTIT2 (fdim, -6.0, -5.0, 0.0); /* fdim(-6,-5) == 0 */
|
306 |
|
|
|
307 |
|
|
TESTIT2 (fmin, 5.0, 6.0, 5.0); /* fmin(5,6) == 5 */
|
308 |
|
|
TESTIT2 (fmin, 6.0, 5.0, 5.0); /* fmin(6,5) == 5 */
|
309 |
|
|
TESTIT2 (fmin, -5.0, -6.0, -6.0); /* fmin(-5,-6) == -6 */
|
310 |
|
|
TESTIT2 (fmin, -6.0, -5.0, -6.0); /* fmin(-6,-5) == -6 */
|
311 |
|
|
TESTIT2 (fmin, -0.0, 0.0, -0.0); /* fmin(-0,0) == -0 */
|
312 |
|
|
TESTIT2 (fmin, 0.0, -0.0, -0.0); /* fmin(-0,0) == -0 */
|
313 |
|
|
|
314 |
|
|
TESTIT2 (fmax, 5.0, 6.0, 6.0); /* fmax(5,6) == 6 */
|
315 |
|
|
TESTIT2 (fmax, 6.0, 5.0, 6.0); /* fmax(6,5) == 6 */
|
316 |
|
|
TESTIT2 (fmax, -5.0, -6.0, -5.0); /* fmax(-5,-6) == -5 */
|
317 |
|
|
TESTIT2 (fmax, -6.0, -5.0, -5.0); /* fmax(-6,-5) == -5 */
|
318 |
|
|
TESTIT2 (fmax, -0.0, 0.0, 0.0); /* fmax(-0,0) == 0 */
|
319 |
|
|
TESTIT2 (fmax, 0.0, -0.0, 0.0); /* fmax(-0,0) == 0 */
|
320 |
|
|
|
321 |
|
|
TESTIT3 (fma, 2.0, 3.0, 4.0, 10.0); /* fma(2,3,4) == 10 */
|
322 |
|
|
TESTIT3 (fma, 2.0, -3.0, 4.0, -2.0); /* fma(2,-3,4) == -2 */
|
323 |
|
|
TESTIT3 (fma, 2.0, 3.0, -4.0, 2.0); /* fma(2,3,-4) == 2 */
|
324 |
|
|
TESTIT3 (fma, 2.0, -3.0, -4.0, -10.0); /* fma(2,-3,-4) == -10 */
|
325 |
|
|
TESTIT3 (fma, -2.0, -3.0, -4.0, 2.0); /* fma(-2,-3,-4) == 2 */
|
326 |
|
|
TESTIT3 (fma, 6.0, -0.0, 0.0, 0.0); /* fma(6,-0,0) == 0 */
|
327 |
|
|
TESTIT3 (fma, -0.0, 6.0, 0.0, 0.0); /* fma(-0,6,0) == 0 */
|
328 |
|
|
TESTIT3 (fma, 6.0, -0.0, -0.0, -0.0); /* fma(6,-0,-0) == -0 */
|
329 |
|
|
TESTIT3 (fma, -0.0, 6.0, -0.0, -0.0); /* fma(-0,6,-0) == -0 */
|
330 |
|
|
TESTIT3 (fma, 0.0, 0.0, 0.0, 0.0); /* fma(0,0,0) == 0 */
|
331 |
|
|
TESTIT3 (fma, -0.0, 0.0, 0.0, 0.0); /* fma(-0,0,0) == 0 */
|
332 |
|
|
TESTIT3 (fma, 0.0, -0.0, 0.0, 0.0); /* fma(0,-0,0) == 0 */
|
333 |
|
|
TESTIT3 (fma, -0.0, -0.0, 0.0, 0.0); /* fma(-0,-0,0) == 0 */
|
334 |
|
|
TESTIT3 (fma, 0.0, 0.0, -0.0, 0.0); /* fma(0,0,-0) == 0 */
|
335 |
|
|
TESTIT3 (fma, -0.0, 0.0, -0.0, -0.0); /* fma(-0,0,-0) == -0 */
|
336 |
|
|
TESTIT3 (fma, 0.0, -0.0, -0.0, -0.0); /* fma(0,-0,-0) == -0 */
|
337 |
|
|
TESTIT3 (fma, -0.0, -0.0, -0.0, 0.0); /* fma(-0,-0,-0) == 0 */
|
338 |
|
|
|
339 |
|
|
if (__builtin_fmaf(__FLT_MAX__, 2.0F, -__FLT_MAX__) != __FLT_MAX__)
|
340 |
|
|
link_error (__LINE__);
|
341 |
|
|
if (__builtin_fmaf(2.0F,__FLT_MAX__, -__FLT_MAX__) != __FLT_MAX__)
|
342 |
|
|
link_error (__LINE__);
|
343 |
|
|
if (__builtin_fmaf(__FLT_MIN__, 0.5F, __FLT_MIN__) != __FLT_MIN__*1.5F)
|
344 |
|
|
link_error (__LINE__);
|
345 |
|
|
if (__builtin_fmaf(0.5F,__FLT_MIN__, __FLT_MIN__) != __FLT_MIN__*1.5F)
|
346 |
|
|
link_error (__LINE__);
|
347 |
|
|
|
348 |
|
|
if (__builtin_fma(__DBL_MAX__, 2.0, -__DBL_MAX__) != __DBL_MAX__)
|
349 |
|
|
link_error (__LINE__);
|
350 |
|
|
if (__builtin_fma(2.0,__DBL_MAX__, -__DBL_MAX__) != __DBL_MAX__)
|
351 |
|
|
link_error (__LINE__);
|
352 |
|
|
if (__builtin_fma(__DBL_MIN__, 0.5, __DBL_MIN__) != __DBL_MIN__*1.5)
|
353 |
|
|
link_error (__LINE__);
|
354 |
|
|
if (__builtin_fma(0.5,__DBL_MIN__, __DBL_MIN__) != __DBL_MIN__*1.5)
|
355 |
|
|
link_error (__LINE__);
|
356 |
|
|
|
357 |
|
|
if (__builtin_fmal(__LDBL_MAX__, 2.0L, -__LDBL_MAX__) != __LDBL_MAX__)
|
358 |
|
|
link_error (__LINE__);
|
359 |
|
|
if (__builtin_fmal(2.0L,__LDBL_MAX__, -__LDBL_MAX__) != __LDBL_MAX__)
|
360 |
|
|
link_error (__LINE__);
|
361 |
|
|
if (__builtin_fmal(__LDBL_MIN__, 0.5L, __LDBL_MIN__) != __LDBL_MIN__*1.5L)
|
362 |
|
|
link_error (__LINE__);
|
363 |
|
|
if (__builtin_fmal(0.5L,__LDBL_MIN__, __LDBL_MIN__) != __LDBL_MIN__*1.5L)
|
364 |
|
|
link_error (__LINE__);
|
365 |
|
|
|
366 |
|
|
TESTIT (sqrt, -0.0, -0.0); /* sqrt(-0) == -0 */
|
367 |
|
|
TESTIT (sqrt, 0.0, 0.0); /* sqrt(0) == 0 */
|
368 |
|
|
TESTIT (sqrt, 1.0, 1.0); /* sqrt(1) == 1 */
|
369 |
|
|
TESTIT (sqrt, 4.0, 2.0); /* sqrt(4) == 2 */
|
370 |
|
|
TESTIT_R (sqrt, 1.5, 1.22, 1.23); /* sqrt(1.5) == 1.224... */
|
371 |
|
|
TESTIT_R (sqrt, 2.0, 1.41, 1.42); /* sqrt(2) == 1.414... */
|
372 |
|
|
|
373 |
|
|
return 0;
|
374 |
|
|
}
|