Updated to fedora-glibc-20050523T1354
[kopensolaris-gnu/glibc.git] / math / basic-test.c
1 /* Copyright (C) 1999 Free Software Foundation, Inc.
2    This file is part of the GNU C Library.
3    Contributed by Andreas Jaeger <aj@suse.de>, 1999.
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, write to the Free
17    Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
18    02111-1307 USA.  */
19
20 #include <math.h>
21 #include <float.h>
22 #include <stdio.h>
23
24 static int errors = 0;
25
26
27 static void
28 check (const char *testname, int result)
29 {
30   if (!result) {
31     printf ("Failure: %s\n", testname);
32     errors++;
33   }
34 }
35
36 #define TEST_FUNC(NAME, FLOAT, NANFUNC, EPSILON, HUGEVAL) \
37 static void                                                                   \
38 NAME (void)                                                                   \
39 {                                                                             \
40   /* Variables are declared volatile to forbid some compiler                  \
41      optimizations.  */                                                       \
42   volatile FLOAT Inf_var, NaN_var, zero_var, one_var;                         \
43   FLOAT x1, x2;                                                               \
44                                                                               \
45   zero_var = 0.0;                                                             \
46   one_var = 1.0;                                                              \
47   NaN_var = zero_var/zero_var;                                                \
48   Inf_var = one_var / zero_var;                                               \
49                                                                               \
50   (void) &zero_var;                                                           \
51   (void) &one_var;                                                            \
52   (void) &NaN_var;                                                            \
53   (void) &Inf_var;                                                            \
54                                                                               \
55                                                                               \
56   check (#FLOAT " isinf (inf) == 1", isinf (Inf_var) == 1);                   \
57   check (#FLOAT " isinf (-inf) == -1", isinf (-Inf_var) == -1);               \
58   check (#FLOAT " !isinf (1)", !(isinf (one_var)));                           \
59   check (#FLOAT " !isinf (NaN)", !(isinf (NaN_var)));                         \
60                                                                               \
61   check (#FLOAT " isnan (NaN)", isnan (NaN_var));                             \
62   check (#FLOAT " isnan (-NaN)", isnan (-NaN_var));                           \
63   check (#FLOAT " !isnan (1)", !(isnan (one_var)));                           \
64   check (#FLOAT " !isnan (inf)", !(isnan (Inf_var)));                         \
65                                                                               \
66   check (#FLOAT " inf == inf", Inf_var == Inf_var);                           \
67   check (#FLOAT " -inf == -inf", -Inf_var == -Inf_var);                       \
68   check (#FLOAT " inf != -inf", Inf_var != -Inf_var);                         \
69   check (#FLOAT " NaN != NaN", NaN_var != NaN_var);                           \
70                                                                               \
71   /*                                                                          \
72      the same tests but this time with NAN from <bits/nan.h>                  \
73      NAN is a double const                                                    \
74    */                                                                         \
75   check (#FLOAT " isnan (NAN)", isnan (NAN));                                 \
76   check (#FLOAT " isnan (-NAN)", isnan (-NAN));                               \
77   check (#FLOAT " !isinf (NAN)", !(isinf (NAN)));                             \
78   check (#FLOAT " !isinf (-NAN)", !(isinf (-NAN)));                           \
79   check (#FLOAT " NAN != NAN", NAN != NAN);                                   \
80                                                                               \
81   /*                                                                          \
82      And again with the value returned by the `nan' function.                 \
83    */                                                                         \
84   check (#FLOAT " isnan (NAN)", isnan (NANFUNC ("")));                        \
85   check (#FLOAT " isnan (-NAN)", isnan (-NANFUNC ("")));                      \
86   check (#FLOAT " !isinf (NAN)", !(isinf (NANFUNC (""))));                    \
87   check (#FLOAT " !isinf (-NAN)", !(isinf (-NANFUNC (""))));                  \
88   check (#FLOAT " NAN != NAN", NANFUNC ("") != NANFUNC (""));                 \
89                                                                               \
90   /* test if EPSILON is ok */                                                 \
91   x1 = 1.0;                                                                   \
92   x2 = x1 + EPSILON;                                                          \
93   check (#FLOAT " 1 != 1+EPSILON", x1 != x2);                                 \
94                                                                               \
95   x1 = 1.0;                                                                   \
96   x2 = x1 - EPSILON;                                                          \
97   check (#FLOAT " 1 != 1-EPSILON", x1 != x2);                                 \
98                                                                               \
99   /* test if HUGE_VALx is ok */                                               \
100   x1 = HUGEVAL;                                                               \
101   check (#FLOAT " isinf (HUGE_VALx) == +1", isinf (x1) == +1);                \
102   x1 = - HUGEVAL;                                                             \
103   check (#FLOAT " isinf (-HUGE_VALx) == -1", isinf (x1) == -1);               \
104 }
105
106 TEST_FUNC (float_test, float, nanf, FLT_EPSILON, HUGE_VALF)
107 TEST_FUNC (double_test, double, nan, DBL_EPSILON, HUGE_VAL)
108 #ifndef NO_LONG_DOUBLE
109 TEST_FUNC (ldouble_test, long double, nanl, LDBL_EPSILON, HUGE_VALL)
110 #endif
111
112 int
113 main (void)
114 {
115   float_test ();
116   double_test ();
117
118 #ifndef NO_LONG_DOUBLE
119   ldouble_test ();
120 #endif
121
122   return errors != 0;
123 }