Test for module 'ceill'.
[pspp] / tests / test-ceilf.c
1 /* Test of rounding towards positive infinity.
2    Copyright (C) 2007 Free Software Foundation, Inc.
3
4    This program is free software; you can redistribute it and/or modify
5    it under the terms of the GNU General Public License as published by
6    the Free Software Foundation; either version 2, or (at your option)
7    any later version.
8
9    This program is distributed in the hope that it will be useful,
10    but WITHOUT ANY WARRANTY; without even the implied warranty of
11    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
12    GNU General Public License for more details.
13
14    You should have received a copy of the GNU General Public License
15    along with this program; if not, write to the Free Software Foundation,
16    Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.  */
17
18 /* Written by Bruno Haible <bruno@clisp.org>, 2007.  */
19
20 #include <config.h>
21
22 #include <math.h>
23
24 #include <stdio.h>
25 #include <stdlib.h>
26
27 #define ASSERT(expr) \
28   do                                                                         \
29     {                                                                        \
30       if (!(expr))                                                           \
31         {                                                                    \
32           fprintf (stderr, "%s:%d: assertion failed\n", __FILE__, __LINE__); \
33           abort ();                                                          \
34         }                                                                    \
35     }                                                                        \
36   while (0)
37
38 /* The Compaq (ex-DEC) C 6.4 compiler chokes on the expression 0.0 / 0.0.  */
39 #ifdef __DECC
40 static float
41 NaN ()
42 {
43   static float zero = 0.0f;
44   return zero / zero;
45 }
46 #else
47 # define NaN() (0.0f / 0.0f)
48 #endif
49
50 int
51 main ()
52 {
53   /* Zero.  */
54   ASSERT (ceilf (0.0f) == 0.0f);
55   ASSERT (ceilf (-0.0f) == 0.0f);
56   /* Positive numbers.  */
57   ASSERT (ceilf (0.3f) == 1.0f);
58   ASSERT (ceilf (0.7f) == 1.0f);
59   ASSERT (ceilf (1.0f) == 1.0f);
60   ASSERT (ceilf (1.001f) == 2.0f);
61   ASSERT (ceilf (1.5f) == 2.0f);
62   ASSERT (ceilf (1.999f) == 2.0f);
63   ASSERT (ceilf (2.0f) == 2.0f);
64   ASSERT (ceilf (65535.99f) == 65536.0f);
65   ASSERT (ceilf (65536.0f) == 65536.0f);
66   ASSERT (ceilf (2.341e31f) == 2.341e31f);
67   /* Negative numbers.  */
68   ASSERT (ceilf (-0.3f) == 0.0f);
69   ASSERT (ceilf (-0.7f) == 0.0f);
70   ASSERT (ceilf (-1.0f) == -1.0f);
71   ASSERT (ceilf (-1.5f) == -1.0f);
72   ASSERT (ceilf (-1.999f) == -1.0f);
73   ASSERT (ceilf (-2.0f) == -2.0f);
74   ASSERT (ceilf (-65535.99f) == -65535.0f);
75   ASSERT (ceilf (-65536.0f) == -65536.0f);
76   ASSERT (ceilf (-2.341e31f) == -2.341e31f);
77   /* Infinite numbers.  */
78   ASSERT (ceilf (1.0f / 0.0f) == 1.0f / 0.0f);
79   ASSERT (ceilf (-1.0f / 0.0f) == -1.0f / 0.0f);
80   /* NaNs.  */
81   ASSERT (isnan (ceilf (NaN ())));
82
83   return 0;
84 }