Google Git
Sign in
chromium / external / github.com / llvm / llvm-project.git / 412dbebe1b8adb9d6c3365e3b076139c675d0e03 / . / libcxx / test / numerics / complex.number / cases.h
blob: cd6bd165a9d160ab0555d61e59e8a2d84984acde [file] [log] [blame]
Howard Hinnant3e519522010-05-11 19:42:16[diff] [blame]1//===----------------------------------------------------------------------===//
2//
Howard Hinnant5b08a8a2010-05-11 21:36:01[diff] [blame]3// The LLVM Compiler Infrastructure
Howard Hinnant3e519522010-05-11 19:42:16[diff] [blame]4//
Howard Hinnant412dbeb2010-11-16 22:09:02[diff] [blame^]5// This file is dual licensed under the MIT and the University of Illinois Open
6// Source Licenses. See LICENSE.TXT for details.
Howard Hinnant3e519522010-05-11 19:42:16[diff] [blame]7//
8//===----------------------------------------------------------------------===//
9
10// <complex>
11
12// test cases
13
14#ifndef CASES_H
15#define CASES_H
16
17#include <complex>
18#include <cassert>
19
20std::complex<double> x[] =
21{
22 std::complex<double>( 1.e-6, 1.e-6),
23 std::complex<double>(-1.e-6, 1.e-6),
24 std::complex<double>(-1.e-6, -1.e-6),
25 std::complex<double>( 1.e-6, -1.e-6),
26
27 std::complex<double>( 1.e+6, 1.e-6),
28 std::complex<double>(-1.e+6, 1.e-6),
29 std::complex<double>(-1.e+6, -1.e-6),
30 std::complex<double>( 1.e+6, -1.e-6),
31
32 std::complex<double>( 1.e-6, 1.e+6),
33 std::complex<double>(-1.e-6, 1.e+6),
34 std::complex<double>(-1.e-6, -1.e+6),
35 std::complex<double>( 1.e-6, -1.e+6),
36
37 std::complex<double>( 1.e+6, 1.e+6),
38 std::complex<double>(-1.e+6, 1.e+6),
39 std::complex<double>(-1.e+6, -1.e+6),
40 std::complex<double>( 1.e+6, -1.e+6),
41
42 std::complex<double>(NAN, NAN),
43 std::complex<double>(-INFINITY, NAN),
44 std::complex<double>(-2, NAN),
45 std::complex<double>(-1, NAN),
46 std::complex<double>(-0.5, NAN),
47 std::complex<double>(-0., NAN),
48 std::complex<double>(+0., NAN),
49 std::complex<double>(0.5, NAN),
50 std::complex<double>(1, NAN),
51 std::complex<double>(2, NAN),
52 std::complex<double>(INFINITY, NAN),
53
54 std::complex<double>(NAN, -INFINITY),
55 std::complex<double>(-INFINITY, -INFINITY),
56 std::complex<double>(-2, -INFINITY),
57 std::complex<double>(-1, -INFINITY),
58 std::complex<double>(-0.5, -INFINITY),
59 std::complex<double>(-0., -INFINITY),
60 std::complex<double>(+0., -INFINITY),
61 std::complex<double>(0.5, -INFINITY),
62 std::complex<double>(1, -INFINITY),
63 std::complex<double>(2, -INFINITY),
64 std::complex<double>(INFINITY, -INFINITY),
65
66 std::complex<double>(NAN, -2),
67 std::complex<double>(-INFINITY, -2),
68 std::complex<double>(-2, -2),
69 std::complex<double>(-1, -2),
70 std::complex<double>(-0.5, -2),
71 std::complex<double>(-0., -2),
72 std::complex<double>(+0., -2),
73 std::complex<double>(0.5, -2),
74 std::complex<double>(1, -2),
75 std::complex<double>(2, -2),
76 std::complex<double>(INFINITY, -2),
77
78 std::complex<double>(NAN, -1),
79 std::complex<double>(-INFINITY, -1),
80 std::complex<double>(-2, -1),
81 std::complex<double>(-1, -1),
82 std::complex<double>(-0.5, -1),
83 std::complex<double>(-0., -1),
84 std::complex<double>(+0., -1),
85 std::complex<double>(0.5, -1),
86 std::complex<double>(1, -1),
87 std::complex<double>(2, -1),
88 std::complex<double>(INFINITY, -1),
89
90 std::complex<double>(NAN, -0.5),
91 std::complex<double>(-INFINITY, -0.5),
92 std::complex<double>(-2, -0.5),
93 std::complex<double>(-1, -0.5),
94 std::complex<double>(-0.5, -0.5),
95 std::complex<double>(-0., -0.5),
96 std::complex<double>(+0., -0.5),
97 std::complex<double>(0.5, -0.5),
98 std::complex<double>(1, -0.5),
99 std::complex<double>(2, -0.5),
100 std::complex<double>(INFINITY, -0.5),
101
102 std::complex<double>(NAN, -0.),
103 std::complex<double>(-INFINITY, -0.),
104 std::complex<double>(-2, -0.),
105 std::complex<double>(-1, -0.),
106 std::complex<double>(-0.5, -0.),
107 std::complex<double>(-0., -0.),
108 std::complex<double>(+0., -0.),
109 std::complex<double>(0.5, -0.),
110 std::complex<double>(1, -0.),
111 std::complex<double>(2, -0.),
112 std::complex<double>(INFINITY, -0.),
113
114 std::complex<double>(NAN, +0.),
115 std::complex<double>(-INFINITY, +0.),
116 std::complex<double>(-2, +0.),
117 std::complex<double>(-1, +0.),
118 std::complex<double>(-0.5, +0.),
119 std::complex<double>(-0., +0.),
120 std::complex<double>(+0., +0.),
121 std::complex<double>(0.5, +0.),
122 std::complex<double>(1, +0.),
123 std::complex<double>(2, +0.),
124 std::complex<double>(INFINITY, +0.),
125
126 std::complex<double>(NAN, 0.5),
127 std::complex<double>(-INFINITY, 0.5),
128 std::complex<double>(-2, 0.5),
129 std::complex<double>(-1, 0.5),
130 std::complex<double>(-0.5, 0.5),
131 std::complex<double>(-0., 0.5),
132 std::complex<double>(+0., 0.5),
133 std::complex<double>(0.5, 0.5),
134 std::complex<double>(1, 0.5),
135 std::complex<double>(2, 0.5),
136 std::complex<double>(INFINITY, 0.5),
137
138 std::complex<double>(NAN, 1),
139 std::complex<double>(-INFINITY, 1),
140 std::complex<double>(-2, 1),
141 std::complex<double>(-1, 1),
142 std::complex<double>(-0.5, 1),
143 std::complex<double>(-0., 1),
144 std::complex<double>(+0., 1),
145 std::complex<double>(0.5, 1),
146 std::complex<double>(1, 1),
147 std::complex<double>(2, 1),
148 std::complex<double>(INFINITY, 1),
149
150 std::complex<double>(NAN, 2),
151 std::complex<double>(-INFINITY, 2),
152 std::complex<double>(-2, 2),
153 std::complex<double>(-1, 2),
154 std::complex<double>(-0.5, 2),
155 std::complex<double>(-0., 2),
156 std::complex<double>(+0., 2),
157 std::complex<double>(0.5, 2),
158 std::complex<double>(1, 2),
159 std::complex<double>(2, 2),
160 std::complex<double>(INFINITY, 2),
161
162 std::complex<double>(NAN, INFINITY),
163 std::complex<double>(-INFINITY, INFINITY),
164 std::complex<double>(-2, INFINITY),
165 std::complex<double>(-1, INFINITY),
166 std::complex<double>(-0.5, INFINITY),
167 std::complex<double>(-0., INFINITY),
168 std::complex<double>(+0., INFINITY),
169 std::complex<double>(0.5, INFINITY),
170 std::complex<double>(1, INFINITY),
171 std::complex<double>(2, INFINITY),
172 std::complex<double>(INFINITY, INFINITY)
173};
174
175enum {zero, non_zero, inf, NaN, non_zero_nan};
176
177template <class T>
178int
179classify(const std::complex<T>& x)
180{
181 if (x == std::complex<T>())
182 return zero;
183 if (std::isinf(x.real()) || std::isinf(x.imag()))
184 return inf;
185 if (std::isnan(x.real()) && std::isnan(x.imag()))
186 return NaN;
187 if (std::isnan(x.real()))
188 {
189 if (x.imag() == T(0))
190 return NaN;
191 return non_zero_nan;
192 }
193 if (std::isnan(x.imag()))
194 {
195 if (x.real() == T(0))
196 return NaN;
197 return non_zero_nan;
198 }
199 return non_zero;
200}
201
202inline
203int
204classify(double x)
205{
206 if (x == 0)
207 return zero;
208 if (std::isinf(x))
209 return inf;
210 if (std::isnan(x))
211 return NaN;
212 return non_zero;
213}
214
215void is_about(float x, float y)
216{
217 assert(std::abs((x-y)/(x+y)) < 1.e-6);
218}
219
220void is_about(double x, double y)
221{
222 assert(std::abs((x-y)/(x+y)) < 1.e-14);
223}
224
225void is_about(long double x, long double y)
226{
227 assert(std::abs((x-y)/(x+y)) < 1.e-14);
228}
229
Howard Hinnanta0eaf602010-08-22 00:42:21[diff] [blame]230#endif // CASES_H