try to preserve x bits
[pspp] / tests / test-trunc1.c
1 /* Test of rounding towards zero.
2    Copyright (C) 2007-2008 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 3 of the License, or
7    (at your option) 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, see <http://www.gnu.org/licenses/>.  */
16
17 /* Written by Bruno Haible <bruno@clisp.org>, 2007.  */
18
19 #include <config.h>
20
21 #include <math.h>
22
23 #include <stdio.h>
24 #include <stdlib.h>
25
26 #include "isnand.h"
27
28 #define ASSERT(expr) \
29   do                                                                         \
30     {                                                                        \
31       if (!(expr))                                                           \
32         {                                                                    \
33           fprintf (stderr, "%s:%d: assertion failed\n", __FILE__, __LINE__); \
34           abort ();                                                          \
35         }                                                                    \
36     }                                                                        \
37   while (0)
38
39 /* The Compaq (ex-DEC) C 6.4 compiler chokes on the expression 0.0 / 0.0.  */
40 #ifdef __DECC
41 static double
42 NaN ()
43 {
44   static double zero = 0.0;
45   return zero / zero;
46 }
47 #else
48 # define NaN() (0.0 / 0.0)
49 #endif
50
51 int
52 main ()
53 {
54   /* Zero.  */
55   ASSERT (trunc (0.0) == 0.0);
56   ASSERT (trunc (-0.0) == 0.0);
57   /* Positive numbers.  */
58   ASSERT (trunc (0.3) == 0.0);
59   ASSERT (trunc (0.7) == 0.0);
60   ASSERT (trunc (1.0) == 1.0);
61   ASSERT (trunc (1.5) == 1.0);
62   ASSERT (trunc (1.999) == 1.0);
63   ASSERT (trunc (2.0) == 2.0);
64   ASSERT (trunc (65535.999) == 65535.0);
65   ASSERT (trunc (65536.0) == 65536.0);
66   ASSERT (trunc (2.341e31) == 2.341e31);
67   /* Negative numbers.  */
68   ASSERT (trunc (-0.3) == 0.0);
69   ASSERT (trunc (-0.7) == 0.0);
70   ASSERT (trunc (-1.0) == -1.0);
71   ASSERT (trunc (-1.5) == -1.0);
72   ASSERT (trunc (-1.999) == -1.0);
73   ASSERT (trunc (-2.0) == -2.0);
74   ASSERT (trunc (-65535.999) == -65535.0);
75   ASSERT (trunc (-65536.0) == -65536.0);
76   ASSERT (trunc (-2.341e31) == -2.341e31);
77   /* Infinite numbers.  */
78   ASSERT (trunc (1.0 / 0.0) == 1.0 / 0.0);
79   ASSERT (trunc (-1.0 / 0.0) == -1.0 / 0.0);
80   /* NaNs.  */
81   ASSERT (isnand (trunc (NaN ())));
82
83   return 0;
84 }