(obstack_free): Explicitly convert __obj to
[kopensolaris-gnu/glibc.git] / math / test-fenv.c
1 /* Copyright (C) 1997, 1998 Free Software Foundation, Inc.
2    This file is part of the GNU C Library.
3    Contributed by Andreas Jaeger <aj@arthur.rhein-neckar.de> and
4    Ulrich Drepper <drepper@cygnus.com>, 1997.
5
6    The GNU C Library is free software; you can redistribute it and/or
7    modify it under the terms of the GNU Library General Public License as
8    published by the Free Software Foundation; either version 2 of the
9    License, or (at your option) any later version.
10
11    The GNU C Library is distributed in the hope that it will be useful,
12    but WITHOUT ANY WARRANTY; without even the implied warranty of
13    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14    Library General Public License for more details.
15
16    You should have received a copy of the GNU Library General Public
17    License along with the GNU C Library; see the file COPYING.LIB.  If not,
18    write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
19    Boston, MA 02111-1307, USA.  */
20
21 /* Tests for ISO C 9X 7.6: Floating-point environment  */
22
23 #ifndef _GNU_SOURCE
24 # define _GNU_SOURCE
25 #endif
26
27 #include <complex.h>
28 #include <math.h>
29 #include <float.h>
30 #include <fenv.h>
31
32 #include <errno.h>
33 #include <signal.h>
34 #include <stdlib.h>
35 #include <stdio.h>
36 #include <string.h>
37 #include <unistd.h>
38 #include <sys/wait.h>
39 #include <sys/resource.h>
40
41 /*
42   Since not all architectures might define all exceptions, we define
43   a private set and map accordingly.
44 */
45 #define NO_EXC 0
46 #define INEXACT_EXC 0x1
47 #define DIVBYZERO_EXC 0x2
48 #define UNDERFLOW_EXC 0x04
49 #define OVERFLOW_EXC 0x08
50 #define INVALID_EXC 0x10
51 #define ALL_EXC \
52         (INEXACT_EXC | DIVBYZERO_EXC | UNDERFLOW_EXC | OVERFLOW_EXC | \
53          INVALID_EXC)
54
55 static int count_errors;
56
57 /* Test whether a given exception was raised.  */
58 static void
59 test_single_exception (short int exception,
60                        short int exc_flag,
61                        fexcept_t fe_flag,
62                        const char *flag_name)
63 {
64   if (exception & exc_flag)
65     {
66       if (fetestexcept (fe_flag))
67         printf ("  Pass: Exception \"%s\" is set\n", flag_name);
68       else
69         {
70           printf ("  Fail: Exception \"%s\" is not set\n", flag_name);
71           ++count_errors;
72         }
73     }
74   else
75     {
76       if (fetestexcept (fe_flag))
77         {
78           printf ("  Fail: Exception \"%s\" is set\n", flag_name);
79           ++count_errors;
80         }
81       else
82         {
83           printf ("  Pass: Exception \"%s\" is not set\n", flag_name);
84         }
85     }
86 }
87
88 static void
89 test_exceptions (const char *test_name, short int exception,
90                  int ignore_inexact)
91 {
92   printf ("Test: %s\n", test_name);
93 #ifdef FE_DIVBYZERO
94   test_single_exception (exception, DIVBYZERO_EXC, FE_DIVBYZERO,
95                          "DIVBYZERO");
96 #endif
97 #ifdef FE_INVALID
98   test_single_exception (exception, INVALID_EXC, FE_INVALID,
99                          "INVALID");
100 #endif
101 #ifdef FE_INEXACT
102   if (!ignore_inexact)
103     test_single_exception (exception, INEXACT_EXC, FE_INEXACT,
104                            "INEXACT");
105 #endif
106 #ifdef FE_UNDERFLOW
107   test_single_exception (exception, UNDERFLOW_EXC, FE_UNDERFLOW,
108                          "UNDERFLOW");
109 #endif
110 #ifdef FE_OVERFLOW
111   test_single_exception (exception, OVERFLOW_EXC, FE_OVERFLOW,
112                          "OVERFLOW");
113 #endif
114 }
115
116 static void
117 print_rounding (int rounding)
118 {
119
120   switch (rounding)
121     {
122 #ifdef FE_TONEAREST
123     case FE_TONEAREST:
124       printf ("TONEAREST");
125       break;
126 #endif
127 #ifdef FE_UPWARD
128     case FE_UPWARD:
129       printf ("UPWARD");
130       break;
131 #endif
132 #ifdef FE_DOWNWARD
133     case FE_DOWNWARD:
134       printf ("DOWNWARD");
135       break;
136 #endif
137 #ifdef FE_TOWARDZERO
138     case FE_TOWARDZERO:
139       printf ("TOWARDZERO");
140       break;
141 #endif
142     }
143   printf (".\n");
144 }
145
146
147 static void
148 test_rounding (const char *test_name, int rounding_mode)
149 {
150   int curr_rounding = fegetround ();
151
152   printf ("Test: %s\n", test_name);
153   if (curr_rounding == rounding_mode)
154     {
155       printf ("  Pass: Rounding mode is ");
156       print_rounding (curr_rounding);
157     }
158   else
159     {
160       ++count_errors;
161       printf ("  Fail: Rounding mode is ");
162       print_rounding (curr_rounding);
163     }
164 }
165
166
167 static void
168 set_single_exc (const char *test_name, int fe_exc, fexcept_t exception)
169 {
170   char str[200];
171   /* The standard allows the inexact exception to be set together with the
172      underflow and overflow exceptions.  So ignore the inexact flag if the
173      others are raised.  */
174   int ignore_inexact = (fe_exc & (UNDERFLOW_EXC | OVERFLOW_EXC)) != 0;
175
176   strcpy (str, test_name);
177   strcat (str, ": set flag, with rest not set");
178   feclearexcept (FE_ALL_EXCEPT);
179   feraiseexcept (exception);
180   test_exceptions (str, fe_exc, ignore_inexact);
181
182   strcpy (str, test_name);
183   strcat (str, ": clear flag, rest also unset");
184   feclearexcept (exception);
185   test_exceptions (str, NO_EXC, ignore_inexact);
186
187   strcpy (str, test_name);
188   strcat (str, ": set flag, with rest set");
189   feraiseexcept (FE_ALL_EXCEPT ^ exception);
190   feraiseexcept (exception);
191   test_exceptions (str, ALL_EXC, 0);
192
193   strcpy (str, test_name);
194   strcat (str, ": clear flag, leave rest set");
195   feclearexcept (exception);
196   test_exceptions (str, ALL_EXC ^ fe_exc, 0);
197 }
198
199 static void
200 fe_tests (void)
201 {
202   /* clear all exceptions and test if all are cleared */
203   feclearexcept (FE_ALL_EXCEPT);
204   test_exceptions ("feclearexcept (FE_ALL_EXCEPT) clears all exceptions",
205                    NO_EXC, 0);
206
207   /* raise all exceptions and test if all are raised */
208   feraiseexcept (FE_ALL_EXCEPT);
209   test_exceptions ("feraiseexcept (FE_ALL_EXCEPT) raises all exceptions",
210                    ALL_EXC, 0);
211   feclearexcept (FE_ALL_EXCEPT);
212
213 #ifdef FE_DIVBYZERO
214   set_single_exc ("Set/Clear FE_DIVBYZERO", DIVBYZERO_EXC, FE_DIVBYZERO);
215 #endif
216 #ifdef FE_INVALID
217   set_single_exc ("Set/Clear FE_INVALID", INVALID_EXC, FE_INVALID);
218 #endif
219 #ifdef FE_INEXACT
220   set_single_exc ("Set/Clear FE_INEXACT", INEXACT_EXC, FE_INEXACT);
221 #endif
222 #ifdef FE_UNDERFLOW
223   set_single_exc ("Set/Clear FE_UNDERFLOW", UNDERFLOW_EXC, FE_UNDERFLOW);
224 #endif
225 #ifdef FE_OVERFLOW
226   set_single_exc ("Set/Clear FE_OVERFLOW", OVERFLOW_EXC, FE_OVERFLOW);
227 #endif
228 }
229
230 /* Test that program aborts with no masked interrupts */
231 static void
232 feenv_nomask_test (const char *flag_name, int fe_exc)
233 {
234 #if defined FE_NOMASK_ENV
235   int status;
236   pid_t pid;
237   fenv_t saved;
238
239   fegetenv (&saved);
240   errno = 0;
241   fesetenv (FE_NOMASK_ENV);
242   status = errno;
243   fesetenv (&saved);
244   if (status == ENOSYS)
245     {
246       printf ("Test: not testing FE_NOMASK_ENV, it isn't implemented.\n");
247       return;
248     }
249
250   printf ("Test: after fesetenv (FE_NOMASK_ENV) processes will abort\n");
251   printf ("      when feraiseexcept (%s) is called.\n", flag_name);
252   pid = fork ();
253   if (pid == 0)
254     {
255 #ifdef RLIMIT_CORE
256       /* Try to avoid dumping core.  */
257       struct rlimit core_limit;
258       core_limit.rlim_cur = 0;
259       core_limit.rlim_max = 0;
260       setrlimit (RLIMIT_CORE, &core_limit);
261 #endif
262
263       fesetenv (FE_NOMASK_ENV);
264       feraiseexcept (fe_exc);
265       exit (2);
266     }
267   else if (pid < 0)
268     {
269       if (errno != ENOSYS)
270         {
271           printf ("  Fail: Could not fork.\n");
272           ++count_errors;
273         }
274       else
275         printf ("  `fork' not implemented, test ignored.\n");
276     }
277   else {
278     if (waitpid (pid, &status, 0) != pid)
279       {
280         printf ("  Fail: waitpid call failed.\n");
281         ++count_errors;
282       }
283     else if (WIFSIGNALED (status) && WTERMSIG (status) == SIGFPE)
284       printf ("  Pass: Process received SIGFPE.\n");
285     else
286       {
287         printf ("  Fail: Process didn't receive signal and exited with status %d.\n",
288                 status);
289         ++count_errors;
290       }
291   }
292 #endif
293 }
294
295 /* Test that program doesn't abort with default environment */
296 static void
297 feenv_mask_test (const char *flag_name, int fe_exc)
298 {
299   int status;
300   pid_t pid;
301
302   printf ("Test: after fesetenv (FE_DFL_ENV) processes will not abort\n");
303   printf ("      when feraiseexcept (%s) is called.\n", flag_name);
304   pid = fork ();
305   if (pid == 0)
306     {
307 #ifdef RLIMIT_CORE
308       /* Try to avoid dumping core.  */
309       struct rlimit core_limit;
310       core_limit.rlim_cur = 0;
311       core_limit.rlim_max = 0;
312       setrlimit (RLIMIT_CORE, &core_limit);
313 #endif
314
315       fesetenv (FE_DFL_ENV);
316       feraiseexcept (fe_exc);
317       exit (2);
318     }
319   else if (pid < 0)
320     {
321       if (errno != ENOSYS)
322         {
323           printf ("  Fail: Could not fork.\n");
324           ++count_errors;
325         }
326       else
327         printf ("  `fork' not implemented, test ignored.\n");
328     }
329   else {
330     if (waitpid (pid, &status, 0) != pid)
331       {
332         printf ("  Fail: waitpid call failed.\n");
333         ++count_errors;
334       }
335     else if (WIFEXITED (status) && WEXITSTATUS (status) == 2)
336       printf ("  Pass: Process exited normally.\n");
337     else
338       {
339         printf ("  Fail: Process exited abnormally with status %d.\n",
340                 status);
341         ++count_errors;
342       }
343   }
344 }
345
346
347
348 static void
349 feenv_tests (void)
350 {
351
352 #ifdef FE_DIVBYZERO
353   feenv_nomask_test ("FE_DIVBYZERO", FE_DIVBYZERO);
354   feenv_mask_test ("FE_DIVBYZERO", FE_DIVBYZERO);
355 #endif
356 #ifdef FE_INVALID
357   feenv_nomask_test ("FE_INVALID", FE_INVALID);
358   feenv_mask_test ("FE_INVALID", FE_INVALID);
359 #endif
360 #ifdef FE_INEXACT
361   feenv_nomask_test ("FE_INEXACT", FE_INEXACT);
362   feenv_mask_test ("FE_INEXACT", FE_INEXACT);
363 #endif
364 #ifdef FE_UNDERFLOW
365   feenv_nomask_test ("FE_UNDERFLOW", FE_UNDERFLOW);
366   feenv_mask_test ("FE_UNDERFLOW", FE_UNDERFLOW);
367 #endif
368 #ifdef FE_OVERFLOW
369   feenv_nomask_test ("FE_OVERFLOW", FE_OVERFLOW);
370   feenv_mask_test ("FE_OVERFLOW", FE_OVERFLOW);
371 #endif
372   fesetenv (FE_DFL_ENV);
373 }
374
375
376 /* IEC 559 and ISO C 9X define a default startup environment */
377 static void
378 initial_tests (void)
379 {
380   test_exceptions ("Initially all exceptions should be cleared",
381                    NO_EXC, 0);
382   test_rounding ("Rounding direction should be initalized to nearest",
383                  FE_TONEAREST);
384 }
385
386 int
387 main (void)
388 {
389   initial_tests ();
390   fe_tests ();
391   feenv_tests ();
392
393   if (count_errors)
394     {
395       printf ("\n%d errors occured.\n", count_errors);
396       exit (1);
397     }
398   printf ("\n All tests passed successfully.\n");
399   exit (0);
400 }