aboutsummaryrefslogtreecommitdiff
path: root/test/test-fp_Infcases.c
blob: 88a89eadb01e8160cfc080f2458e00074f0ff5ee (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
/* test cases for floating point comparison transformations
 * compile with -DFLOAT_TYPE=float
 *          or  -DFLOAT_TYPE=double
 *          or  -DFLOAT_TYPE="long double"
 */

#include <assert.h>
#define _GNU_SOURCE
#include <math.h> /* for NaNs and infinity values */

int main() {

  volatile FLOAT_TYPE a, b;

#ifdef INFINITY
  FLOAT_TYPE inf = (FLOAT_TYPE) INFINITY;
#else
  FLOAT_TYPE inf = 1.0 / 0.0; /* produces infinity */
#endif
  FLOAT_TYPE negZero = 1.0 / -inf;
  FLOAT_TYPE posZero = 0.0;

  /* plus infinity */
  a = (1.0 / 0.0);	/* positive infinity */
  b = (1.0 / 0.0);	/* positive infinity */
  assert(!(a < b));
  assert((a <= b));
  assert(!(a > b));
  assert((a >= b));
  assert(!(a != b));
  assert((a == b));

  b = -(1.0 / 0.0);	/* negative infinity */
  assert(!(a < b));
  assert(!(a <= b));
  assert((a > b));
  assert((a >= b));
  assert((a != b));
  assert(!(a == b));

  b = 1.0 / -(1.0 / 0.0);	/* negative 0 */
  assert(!(a < b));
  assert(!(a <= b));
  assert((a > b));
  assert((a >= b));
  assert((a != b));
  assert(!(a == b));

  b = 0.0;	/* positive 0 */
  assert(!(a < b));
  assert(!(a <= b));
  assert((a > b));
  assert((a >= b));
  assert((a != b));
  assert(!(a == b));

  b = -42.0;
  assert(!(a < b));
  assert(!(a <= b));
  assert((a > b));
  assert((a >= b));
  assert((a != b));
  assert(!(a == b));

  b = 42.0;
  assert(!(a < b));
  assert(!(a <= b));
  assert((a > b));
  assert((a >= b));
  assert((a != b));
  assert(!(a == b));

  /* negative infinity */
  a = -(1.0 / 0.0);	
  b = (1.0 / 0.0);	/* positive infinity */
  assert((a < b));
  assert((a <= b));
  assert(!(a > b));
  assert(!(a >= b));
  assert((a != b));
  assert(!(a == b));

  b = -(1.0 / 0.0);	/* negative infinity */
  assert(!(a < b));
  assert((a <= b));
  assert(!(a > b));
  assert((a >= b));
  assert(!(a != b));
  assert((a == b));

  b = 1.0 / -(1.0 / 0.0);	/* negative 0 */
  assert((a < b));
  assert((a <= b));
  assert(!(a > b));
  assert(!(a >= b));
  assert((a != b));
  assert(!(a == b));

  b = 0.0;	/* positive 0 */
  assert((a < b));
  assert((a <= b));
  assert(!(a > b));
  assert(!(a >= b));
  assert((a != b));
  assert(!(a == b));

  b = -42.0;
  assert((a < b));
  assert((a <= b));
  assert(!(a > b));
  assert(!(a >= b));
  assert((a != b));
  assert(!(a == b));

  b = 42.0;
  assert((a < b));
  assert((a <= b));
  assert(!(a > b));
  assert(!(a >= b));
  assert((a != b));
  assert(!(a == b));

}