1 /* 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 }