1/* Test femode_t functions.
2 Copyright (C) 2016-2022 Free Software Foundation, Inc.
3 This file is part of the GNU C Library.
4
5 The GNU C Library is free software; you can redistribute it and/or
6 modify it under the terms of the GNU Lesser General Public
7 License as published by the Free Software Foundation; either
8 version 2.1 of the License, or (at your option) any later version.
9
10 The GNU C Library is distributed in the hope that it will be useful,
11 but WITHOUT ANY WARRANTY; without even the implied warranty of
12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 Lesser General Public License for more details.
14
15 You should have received a copy of the GNU Lesser General Public
16 License along with the GNU C Library; if not, see
17 <https://www.gnu.org/licenses/>. */
18
19#include <fenv.h>
20#include <stdio.h>
21#include <math-tests.h>
22
23static int
24test_mmee (int mode1, int mode2, int exc1, int exc2)
25{
26 int result = 0;
27 printf (format: "testing %x %x %x %x\n", (unsigned int) mode1, (unsigned int) mode2,
28 (unsigned int) exc1, (unsigned int) exc2);
29
30 feclearexcept (FE_ALL_EXCEPT);
31 int ret = fesetround (rounding_direction: mode1);
32 if (ret != 0)
33 {
34 if (ROUNDING_TESTS (float, mode1))
35 {
36 puts (s: "first fesetround failed unexpectedly");
37 result = 1;
38 }
39 else
40 puts (s: "first fesetround failed, cannot test");
41 return result;
42 }
43 ret = fesetexcept (excepts: exc1);
44 if (ret != 0)
45 {
46 if (EXCEPTION_TESTS (float) || exc1 == 0)
47 {
48 puts (s: "first fesetexcept failed unexpectedly");
49 result = 1;
50 }
51 else
52 puts (s: "first fesetexcept failed, cannot test");
53 return result;
54 }
55 femode_t saved;
56 ret = fegetmode (modep: &saved);
57 if (ret != 0)
58 {
59 puts (s: "fegetmode failed");
60 result = 1;
61 return result;
62 }
63 feclearexcept (FE_ALL_EXCEPT);
64 ret = fesetround (rounding_direction: mode2);
65 if (ret != 0)
66 {
67 if (ROUNDING_TESTS (float, mode2))
68 {
69 puts (s: "second fesetround failed unexpectedly");
70 result = 1;
71 }
72 else
73 puts (s: "second fesetround failed, cannot test");
74 return result;
75 }
76 ret = fesetexcept (excepts: exc2);
77 if (ret != 0)
78 {
79 if (EXCEPTION_TESTS (float) || exc2 == 0)
80 {
81 puts (s: "second fesetexcept failed unexpectedly");
82 result = 1;
83 }
84 else
85 puts (s: "second fesetexcept failed, cannot test");
86 return result;
87 }
88 ret = fesetmode (modep: &saved);
89 if (ret != 0)
90 {
91 puts (s: "fesetmode failed");
92 result = 1;
93 return result;
94 }
95 /* Verify that the rounding mode was restored but the exception
96 flags remain unchanged. */
97 ret = fegetround ();
98 if (ret != mode1)
99 {
100 printf (format: "restored rounding mode %x not %x\n", (unsigned int) ret,
101 (unsigned int) mode1);
102 result = 1;
103 }
104 ret = fetestexcept (FE_ALL_EXCEPT);
105 if (ret != exc2)
106 {
107 printf (format: "exceptions %x not %x\n", (unsigned int) ret,
108 (unsigned int) exc2);
109 result = 1;
110 }
111 /* Likewise, with default modes. */
112 ret = fesetmode (FE_DFL_MODE);
113 if (ret != 0)
114 {
115 puts (s: "fesetmode (FE_DFL_MODE) failed");
116 result = 1;
117 return result;
118 }
119 ret = fegetround ();
120 if (ret != FE_TONEAREST)
121 {
122 printf (format: "FE_DFL_MODE rounding mode %x not %x\n", (unsigned int) ret,
123 (unsigned int) FE_TONEAREST);
124 result = 1;
125 }
126 ret = fetestexcept (FE_ALL_EXCEPT);
127 if (ret != exc2)
128 {
129 printf (format: "FE_DFL_MODE exceptions %x not %x\n", (unsigned int) ret,
130 (unsigned int) exc2);
131 result = 1;
132 }
133 return result;
134}
135
136static int
137test_mme (int mode1, int mode2, int exc1)
138{
139 int result = 0;
140
141 result |= test_mmee (mode1, mode2, exc1, exc2: 0);
142 result |= test_mmee (mode1, mode2, exc1, FE_ALL_EXCEPT);
143#ifdef FE_DIVBYZERO
144 result |= test_mmee (mode1, mode2, exc1, FE_DIVBYZERO);
145#endif
146#ifdef FE_INEXACT
147 result |= test_mmee (mode1, mode2, exc1, FE_INEXACT);
148#endif
149#ifdef FE_INVALID
150 result |= test_mmee (mode1, mode2, exc1, FE_INVALID);
151#endif
152#ifdef FE_OVERFLOW
153 result |= test_mmee (mode1, mode2, exc1, FE_OVERFLOW);
154#endif
155#ifdef FE_UNDERFLOW
156 result |= test_mmee (mode1, mode2, exc1, FE_UNDERFLOW);
157#endif
158
159 return result;
160}
161
162static int
163test_mm (int mode1, int mode2)
164{
165 int result = 0;
166
167 result |= test_mme (mode1, mode2, exc1: 0);
168 result |= test_mme (mode1, mode2, FE_ALL_EXCEPT);
169#ifdef FE_DIVBYZERO
170 result |= test_mme (mode1, mode2, FE_DIVBYZERO);
171#endif
172#ifdef FE_INEXACT
173 result |= test_mme (mode1, mode2, FE_INEXACT);
174#endif
175#ifdef FE_INVALID
176 result |= test_mme (mode1, mode2, FE_INVALID);
177#endif
178#ifdef FE_OVERFLOW
179 result |= test_mme (mode1, mode2, FE_OVERFLOW);
180#endif
181#ifdef FE_UNDERFLOW
182 result |= test_mme (mode1, mode2, FE_UNDERFLOW);
183#endif
184
185 return result;
186}
187
188static int
189test_m (int mode1)
190{
191 int result = 0;
192
193#ifdef FE_DOWNWARD
194 result |= test_mm (mode1, FE_DOWNWARD);
195#endif
196#ifdef FE_TONEAREST
197 result |= test_mm (mode1, FE_TONEAREST);
198#endif
199#ifdef FE_TOWARDZERO
200 result |= test_mm (mode1, FE_TOWARDZERO);
201#endif
202#ifdef FE_UPWARD
203 result |= test_mm (mode1, FE_UPWARD);
204#endif
205
206 return result;
207}
208
209static int
210do_test (void)
211{
212 int result = 0;
213
214#ifdef FE_DOWNWARD
215 result |= test_m (FE_DOWNWARD);
216#endif
217#ifdef FE_TONEAREST
218 result |= test_m (FE_TONEAREST);
219#endif
220#ifdef FE_TOWARDZERO
221 result |= test_m (FE_TOWARDZERO);
222#endif
223#ifdef FE_UPWARD
224 result |= test_m (FE_UPWARD);
225#endif
226
227 return result;
228}
229
230#define TEST_FUNCTION do_test ()
231#include "../test-skeleton.c"
232

source code of glibc/math/test-femode.c