blob: 9bcd63dd2302cb84ef703cc28bbd51c09824072c [file] [log] [blame]
Michael W. Hudsonba283e22005-05-27 15:23:20 +00001
2import unittest, struct
Christian Heimes827b35c2007-12-10 22:19:17 +00003import os
Eric Smith0923d1d2009-04-16 20:16:10 +00004import sys
Benjamin Petersonee8712c2008-05-20 21:35:26 +00005from test import support
Christian Heimes53876d92008-04-19 00:31:39 +00006import math
Mark Dickinson65fe25e2008-07-16 11:30:51 +00007from math import isinf, isnan, copysign, ldexp
Christian Heimes53876d92008-04-19 00:31:39 +00008import operator
Amaury Forgeot d'Arc7e958d12008-09-06 21:03:22 +00009import random, fractions
Michael W. Hudsonba283e22005-05-27 15:23:20 +000010
Christian Heimes53876d92008-04-19 00:31:39 +000011INF = float("inf")
12NAN = float("nan")
Christian Heimes99170a52007-12-19 02:07:34 +000013
Benjamin Petersone401c682010-07-02 23:25:44 +000014have_getformat = hasattr(float, "__getformat__")
15requires_getformat = unittest.skipUnless(have_getformat,
16 "requires __getformat__")
17requires_setformat = unittest.skipUnless(hasattr(float, "__setformat__"),
18 "requires __setformat__")
Mark Dickinson9ab44b52009-12-30 16:22:49 +000019# decorator for skipping tests on non-IEEE 754 platforms
Benjamin Petersone401c682010-07-02 23:25:44 +000020requires_IEEE_754 = unittest.skipUnless(have_getformat and
Mark Dickinson9ab44b52009-12-30 16:22:49 +000021 float.__getformat__("double").startswith("IEEE"),
22 "test requires IEEE 754 doubles")
23
Eric Smith0923d1d2009-04-16 20:16:10 +000024#locate file with float format test values
25test_dir = os.path.dirname(__file__) or os.curdir
26format_testfile = os.path.join(test_dir, 'formatfloat_testcases.txt')
27
Christian Heimes81ee3ef2008-05-04 22:42:01 +000028class GeneralFloatCases(unittest.TestCase):
29
30 def test_float(self):
31 self.assertEqual(float(3.14), 3.14)
32 self.assertEqual(float(314), 314.0)
33 self.assertEqual(float(" 3.14 "), 3.14)
Amaury Forgeot d'Arc7e958d12008-09-06 21:03:22 +000034 self.assertEqual(float(b" 3.14 "), 3.14)
Christian Heimes81ee3ef2008-05-04 22:42:01 +000035 self.assertRaises(ValueError, float, " 0x3.1 ")
36 self.assertRaises(ValueError, float, " -0x3.p-1 ")
37 self.assertRaises(ValueError, float, " +0x3.p-1 ")
38 self.assertRaises(ValueError, float, "++3.14")
39 self.assertRaises(ValueError, float, "+-3.14")
40 self.assertRaises(ValueError, float, "-+3.14")
41 self.assertRaises(ValueError, float, "--3.14")
Eric Smith0923d1d2009-04-16 20:16:10 +000042 self.assertRaises(ValueError, float, ".nan")
43 self.assertRaises(ValueError, float, "+.inf")
44 self.assertRaises(ValueError, float, ".")
45 self.assertRaises(ValueError, float, "-.")
Alexander Belopolsky942af5a2010-12-04 03:38:46 +000046 self.assertRaises(ValueError, float, b"-")
47 self.assertRaises(TypeError, float, {})
48 # Lone surrogate
49 self.assertRaises(UnicodeEncodeError, float, '\uD8F0')
Mark Dickinsonc2d86892010-02-12 21:18:34 +000050 # check that we don't accept alternate exponent markers
51 self.assertRaises(ValueError, float, "-1.7d29")
52 self.assertRaises(ValueError, float, "3D-14")
Alexander Belopolsky942af5a2010-12-04 03:38:46 +000053 self.assertEqual(float(" \u0663.\u0661\u0664 "), 3.14)
54 self.assertEqual(float("\N{EM SPACE}3.14\N{EN SPACE}"), 3.14)
Mark Dickinsona9023be2009-10-27 22:12:20 +000055 # extra long strings should not be a problem
56 float(b'.' + b'1'*1000)
57 float('.' + '1'*1000)
Christian Heimes81ee3ef2008-05-04 22:42:01 +000058
Alexander Belopolsky942af5a2010-12-04 03:38:46 +000059 def test_error_message(self):
60 testlist = ('\xbd', '123\xbd', ' 123 456 ')
61 for s in testlist:
62 try:
63 float(s)
64 except ValueError as e:
65 self.assertIn(s.strip(), e.args[0])
66 else:
67 self.fail("Expected int(%r) to raise a ValueError", s)
68
69
Benjamin Petersonee8712c2008-05-20 21:35:26 +000070 @support.run_with_locale('LC_NUMERIC', 'fr_FR', 'de_DE')
Christian Heimes81ee3ef2008-05-04 22:42:01 +000071 def test_float_with_comma(self):
72 # set locale to something that doesn't use '.' for the decimal point
73 # float must not accept the locale specific decimal point but
74 # it still has to accept the normal python syntac
75 import locale
76 if not locale.localeconv()['decimal_point'] == ',':
77 return
78
79 self.assertEqual(float(" 3.14 "), 3.14)
80 self.assertEqual(float("+3.14 "), 3.14)
81 self.assertEqual(float("-3.14 "), -3.14)
82 self.assertEqual(float(".14 "), .14)
83 self.assertEqual(float("3. "), 3.0)
84 self.assertEqual(float("3.e3 "), 3000.0)
85 self.assertEqual(float("3.2e3 "), 3200.0)
86 self.assertEqual(float("2.5e-1 "), 0.25)
87 self.assertEqual(float("5e-1"), 0.5)
88 self.assertRaises(ValueError, float, " 3,14 ")
89 self.assertRaises(ValueError, float, " +3,14 ")
90 self.assertRaises(ValueError, float, " -3,14 ")
91 self.assertRaises(ValueError, float, " 0x3.1 ")
92 self.assertRaises(ValueError, float, " -0x3.p-1 ")
93 self.assertRaises(ValueError, float, " +0x3.p-1 ")
94 self.assertEqual(float(" 25.e-1 "), 2.5)
Benjamin Petersond43029b2008-09-06 23:33:21 +000095 self.assertEqual(support.fcmp(float(" .25e-1 "), .025), 0)
Christian Heimes81ee3ef2008-05-04 22:42:01 +000096
97 def test_floatconversion(self):
98 # Make sure that calls to __float__() work properly
99 class Foo0:
100 def __float__(self):
101 return 42.
102
103 class Foo1(object):
104 def __float__(self):
105 return 42.
106
107 class Foo2(float):
108 def __float__(self):
109 return 42.
110
111 class Foo3(float):
112 def __new__(cls, value=0.):
113 return float.__new__(cls, 2*value)
114
115 def __float__(self):
116 return self
117
118 class Foo4(float):
119 def __float__(self):
120 return 42
121
Benjamin Peterson2808d3c2009-04-15 21:34:27 +0000122 # Issue 5759: __float__ not called on str subclasses (though it is on
123 # unicode subclasses).
124 class FooStr(str):
125 def __float__(self):
126 return float(str(self)) + 1
127
Christian Heimes81ee3ef2008-05-04 22:42:01 +0000128 self.assertAlmostEqual(float(Foo0()), 42.)
129 self.assertAlmostEqual(float(Foo1()), 42.)
130 self.assertAlmostEqual(float(Foo2()), 42.)
131 self.assertAlmostEqual(float(Foo3(21)), 42.)
132 self.assertRaises(TypeError, float, Foo4(42))
Benjamin Peterson2808d3c2009-04-15 21:34:27 +0000133 self.assertAlmostEqual(float(FooStr('8')), 9.)
Christian Heimes81ee3ef2008-05-04 22:42:01 +0000134
135 def test_floatasratio(self):
136 for f, ratio in [
137 (0.875, (7, 8)),
138 (-0.875, (-7, 8)),
139 (0.0, (0, 1)),
140 (11.5, (23, 2)),
141 ]:
142 self.assertEqual(f.as_integer_ratio(), ratio)
143
144 for i in range(10000):
145 f = random.random()
146 f *= 10 ** random.randint(-100, 100)
147 n, d = f.as_integer_ratio()
148 self.assertEqual(float(n).__truediv__(d), f)
149
150 R = fractions.Fraction
151 self.assertEqual(R(0, 1),
152 R(*float(0.0).as_integer_ratio()))
153 self.assertEqual(R(5, 2),
154 R(*float(2.5).as_integer_ratio()))
155 self.assertEqual(R(1, 2),
156 R(*float(0.5).as_integer_ratio()))
157 self.assertEqual(R(4728779608739021, 2251799813685248),
158 R(*float(2.1).as_integer_ratio()))
159 self.assertEqual(R(-4728779608739021, 2251799813685248),
160 R(*float(-2.1).as_integer_ratio()))
161 self.assertEqual(R(-2100, 1),
162 R(*float(-2100.0).as_integer_ratio()))
163
164 self.assertRaises(OverflowError, float('inf').as_integer_ratio)
165 self.assertRaises(OverflowError, float('-inf').as_integer_ratio)
166 self.assertRaises(ValueError, float('nan').as_integer_ratio)
167
Mark Dickinson4a1f5932008-11-12 23:23:36 +0000168 def test_float_containment(self):
169 floats = (INF, -INF, 0.0, 1.0, NAN)
170 for f in floats:
Benjamin Peterson577473f2010-01-19 00:09:57 +0000171 self.assertIn(f, [f])
Benjamin Peterson577473f2010-01-19 00:09:57 +0000172 self.assertIn(f, (f,))
Benjamin Peterson577473f2010-01-19 00:09:57 +0000173 self.assertIn(f, {f})
Benjamin Peterson577473f2010-01-19 00:09:57 +0000174 self.assertIn(f, {f: None})
Mark Dickinson4a1f5932008-11-12 23:23:36 +0000175 self.assertEqual([f].count(f), 1, "[].count('%r') != 1" % f)
Benjamin Peterson577473f2010-01-19 00:09:57 +0000176 self.assertIn(f, floats)
Mark Dickinson4a1f5932008-11-12 23:23:36 +0000177
178 for f in floats:
179 # nonidentical containers, same type, same contents
Benjamin Petersonc9c0f202009-06-30 23:06:06 +0000180 self.assertTrue([f] == [f], "[%r] != [%r]" % (f, f))
181 self.assertTrue((f,) == (f,), "(%r,) != (%r,)" % (f, f))
182 self.assertTrue({f} == {f}, "{%r} != {%r}" % (f, f))
183 self.assertTrue({f : None} == {f: None}, "{%r : None} != "
Mark Dickinson4a1f5932008-11-12 23:23:36 +0000184 "{%r : None}" % (f, f))
185
186 # identical containers
187 l, t, s, d = [f], (f,), {f}, {f: None}
Benjamin Petersonc9c0f202009-06-30 23:06:06 +0000188 self.assertTrue(l == l, "[%r] not equal to itself" % f)
189 self.assertTrue(t == t, "(%r,) not equal to itself" % f)
190 self.assertTrue(s == s, "{%r} not equal to itself" % f)
191 self.assertTrue(d == d, "{%r : None} not equal to itself" % f)
Mark Dickinson4a1f5932008-11-12 23:23:36 +0000192
Mark Dickinson9ab44b52009-12-30 16:22:49 +0000193 def assertEqualAndEqualSign(self, a, b):
194 # fail unless a == b and a and b have the same sign bit;
195 # the only difference from assertEqual is that this test
196 # distingishes -0.0 and 0.0.
197 self.assertEqual((a, copysign(1.0, a)), (b, copysign(1.0, b)))
198
199 @requires_IEEE_754
200 def test_float_pow(self):
201 # test builtin pow and ** operator for IEEE 754 special cases.
202 # Special cases taken from section F.9.4.4 of the C99 specification
203
204 for pow_op in pow, operator.pow:
205 # x**NAN is NAN for any x except 1
206 self.assertTrue(isnan(pow_op(-INF, NAN)))
207 self.assertTrue(isnan(pow_op(-2.0, NAN)))
208 self.assertTrue(isnan(pow_op(-1.0, NAN)))
209 self.assertTrue(isnan(pow_op(-0.5, NAN)))
210 self.assertTrue(isnan(pow_op(-0.0, NAN)))
211 self.assertTrue(isnan(pow_op(0.0, NAN)))
212 self.assertTrue(isnan(pow_op(0.5, NAN)))
213 self.assertTrue(isnan(pow_op(2.0, NAN)))
214 self.assertTrue(isnan(pow_op(INF, NAN)))
215 self.assertTrue(isnan(pow_op(NAN, NAN)))
216
217 # NAN**y is NAN for any y except +-0
218 self.assertTrue(isnan(pow_op(NAN, -INF)))
219 self.assertTrue(isnan(pow_op(NAN, -2.0)))
220 self.assertTrue(isnan(pow_op(NAN, -1.0)))
221 self.assertTrue(isnan(pow_op(NAN, -0.5)))
222 self.assertTrue(isnan(pow_op(NAN, 0.5)))
223 self.assertTrue(isnan(pow_op(NAN, 1.0)))
224 self.assertTrue(isnan(pow_op(NAN, 2.0)))
225 self.assertTrue(isnan(pow_op(NAN, INF)))
226
227 # (+-0)**y raises ZeroDivisionError for y a negative odd integer
228 self.assertRaises(ZeroDivisionError, pow_op, -0.0, -1.0)
229 self.assertRaises(ZeroDivisionError, pow_op, 0.0, -1.0)
230
231 # (+-0)**y raises ZeroDivisionError for y finite and negative
232 # but not an odd integer
233 self.assertRaises(ZeroDivisionError, pow_op, -0.0, -2.0)
234 self.assertRaises(ZeroDivisionError, pow_op, -0.0, -0.5)
235 self.assertRaises(ZeroDivisionError, pow_op, 0.0, -2.0)
236 self.assertRaises(ZeroDivisionError, pow_op, 0.0, -0.5)
237
238 # (+-0)**y is +-0 for y a positive odd integer
239 self.assertEqualAndEqualSign(pow_op(-0.0, 1.0), -0.0)
240 self.assertEqualAndEqualSign(pow_op(0.0, 1.0), 0.0)
241
242 # (+-0)**y is 0 for y finite and positive but not an odd integer
243 self.assertEqualAndEqualSign(pow_op(-0.0, 0.5), 0.0)
244 self.assertEqualAndEqualSign(pow_op(-0.0, 2.0), 0.0)
245 self.assertEqualAndEqualSign(pow_op(0.0, 0.5), 0.0)
246 self.assertEqualAndEqualSign(pow_op(0.0, 2.0), 0.0)
247
248 # (-1)**+-inf is 1
249 self.assertEqualAndEqualSign(pow_op(-1.0, -INF), 1.0)
250 self.assertEqualAndEqualSign(pow_op(-1.0, INF), 1.0)
251
252 # 1**y is 1 for any y, even if y is an infinity or nan
253 self.assertEqualAndEqualSign(pow_op(1.0, -INF), 1.0)
254 self.assertEqualAndEqualSign(pow_op(1.0, -2.0), 1.0)
255 self.assertEqualAndEqualSign(pow_op(1.0, -1.0), 1.0)
256 self.assertEqualAndEqualSign(pow_op(1.0, -0.5), 1.0)
257 self.assertEqualAndEqualSign(pow_op(1.0, -0.0), 1.0)
258 self.assertEqualAndEqualSign(pow_op(1.0, 0.0), 1.0)
259 self.assertEqualAndEqualSign(pow_op(1.0, 0.5), 1.0)
260 self.assertEqualAndEqualSign(pow_op(1.0, 1.0), 1.0)
261 self.assertEqualAndEqualSign(pow_op(1.0, 2.0), 1.0)
262 self.assertEqualAndEqualSign(pow_op(1.0, INF), 1.0)
263 self.assertEqualAndEqualSign(pow_op(1.0, NAN), 1.0)
264
265 # x**+-0 is 1 for any x, even if x is a zero, infinity, or nan
266 self.assertEqualAndEqualSign(pow_op(-INF, 0.0), 1.0)
267 self.assertEqualAndEqualSign(pow_op(-2.0, 0.0), 1.0)
268 self.assertEqualAndEqualSign(pow_op(-1.0, 0.0), 1.0)
269 self.assertEqualAndEqualSign(pow_op(-0.5, 0.0), 1.0)
270 self.assertEqualAndEqualSign(pow_op(-0.0, 0.0), 1.0)
271 self.assertEqualAndEqualSign(pow_op(0.0, 0.0), 1.0)
272 self.assertEqualAndEqualSign(pow_op(0.5, 0.0), 1.0)
273 self.assertEqualAndEqualSign(pow_op(1.0, 0.0), 1.0)
274 self.assertEqualAndEqualSign(pow_op(2.0, 0.0), 1.0)
275 self.assertEqualAndEqualSign(pow_op(INF, 0.0), 1.0)
276 self.assertEqualAndEqualSign(pow_op(NAN, 0.0), 1.0)
277 self.assertEqualAndEqualSign(pow_op(-INF, -0.0), 1.0)
278 self.assertEqualAndEqualSign(pow_op(-2.0, -0.0), 1.0)
279 self.assertEqualAndEqualSign(pow_op(-1.0, -0.0), 1.0)
280 self.assertEqualAndEqualSign(pow_op(-0.5, -0.0), 1.0)
281 self.assertEqualAndEqualSign(pow_op(-0.0, -0.0), 1.0)
282 self.assertEqualAndEqualSign(pow_op(0.0, -0.0), 1.0)
283 self.assertEqualAndEqualSign(pow_op(0.5, -0.0), 1.0)
284 self.assertEqualAndEqualSign(pow_op(1.0, -0.0), 1.0)
285 self.assertEqualAndEqualSign(pow_op(2.0, -0.0), 1.0)
286 self.assertEqualAndEqualSign(pow_op(INF, -0.0), 1.0)
287 self.assertEqualAndEqualSign(pow_op(NAN, -0.0), 1.0)
288
289 # x**y defers to complex pow for finite negative x and
290 # non-integral y.
291 self.assertEqual(type(pow_op(-2.0, -0.5)), complex)
292 self.assertEqual(type(pow_op(-2.0, 0.5)), complex)
293 self.assertEqual(type(pow_op(-1.0, -0.5)), complex)
294 self.assertEqual(type(pow_op(-1.0, 0.5)), complex)
295 self.assertEqual(type(pow_op(-0.5, -0.5)), complex)
296 self.assertEqual(type(pow_op(-0.5, 0.5)), complex)
297
298 # x**-INF is INF for abs(x) < 1
299 self.assertEqualAndEqualSign(pow_op(-0.5, -INF), INF)
300 self.assertEqualAndEqualSign(pow_op(-0.0, -INF), INF)
301 self.assertEqualAndEqualSign(pow_op(0.0, -INF), INF)
302 self.assertEqualAndEqualSign(pow_op(0.5, -INF), INF)
303
304 # x**-INF is 0 for abs(x) > 1
305 self.assertEqualAndEqualSign(pow_op(-INF, -INF), 0.0)
306 self.assertEqualAndEqualSign(pow_op(-2.0, -INF), 0.0)
307 self.assertEqualAndEqualSign(pow_op(2.0, -INF), 0.0)
308 self.assertEqualAndEqualSign(pow_op(INF, -INF), 0.0)
309
310 # x**INF is 0 for abs(x) < 1
311 self.assertEqualAndEqualSign(pow_op(-0.5, INF), 0.0)
312 self.assertEqualAndEqualSign(pow_op(-0.0, INF), 0.0)
313 self.assertEqualAndEqualSign(pow_op(0.0, INF), 0.0)
314 self.assertEqualAndEqualSign(pow_op(0.5, INF), 0.0)
315
316 # x**INF is INF for abs(x) > 1
317 self.assertEqualAndEqualSign(pow_op(-INF, INF), INF)
318 self.assertEqualAndEqualSign(pow_op(-2.0, INF), INF)
319 self.assertEqualAndEqualSign(pow_op(2.0, INF), INF)
320 self.assertEqualAndEqualSign(pow_op(INF, INF), INF)
321
322 # (-INF)**y is -0.0 for y a negative odd integer
323 self.assertEqualAndEqualSign(pow_op(-INF, -1.0), -0.0)
324
325 # (-INF)**y is 0.0 for y negative but not an odd integer
326 self.assertEqualAndEqualSign(pow_op(-INF, -0.5), 0.0)
327 self.assertEqualAndEqualSign(pow_op(-INF, -2.0), 0.0)
328
329 # (-INF)**y is -INF for y a positive odd integer
330 self.assertEqualAndEqualSign(pow_op(-INF, 1.0), -INF)
331
332 # (-INF)**y is INF for y positive but not an odd integer
333 self.assertEqualAndEqualSign(pow_op(-INF, 0.5), INF)
334 self.assertEqualAndEqualSign(pow_op(-INF, 2.0), INF)
335
336 # INF**y is INF for y positive
337 self.assertEqualAndEqualSign(pow_op(INF, 0.5), INF)
338 self.assertEqualAndEqualSign(pow_op(INF, 1.0), INF)
339 self.assertEqualAndEqualSign(pow_op(INF, 2.0), INF)
340
341 # INF**y is 0.0 for y negative
342 self.assertEqualAndEqualSign(pow_op(INF, -2.0), 0.0)
343 self.assertEqualAndEqualSign(pow_op(INF, -1.0), 0.0)
344 self.assertEqualAndEqualSign(pow_op(INF, -0.5), 0.0)
345
346 # basic checks not covered by the special cases above
347 self.assertEqualAndEqualSign(pow_op(-2.0, -2.0), 0.25)
348 self.assertEqualAndEqualSign(pow_op(-2.0, -1.0), -0.5)
349 self.assertEqualAndEqualSign(pow_op(-2.0, -0.0), 1.0)
350 self.assertEqualAndEqualSign(pow_op(-2.0, 0.0), 1.0)
351 self.assertEqualAndEqualSign(pow_op(-2.0, 1.0), -2.0)
352 self.assertEqualAndEqualSign(pow_op(-2.0, 2.0), 4.0)
353 self.assertEqualAndEqualSign(pow_op(-1.0, -2.0), 1.0)
354 self.assertEqualAndEqualSign(pow_op(-1.0, -1.0), -1.0)
355 self.assertEqualAndEqualSign(pow_op(-1.0, -0.0), 1.0)
356 self.assertEqualAndEqualSign(pow_op(-1.0, 0.0), 1.0)
357 self.assertEqualAndEqualSign(pow_op(-1.0, 1.0), -1.0)
358 self.assertEqualAndEqualSign(pow_op(-1.0, 2.0), 1.0)
359 self.assertEqualAndEqualSign(pow_op(2.0, -2.0), 0.25)
360 self.assertEqualAndEqualSign(pow_op(2.0, -1.0), 0.5)
361 self.assertEqualAndEqualSign(pow_op(2.0, -0.0), 1.0)
362 self.assertEqualAndEqualSign(pow_op(2.0, 0.0), 1.0)
363 self.assertEqualAndEqualSign(pow_op(2.0, 1.0), 2.0)
364 self.assertEqualAndEqualSign(pow_op(2.0, 2.0), 4.0)
365
366 # 1 ** large and -1 ** large; some libms apparently
367 # have problems with these
368 self.assertEqualAndEqualSign(pow_op(1.0, -1e100), 1.0)
369 self.assertEqualAndEqualSign(pow_op(1.0, 1e100), 1.0)
370 self.assertEqualAndEqualSign(pow_op(-1.0, -1e100), 1.0)
371 self.assertEqualAndEqualSign(pow_op(-1.0, 1e100), 1.0)
372
373 # check sign for results that underflow to 0
374 self.assertEqualAndEqualSign(pow_op(-2.0, -2000.0), 0.0)
375 self.assertEqual(type(pow_op(-2.0, -2000.5)), complex)
376 self.assertEqualAndEqualSign(pow_op(-2.0, -2001.0), -0.0)
377 self.assertEqualAndEqualSign(pow_op(2.0, -2000.0), 0.0)
378 self.assertEqualAndEqualSign(pow_op(2.0, -2000.5), 0.0)
379 self.assertEqualAndEqualSign(pow_op(2.0, -2001.0), 0.0)
380 self.assertEqualAndEqualSign(pow_op(-0.5, 2000.0), 0.0)
381 self.assertEqual(type(pow_op(-0.5, 2000.5)), complex)
382 self.assertEqualAndEqualSign(pow_op(-0.5, 2001.0), -0.0)
383 self.assertEqualAndEqualSign(pow_op(0.5, 2000.0), 0.0)
384 self.assertEqualAndEqualSign(pow_op(0.5, 2000.5), 0.0)
385 self.assertEqualAndEqualSign(pow_op(0.5, 2001.0), 0.0)
386
387 # check we don't raise an exception for subnormal results,
388 # and validate signs. Tests currently disabled, since
389 # they fail on systems where a subnormal result from pow
390 # is flushed to zero (e.g. Debian/ia64.)
391 #self.assertTrue(0.0 < pow_op(0.5, 1048) < 1e-315)
392 #self.assertTrue(0.0 < pow_op(-0.5, 1048) < 1e-315)
393 #self.assertTrue(0.0 < pow_op(0.5, 1047) < 1e-315)
394 #self.assertTrue(0.0 > pow_op(-0.5, 1047) > -1e-315)
395 #self.assertTrue(0.0 < pow_op(2.0, -1048) < 1e-315)
396 #self.assertTrue(0.0 < pow_op(-2.0, -1048) < 1e-315)
397 #self.assertTrue(0.0 < pow_op(2.0, -1047) < 1e-315)
398 #self.assertTrue(0.0 > pow_op(-2.0, -1047) > -1e-315)
Mark Dickinson4a1f5932008-11-12 23:23:36 +0000399
400
Benjamin Petersone401c682010-07-02 23:25:44 +0000401@requires_setformat
Michael W. Hudsonba283e22005-05-27 15:23:20 +0000402class FormatFunctionsTestCase(unittest.TestCase):
403
404 def setUp(self):
405 self.save_formats = {'double':float.__getformat__('double'),
406 'float':float.__getformat__('float')}
407
408 def tearDown(self):
409 float.__setformat__('double', self.save_formats['double'])
410 float.__setformat__('float', self.save_formats['float'])
411
412 def test_getformat(self):
Ezio Melottib58e0bd2010-01-23 15:40:09 +0000413 self.assertIn(float.__getformat__('double'),
414 ['unknown', 'IEEE, big-endian', 'IEEE, little-endian'])
415 self.assertIn(float.__getformat__('float'),
416 ['unknown', 'IEEE, big-endian', 'IEEE, little-endian'])
Michael W. Hudsonba283e22005-05-27 15:23:20 +0000417 self.assertRaises(ValueError, float.__getformat__, 'chicken')
418 self.assertRaises(TypeError, float.__getformat__, 1)
419
420 def test_setformat(self):
421 for t in 'double', 'float':
422 float.__setformat__(t, 'unknown')
423 if self.save_formats[t] == 'IEEE, big-endian':
424 self.assertRaises(ValueError, float.__setformat__,
425 t, 'IEEE, little-endian')
426 elif self.save_formats[t] == 'IEEE, little-endian':
427 self.assertRaises(ValueError, float.__setformat__,
428 t, 'IEEE, big-endian')
429 else:
430 self.assertRaises(ValueError, float.__setformat__,
431 t, 'IEEE, big-endian')
432 self.assertRaises(ValueError, float.__setformat__,
433 t, 'IEEE, little-endian')
434 self.assertRaises(ValueError, float.__setformat__,
435 t, 'chicken')
436 self.assertRaises(ValueError, float.__setformat__,
437 'chicken', 'unknown')
438
Guido van Rossum2be161d2007-05-15 20:43:51 +0000439BE_DOUBLE_INF = b'\x7f\xf0\x00\x00\x00\x00\x00\x00'
Guido van Rossum254348e2007-11-21 19:29:53 +0000440LE_DOUBLE_INF = bytes(reversed(BE_DOUBLE_INF))
Guido van Rossum2be161d2007-05-15 20:43:51 +0000441BE_DOUBLE_NAN = b'\x7f\xf8\x00\x00\x00\x00\x00\x00'
Guido van Rossum254348e2007-11-21 19:29:53 +0000442LE_DOUBLE_NAN = bytes(reversed(BE_DOUBLE_NAN))
Michael W. Hudsonba283e22005-05-27 15:23:20 +0000443
Guido van Rossum2be161d2007-05-15 20:43:51 +0000444BE_FLOAT_INF = b'\x7f\x80\x00\x00'
Guido van Rossum254348e2007-11-21 19:29:53 +0000445LE_FLOAT_INF = bytes(reversed(BE_FLOAT_INF))
Guido van Rossum2be161d2007-05-15 20:43:51 +0000446BE_FLOAT_NAN = b'\x7f\xc0\x00\x00'
Guido van Rossum254348e2007-11-21 19:29:53 +0000447LE_FLOAT_NAN = bytes(reversed(BE_FLOAT_NAN))
Michael W. Hudsonba283e22005-05-27 15:23:20 +0000448
449# on non-IEEE platforms, attempting to unpack a bit pattern
450# representing an infinity or a NaN should raise an exception.
451
Benjamin Petersone401c682010-07-02 23:25:44 +0000452@requires_setformat
Michael W. Hudsonba283e22005-05-27 15:23:20 +0000453class UnknownFormatTestCase(unittest.TestCase):
454 def setUp(self):
455 self.save_formats = {'double':float.__getformat__('double'),
456 'float':float.__getformat__('float')}
457 float.__setformat__('double', 'unknown')
458 float.__setformat__('float', 'unknown')
Tim Peters5d36a552005-06-03 22:40:27 +0000459
Michael W. Hudsonba283e22005-05-27 15:23:20 +0000460 def tearDown(self):
461 float.__setformat__('double', self.save_formats['double'])
462 float.__setformat__('float', self.save_formats['float'])
463
464 def test_double_specials_dont_unpack(self):
465 for fmt, data in [('>d', BE_DOUBLE_INF),
466 ('>d', BE_DOUBLE_NAN),
467 ('<d', LE_DOUBLE_INF),
468 ('<d', LE_DOUBLE_NAN)]:
469 self.assertRaises(ValueError, struct.unpack, fmt, data)
470
471 def test_float_specials_dont_unpack(self):
472 for fmt, data in [('>f', BE_FLOAT_INF),
473 ('>f', BE_FLOAT_NAN),
474 ('<f', LE_FLOAT_INF),
475 ('<f', LE_FLOAT_NAN)]:
476 self.assertRaises(ValueError, struct.unpack, fmt, data)
477
478
479# on an IEEE platform, all we guarantee is that bit patterns
480# representing infinities or NaNs do not raise an exception; all else
481# is accident (today).
Guido van Rossum04110fb2007-08-24 16:32:05 +0000482# let's also try to guarantee that -0.0 and 0.0 don't get confused.
Michael W. Hudsonba283e22005-05-27 15:23:20 +0000483
484class IEEEFormatTestCase(unittest.TestCase):
Michael W. Hudsonba283e22005-05-27 15:23:20 +0000485
Benjamin Petersone401c682010-07-02 23:25:44 +0000486 @requires_IEEE_754
487 def test_double_specials_do_unpack(self):
488 for fmt, data in [('>d', BE_DOUBLE_INF),
489 ('>d', BE_DOUBLE_NAN),
490 ('<d', LE_DOUBLE_INF),
491 ('<d', LE_DOUBLE_NAN)]:
492 struct.unpack(fmt, data)
Michael W. Hudsonba283e22005-05-27 15:23:20 +0000493
Benjamin Petersone401c682010-07-02 23:25:44 +0000494 @requires_IEEE_754
495 def test_float_specials_do_unpack(self):
496 for fmt, data in [('>f', BE_FLOAT_INF),
497 ('>f', BE_FLOAT_NAN),
498 ('<f', LE_FLOAT_INF),
499 ('<f', LE_FLOAT_NAN)]:
500 struct.unpack(fmt, data)
Guido van Rossum04110fb2007-08-24 16:32:05 +0000501
Eric Smith8c663262007-08-25 02:26:07 +0000502class FormatTestCase(unittest.TestCase):
Benjamin Petersone401c682010-07-02 23:25:44 +0000503
Eric Smith11fe3e02007-08-31 01:33:06 +0000504 def test_format(self):
Eric Smith8c663262007-08-25 02:26:07 +0000505 # these should be rewritten to use both format(x, spec) and
506 # x.__format__(spec)
507
508 self.assertEqual(format(0.0, 'f'), '0.000000')
509
510 # the default is 'g', except for empty format spec
511 self.assertEqual(format(0.0, ''), '0.0')
512 self.assertEqual(format(0.01, ''), '0.01')
513 self.assertEqual(format(0.01, 'g'), '0.01')
514
Eric Smith63376222009-05-05 14:04:18 +0000515 # empty presentation type should format in the same way as str
516 # (issue 5920)
517 x = 100/7.
518 self.assertEqual(format(x, ''), str(x))
519 self.assertEqual(format(x, '-'), str(x))
520 self.assertEqual(format(x, '>'), str(x))
521 self.assertEqual(format(x, '2'), str(x))
Eric Smith8c663262007-08-25 02:26:07 +0000522
523 self.assertEqual(format(1.0, 'f'), '1.000000')
Eric Smith8c663262007-08-25 02:26:07 +0000524
525 self.assertEqual(format(-1.0, 'f'), '-1.000000')
Eric Smith8c663262007-08-25 02:26:07 +0000526
527 self.assertEqual(format( 1.0, ' f'), ' 1.000000')
528 self.assertEqual(format(-1.0, ' f'), '-1.000000')
529 self.assertEqual(format( 1.0, '+f'), '+1.000000')
530 self.assertEqual(format(-1.0, '+f'), '-1.000000')
531
532 # % formatting
533 self.assertEqual(format(-1.0, '%'), '-100.000000%')
534
535 # conversion to string should fail
536 self.assertRaises(ValueError, format, 3.0, "s")
537
Eric Smith7b69c6c2008-01-27 21:07:59 +0000538 # other format specifiers shouldn't work on floats,
539 # in particular int specifiers
540 for format_spec in ([chr(x) for x in range(ord('a'), ord('z')+1)] +
541 [chr(x) for x in range(ord('A'), ord('Z')+1)]):
542 if not format_spec in 'eEfFgGn%':
543 self.assertRaises(ValueError, format, 0.0, format_spec)
544 self.assertRaises(ValueError, format, 1.0, format_spec)
545 self.assertRaises(ValueError, format, -1.0, format_spec)
546 self.assertRaises(ValueError, format, 1e100, format_spec)
547 self.assertRaises(ValueError, format, -1e100, format_spec)
548 self.assertRaises(ValueError, format, 1e-100, format_spec)
549 self.assertRaises(ValueError, format, -1e-100, format_spec)
Michael W. Hudsonba283e22005-05-27 15:23:20 +0000550
Eric Smith741191f2009-05-06 13:08:15 +0000551 # issue 3382
552 self.assertEqual(format(NAN, 'f'), 'nan')
553 self.assertEqual(format(NAN, 'F'), 'NAN')
554 self.assertEqual(format(INF, 'f'), 'inf')
555 self.assertEqual(format(INF, 'F'), 'INF')
556
Benjamin Petersone401c682010-07-02 23:25:44 +0000557 @requires_IEEE_754
Eric Smith0923d1d2009-04-16 20:16:10 +0000558 def test_format_testfile(self):
Brian Curtin076623b2010-10-31 00:03:45 +0000559 with open(format_testfile) as testfile:
560 for line in testfile:
561 if line.startswith('--'):
562 continue
563 line = line.strip()
564 if not line:
565 continue
Eric Smith0923d1d2009-04-16 20:16:10 +0000566
Brian Curtin076623b2010-10-31 00:03:45 +0000567 lhs, rhs = map(str.strip, line.split('->'))
568 fmt, arg = lhs.split()
569 self.assertEqual(fmt % float(arg), rhs)
570 self.assertEqual(fmt % -float(arg), '-' + rhs)
Eric Smith0923d1d2009-04-16 20:16:10 +0000571
Mark Dickinsond3ca5572009-04-29 18:47:07 +0000572 def test_issue5864(self):
Ezio Melottib3aedd42010-11-20 19:04:17 +0000573 self.assertEqual(format(123.456, '.4'), '123.5')
574 self.assertEqual(format(1234.56, '.4'), '1.235e+03')
575 self.assertEqual(format(12345.6, '.4'), '1.235e+04')
Mark Dickinsond3ca5572009-04-29 18:47:07 +0000576
Mark Dickinson7efad9e2009-04-17 20:59:58 +0000577class ReprTestCase(unittest.TestCase):
578 def test_repr(self):
579 floats_file = open(os.path.join(os.path.split(__file__)[0],
580 'floating_points.txt'))
581 for line in floats_file:
582 line = line.strip()
583 if not line or line.startswith('#'):
584 continue
585 v = eval(line)
586 self.assertEqual(v, eval(repr(v)))
587 floats_file.close()
588
Eric Smith0923d1d2009-04-16 20:16:10 +0000589 @unittest.skipUnless(getattr(sys, 'float_repr_style', '') == 'short',
590 "applies only when using short float repr style")
591 def test_short_repr(self):
592 # test short float repr introduced in Python 3.1. One aspect
593 # of this repr is that we get some degree of str -> float ->
594 # str roundtripping. In particular, for any numeric string
595 # containing 15 or fewer significant digits, those exact same
596 # digits (modulo trailing zeros) should appear in the output.
597 # No more repr(0.03) -> "0.029999999999999999"!
598
599 test_strings = [
600 # output always includes *either* a decimal point and at
601 # least one digit after that point, or an exponent.
602 '0.0',
603 '1.0',
604 '0.01',
605 '0.02',
606 '0.03',
607 '0.04',
608 '0.05',
609 '1.23456789',
610 '10.0',
611 '100.0',
612 # values >= 1e16 get an exponent...
613 '1000000000000000.0',
614 '9999999999999990.0',
615 '1e+16',
616 '1e+17',
617 # ... and so do values < 1e-4
618 '0.001',
619 '0.001001',
620 '0.00010000000000001',
621 '0.0001',
622 '9.999999999999e-05',
623 '1e-05',
624 # values designed to provoke failure if the FPU rounding
625 # precision isn't set correctly
626 '8.72293771110361e+25',
627 '7.47005307342313e+26',
628 '2.86438000439698e+28',
629 '8.89142905246179e+28',
630 '3.08578087079232e+35',
631 ]
632
633 for s in test_strings:
634 negs = '-'+s
635 self.assertEqual(s, repr(float(s)))
636 self.assertEqual(negs, repr(float(negs)))
Mark Dickinson388122d2010-08-04 20:56:28 +0000637 # Since Python 3.2, repr and str are identical
638 self.assertEqual(repr(float(s)), str(float(s)))
639 self.assertEqual(repr(float(negs)), str(float(negs)))
Benjamin Petersone401c682010-07-02 23:25:44 +0000640
641@requires_IEEE_754
Mark Dickinsone6a076d2009-04-18 11:48:33 +0000642class RoundTestCase(unittest.TestCase):
Benjamin Petersone401c682010-07-02 23:25:44 +0000643
Mark Dickinsone6a076d2009-04-18 11:48:33 +0000644 def test_inf_nan(self):
645 self.assertRaises(OverflowError, round, INF)
646 self.assertRaises(OverflowError, round, -INF)
647 self.assertRaises(ValueError, round, NAN)
Mark Dickinson4ca33d12009-11-24 10:59:34 +0000648 self.assertRaises(TypeError, round, INF, 0.0)
649 self.assertRaises(TypeError, round, -INF, 1.0)
650 self.assertRaises(TypeError, round, NAN, "ceci n'est pas un integer")
651 self.assertRaises(TypeError, round, -0.0, 1j)
Mark Dickinsone6a076d2009-04-18 11:48:33 +0000652
Mark Dickinsone6a076d2009-04-18 11:48:33 +0000653 def test_large_n(self):
654 for n in [324, 325, 400, 2**31-1, 2**31, 2**32, 2**100]:
655 self.assertEqual(round(123.456, n), 123.456)
656 self.assertEqual(round(-123.456, n), -123.456)
657 self.assertEqual(round(1e300, n), 1e300)
658 self.assertEqual(round(1e-320, n), 1e-320)
659 self.assertEqual(round(1e150, 300), 1e150)
660 self.assertEqual(round(1e300, 307), 1e300)
661 self.assertEqual(round(-3.1415, 308), -3.1415)
662 self.assertEqual(round(1e150, 309), 1e150)
663 self.assertEqual(round(1.4e-315, 315), 1e-315)
664
Mark Dickinsone6a076d2009-04-18 11:48:33 +0000665 def test_small_n(self):
666 for n in [-308, -309, -400, 1-2**31, -2**31, -2**31-1, -2**100]:
667 self.assertEqual(round(123.456, n), 0.0)
668 self.assertEqual(round(-123.456, n), -0.0)
669 self.assertEqual(round(1e300, n), 0.0)
670 self.assertEqual(round(1e-320, n), 0.0)
671
Mark Dickinsone6a076d2009-04-18 11:48:33 +0000672 def test_overflow(self):
673 self.assertRaises(OverflowError, round, 1.6e308, -308)
674 self.assertRaises(OverflowError, round, -1.7e308, -308)
675
676 @unittest.skipUnless(getattr(sys, 'float_repr_style', '') == 'short',
677 "applies only when using short float repr style")
678 def test_previous_round_bugs(self):
679 # particular cases that have occurred in bug reports
680 self.assertEqual(round(562949953421312.5, 1),
681 562949953421312.5)
682 self.assertEqual(round(56294995342131.5, 3),
683 56294995342131.5)
684 # round-half-even
685 self.assertEqual(round(25.0, -1), 20.0)
686 self.assertEqual(round(35.0, -1), 40.0)
687 self.assertEqual(round(45.0, -1), 40.0)
688 self.assertEqual(round(55.0, -1), 60.0)
689 self.assertEqual(round(65.0, -1), 60.0)
690 self.assertEqual(round(75.0, -1), 80.0)
691 self.assertEqual(round(85.0, -1), 80.0)
692 self.assertEqual(round(95.0, -1), 100.0)
693
694 @unittest.skipUnless(getattr(sys, 'float_repr_style', '') == 'short',
695 "applies only when using short float repr style")
696 def test_matches_float_format(self):
697 # round should give the same results as float formatting
698 for i in range(500):
699 x = i/1000.
700 self.assertEqual(float(format(x, '.0f')), round(x, 0))
701 self.assertEqual(float(format(x, '.1f')), round(x, 1))
702 self.assertEqual(float(format(x, '.2f')), round(x, 2))
703 self.assertEqual(float(format(x, '.3f')), round(x, 3))
704
705 for i in range(5, 5000, 10):
706 x = i/1000.
707 self.assertEqual(float(format(x, '.0f')), round(x, 0))
708 self.assertEqual(float(format(x, '.1f')), round(x, 1))
709 self.assertEqual(float(format(x, '.2f')), round(x, 2))
710 self.assertEqual(float(format(x, '.3f')), round(x, 3))
711
712 for i in range(500):
713 x = random.random()
714 self.assertEqual(float(format(x, '.0f')), round(x, 0))
715 self.assertEqual(float(format(x, '.1f')), round(x, 1))
716 self.assertEqual(float(format(x, '.2f')), round(x, 2))
717 self.assertEqual(float(format(x, '.3f')), round(x, 3))
718
Eric Smith8a10ecc2009-12-02 17:58:24 +0000719 def test_format_specials(self):
720 # Test formatting of nans and infs.
721
722 def test(fmt, value, expected):
723 # Test with both % and format().
724 self.assertEqual(fmt % value, expected, fmt)
Eric Smith984bb582010-11-25 16:08:06 +0000725 fmt = fmt[1:] # strip off the %
726 self.assertEqual(format(value, fmt), expected, fmt)
Eric Smith8a10ecc2009-12-02 17:58:24 +0000727
728 for fmt in ['%e', '%f', '%g', '%.0e', '%.6f', '%.20g',
729 '%#e', '%#f', '%#g', '%#.20e', '%#.15f', '%#.3g']:
730 pfmt = '%+' + fmt[1:]
731 sfmt = '% ' + fmt[1:]
732 test(fmt, INF, 'inf')
733 test(fmt, -INF, '-inf')
734 test(fmt, NAN, 'nan')
735 test(fmt, -NAN, 'nan')
736 # When asking for a sign, it's always provided. nans are
737 # always positive.
738 test(pfmt, INF, '+inf')
739 test(pfmt, -INF, '-inf')
740 test(pfmt, NAN, '+nan')
741 test(pfmt, -NAN, '+nan')
742 # When using ' ' for a sign code, only infs can be negative.
743 # Others have a space.
744 test(sfmt, INF, ' inf')
745 test(sfmt, -INF, '-inf')
746 test(sfmt, NAN, ' nan')
747 test(sfmt, -NAN, ' nan')
748
Mark Dickinsone6a076d2009-04-18 11:48:33 +0000749
Christian Heimes99170a52007-12-19 02:07:34 +0000750# Beginning with Python 2.6 float has cross platform compatible
Georg Brandl2ee470f2008-07-16 12:55:28 +0000751# ways to create and represent inf and nan
Christian Heimes99170a52007-12-19 02:07:34 +0000752class InfNanTest(unittest.TestCase):
753 def test_inf_from_str(self):
Benjamin Petersonc9c0f202009-06-30 23:06:06 +0000754 self.assertTrue(isinf(float("inf")))
755 self.assertTrue(isinf(float("+inf")))
756 self.assertTrue(isinf(float("-inf")))
757 self.assertTrue(isinf(float("infinity")))
758 self.assertTrue(isinf(float("+infinity")))
759 self.assertTrue(isinf(float("-infinity")))
Christian Heimes99170a52007-12-19 02:07:34 +0000760
761 self.assertEqual(repr(float("inf")), "inf")
762 self.assertEqual(repr(float("+inf")), "inf")
763 self.assertEqual(repr(float("-inf")), "-inf")
Georg Brandl2ee470f2008-07-16 12:55:28 +0000764 self.assertEqual(repr(float("infinity")), "inf")
765 self.assertEqual(repr(float("+infinity")), "inf")
766 self.assertEqual(repr(float("-infinity")), "-inf")
Christian Heimes99170a52007-12-19 02:07:34 +0000767
768 self.assertEqual(repr(float("INF")), "inf")
769 self.assertEqual(repr(float("+Inf")), "inf")
770 self.assertEqual(repr(float("-iNF")), "-inf")
Georg Brandl2ee470f2008-07-16 12:55:28 +0000771 self.assertEqual(repr(float("Infinity")), "inf")
772 self.assertEqual(repr(float("+iNfInItY")), "inf")
773 self.assertEqual(repr(float("-INFINITY")), "-inf")
Christian Heimes99170a52007-12-19 02:07:34 +0000774
775 self.assertEqual(str(float("inf")), "inf")
776 self.assertEqual(str(float("+inf")), "inf")
777 self.assertEqual(str(float("-inf")), "-inf")
Georg Brandl2ee470f2008-07-16 12:55:28 +0000778 self.assertEqual(str(float("infinity")), "inf")
779 self.assertEqual(str(float("+infinity")), "inf")
780 self.assertEqual(str(float("-infinity")), "-inf")
Christian Heimes99170a52007-12-19 02:07:34 +0000781
782 self.assertRaises(ValueError, float, "info")
783 self.assertRaises(ValueError, float, "+info")
784 self.assertRaises(ValueError, float, "-info")
785 self.assertRaises(ValueError, float, "in")
786 self.assertRaises(ValueError, float, "+in")
787 self.assertRaises(ValueError, float, "-in")
Georg Brandl2ee470f2008-07-16 12:55:28 +0000788 self.assertRaises(ValueError, float, "infinit")
789 self.assertRaises(ValueError, float, "+Infin")
790 self.assertRaises(ValueError, float, "-INFI")
791 self.assertRaises(ValueError, float, "infinitys")
Christian Heimes99170a52007-12-19 02:07:34 +0000792
Mark Dickinsonbd16edd2009-05-20 22:05:25 +0000793 self.assertRaises(ValueError, float, "++Inf")
794 self.assertRaises(ValueError, float, "-+inf")
795 self.assertRaises(ValueError, float, "+-infinity")
796 self.assertRaises(ValueError, float, "--Infinity")
797
Christian Heimes99170a52007-12-19 02:07:34 +0000798 def test_inf_as_str(self):
799 self.assertEqual(repr(1e300 * 1e300), "inf")
800 self.assertEqual(repr(-1e300 * 1e300), "-inf")
801
802 self.assertEqual(str(1e300 * 1e300), "inf")
803 self.assertEqual(str(-1e300 * 1e300), "-inf")
804
805 def test_nan_from_str(self):
Benjamin Petersonc9c0f202009-06-30 23:06:06 +0000806 self.assertTrue(isnan(float("nan")))
807 self.assertTrue(isnan(float("+nan")))
808 self.assertTrue(isnan(float("-nan")))
Christian Heimes99170a52007-12-19 02:07:34 +0000809
810 self.assertEqual(repr(float("nan")), "nan")
811 self.assertEqual(repr(float("+nan")), "nan")
812 self.assertEqual(repr(float("-nan")), "nan")
813
814 self.assertEqual(repr(float("NAN")), "nan")
815 self.assertEqual(repr(float("+NAn")), "nan")
816 self.assertEqual(repr(float("-NaN")), "nan")
817
818 self.assertEqual(str(float("nan")), "nan")
819 self.assertEqual(str(float("+nan")), "nan")
820 self.assertEqual(str(float("-nan")), "nan")
821
822 self.assertRaises(ValueError, float, "nana")
823 self.assertRaises(ValueError, float, "+nana")
824 self.assertRaises(ValueError, float, "-nana")
825 self.assertRaises(ValueError, float, "na")
826 self.assertRaises(ValueError, float, "+na")
827 self.assertRaises(ValueError, float, "-na")
828
Mark Dickinsonbd16edd2009-05-20 22:05:25 +0000829 self.assertRaises(ValueError, float, "++nan")
830 self.assertRaises(ValueError, float, "-+NAN")
831 self.assertRaises(ValueError, float, "+-NaN")
832 self.assertRaises(ValueError, float, "--nAn")
833
Christian Heimes99170a52007-12-19 02:07:34 +0000834 def test_nan_as_str(self):
835 self.assertEqual(repr(1e300 * 1e300 * 0), "nan")
836 self.assertEqual(repr(-1e300 * 1e300 * 0), "nan")
837
838 self.assertEqual(str(1e300 * 1e300 * 0), "nan")
839 self.assertEqual(str(-1e300 * 1e300 * 0), "nan")
Christian Heimes827b35c2007-12-10 22:19:17 +0000840
Christian Heimes53876d92008-04-19 00:31:39 +0000841 def notest_float_nan(self):
Benjamin Petersonc9c0f202009-06-30 23:06:06 +0000842 self.assertTrue(NAN.is_nan())
843 self.assertFalse(INF.is_nan())
844 self.assertFalse((0.).is_nan())
Christian Heimes53876d92008-04-19 00:31:39 +0000845
846 def notest_float_inf(self):
Benjamin Petersonc9c0f202009-06-30 23:06:06 +0000847 self.assertTrue(INF.is_inf())
848 self.assertFalse(NAN.is_inf())
849 self.assertFalse((0.).is_inf())
Christian Heimes53876d92008-04-19 00:31:39 +0000850
Mark Dickinson65fe25e2008-07-16 11:30:51 +0000851fromHex = float.fromhex
852toHex = float.hex
853class HexFloatTestCase(unittest.TestCase):
854 MAX = fromHex('0x.fffffffffffff8p+1024') # max normal
855 MIN = fromHex('0x1p-1022') # min normal
856 TINY = fromHex('0x0.0000000000001p-1022') # min subnormal
857 EPS = fromHex('0x0.0000000000001p0') # diff between 1.0 and next float up
858
859 def identical(self, x, y):
860 # check that floats x and y are identical, or that both
861 # are NaNs
862 if isnan(x) or isnan(y):
863 if isnan(x) == isnan(y):
864 return
865 elif x == y and (x != 0.0 or copysign(1.0, x) == copysign(1.0, y)):
866 return
867 self.fail('%r not identical to %r' % (x, y))
868
869 def test_ends(self):
Mark Dickinson38bbc482008-07-16 11:32:23 +0000870 self.identical(self.MIN, ldexp(1.0, -1022))
871 self.identical(self.TINY, ldexp(1.0, -1074))
872 self.identical(self.EPS, ldexp(1.0, -52))
873 self.identical(self.MAX, 2.*(ldexp(1.0, 1023) - ldexp(1.0, 970)))
Mark Dickinson65fe25e2008-07-16 11:30:51 +0000874
875 def test_invalid_inputs(self):
876 invalid_inputs = [
877 'infi', # misspelt infinities and nans
878 '-Infinit',
879 '++inf',
880 '-+Inf',
881 '--nan',
882 '+-NaN',
883 'snan',
884 'NaNs',
885 'nna',
Mark Dickinsond1ec8b22009-05-11 15:45:15 +0000886 'an',
887 'nf',
888 'nfinity',
889 'inity',
890 'iinity',
Mark Dickinson65fe25e2008-07-16 11:30:51 +0000891 '0xnan',
892 '',
893 ' ',
894 'x1.0p0',
895 '0xX1.0p0',
896 '+ 0x1.0p0', # internal whitespace
897 '- 0x1.0p0',
898 '0 x1.0p0',
899 '0x 1.0p0',
900 '0x1 2.0p0',
901 '+0x1 .0p0',
902 '0x1. 0p0',
903 '-0x1.0 1p0',
904 '-0x1.0 p0',
905 '+0x1.0p +0',
906 '0x1.0p -0',
907 '0x1.0p 0',
908 '+0x1.0p+ 0',
909 '-0x1.0p- 0',
910 '++0x1.0p-0', # double signs
911 '--0x1.0p0',
912 '+-0x1.0p+0',
913 '-+0x1.0p0',
914 '0x1.0p++0',
915 '+0x1.0p+-0',
916 '-0x1.0p-+0',
917 '0x1.0p--0',
918 '0x1.0.p0',
919 '0x.p0', # no hex digits before or after point
920 '0x1,p0', # wrong decimal point character
921 '0x1pa',
922 '0x1p\uff10', # fullwidth Unicode digits
923 '\uff10x1p0',
924 '0x\uff11p0',
925 '0x1.\uff10p0',
926 '0x1p0 \n 0x2p0',
927 '0x1p0\0 0x1p0', # embedded null byte is not end of string
928 ]
929 for x in invalid_inputs:
Mark Dickinson589b7952008-08-21 20:05:56 +0000930 try:
931 result = fromHex(x)
932 except ValueError:
933 pass
934 else:
935 self.fail('Expected float.fromhex(%r) to raise ValueError; '
936 'got %r instead' % (x, result))
Mark Dickinson65fe25e2008-07-16 11:30:51 +0000937
938
Mark Dickinsond1ec8b22009-05-11 15:45:15 +0000939 def test_whitespace(self):
940 value_pairs = [
941 ('inf', INF),
942 ('-Infinity', -INF),
943 ('nan', NAN),
944 ('1.0', 1.0),
945 ('-0x.2', -0.125),
946 ('-0.0', -0.0)
947 ]
948 whitespace = [
949 '',
950 ' ',
951 '\t',
952 '\n',
953 '\n \t',
954 '\f',
955 '\v',
956 '\r'
957 ]
958 for inp, expected in value_pairs:
959 for lead in whitespace:
960 for trail in whitespace:
961 got = fromHex(lead + inp + trail)
962 self.identical(got, expected)
963
964
Mark Dickinson65fe25e2008-07-16 11:30:51 +0000965 def test_from_hex(self):
966 MIN = self.MIN;
967 MAX = self.MAX;
968 TINY = self.TINY;
969 EPS = self.EPS;
970
971 # two spellings of infinity, with optional signs; case-insensitive
972 self.identical(fromHex('inf'), INF)
973 self.identical(fromHex('+Inf'), INF)
974 self.identical(fromHex('-INF'), -INF)
975 self.identical(fromHex('iNf'), INF)
976 self.identical(fromHex('Infinity'), INF)
977 self.identical(fromHex('+INFINITY'), INF)
978 self.identical(fromHex('-infinity'), -INF)
979 self.identical(fromHex('-iNFiNitY'), -INF)
980
981 # nans with optional sign; case insensitive
982 self.identical(fromHex('nan'), NAN)
983 self.identical(fromHex('+NaN'), NAN)
984 self.identical(fromHex('-NaN'), NAN)
985 self.identical(fromHex('-nAN'), NAN)
986
987 # variations in input format
988 self.identical(fromHex('1'), 1.0)
989 self.identical(fromHex('+1'), 1.0)
990 self.identical(fromHex('1.'), 1.0)
991 self.identical(fromHex('1.0'), 1.0)
992 self.identical(fromHex('1.0p0'), 1.0)
993 self.identical(fromHex('01'), 1.0)
994 self.identical(fromHex('01.'), 1.0)
995 self.identical(fromHex('0x1'), 1.0)
996 self.identical(fromHex('0x1.'), 1.0)
997 self.identical(fromHex('0x1.0'), 1.0)
998 self.identical(fromHex('+0x1.0'), 1.0)
999 self.identical(fromHex('0x1p0'), 1.0)
1000 self.identical(fromHex('0X1p0'), 1.0)
1001 self.identical(fromHex('0X1P0'), 1.0)
1002 self.identical(fromHex('0x1P0'), 1.0)
1003 self.identical(fromHex('0x1.p0'), 1.0)
1004 self.identical(fromHex('0x1.0p0'), 1.0)
1005 self.identical(fromHex('0x.1p4'), 1.0)
1006 self.identical(fromHex('0x.1p04'), 1.0)
1007 self.identical(fromHex('0x.1p004'), 1.0)
1008 self.identical(fromHex('0x1p+0'), 1.0)
1009 self.identical(fromHex('0x1P-0'), 1.0)
1010 self.identical(fromHex('+0x1p0'), 1.0)
1011 self.identical(fromHex('0x01p0'), 1.0)
1012 self.identical(fromHex('0x1p00'), 1.0)
1013 self.identical(fromHex(' 0x1p0 '), 1.0)
1014 self.identical(fromHex('\n 0x1p0'), 1.0)
1015 self.identical(fromHex('0x1p0 \t'), 1.0)
1016 self.identical(fromHex('0xap0'), 10.0)
1017 self.identical(fromHex('0xAp0'), 10.0)
1018 self.identical(fromHex('0xaP0'), 10.0)
1019 self.identical(fromHex('0xAP0'), 10.0)
1020 self.identical(fromHex('0xbep0'), 190.0)
1021 self.identical(fromHex('0xBep0'), 190.0)
1022 self.identical(fromHex('0xbEp0'), 190.0)
1023 self.identical(fromHex('0XBE0P-4'), 190.0)
1024 self.identical(fromHex('0xBEp0'), 190.0)
1025 self.identical(fromHex('0xB.Ep4'), 190.0)
1026 self.identical(fromHex('0x.BEp8'), 190.0)
1027 self.identical(fromHex('0x.0BEp12'), 190.0)
1028
1029 # moving the point around
1030 pi = fromHex('0x1.921fb54442d18p1')
1031 self.identical(fromHex('0x.006487ed5110b46p11'), pi)
1032 self.identical(fromHex('0x.00c90fdaa22168cp10'), pi)
1033 self.identical(fromHex('0x.01921fb54442d18p9'), pi)
1034 self.identical(fromHex('0x.03243f6a8885a3p8'), pi)
1035 self.identical(fromHex('0x.06487ed5110b46p7'), pi)
1036 self.identical(fromHex('0x.0c90fdaa22168cp6'), pi)
1037 self.identical(fromHex('0x.1921fb54442d18p5'), pi)
1038 self.identical(fromHex('0x.3243f6a8885a3p4'), pi)
1039 self.identical(fromHex('0x.6487ed5110b46p3'), pi)
1040 self.identical(fromHex('0x.c90fdaa22168cp2'), pi)
1041 self.identical(fromHex('0x1.921fb54442d18p1'), pi)
1042 self.identical(fromHex('0x3.243f6a8885a3p0'), pi)
1043 self.identical(fromHex('0x6.487ed5110b46p-1'), pi)
1044 self.identical(fromHex('0xc.90fdaa22168cp-2'), pi)
1045 self.identical(fromHex('0x19.21fb54442d18p-3'), pi)
1046 self.identical(fromHex('0x32.43f6a8885a3p-4'), pi)
1047 self.identical(fromHex('0x64.87ed5110b46p-5'), pi)
1048 self.identical(fromHex('0xc9.0fdaa22168cp-6'), pi)
1049 self.identical(fromHex('0x192.1fb54442d18p-7'), pi)
1050 self.identical(fromHex('0x324.3f6a8885a3p-8'), pi)
1051 self.identical(fromHex('0x648.7ed5110b46p-9'), pi)
1052 self.identical(fromHex('0xc90.fdaa22168cp-10'), pi)
1053 self.identical(fromHex('0x1921.fb54442d18p-11'), pi)
1054 # ...
1055 self.identical(fromHex('0x1921fb54442d1.8p-47'), pi)
1056 self.identical(fromHex('0x3243f6a8885a3p-48'), pi)
1057 self.identical(fromHex('0x6487ed5110b46p-49'), pi)
1058 self.identical(fromHex('0xc90fdaa22168cp-50'), pi)
1059 self.identical(fromHex('0x1921fb54442d18p-51'), pi)
1060 self.identical(fromHex('0x3243f6a8885a30p-52'), pi)
1061 self.identical(fromHex('0x6487ed5110b460p-53'), pi)
1062 self.identical(fromHex('0xc90fdaa22168c0p-54'), pi)
1063 self.identical(fromHex('0x1921fb54442d180p-55'), pi)
1064
1065
1066 # results that should overflow...
1067 self.assertRaises(OverflowError, fromHex, '-0x1p1024')
1068 self.assertRaises(OverflowError, fromHex, '0x1p+1025')
1069 self.assertRaises(OverflowError, fromHex, '+0X1p1030')
1070 self.assertRaises(OverflowError, fromHex, '-0x1p+1100')
1071 self.assertRaises(OverflowError, fromHex, '0X1p123456789123456789')
1072 self.assertRaises(OverflowError, fromHex, '+0X.8p+1025')
1073 self.assertRaises(OverflowError, fromHex, '+0x0.8p1025')
1074 self.assertRaises(OverflowError, fromHex, '-0x0.4p1026')
1075 self.assertRaises(OverflowError, fromHex, '0X2p+1023')
1076 self.assertRaises(OverflowError, fromHex, '0x2.p1023')
1077 self.assertRaises(OverflowError, fromHex, '-0x2.0p+1023')
1078 self.assertRaises(OverflowError, fromHex, '+0X4p+1022')
1079 self.assertRaises(OverflowError, fromHex, '0x1.ffffffffffffffp+1023')
1080 self.assertRaises(OverflowError, fromHex, '-0X1.fffffffffffff9p1023')
1081 self.assertRaises(OverflowError, fromHex, '0X1.fffffffffffff8p1023')
1082 self.assertRaises(OverflowError, fromHex, '+0x3.fffffffffffffp1022')
1083 self.assertRaises(OverflowError, fromHex, '0x3fffffffffffffp+970')
1084 self.assertRaises(OverflowError, fromHex, '0x10000000000000000p960')
1085 self.assertRaises(OverflowError, fromHex, '-0Xffffffffffffffffp960')
1086
1087 # ...and those that round to +-max float
1088 self.identical(fromHex('+0x1.fffffffffffffp+1023'), MAX)
1089 self.identical(fromHex('-0X1.fffffffffffff7p1023'), -MAX)
1090 self.identical(fromHex('0X1.fffffffffffff7fffffffffffffp1023'), MAX)
1091
1092 # zeros
1093 self.identical(fromHex('0x0p0'), 0.0)
1094 self.identical(fromHex('0x0p1000'), 0.0)
1095 self.identical(fromHex('-0x0p1023'), -0.0)
1096 self.identical(fromHex('0X0p1024'), 0.0)
1097 self.identical(fromHex('-0x0p1025'), -0.0)
1098 self.identical(fromHex('0X0p2000'), 0.0)
1099 self.identical(fromHex('0x0p123456789123456789'), 0.0)
1100 self.identical(fromHex('-0X0p-0'), -0.0)
1101 self.identical(fromHex('-0X0p-1000'), -0.0)
1102 self.identical(fromHex('0x0p-1023'), 0.0)
1103 self.identical(fromHex('-0X0p-1024'), -0.0)
1104 self.identical(fromHex('-0x0p-1025'), -0.0)
1105 self.identical(fromHex('-0x0p-1072'), -0.0)
1106 self.identical(fromHex('0X0p-1073'), 0.0)
1107 self.identical(fromHex('-0x0p-1074'), -0.0)
1108 self.identical(fromHex('0x0p-1075'), 0.0)
1109 self.identical(fromHex('0X0p-1076'), 0.0)
1110 self.identical(fromHex('-0X0p-2000'), -0.0)
1111 self.identical(fromHex('-0x0p-123456789123456789'), -0.0)
1112
1113 # values that should underflow to 0
1114 self.identical(fromHex('0X1p-1075'), 0.0)
1115 self.identical(fromHex('-0X1p-1075'), -0.0)
1116 self.identical(fromHex('-0x1p-123456789123456789'), -0.0)
1117 self.identical(fromHex('0x1.00000000000000001p-1075'), TINY)
1118 self.identical(fromHex('-0x1.1p-1075'), -TINY)
1119 self.identical(fromHex('0x1.fffffffffffffffffp-1075'), TINY)
1120
1121 # check round-half-even is working correctly near 0 ...
1122 self.identical(fromHex('0x1p-1076'), 0.0)
1123 self.identical(fromHex('0X2p-1076'), 0.0)
1124 self.identical(fromHex('0X3p-1076'), TINY)
1125 self.identical(fromHex('0x4p-1076'), TINY)
1126 self.identical(fromHex('0X5p-1076'), TINY)
1127 self.identical(fromHex('0X6p-1076'), 2*TINY)
1128 self.identical(fromHex('0x7p-1076'), 2*TINY)
1129 self.identical(fromHex('0X8p-1076'), 2*TINY)
1130 self.identical(fromHex('0X9p-1076'), 2*TINY)
1131 self.identical(fromHex('0xap-1076'), 2*TINY)
1132 self.identical(fromHex('0Xbp-1076'), 3*TINY)
1133 self.identical(fromHex('0xcp-1076'), 3*TINY)
1134 self.identical(fromHex('0Xdp-1076'), 3*TINY)
1135 self.identical(fromHex('0Xep-1076'), 4*TINY)
1136 self.identical(fromHex('0xfp-1076'), 4*TINY)
1137 self.identical(fromHex('0x10p-1076'), 4*TINY)
1138 self.identical(fromHex('-0x1p-1076'), -0.0)
1139 self.identical(fromHex('-0X2p-1076'), -0.0)
1140 self.identical(fromHex('-0x3p-1076'), -TINY)
1141 self.identical(fromHex('-0X4p-1076'), -TINY)
1142 self.identical(fromHex('-0x5p-1076'), -TINY)
1143 self.identical(fromHex('-0x6p-1076'), -2*TINY)
1144 self.identical(fromHex('-0X7p-1076'), -2*TINY)
1145 self.identical(fromHex('-0X8p-1076'), -2*TINY)
1146 self.identical(fromHex('-0X9p-1076'), -2*TINY)
1147 self.identical(fromHex('-0Xap-1076'), -2*TINY)
1148 self.identical(fromHex('-0xbp-1076'), -3*TINY)
1149 self.identical(fromHex('-0xcp-1076'), -3*TINY)
1150 self.identical(fromHex('-0Xdp-1076'), -3*TINY)
1151 self.identical(fromHex('-0xep-1076'), -4*TINY)
1152 self.identical(fromHex('-0Xfp-1076'), -4*TINY)
1153 self.identical(fromHex('-0X10p-1076'), -4*TINY)
1154
1155 # ... and near MIN ...
1156 self.identical(fromHex('0x0.ffffffffffffd6p-1022'), MIN-3*TINY)
1157 self.identical(fromHex('0x0.ffffffffffffd8p-1022'), MIN-2*TINY)
1158 self.identical(fromHex('0x0.ffffffffffffdap-1022'), MIN-2*TINY)
1159 self.identical(fromHex('0x0.ffffffffffffdcp-1022'), MIN-2*TINY)
1160 self.identical(fromHex('0x0.ffffffffffffdep-1022'), MIN-2*TINY)
1161 self.identical(fromHex('0x0.ffffffffffffe0p-1022'), MIN-2*TINY)
1162 self.identical(fromHex('0x0.ffffffffffffe2p-1022'), MIN-2*TINY)
1163 self.identical(fromHex('0x0.ffffffffffffe4p-1022'), MIN-2*TINY)
1164 self.identical(fromHex('0x0.ffffffffffffe6p-1022'), MIN-2*TINY)
1165 self.identical(fromHex('0x0.ffffffffffffe8p-1022'), MIN-2*TINY)
1166 self.identical(fromHex('0x0.ffffffffffffeap-1022'), MIN-TINY)
1167 self.identical(fromHex('0x0.ffffffffffffecp-1022'), MIN-TINY)
1168 self.identical(fromHex('0x0.ffffffffffffeep-1022'), MIN-TINY)
1169 self.identical(fromHex('0x0.fffffffffffff0p-1022'), MIN-TINY)
1170 self.identical(fromHex('0x0.fffffffffffff2p-1022'), MIN-TINY)
1171 self.identical(fromHex('0x0.fffffffffffff4p-1022'), MIN-TINY)
1172 self.identical(fromHex('0x0.fffffffffffff6p-1022'), MIN-TINY)
1173 self.identical(fromHex('0x0.fffffffffffff8p-1022'), MIN)
1174 self.identical(fromHex('0x0.fffffffffffffap-1022'), MIN)
1175 self.identical(fromHex('0x0.fffffffffffffcp-1022'), MIN)
1176 self.identical(fromHex('0x0.fffffffffffffep-1022'), MIN)
1177 self.identical(fromHex('0x1.00000000000000p-1022'), MIN)
1178 self.identical(fromHex('0x1.00000000000002p-1022'), MIN)
1179 self.identical(fromHex('0x1.00000000000004p-1022'), MIN)
1180 self.identical(fromHex('0x1.00000000000006p-1022'), MIN)
1181 self.identical(fromHex('0x1.00000000000008p-1022'), MIN)
1182 self.identical(fromHex('0x1.0000000000000ap-1022'), MIN+TINY)
1183 self.identical(fromHex('0x1.0000000000000cp-1022'), MIN+TINY)
1184 self.identical(fromHex('0x1.0000000000000ep-1022'), MIN+TINY)
1185 self.identical(fromHex('0x1.00000000000010p-1022'), MIN+TINY)
1186 self.identical(fromHex('0x1.00000000000012p-1022'), MIN+TINY)
1187 self.identical(fromHex('0x1.00000000000014p-1022'), MIN+TINY)
1188 self.identical(fromHex('0x1.00000000000016p-1022'), MIN+TINY)
1189 self.identical(fromHex('0x1.00000000000018p-1022'), MIN+2*TINY)
1190
1191 # ... and near 1.0.
1192 self.identical(fromHex('0x0.fffffffffffff0p0'), 1.0-EPS)
1193 self.identical(fromHex('0x0.fffffffffffff1p0'), 1.0-EPS)
1194 self.identical(fromHex('0X0.fffffffffffff2p0'), 1.0-EPS)
1195 self.identical(fromHex('0x0.fffffffffffff3p0'), 1.0-EPS)
1196 self.identical(fromHex('0X0.fffffffffffff4p0'), 1.0-EPS)
1197 self.identical(fromHex('0X0.fffffffffffff5p0'), 1.0-EPS/2)
1198 self.identical(fromHex('0X0.fffffffffffff6p0'), 1.0-EPS/2)
1199 self.identical(fromHex('0x0.fffffffffffff7p0'), 1.0-EPS/2)
1200 self.identical(fromHex('0x0.fffffffffffff8p0'), 1.0-EPS/2)
1201 self.identical(fromHex('0X0.fffffffffffff9p0'), 1.0-EPS/2)
1202 self.identical(fromHex('0X0.fffffffffffffap0'), 1.0-EPS/2)
1203 self.identical(fromHex('0x0.fffffffffffffbp0'), 1.0-EPS/2)
1204 self.identical(fromHex('0X0.fffffffffffffcp0'), 1.0)
1205 self.identical(fromHex('0x0.fffffffffffffdp0'), 1.0)
1206 self.identical(fromHex('0X0.fffffffffffffep0'), 1.0)
1207 self.identical(fromHex('0x0.ffffffffffffffp0'), 1.0)
1208 self.identical(fromHex('0X1.00000000000000p0'), 1.0)
1209 self.identical(fromHex('0X1.00000000000001p0'), 1.0)
1210 self.identical(fromHex('0x1.00000000000002p0'), 1.0)
1211 self.identical(fromHex('0X1.00000000000003p0'), 1.0)
1212 self.identical(fromHex('0x1.00000000000004p0'), 1.0)
1213 self.identical(fromHex('0X1.00000000000005p0'), 1.0)
1214 self.identical(fromHex('0X1.00000000000006p0'), 1.0)
1215 self.identical(fromHex('0X1.00000000000007p0'), 1.0)
1216 self.identical(fromHex('0x1.00000000000007ffffffffffffffffffffp0'),
1217 1.0)
1218 self.identical(fromHex('0x1.00000000000008p0'), 1.0)
1219 self.identical(fromHex('0x1.00000000000008000000000000000001p0'),
1220 1+EPS)
1221 self.identical(fromHex('0X1.00000000000009p0'), 1.0+EPS)
1222 self.identical(fromHex('0x1.0000000000000ap0'), 1.0+EPS)
1223 self.identical(fromHex('0x1.0000000000000bp0'), 1.0+EPS)
1224 self.identical(fromHex('0X1.0000000000000cp0'), 1.0+EPS)
1225 self.identical(fromHex('0x1.0000000000000dp0'), 1.0+EPS)
1226 self.identical(fromHex('0x1.0000000000000ep0'), 1.0+EPS)
1227 self.identical(fromHex('0X1.0000000000000fp0'), 1.0+EPS)
1228 self.identical(fromHex('0x1.00000000000010p0'), 1.0+EPS)
1229 self.identical(fromHex('0X1.00000000000011p0'), 1.0+EPS)
1230 self.identical(fromHex('0x1.00000000000012p0'), 1.0+EPS)
1231 self.identical(fromHex('0X1.00000000000013p0'), 1.0+EPS)
1232 self.identical(fromHex('0X1.00000000000014p0'), 1.0+EPS)
1233 self.identical(fromHex('0x1.00000000000015p0'), 1.0+EPS)
1234 self.identical(fromHex('0x1.00000000000016p0'), 1.0+EPS)
1235 self.identical(fromHex('0X1.00000000000017p0'), 1.0+EPS)
1236 self.identical(fromHex('0x1.00000000000017ffffffffffffffffffffp0'),
1237 1.0+EPS)
1238 self.identical(fromHex('0x1.00000000000018p0'), 1.0+2*EPS)
1239 self.identical(fromHex('0X1.00000000000018000000000000000001p0'),
1240 1.0+2*EPS)
1241 self.identical(fromHex('0x1.00000000000019p0'), 1.0+2*EPS)
1242 self.identical(fromHex('0X1.0000000000001ap0'), 1.0+2*EPS)
1243 self.identical(fromHex('0X1.0000000000001bp0'), 1.0+2*EPS)
1244 self.identical(fromHex('0x1.0000000000001cp0'), 1.0+2*EPS)
1245 self.identical(fromHex('0x1.0000000000001dp0'), 1.0+2*EPS)
1246 self.identical(fromHex('0x1.0000000000001ep0'), 1.0+2*EPS)
1247 self.identical(fromHex('0X1.0000000000001fp0'), 1.0+2*EPS)
1248 self.identical(fromHex('0x1.00000000000020p0'), 1.0+2*EPS)
1249
1250 def test_roundtrip(self):
1251 def roundtrip(x):
1252 return fromHex(toHex(x))
1253
1254 for x in [NAN, INF, self.MAX, self.MIN, self.MIN-self.TINY, self.TINY, 0.0]:
1255 self.identical(x, roundtrip(x))
1256 self.identical(-x, roundtrip(-x))
1257
1258 # fromHex(toHex(x)) should exactly recover x, for any non-NaN float x.
1259 import random
1260 for i in range(10000):
1261 e = random.randrange(-1200, 1200)
1262 m = random.random()
1263 s = random.choice([1.0, -1.0])
1264 try:
1265 x = s*ldexp(m, e)
1266 except OverflowError:
1267 pass
1268 else:
1269 self.identical(x, fromHex(toHex(x)))
1270
Christian Heimes53876d92008-04-19 00:31:39 +00001271
Michael W. Hudsonba283e22005-05-27 15:23:20 +00001272def test_main():
Benjamin Petersonee8712c2008-05-20 21:35:26 +00001273 support.run_unittest(
Amaury Forgeot d'Arc7e958d12008-09-06 21:03:22 +00001274 GeneralFloatCases,
Michael W. Hudsonba283e22005-05-27 15:23:20 +00001275 FormatFunctionsTestCase,
1276 UnknownFormatTestCase,
Eric Smith8c663262007-08-25 02:26:07 +00001277 IEEEFormatTestCase,
Christian Heimes827b35c2007-12-10 22:19:17 +00001278 FormatTestCase,
Christian Heimes99170a52007-12-19 02:07:34 +00001279 ReprTestCase,
Mark Dickinsone6a076d2009-04-18 11:48:33 +00001280 RoundTestCase,
Christian Heimes99170a52007-12-19 02:07:34 +00001281 InfNanTest,
Mark Dickinson65fe25e2008-07-16 11:30:51 +00001282 HexFloatTestCase,
Christian Heimesb76922a2007-12-11 01:06:40 +00001283 )
Michael W. Hudsonba283e22005-05-27 15:23:20 +00001284
1285if __name__ == '__main__':
1286 test_main()