test_long.py 62 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283848586878889909192939495969798991001011021031041051061071081091101111121131141151161171181191201211221231241251261271281291301311321331341351361371381391401411421431441451461471481491501511521531541551561571581591601611621631641651661671681691701711721731741751761771781791801811821831841851861871881891901911921931941951961971981992002012022032042052062072082092102112122132142152162172182192202212222232242252262272282292302312322332342352362372382392402412422432442452462472482492502512522532542552562572582592602612622632642652662672682692702712722732742752762772782792802812822832842852862872882892902912922932942952962972982993003013023033043053063073083093103113123133143153163173183193203213223233243253263273283293303313323333343353363373383393403413423433443453463473483493503513523533543553563573583593603613623633643653663673683693703713723733743753763773783793803813823833843853863873883893903913923933943953963973983994004014024034044054064074084094104114124134144154164174184194204214224234244254264274284294304314324334344354364374384394404414424434444454464474484494504514524534544554564574584594604614624634644654664674684694704714724734744754764774784794804814824834844854864874884894904914924934944954964974984995005015025035045055065075085095105115125135145155165175185195205215225235245255265275285295305315325335345355365375385395405415425435445455465475485495505515525535545555565575585595605615625635645655665675685695705715725735745755765775785795805815825835845855865875885895905915925935945955965975985996006016026036046056066076086096106116126136146156166176186196206216226236246256266276286296306316326336346356366376386396406416426436446456466476486496506516526536546556566576586596606616626636646656666676686696706716726736746756766776786796806816826836846856866876886896906916926936946956966976986997007017027037047057067077087097107117127137147157167177187197207217227237247257267277287297307317327337347357367377387397407417427437447457467477487497507517527537547557567577587597607617627637647657667677687697707717727737747757767777787797807817827837847857867877887897907917927937947957967977987998008018028038048058068078088098108118128138148158168178188198208218228238248258268278288298308318328338348358368378388398408418428438448458468478488498508518528538548558568578588598608618628638648658668678688698708718728738748758768778788798808818828838848858868878888898908918928938948958968978988999009019029039049059069079089099109119129139149159169179189199209219229239249259269279289299309319329339349359369379389399409419429439449459469479489499509519529539549559569579589599609619629639649659669679689699709719729739749759769779789799809819829839849859869879889899909919929939949959969979989991000100110021003100410051006100710081009101010111012101310141015101610171018101910201021102210231024102510261027102810291030103110321033103410351036103710381039104010411042104310441045104610471048104910501051105210531054105510561057105810591060106110621063106410651066106710681069107010711072107310741075107610771078107910801081108210831084108510861087108810891090109110921093109410951096109710981099110011011102110311041105110611071108110911101111111211131114111511161117111811191120112111221123112411251126112711281129113011311132113311341135113611371138113911401141114211431144114511461147114811491150115111521153115411551156115711581159116011611162116311641165116611671168116911701171117211731174117511761177117811791180118111821183118411851186118711881189119011911192119311941195119611971198119912001201120212031204120512061207120812091210121112121213121412151216121712181219122012211222122312241225122612271228122912301231123212331234123512361237123812391240124112421243124412451246124712481249125012511252125312541255125612571258125912601261126212631264126512661267126812691270127112721273127412751276127712781279128012811282128312841285128612871288128912901291129212931294129512961297129812991300130113021303130413051306130713081309131013111312131313141315131613171318131913201321132213231324132513261327132813291330133113321333133413351336133713381339134013411342134313441345134613471348134913501351135213531354135513561357135813591360136113621363136413651366136713681369137013711372137313741375137613771378137913801381138213831384138513861387138813891390139113921393139413951396139713981399140014011402140314041405140614071408140914101411141214131414141514161417141814191420142114221423142414251426142714281429143014311432143314341435143614371438143914401441144214431444144514461447144814491450145114521453145414551456145714581459146014611462146314641465146614671468146914701471147214731474147514761477147814791480148114821483148414851486148714881489149014911492149314941495149614971498149915001501150215031504150515061507150815091510151115121513151415151516151715181519152015211522152315241525152615271528152915301531153215331534153515361537153815391540154115421543154415451546154715481549155015511552155315541555155615571558155915601561156215631564156515661567156815691570157115721573157415751576157715781579158015811582158315841585158615871588158915901591159215931594159515961597159815991600
  1. import unittest
  2. from test import support
  3. import sys
  4. import random
  5. import math
  6. import array
  7. # SHIFT should match the value in longintrepr.h for best testing.
  8. SHIFT = sys.int_info.bits_per_digit
  9. BASE = 2 ** SHIFT
  10. MASK = BASE - 1
  11. KARATSUBA_CUTOFF = 70 # from longobject.c
  12. # Max number of base BASE digits to use in test cases. Doubling
  13. # this will more than double the runtime.
  14. MAXDIGITS = 15
  15. # build some special values
  16. special = [0, 1, 2, BASE, BASE >> 1, 0x5555555555555555, 0xaaaaaaaaaaaaaaaa]
  17. # some solid strings of one bits
  18. p2 = 4 # 0 and 1 already added
  19. for i in range(2*SHIFT):
  20. special.append(p2 - 1)
  21. p2 = p2 << 1
  22. del p2
  23. # add complements & negations
  24. special += [~x for x in special] + [-x for x in special]
  25. DBL_MAX = sys.float_info.max
  26. DBL_MAX_EXP = sys.float_info.max_exp
  27. DBL_MIN_EXP = sys.float_info.min_exp
  28. DBL_MANT_DIG = sys.float_info.mant_dig
  29. DBL_MIN_OVERFLOW = 2**DBL_MAX_EXP - 2**(DBL_MAX_EXP - DBL_MANT_DIG - 1)
  30. # Pure Python version of correctly-rounded integer-to-float conversion.
  31. def int_to_float(n):
  32. """
  33. Correctly-rounded integer-to-float conversion.
  34. """
  35. # Constants, depending only on the floating-point format in use.
  36. # We use an extra 2 bits of precision for rounding purposes.
  37. PRECISION = sys.float_info.mant_dig + 2
  38. SHIFT_MAX = sys.float_info.max_exp - PRECISION
  39. Q_MAX = 1 << PRECISION
  40. ROUND_HALF_TO_EVEN_CORRECTION = [0, -1, -2, 1, 0, -1, 2, 1]
  41. # Reduce to the case where n is positive.
  42. if n == 0:
  43. return 0.0
  44. elif n < 0:
  45. return -int_to_float(-n)
  46. # Convert n to a 'floating-point' number q * 2**shift, where q is an
  47. # integer with 'PRECISION' significant bits. When shifting n to create q,
  48. # the least significant bit of q is treated as 'sticky'. That is, the
  49. # least significant bit of q is set if either the corresponding bit of n
  50. # was already set, or any one of the bits of n lost in the shift was set.
  51. shift = n.bit_length() - PRECISION
  52. q = n << -shift if shift < 0 else (n >> shift) | bool(n & ~(-1 << shift))
  53. # Round half to even (actually rounds to the nearest multiple of 4,
  54. # rounding ties to a multiple of 8).
  55. q += ROUND_HALF_TO_EVEN_CORRECTION[q & 7]
  56. # Detect overflow.
  57. if shift + (q == Q_MAX) > SHIFT_MAX:
  58. raise OverflowError("integer too large to convert to float")
  59. # Checks: q is exactly representable, and q**2**shift doesn't overflow.
  60. assert q % 4 == 0 and q // 4 <= 2**(sys.float_info.mant_dig)
  61. assert q * 2**shift <= sys.float_info.max
  62. # Some circularity here, since float(q) is doing an int-to-float
  63. # conversion. But here q is of bounded size, and is exactly representable
  64. # as a float. In a low-level C-like language, this operation would be a
  65. # simple cast (e.g., from unsigned long long to double).
  66. return math.ldexp(float(q), shift)
  67. # pure Python version of correctly-rounded true division
  68. def truediv(a, b):
  69. """Correctly-rounded true division for integers."""
  70. negative = a^b < 0
  71. a, b = abs(a), abs(b)
  72. # exceptions: division by zero, overflow
  73. if not b:
  74. raise ZeroDivisionError("division by zero")
  75. if a >= DBL_MIN_OVERFLOW * b:
  76. raise OverflowError("int/int too large to represent as a float")
  77. # find integer d satisfying 2**(d - 1) <= a/b < 2**d
  78. d = a.bit_length() - b.bit_length()
  79. if d >= 0 and a >= 2**d * b or d < 0 and a * 2**-d >= b:
  80. d += 1
  81. # compute 2**-exp * a / b for suitable exp
  82. exp = max(d, DBL_MIN_EXP) - DBL_MANT_DIG
  83. a, b = a << max(-exp, 0), b << max(exp, 0)
  84. q, r = divmod(a, b)
  85. # round-half-to-even: fractional part is r/b, which is > 0.5 iff
  86. # 2*r > b, and == 0.5 iff 2*r == b.
  87. if 2*r > b or 2*r == b and q % 2 == 1:
  88. q += 1
  89. result = math.ldexp(q, exp)
  90. return -result if negative else result
  91. class LongTest(unittest.TestCase):
  92. # Get quasi-random long consisting of ndigits digits (in base BASE).
  93. # quasi == the most-significant digit will not be 0, and the number
  94. # is constructed to contain long strings of 0 and 1 bits. These are
  95. # more likely than random bits to provoke digit-boundary errors.
  96. # The sign of the number is also random.
  97. def getran(self, ndigits):
  98. self.assertGreater(ndigits, 0)
  99. nbits_hi = ndigits * SHIFT
  100. nbits_lo = nbits_hi - SHIFT + 1
  101. answer = 0
  102. nbits = 0
  103. r = int(random.random() * (SHIFT * 2)) | 1 # force 1 bits to start
  104. while nbits < nbits_lo:
  105. bits = (r >> 1) + 1
  106. bits = min(bits, nbits_hi - nbits)
  107. self.assertTrue(1 <= bits <= SHIFT)
  108. nbits = nbits + bits
  109. answer = answer << bits
  110. if r & 1:
  111. answer = answer | ((1 << bits) - 1)
  112. r = int(random.random() * (SHIFT * 2))
  113. self.assertTrue(nbits_lo <= nbits <= nbits_hi)
  114. if random.random() < 0.5:
  115. answer = -answer
  116. return answer
  117. # Get random long consisting of ndigits random digits (relative to base
  118. # BASE). The sign bit is also random.
  119. def getran2(ndigits):
  120. answer = 0
  121. for i in range(ndigits):
  122. answer = (answer << SHIFT) | random.randint(0, MASK)
  123. if random.random() < 0.5:
  124. answer = -answer
  125. return answer
  126. def check_division(self, x, y):
  127. eq = self.assertEqual
  128. with self.subTest(x=x, y=y):
  129. q, r = divmod(x, y)
  130. q2, r2 = x//y, x%y
  131. pab, pba = x*y, y*x
  132. eq(pab, pba, "multiplication does not commute")
  133. eq(q, q2, "divmod returns different quotient than /")
  134. eq(r, r2, "divmod returns different mod than %")
  135. eq(x, q*y + r, "x != q*y + r after divmod")
  136. if y > 0:
  137. self.assertTrue(0 <= r < y, "bad mod from divmod")
  138. else:
  139. self.assertTrue(y < r <= 0, "bad mod from divmod")
  140. def test_division(self):
  141. digits = list(range(1, MAXDIGITS+1)) + list(range(KARATSUBA_CUTOFF,
  142. KARATSUBA_CUTOFF + 14))
  143. digits.append(KARATSUBA_CUTOFF * 3)
  144. for lenx in digits:
  145. x = self.getran(lenx)
  146. for leny in digits:
  147. y = self.getran(leny) or 1
  148. self.check_division(x, y)
  149. # specific numbers chosen to exercise corner cases of the
  150. # current long division implementation
  151. # 30-bit cases involving a quotient digit estimate of BASE+1
  152. self.check_division(1231948412290879395966702881,
  153. 1147341367131428698)
  154. self.check_division(815427756481275430342312021515587883,
  155. 707270836069027745)
  156. self.check_division(627976073697012820849443363563599041,
  157. 643588798496057020)
  158. self.check_division(1115141373653752303710932756325578065,
  159. 1038556335171453937726882627)
  160. # 30-bit cases that require the post-subtraction correction step
  161. self.check_division(922498905405436751940989320930368494,
  162. 949985870686786135626943396)
  163. self.check_division(768235853328091167204009652174031844,
  164. 1091555541180371554426545266)
  165. # 15-bit cases involving a quotient digit estimate of BASE+1
  166. self.check_division(20172188947443, 615611397)
  167. self.check_division(1020908530270155025, 950795710)
  168. self.check_division(128589565723112408, 736393718)
  169. self.check_division(609919780285761575, 18613274546784)
  170. # 15-bit cases that require the post-subtraction correction step
  171. self.check_division(710031681576388032, 26769404391308)
  172. self.check_division(1933622614268221, 30212853348836)
  173. def test_karatsuba(self):
  174. digits = list(range(1, 5)) + list(range(KARATSUBA_CUTOFF,
  175. KARATSUBA_CUTOFF + 10))
  176. digits.extend([KARATSUBA_CUTOFF * 10, KARATSUBA_CUTOFF * 100])
  177. bits = [digit * SHIFT for digit in digits]
  178. # Test products of long strings of 1 bits -- (2**x-1)*(2**y-1) ==
  179. # 2**(x+y) - 2**x - 2**y + 1, so the proper result is easy to check.
  180. for abits in bits:
  181. a = (1 << abits) - 1
  182. for bbits in bits:
  183. if bbits < abits:
  184. continue
  185. with self.subTest(abits=abits, bbits=bbits):
  186. b = (1 << bbits) - 1
  187. x = a * b
  188. y = ((1 << (abits + bbits)) -
  189. (1 << abits) -
  190. (1 << bbits) +
  191. 1)
  192. self.assertEqual(x, y)
  193. def check_bitop_identities_1(self, x):
  194. eq = self.assertEqual
  195. with self.subTest(x=x):
  196. eq(x & 0, 0)
  197. eq(x | 0, x)
  198. eq(x ^ 0, x)
  199. eq(x & -1, x)
  200. eq(x | -1, -1)
  201. eq(x ^ -1, ~x)
  202. eq(x, ~~x)
  203. eq(x & x, x)
  204. eq(x | x, x)
  205. eq(x ^ x, 0)
  206. eq(x & ~x, 0)
  207. eq(x | ~x, -1)
  208. eq(x ^ ~x, -1)
  209. eq(-x, 1 + ~x)
  210. eq(-x, ~(x-1))
  211. for n in range(2*SHIFT):
  212. p2 = 2 ** n
  213. with self.subTest(x=x, n=n, p2=p2):
  214. eq(x << n >> n, x)
  215. eq(x // p2, x >> n)
  216. eq(x * p2, x << n)
  217. eq(x & -p2, x >> n << n)
  218. eq(x & -p2, x & ~(p2 - 1))
  219. def check_bitop_identities_2(self, x, y):
  220. eq = self.assertEqual
  221. with self.subTest(x=x, y=y):
  222. eq(x & y, y & x)
  223. eq(x | y, y | x)
  224. eq(x ^ y, y ^ x)
  225. eq(x ^ y ^ x, y)
  226. eq(x & y, ~(~x | ~y))
  227. eq(x | y, ~(~x & ~y))
  228. eq(x ^ y, (x | y) & ~(x & y))
  229. eq(x ^ y, (x & ~y) | (~x & y))
  230. eq(x ^ y, (x | y) & (~x | ~y))
  231. def check_bitop_identities_3(self, x, y, z):
  232. eq = self.assertEqual
  233. with self.subTest(x=x, y=y, z=z):
  234. eq((x & y) & z, x & (y & z))
  235. eq((x | y) | z, x | (y | z))
  236. eq((x ^ y) ^ z, x ^ (y ^ z))
  237. eq(x & (y | z), (x & y) | (x & z))
  238. eq(x | (y & z), (x | y) & (x | z))
  239. def test_bitop_identities(self):
  240. for x in special:
  241. self.check_bitop_identities_1(x)
  242. digits = range(1, MAXDIGITS+1)
  243. for lenx in digits:
  244. x = self.getran(lenx)
  245. self.check_bitop_identities_1(x)
  246. for leny in digits:
  247. y = self.getran(leny)
  248. self.check_bitop_identities_2(x, y)
  249. self.check_bitop_identities_3(x, y, self.getran((lenx + leny)//2))
  250. def slow_format(self, x, base):
  251. digits = []
  252. sign = 0
  253. if x < 0:
  254. sign, x = 1, -x
  255. while x:
  256. x, r = divmod(x, base)
  257. digits.append(int(r))
  258. digits.reverse()
  259. digits = digits or [0]
  260. return '-'[:sign] + \
  261. {2: '0b', 8: '0o', 10: '', 16: '0x'}[base] + \
  262. "".join("0123456789abcdef"[i] for i in digits)
  263. def check_format_1(self, x):
  264. for base, mapper in (2, bin), (8, oct), (10, str), (10, repr), (16, hex):
  265. got = mapper(x)
  266. with self.subTest(x=x, mapper=mapper.__name__):
  267. expected = self.slow_format(x, base)
  268. self.assertEqual(got, expected)
  269. with self.subTest(got=got):
  270. self.assertEqual(int(got, 0), x)
  271. def test_format(self):
  272. for x in special:
  273. self.check_format_1(x)
  274. for i in range(10):
  275. for lenx in range(1, MAXDIGITS+1):
  276. x = self.getran(lenx)
  277. self.check_format_1(x)
  278. def test_long(self):
  279. # Check conversions from string
  280. LL = [
  281. ('1' + '0'*20, 10**20),
  282. ('1' + '0'*100, 10**100)
  283. ]
  284. for s, v in LL:
  285. for sign in "", "+", "-":
  286. for prefix in "", " ", "\t", " \t\t ":
  287. ss = prefix + sign + s
  288. vv = v
  289. if sign == "-" and v is not ValueError:
  290. vv = -v
  291. try:
  292. self.assertEqual(int(ss), vv)
  293. except ValueError:
  294. pass
  295. # trailing L should no longer be accepted...
  296. self.assertRaises(ValueError, int, '123L')
  297. self.assertRaises(ValueError, int, '123l')
  298. self.assertRaises(ValueError, int, '0L')
  299. self.assertRaises(ValueError, int, '-37L')
  300. self.assertRaises(ValueError, int, '0x32L', 16)
  301. self.assertRaises(ValueError, int, '1L', 21)
  302. # ... but it's just a normal digit if base >= 22
  303. self.assertEqual(int('1L', 22), 43)
  304. # tests with base 0
  305. self.assertEqual(int('000', 0), 0)
  306. self.assertEqual(int('0o123', 0), 83)
  307. self.assertEqual(int('0x123', 0), 291)
  308. self.assertEqual(int('0b100', 0), 4)
  309. self.assertEqual(int(' 0O123 ', 0), 83)
  310. self.assertEqual(int(' 0X123 ', 0), 291)
  311. self.assertEqual(int(' 0B100 ', 0), 4)
  312. self.assertEqual(int('0', 0), 0)
  313. self.assertEqual(int('+0', 0), 0)
  314. self.assertEqual(int('-0', 0), 0)
  315. self.assertEqual(int('00', 0), 0)
  316. self.assertRaises(ValueError, int, '08', 0)
  317. self.assertRaises(ValueError, int, '-012395', 0)
  318. # invalid bases
  319. invalid_bases = [-909,
  320. 2**31-1, 2**31, -2**31, -2**31-1,
  321. 2**63-1, 2**63, -2**63, -2**63-1,
  322. 2**100, -2**100,
  323. ]
  324. for base in invalid_bases:
  325. self.assertRaises(ValueError, int, '42', base)
  326. # Invalid unicode string
  327. # See bpo-34087
  328. self.assertRaises(ValueError, int, '\u3053\u3093\u306b\u3061\u306f')
  329. def test_conversion(self):
  330. class JustLong:
  331. # test that __long__ no longer used in 3.x
  332. def __long__(self):
  333. return 42
  334. self.assertRaises(TypeError, int, JustLong())
  335. class LongTrunc:
  336. # __long__ should be ignored in 3.x
  337. def __long__(self):
  338. return 42
  339. def __trunc__(self):
  340. return 1729
  341. with self.assertWarns(DeprecationWarning):
  342. self.assertEqual(int(LongTrunc()), 1729)
  343. def check_float_conversion(self, n):
  344. # Check that int -> float conversion behaviour matches
  345. # that of the pure Python version above.
  346. try:
  347. actual = float(n)
  348. except OverflowError:
  349. actual = 'overflow'
  350. try:
  351. expected = int_to_float(n)
  352. except OverflowError:
  353. expected = 'overflow'
  354. msg = ("Error in conversion of integer {} to float. "
  355. "Got {}, expected {}.".format(n, actual, expected))
  356. self.assertEqual(actual, expected, msg)
  357. @support.requires_IEEE_754
  358. def test_float_conversion(self):
  359. exact_values = [0, 1, 2,
  360. 2**53-3,
  361. 2**53-2,
  362. 2**53-1,
  363. 2**53,
  364. 2**53+2,
  365. 2**54-4,
  366. 2**54-2,
  367. 2**54,
  368. 2**54+4]
  369. for x in exact_values:
  370. self.assertEqual(float(x), x)
  371. self.assertEqual(float(-x), -x)
  372. # test round-half-even
  373. for x, y in [(1, 0), (2, 2), (3, 4), (4, 4), (5, 4), (6, 6), (7, 8)]:
  374. for p in range(15):
  375. self.assertEqual(int(float(2**p*(2**53+x))), 2**p*(2**53+y))
  376. for x, y in [(0, 0), (1, 0), (2, 0), (3, 4), (4, 4), (5, 4), (6, 8),
  377. (7, 8), (8, 8), (9, 8), (10, 8), (11, 12), (12, 12),
  378. (13, 12), (14, 16), (15, 16)]:
  379. for p in range(15):
  380. self.assertEqual(int(float(2**p*(2**54+x))), 2**p*(2**54+y))
  381. # behaviour near extremes of floating-point range
  382. int_dbl_max = int(DBL_MAX)
  383. top_power = 2**DBL_MAX_EXP
  384. halfway = (int_dbl_max + top_power)//2
  385. self.assertEqual(float(int_dbl_max), DBL_MAX)
  386. self.assertEqual(float(int_dbl_max+1), DBL_MAX)
  387. self.assertEqual(float(halfway-1), DBL_MAX)
  388. self.assertRaises(OverflowError, float, halfway)
  389. self.assertEqual(float(1-halfway), -DBL_MAX)
  390. self.assertRaises(OverflowError, float, -halfway)
  391. self.assertRaises(OverflowError, float, top_power-1)
  392. self.assertRaises(OverflowError, float, top_power)
  393. self.assertRaises(OverflowError, float, top_power+1)
  394. self.assertRaises(OverflowError, float, 2*top_power-1)
  395. self.assertRaises(OverflowError, float, 2*top_power)
  396. self.assertRaises(OverflowError, float, top_power*top_power)
  397. for p in range(100):
  398. x = 2**p * (2**53 + 1) + 1
  399. y = 2**p * (2**53 + 2)
  400. self.assertEqual(int(float(x)), y)
  401. x = 2**p * (2**53 + 1)
  402. y = 2**p * 2**53
  403. self.assertEqual(int(float(x)), y)
  404. # Compare builtin float conversion with pure Python int_to_float
  405. # function above.
  406. test_values = [
  407. int_dbl_max-1, int_dbl_max, int_dbl_max+1,
  408. halfway-1, halfway, halfway + 1,
  409. top_power-1, top_power, top_power+1,
  410. 2*top_power-1, 2*top_power, top_power*top_power,
  411. ]
  412. test_values.extend(exact_values)
  413. for p in range(-4, 8):
  414. for x in range(-128, 128):
  415. test_values.append(2**(p+53) + x)
  416. for value in test_values:
  417. self.check_float_conversion(value)
  418. self.check_float_conversion(-value)
  419. def test_float_overflow(self):
  420. for x in -2.0, -1.0, 0.0, 1.0, 2.0:
  421. self.assertEqual(float(int(x)), x)
  422. shuge = '12345' * 120
  423. huge = 1 << 30000
  424. mhuge = -huge
  425. namespace = {'huge': huge, 'mhuge': mhuge, 'shuge': shuge, 'math': math}
  426. for test in ["float(huge)", "float(mhuge)",
  427. "complex(huge)", "complex(mhuge)",
  428. "complex(huge, 1)", "complex(mhuge, 1)",
  429. "complex(1, huge)", "complex(1, mhuge)",
  430. "1. + huge", "huge + 1.", "1. + mhuge", "mhuge + 1.",
  431. "1. - huge", "huge - 1.", "1. - mhuge", "mhuge - 1.",
  432. "1. * huge", "huge * 1.", "1. * mhuge", "mhuge * 1.",
  433. "1. // huge", "huge // 1.", "1. // mhuge", "mhuge // 1.",
  434. "1. / huge", "huge / 1.", "1. / mhuge", "mhuge / 1.",
  435. "1. ** huge", "huge ** 1.", "1. ** mhuge", "mhuge ** 1.",
  436. "math.sin(huge)", "math.sin(mhuge)",
  437. "math.sqrt(huge)", "math.sqrt(mhuge)", # should do better
  438. # math.floor() of an int returns an int now
  439. ##"math.floor(huge)", "math.floor(mhuge)",
  440. ]:
  441. self.assertRaises(OverflowError, eval, test, namespace)
  442. # XXX Perhaps float(shuge) can raise OverflowError on some box?
  443. # The comparison should not.
  444. self.assertNotEqual(float(shuge), int(shuge),
  445. "float(shuge) should not equal int(shuge)")
  446. def test_logs(self):
  447. LOG10E = math.log10(math.e)
  448. for exp in list(range(10)) + [100, 1000, 10000]:
  449. value = 10 ** exp
  450. log10 = math.log10(value)
  451. self.assertAlmostEqual(log10, exp)
  452. # log10(value) == exp, so log(value) == log10(value)/log10(e) ==
  453. # exp/LOG10E
  454. expected = exp / LOG10E
  455. log = math.log(value)
  456. self.assertAlmostEqual(log, expected)
  457. for bad in -(1 << 10000), -2, 0:
  458. self.assertRaises(ValueError, math.log, bad)
  459. self.assertRaises(ValueError, math.log10, bad)
  460. def test_mixed_compares(self):
  461. eq = self.assertEqual
  462. # We're mostly concerned with that mixing floats and ints does the
  463. # right stuff, even when ints are too large to fit in a float.
  464. # The safest way to check the results is to use an entirely different
  465. # method, which we do here via a skeletal rational class (which
  466. # represents all Python ints and floats exactly).
  467. class Rat:
  468. def __init__(self, value):
  469. if isinstance(value, int):
  470. self.n = value
  471. self.d = 1
  472. elif isinstance(value, float):
  473. # Convert to exact rational equivalent.
  474. f, e = math.frexp(abs(value))
  475. assert f == 0 or 0.5 <= f < 1.0
  476. # |value| = f * 2**e exactly
  477. # Suck up CHUNK bits at a time; 28 is enough so that we suck
  478. # up all bits in 2 iterations for all known binary double-
  479. # precision formats, and small enough to fit in an int.
  480. CHUNK = 28
  481. top = 0
  482. # invariant: |value| = (top + f) * 2**e exactly
  483. while f:
  484. f = math.ldexp(f, CHUNK)
  485. digit = int(f)
  486. assert digit >> CHUNK == 0
  487. top = (top << CHUNK) | digit
  488. f -= digit
  489. assert 0.0 <= f < 1.0
  490. e -= CHUNK
  491. # Now |value| = top * 2**e exactly.
  492. if e >= 0:
  493. n = top << e
  494. d = 1
  495. else:
  496. n = top
  497. d = 1 << -e
  498. if value < 0:
  499. n = -n
  500. self.n = n
  501. self.d = d
  502. assert float(n) / float(d) == value
  503. else:
  504. raise TypeError("can't deal with %r" % value)
  505. def _cmp__(self, other):
  506. if not isinstance(other, Rat):
  507. other = Rat(other)
  508. x, y = self.n * other.d, self.d * other.n
  509. return (x > y) - (x < y)
  510. def __eq__(self, other):
  511. return self._cmp__(other) == 0
  512. def __ge__(self, other):
  513. return self._cmp__(other) >= 0
  514. def __gt__(self, other):
  515. return self._cmp__(other) > 0
  516. def __le__(self, other):
  517. return self._cmp__(other) <= 0
  518. def __lt__(self, other):
  519. return self._cmp__(other) < 0
  520. cases = [0, 0.001, 0.99, 1.0, 1.5, 1e20, 1e200]
  521. # 2**48 is an important boundary in the internals. 2**53 is an
  522. # important boundary for IEEE double precision.
  523. for t in 2.0**48, 2.0**50, 2.0**53:
  524. cases.extend([t - 1.0, t - 0.3, t, t + 0.3, t + 1.0,
  525. int(t-1), int(t), int(t+1)])
  526. cases.extend([0, 1, 2, sys.maxsize, float(sys.maxsize)])
  527. # 1 << 20000 should exceed all double formats. int(1e200) is to
  528. # check that we get equality with 1e200 above.
  529. t = int(1e200)
  530. cases.extend([0, 1, 2, 1 << 20000, t-1, t, t+1])
  531. cases.extend([-x for x in cases])
  532. for x in cases:
  533. Rx = Rat(x)
  534. for y in cases:
  535. Ry = Rat(y)
  536. Rcmp = (Rx > Ry) - (Rx < Ry)
  537. with self.subTest(x=x, y=y, Rcmp=Rcmp):
  538. xycmp = (x > y) - (x < y)
  539. eq(Rcmp, xycmp)
  540. eq(x == y, Rcmp == 0)
  541. eq(x != y, Rcmp != 0)
  542. eq(x < y, Rcmp < 0)
  543. eq(x <= y, Rcmp <= 0)
  544. eq(x > y, Rcmp > 0)
  545. eq(x >= y, Rcmp >= 0)
  546. def test__format__(self):
  547. self.assertEqual(format(123456789, 'd'), '123456789')
  548. self.assertEqual(format(123456789, 'd'), '123456789')
  549. self.assertEqual(format(123456789, ','), '123,456,789')
  550. self.assertEqual(format(123456789, '_'), '123_456_789')
  551. # sign and aligning are interdependent
  552. self.assertEqual(format(1, "-"), '1')
  553. self.assertEqual(format(-1, "-"), '-1')
  554. self.assertEqual(format(1, "-3"), ' 1')
  555. self.assertEqual(format(-1, "-3"), ' -1')
  556. self.assertEqual(format(1, "+3"), ' +1')
  557. self.assertEqual(format(-1, "+3"), ' -1')
  558. self.assertEqual(format(1, " 3"), ' 1')
  559. self.assertEqual(format(-1, " 3"), ' -1')
  560. self.assertEqual(format(1, " "), ' 1')
  561. self.assertEqual(format(-1, " "), '-1')
  562. # hex
  563. self.assertEqual(format(3, "x"), "3")
  564. self.assertEqual(format(3, "X"), "3")
  565. self.assertEqual(format(1234, "x"), "4d2")
  566. self.assertEqual(format(-1234, "x"), "-4d2")
  567. self.assertEqual(format(1234, "8x"), " 4d2")
  568. self.assertEqual(format(-1234, "8x"), " -4d2")
  569. self.assertEqual(format(1234, "x"), "4d2")
  570. self.assertEqual(format(-1234, "x"), "-4d2")
  571. self.assertEqual(format(-3, "x"), "-3")
  572. self.assertEqual(format(-3, "X"), "-3")
  573. self.assertEqual(format(int('be', 16), "x"), "be")
  574. self.assertEqual(format(int('be', 16), "X"), "BE")
  575. self.assertEqual(format(-int('be', 16), "x"), "-be")
  576. self.assertEqual(format(-int('be', 16), "X"), "-BE")
  577. self.assertRaises(ValueError, format, 1234567890, ',x')
  578. self.assertEqual(format(1234567890, '_x'), '4996_02d2')
  579. self.assertEqual(format(1234567890, '_X'), '4996_02D2')
  580. # octal
  581. self.assertEqual(format(3, "o"), "3")
  582. self.assertEqual(format(-3, "o"), "-3")
  583. self.assertEqual(format(1234, "o"), "2322")
  584. self.assertEqual(format(-1234, "o"), "-2322")
  585. self.assertEqual(format(1234, "-o"), "2322")
  586. self.assertEqual(format(-1234, "-o"), "-2322")
  587. self.assertEqual(format(1234, " o"), " 2322")
  588. self.assertEqual(format(-1234, " o"), "-2322")
  589. self.assertEqual(format(1234, "+o"), "+2322")
  590. self.assertEqual(format(-1234, "+o"), "-2322")
  591. self.assertRaises(ValueError, format, 1234567890, ',o')
  592. self.assertEqual(format(1234567890, '_o'), '111_4540_1322')
  593. # binary
  594. self.assertEqual(format(3, "b"), "11")
  595. self.assertEqual(format(-3, "b"), "-11")
  596. self.assertEqual(format(1234, "b"), "10011010010")
  597. self.assertEqual(format(-1234, "b"), "-10011010010")
  598. self.assertEqual(format(1234, "-b"), "10011010010")
  599. self.assertEqual(format(-1234, "-b"), "-10011010010")
  600. self.assertEqual(format(1234, " b"), " 10011010010")
  601. self.assertEqual(format(-1234, " b"), "-10011010010")
  602. self.assertEqual(format(1234, "+b"), "+10011010010")
  603. self.assertEqual(format(-1234, "+b"), "-10011010010")
  604. self.assertRaises(ValueError, format, 1234567890, ',b')
  605. self.assertEqual(format(12345, '_b'), '11_0000_0011_1001')
  606. # make sure these are errors
  607. self.assertRaises(ValueError, format, 3, "1.3") # precision disallowed
  608. self.assertRaises(ValueError, format, 3, "_c") # underscore,
  609. self.assertRaises(ValueError, format, 3, ",c") # comma, and
  610. self.assertRaises(ValueError, format, 3, "+c") # sign not allowed
  611. # with 'c'
  612. self.assertRaisesRegex(ValueError, 'Cannot specify both', format, 3, '_,')
  613. self.assertRaisesRegex(ValueError, 'Cannot specify both', format, 3, ',_')
  614. self.assertRaisesRegex(ValueError, 'Cannot specify both', format, 3, '_,d')
  615. self.assertRaisesRegex(ValueError, 'Cannot specify both', format, 3, ',_d')
  616. self.assertRaisesRegex(ValueError, "Cannot specify ',' with 's'", format, 3, ',s')
  617. self.assertRaisesRegex(ValueError, "Cannot specify '_' with 's'", format, 3, '_s')
  618. # ensure that only int and float type specifiers work
  619. for format_spec in ([chr(x) for x in range(ord('a'), ord('z')+1)] +
  620. [chr(x) for x in range(ord('A'), ord('Z')+1)]):
  621. if not format_spec in 'bcdoxXeEfFgGn%':
  622. self.assertRaises(ValueError, format, 0, format_spec)
  623. self.assertRaises(ValueError, format, 1, format_spec)
  624. self.assertRaises(ValueError, format, -1, format_spec)
  625. self.assertRaises(ValueError, format, 2**100, format_spec)
  626. self.assertRaises(ValueError, format, -(2**100), format_spec)
  627. # ensure that float type specifiers work; format converts
  628. # the int to a float
  629. for format_spec in 'eEfFgG%':
  630. for value in [0, 1, -1, 100, -100, 1234567890, -1234567890]:
  631. self.assertEqual(format(value, format_spec),
  632. format(float(value), format_spec))
  633. def test_nan_inf(self):
  634. self.assertRaises(OverflowError, int, float('inf'))
  635. self.assertRaises(OverflowError, int, float('-inf'))
  636. self.assertRaises(ValueError, int, float('nan'))
  637. def test_mod_division(self):
  638. with self.assertRaises(ZeroDivisionError):
  639. _ = 1 % 0
  640. self.assertEqual(13 % 10, 3)
  641. self.assertEqual(-13 % 10, 7)
  642. self.assertEqual(13 % -10, -7)
  643. self.assertEqual(-13 % -10, -3)
  644. self.assertEqual(12 % 4, 0)
  645. self.assertEqual(-12 % 4, 0)
  646. self.assertEqual(12 % -4, 0)
  647. self.assertEqual(-12 % -4, 0)
  648. def test_true_division(self):
  649. huge = 1 << 40000
  650. mhuge = -huge
  651. self.assertEqual(huge / huge, 1.0)
  652. self.assertEqual(mhuge / mhuge, 1.0)
  653. self.assertEqual(huge / mhuge, -1.0)
  654. self.assertEqual(mhuge / huge, -1.0)
  655. self.assertEqual(1 / huge, 0.0)
  656. self.assertEqual(1 / huge, 0.0)
  657. self.assertEqual(1 / mhuge, 0.0)
  658. self.assertEqual(1 / mhuge, 0.0)
  659. self.assertEqual((666 * huge + (huge >> 1)) / huge, 666.5)
  660. self.assertEqual((666 * mhuge + (mhuge >> 1)) / mhuge, 666.5)
  661. self.assertEqual((666 * huge + (huge >> 1)) / mhuge, -666.5)
  662. self.assertEqual((666 * mhuge + (mhuge >> 1)) / huge, -666.5)
  663. self.assertEqual(huge / (huge << 1), 0.5)
  664. self.assertEqual((1000000 * huge) / huge, 1000000)
  665. namespace = {'huge': huge, 'mhuge': mhuge}
  666. for overflow in ["float(huge)", "float(mhuge)",
  667. "huge / 1", "huge / 2", "huge / -1", "huge / -2",
  668. "mhuge / 100", "mhuge / 200"]:
  669. self.assertRaises(OverflowError, eval, overflow, namespace)
  670. for underflow in ["1 / huge", "2 / huge", "-1 / huge", "-2 / huge",
  671. "100 / mhuge", "200 / mhuge"]:
  672. result = eval(underflow, namespace)
  673. self.assertEqual(result, 0.0,
  674. "expected underflow to 0 from %r" % underflow)
  675. for zero in ["huge / 0", "mhuge / 0"]:
  676. self.assertRaises(ZeroDivisionError, eval, zero, namespace)
  677. def test_floordiv(self):
  678. with self.assertRaises(ZeroDivisionError):
  679. _ = 1 // 0
  680. self.assertEqual(2 // 3, 0)
  681. self.assertEqual(2 // -3, -1)
  682. self.assertEqual(-2 // 3, -1)
  683. self.assertEqual(-2 // -3, 0)
  684. self.assertEqual(-11 // -3, 3)
  685. self.assertEqual(-11 // 3, -4)
  686. self.assertEqual(11 // -3, -4)
  687. self.assertEqual(11 // 3, 3)
  688. self.assertEqual(-12 // -3, 4)
  689. self.assertEqual(-12 // 3, -4)
  690. self.assertEqual(12 // -3, -4)
  691. self.assertEqual(12 // 3, 4)
  692. def check_truediv(self, a, b, skip_small=True):
  693. """Verify that the result of a/b is correctly rounded, by
  694. comparing it with a pure Python implementation of correctly
  695. rounded division. b should be nonzero."""
  696. # skip check for small a and b: in this case, the current
  697. # implementation converts the arguments to float directly and
  698. # then applies a float division. This can give doubly-rounded
  699. # results on x87-using machines (particularly 32-bit Linux).
  700. if skip_small and max(abs(a), abs(b)) < 2**DBL_MANT_DIG:
  701. return
  702. try:
  703. # use repr so that we can distinguish between -0.0 and 0.0
  704. expected = repr(truediv(a, b))
  705. except OverflowError:
  706. expected = 'overflow'
  707. except ZeroDivisionError:
  708. expected = 'zerodivision'
  709. try:
  710. got = repr(a / b)
  711. except OverflowError:
  712. got = 'overflow'
  713. except ZeroDivisionError:
  714. got = 'zerodivision'
  715. self.assertEqual(expected, got, "Incorrectly rounded division {}/{}: "
  716. "expected {}, got {}".format(a, b, expected, got))
  717. @support.requires_IEEE_754
  718. def test_correctly_rounded_true_division(self):
  719. # more stringent tests than those above, checking that the
  720. # result of true division of ints is always correctly rounded.
  721. # This test should probably be considered CPython-specific.
  722. # Exercise all the code paths not involving Gb-sized ints.
  723. # ... divisions involving zero
  724. self.check_truediv(123, 0)
  725. self.check_truediv(-456, 0)
  726. self.check_truediv(0, 3)
  727. self.check_truediv(0, -3)
  728. self.check_truediv(0, 0)
  729. # ... overflow or underflow by large margin
  730. self.check_truediv(671 * 12345 * 2**DBL_MAX_EXP, 12345)
  731. self.check_truediv(12345, 345678 * 2**(DBL_MANT_DIG - DBL_MIN_EXP))
  732. # ... a much larger or smaller than b
  733. self.check_truediv(12345*2**100, 98765)
  734. self.check_truediv(12345*2**30, 98765*7**81)
  735. # ... a / b near a boundary: one of 1, 2**DBL_MANT_DIG, 2**DBL_MIN_EXP,
  736. # 2**DBL_MAX_EXP, 2**(DBL_MIN_EXP-DBL_MANT_DIG)
  737. bases = (0, DBL_MANT_DIG, DBL_MIN_EXP,
  738. DBL_MAX_EXP, DBL_MIN_EXP - DBL_MANT_DIG)
  739. for base in bases:
  740. for exp in range(base - 15, base + 15):
  741. self.check_truediv(75312*2**max(exp, 0), 69187*2**max(-exp, 0))
  742. self.check_truediv(69187*2**max(exp, 0), 75312*2**max(-exp, 0))
  743. # overflow corner case
  744. for m in [1, 2, 7, 17, 12345, 7**100,
  745. -1, -2, -5, -23, -67891, -41**50]:
  746. for n in range(-10, 10):
  747. self.check_truediv(m*DBL_MIN_OVERFLOW + n, m)
  748. self.check_truediv(m*DBL_MIN_OVERFLOW + n, -m)
  749. # check detection of inexactness in shifting stage
  750. for n in range(250):
  751. # (2**DBL_MANT_DIG+1)/(2**DBL_MANT_DIG) lies halfway
  752. # between two representable floats, and would usually be
  753. # rounded down under round-half-to-even. The tiniest of
  754. # additions to the numerator should cause it to be rounded
  755. # up instead.
  756. self.check_truediv((2**DBL_MANT_DIG + 1)*12345*2**200 + 2**n,
  757. 2**DBL_MANT_DIG*12345)
  758. # 1/2731 is one of the smallest division cases that's subject
  759. # to double rounding on IEEE 754 machines working internally with
  760. # 64-bit precision. On such machines, the next check would fail,
  761. # were it not explicitly skipped in check_truediv.
  762. self.check_truediv(1, 2731)
  763. # a particularly bad case for the old algorithm: gives an
  764. # error of close to 3.5 ulps.
  765. self.check_truediv(295147931372582273023, 295147932265116303360)
  766. for i in range(1000):
  767. self.check_truediv(10**(i+1), 10**i)
  768. self.check_truediv(10**i, 10**(i+1))
  769. # test round-half-to-even behaviour, normal result
  770. for m in [1, 2, 4, 7, 8, 16, 17, 32, 12345, 7**100,
  771. -1, -2, -5, -23, -67891, -41**50]:
  772. for n in range(-10, 10):
  773. self.check_truediv(2**DBL_MANT_DIG*m + n, m)
  774. # test round-half-to-even, subnormal result
  775. for n in range(-20, 20):
  776. self.check_truediv(n, 2**1076)
  777. # largeish random divisions: a/b where |a| <= |b| <=
  778. # 2*|a|; |ans| is between 0.5 and 1.0, so error should
  779. # always be bounded by 2**-54 with equality possible only
  780. # if the least significant bit of q=ans*2**53 is zero.
  781. for M in [10**10, 10**100, 10**1000]:
  782. for i in range(1000):
  783. a = random.randrange(1, M)
  784. b = random.randrange(a, 2*a+1)
  785. self.check_truediv(a, b)
  786. self.check_truediv(-a, b)
  787. self.check_truediv(a, -b)
  788. self.check_truediv(-a, -b)
  789. # and some (genuinely) random tests
  790. for _ in range(10000):
  791. a_bits = random.randrange(1000)
  792. b_bits = random.randrange(1, 1000)
  793. x = random.randrange(2**a_bits)
  794. y = random.randrange(1, 2**b_bits)
  795. self.check_truediv(x, y)
  796. self.check_truediv(x, -y)
  797. self.check_truediv(-x, y)
  798. self.check_truediv(-x, -y)
  799. def test_negative_shift_count(self):
  800. with self.assertRaises(ValueError):
  801. 42 << -3
  802. with self.assertRaises(ValueError):
  803. 42 << -(1 << 1000)
  804. with self.assertRaises(ValueError):
  805. 42 >> -3
  806. with self.assertRaises(ValueError):
  807. 42 >> -(1 << 1000)
  808. def test_lshift_of_zero(self):
  809. self.assertEqual(0 << 0, 0)
  810. self.assertEqual(0 << 10, 0)
  811. with self.assertRaises(ValueError):
  812. 0 << -1
  813. self.assertEqual(0 << (1 << 1000), 0)
  814. with self.assertRaises(ValueError):
  815. 0 << -(1 << 1000)
  816. @support.cpython_only
  817. def test_huge_lshift_of_zero(self):
  818. # Shouldn't try to allocate memory for a huge shift. See issue #27870.
  819. # Other implementations may have a different boundary for overflow,
  820. # or not raise at all.
  821. self.assertEqual(0 << sys.maxsize, 0)
  822. self.assertEqual(0 << (sys.maxsize + 1), 0)
  823. @support.cpython_only
  824. @support.bigmemtest(sys.maxsize + 1000, memuse=2/15 * 2, dry_run=False)
  825. def test_huge_lshift(self, size):
  826. self.assertEqual(1 << (sys.maxsize + 1000), 1 << 1000 << sys.maxsize)
  827. def test_huge_rshift(self):
  828. huge_shift = 1 << 1000
  829. self.assertEqual(42 >> huge_shift, 0)
  830. self.assertEqual((-42) >> huge_shift, -1)
  831. self.assertEqual(1123 >> huge_shift, 0)
  832. self.assertEqual((-1123) >> huge_shift, -1)
  833. self.assertEqual(2**128 >> huge_shift, 0)
  834. self.assertEqual(-2**128 >> huge_shift, -1)
  835. @support.cpython_only
  836. @support.bigmemtest(sys.maxsize + 500, memuse=2/15, dry_run=False)
  837. def test_huge_rshift_of_huge(self, size):
  838. huge = ((1 << 500) + 11) << sys.maxsize
  839. self.assertEqual(huge >> (sys.maxsize + 1), (1 << 499) + 5)
  840. self.assertEqual(huge >> (sys.maxsize + 1000), 0)
  841. def test_small_rshift(self):
  842. self.assertEqual(42 >> 1, 21)
  843. self.assertEqual((-42) >> 1, -21)
  844. self.assertEqual(43 >> 1, 21)
  845. self.assertEqual((-43) >> 1, -22)
  846. self.assertEqual(1122 >> 1, 561)
  847. self.assertEqual((-1122) >> 1, -561)
  848. self.assertEqual(1123 >> 1, 561)
  849. self.assertEqual((-1123) >> 1, -562)
  850. self.assertEqual(2**128 >> 1, 2**127)
  851. self.assertEqual(-2**128 >> 1, -2**127)
  852. self.assertEqual((2**128 + 1) >> 1, 2**127)
  853. self.assertEqual(-(2**128 + 1) >> 1, -2**127 - 1)
  854. def test_medium_rshift(self):
  855. self.assertEqual(42 >> 9, 0)
  856. self.assertEqual((-42) >> 9, -1)
  857. self.assertEqual(1122 >> 9, 2)
  858. self.assertEqual((-1122) >> 9, -3)
  859. self.assertEqual(2**128 >> 9, 2**119)
  860. self.assertEqual(-2**128 >> 9, -2**119)
  861. # Exercise corner case of the current algorithm, where the result of
  862. # shifting a two-limb int by the limb size still has two limbs.
  863. self.assertEqual((1 - BASE*BASE) >> SHIFT, -BASE)
  864. self.assertEqual((BASE - 1 - BASE*BASE) >> SHIFT, -BASE)
  865. def test_big_rshift(self):
  866. self.assertEqual(42 >> 32, 0)
  867. self.assertEqual((-42) >> 32, -1)
  868. self.assertEqual(1122 >> 32, 0)
  869. self.assertEqual((-1122) >> 32, -1)
  870. self.assertEqual(2**128 >> 32, 2**96)
  871. self.assertEqual(-2**128 >> 32, -2**96)
  872. def test_small_lshift(self):
  873. self.assertEqual(42 << 1, 84)
  874. self.assertEqual((-42) << 1, -84)
  875. self.assertEqual(561 << 1, 1122)
  876. self.assertEqual((-561) << 1, -1122)
  877. self.assertEqual(2**127 << 1, 2**128)
  878. self.assertEqual(-2**127 << 1, -2**128)
  879. def test_medium_lshift(self):
  880. self.assertEqual(42 << 9, 21504)
  881. self.assertEqual((-42) << 9, -21504)
  882. self.assertEqual(1122 << 9, 574464)
  883. self.assertEqual((-1122) << 9, -574464)
  884. def test_big_lshift(self):
  885. self.assertEqual(42 << 32, 42 * 2**32)
  886. self.assertEqual((-42) << 32, -42 * 2**32)
  887. self.assertEqual(1122 << 32, 1122 * 2**32)
  888. self.assertEqual((-1122) << 32, -1122 * 2**32)
  889. self.assertEqual(2**128 << 32, 2**160)
  890. self.assertEqual(-2**128 << 32, -2**160)
  891. @support.cpython_only
  892. def test_small_ints_in_huge_calculation(self):
  893. a = 2 ** 100
  894. b = -a + 1
  895. c = a + 1
  896. self.assertIs(a + b, 1)
  897. self.assertIs(c - a, 1)
  898. @support.cpython_only
  899. def test_pow_uses_cached_small_ints(self):
  900. self.assertIs(pow(10, 3, 998), 2)
  901. self.assertIs(10 ** 3 % 998, 2)
  902. a, p, m = 10, 3, 998
  903. self.assertIs(a ** p % m, 2)
  904. self.assertIs(pow(2, 31, 2 ** 31 - 1), 1)
  905. self.assertIs(2 ** 31 % (2 ** 31 - 1), 1)
  906. a, p, m = 2, 31, 2 ** 31 - 1
  907. self.assertIs(a ** p % m, 1)
  908. self.assertIs(pow(2, 100, 2**100 - 3), 3)
  909. self.assertIs(2 ** 100 % (2 ** 100 - 3), 3)
  910. a, p, m = 2, 100, 2**100 - 3
  911. self.assertIs(a ** p % m, 3)
  912. @support.cpython_only
  913. def test_divmod_uses_cached_small_ints(self):
  914. big = 10 ** 100
  915. self.assertIs((big + 1) % big, 1)
  916. self.assertIs((big + 1) // big, 1)
  917. self.assertIs(big // (big // 2), 2)
  918. self.assertIs(big // (big // -4), -4)
  919. q, r = divmod(2 * big + 3, big)
  920. self.assertIs(q, 2)
  921. self.assertIs(r, 3)
  922. q, r = divmod(-4 * big + 100, big)
  923. self.assertIs(q, -4)
  924. self.assertIs(r, 100)
  925. q, r = divmod(3 * (-big) - 1, -big)
  926. self.assertIs(q, 3)
  927. self.assertIs(r, -1)
  928. q, r = divmod(3 * big - 1, -big)
  929. self.assertIs(q, -3)
  930. self.assertIs(r, -1)
  931. def test_small_ints(self):
  932. for i in range(-5, 257):
  933. self.assertIs(i, i + 0)
  934. self.assertIs(i, i * 1)
  935. self.assertIs(i, i - 0)
  936. self.assertIs(i, i // 1)
  937. self.assertIs(i, i & -1)
  938. self.assertIs(i, i | 0)
  939. self.assertIs(i, i ^ 0)
  940. self.assertIs(i, ~~i)
  941. self.assertIs(i, i**1)
  942. self.assertIs(i, int(str(i)))
  943. self.assertIs(i, i<<2>>2, str(i))
  944. # corner cases
  945. i = 1 << 70
  946. self.assertIs(i - i, 0)
  947. self.assertIs(0 * i, 0)
  948. def test_bit_length(self):
  949. tiny = 1e-10
  950. for x in range(-65000, 65000):
  951. k = x.bit_length()
  952. # Check equivalence with Python version
  953. self.assertEqual(k, len(bin(x).lstrip('-0b')))
  954. # Behaviour as specified in the docs
  955. if x != 0:
  956. self.assertTrue(2**(k-1) <= abs(x) < 2**k)
  957. else:
  958. self.assertEqual(k, 0)
  959. # Alternative definition: x.bit_length() == 1 + floor(log_2(x))
  960. if x != 0:
  961. # When x is an exact power of 2, numeric errors can
  962. # cause floor(log(x)/log(2)) to be one too small; for
  963. # small x this can be fixed by adding a small quantity
  964. # to the quotient before taking the floor.
  965. self.assertEqual(k, 1 + math.floor(
  966. math.log(abs(x))/math.log(2) + tiny))
  967. self.assertEqual((0).bit_length(), 0)
  968. self.assertEqual((1).bit_length(), 1)
  969. self.assertEqual((-1).bit_length(), 1)
  970. self.assertEqual((2).bit_length(), 2)
  971. self.assertEqual((-2).bit_length(), 2)
  972. for i in [2, 3, 15, 16, 17, 31, 32, 33, 63, 64, 234]:
  973. a = 2**i
  974. self.assertEqual((a-1).bit_length(), i)
  975. self.assertEqual((1-a).bit_length(), i)
  976. self.assertEqual((a).bit_length(), i+1)
  977. self.assertEqual((-a).bit_length(), i+1)
  978. self.assertEqual((a+1).bit_length(), i+1)
  979. self.assertEqual((-a-1).bit_length(), i+1)
  980. def test_bit_count(self):
  981. for a in range(-1000, 1000):
  982. self.assertEqual(a.bit_count(), bin(a).count("1"))
  983. for exp in [10, 17, 63, 64, 65, 1009, 70234, 1234567]:
  984. a = 2**exp
  985. self.assertEqual(a.bit_count(), 1)
  986. self.assertEqual((a - 1).bit_count(), exp)
  987. self.assertEqual((a ^ 63).bit_count(), 7)
  988. self.assertEqual(((a - 1) ^ 510).bit_count(), exp - 8)
  989. def test_round(self):
  990. # check round-half-even algorithm. For round to nearest ten;
  991. # rounding map is invariant under adding multiples of 20
  992. test_dict = {0:0, 1:0, 2:0, 3:0, 4:0, 5:0,
  993. 6:10, 7:10, 8:10, 9:10, 10:10, 11:10, 12:10, 13:10, 14:10,
  994. 15:20, 16:20, 17:20, 18:20, 19:20}
  995. for offset in range(-520, 520, 20):
  996. for k, v in test_dict.items():
  997. got = round(k+offset, -1)
  998. expected = v+offset
  999. self.assertEqual(got, expected)
  1000. self.assertIs(type(got), int)
  1001. # larger second argument
  1002. self.assertEqual(round(-150, -2), -200)
  1003. self.assertEqual(round(-149, -2), -100)
  1004. self.assertEqual(round(-51, -2), -100)
  1005. self.assertEqual(round(-50, -2), 0)
  1006. self.assertEqual(round(-49, -2), 0)
  1007. self.assertEqual(round(-1, -2), 0)
  1008. self.assertEqual(round(0, -2), 0)
  1009. self.assertEqual(round(1, -2), 0)
  1010. self.assertEqual(round(49, -2), 0)
  1011. self.assertEqual(round(50, -2), 0)
  1012. self.assertEqual(round(51, -2), 100)
  1013. self.assertEqual(round(149, -2), 100)
  1014. self.assertEqual(round(150, -2), 200)
  1015. self.assertEqual(round(250, -2), 200)
  1016. self.assertEqual(round(251, -2), 300)
  1017. self.assertEqual(round(172500, -3), 172000)
  1018. self.assertEqual(round(173500, -3), 174000)
  1019. self.assertEqual(round(31415926535, -1), 31415926540)
  1020. self.assertEqual(round(31415926535, -2), 31415926500)
  1021. self.assertEqual(round(31415926535, -3), 31415927000)
  1022. self.assertEqual(round(31415926535, -4), 31415930000)
  1023. self.assertEqual(round(31415926535, -5), 31415900000)
  1024. self.assertEqual(round(31415926535, -6), 31416000000)
  1025. self.assertEqual(round(31415926535, -7), 31420000000)
  1026. self.assertEqual(round(31415926535, -8), 31400000000)
  1027. self.assertEqual(round(31415926535, -9), 31000000000)
  1028. self.assertEqual(round(31415926535, -10), 30000000000)
  1029. self.assertEqual(round(31415926535, -11), 0)
  1030. self.assertEqual(round(31415926535, -12), 0)
  1031. self.assertEqual(round(31415926535, -999), 0)
  1032. # should get correct results even for huge inputs
  1033. for k in range(10, 100):
  1034. got = round(10**k + 324678, -3)
  1035. expect = 10**k + 325000
  1036. self.assertEqual(got, expect)
  1037. self.assertIs(type(got), int)
  1038. # nonnegative second argument: round(x, n) should just return x
  1039. for n in range(5):
  1040. for i in range(100):
  1041. x = random.randrange(-10000, 10000)
  1042. got = round(x, n)
  1043. self.assertEqual(got, x)
  1044. self.assertIs(type(got), int)
  1045. for huge_n in 2**31-1, 2**31, 2**63-1, 2**63, 2**100, 10**100:
  1046. self.assertEqual(round(8979323, huge_n), 8979323)
  1047. # omitted second argument
  1048. for i in range(100):
  1049. x = random.randrange(-10000, 10000)
  1050. got = round(x)
  1051. self.assertEqual(got, x)
  1052. self.assertIs(type(got), int)
  1053. # bad second argument
  1054. bad_exponents = ('brian', 2.0, 0j)
  1055. for e in bad_exponents:
  1056. self.assertRaises(TypeError, round, 3, e)
  1057. def test_to_bytes(self):
  1058. def check(tests, byteorder, signed=False):
  1059. def equivalent_python(n, length, byteorder, signed=False):
  1060. if byteorder == 'little':
  1061. order = range(length)
  1062. elif byteorder == 'big':
  1063. order = reversed(range(length))
  1064. return bytes((n >> i*8) & 0xff for i in order)
  1065. for test, expected in tests.items():
  1066. try:
  1067. self.assertEqual(
  1068. test.to_bytes(len(expected), byteorder, signed=signed),
  1069. expected)
  1070. except Exception as err:
  1071. raise AssertionError(
  1072. "failed to convert {} with byteorder={} and signed={}"
  1073. .format(test, byteorder, signed)) from err
  1074. # Test for all default arguments.
  1075. if len(expected) == 1 and byteorder == 'big' and not signed:
  1076. try:
  1077. self.assertEqual(test.to_bytes(), expected)
  1078. except Exception as err:
  1079. raise AssertionError(
  1080. "failed to convert {} with default arguments"
  1081. .format(test)) from err
  1082. try:
  1083. self.assertEqual(
  1084. equivalent_python(
  1085. test, len(expected), byteorder, signed=signed),
  1086. expected
  1087. )
  1088. except Exception as err:
  1089. raise AssertionError(
  1090. "Code equivalent from docs is not equivalent for "
  1091. "conversion of {0} with byteorder byteorder={1} and "
  1092. "signed={2}".format(test, byteorder, signed)) from err
  1093. # Convert integers to signed big-endian byte arrays.
  1094. tests1 = {
  1095. 0: b'\x00',
  1096. 1: b'\x01',
  1097. -1: b'\xff',
  1098. -127: b'\x81',
  1099. -128: b'\x80',
  1100. -129: b'\xff\x7f',
  1101. 127: b'\x7f',
  1102. 129: b'\x00\x81',
  1103. -255: b'\xff\x01',
  1104. -256: b'\xff\x00',
  1105. 255: b'\x00\xff',
  1106. 256: b'\x01\x00',
  1107. 32767: b'\x7f\xff',
  1108. -32768: b'\xff\x80\x00',
  1109. 65535: b'\x00\xff\xff',
  1110. -65536: b'\xff\x00\x00',
  1111. -8388608: b'\x80\x00\x00'
  1112. }
  1113. check(tests1, 'big', signed=True)
  1114. # Convert integers to signed little-endian byte arrays.
  1115. tests2 = {
  1116. 0: b'\x00',
  1117. 1: b'\x01',
  1118. -1: b'\xff',
  1119. -127: b'\x81',
  1120. -128: b'\x80',
  1121. -129: b'\x7f\xff',
  1122. 127: b'\x7f',
  1123. 129: b'\x81\x00',
  1124. -255: b'\x01\xff',
  1125. -256: b'\x00\xff',
  1126. 255: b'\xff\x00',
  1127. 256: b'\x00\x01',
  1128. 32767: b'\xff\x7f',
  1129. -32768: b'\x00\x80',
  1130. 65535: b'\xff\xff\x00',
  1131. -65536: b'\x00\x00\xff',
  1132. -8388608: b'\x00\x00\x80'
  1133. }
  1134. check(tests2, 'little', signed=True)
  1135. # Convert integers to unsigned big-endian byte arrays.
  1136. tests3 = {
  1137. 0: b'\x00',
  1138. 1: b'\x01',
  1139. 127: b'\x7f',
  1140. 128: b'\x80',
  1141. 255: b'\xff',
  1142. 256: b'\x01\x00',
  1143. 32767: b'\x7f\xff',
  1144. 32768: b'\x80\x00',
  1145. 65535: b'\xff\xff',
  1146. 65536: b'\x01\x00\x00'
  1147. }
  1148. check(tests3, 'big', signed=False)
  1149. # Convert integers to unsigned little-endian byte arrays.
  1150. tests4 = {
  1151. 0: b'\x00',
  1152. 1: b'\x01',
  1153. 127: b'\x7f',
  1154. 128: b'\x80',
  1155. 255: b'\xff',
  1156. 256: b'\x00\x01',
  1157. 32767: b'\xff\x7f',
  1158. 32768: b'\x00\x80',
  1159. 65535: b'\xff\xff',
  1160. 65536: b'\x00\x00\x01'
  1161. }
  1162. check(tests4, 'little', signed=False)
  1163. self.assertRaises(OverflowError, (256).to_bytes, 1, 'big', signed=False)
  1164. self.assertRaises(OverflowError, (256).to_bytes, 1, 'big', signed=True)
  1165. self.assertRaises(OverflowError, (256).to_bytes, 1, 'little', signed=False)
  1166. self.assertRaises(OverflowError, (256).to_bytes, 1, 'little', signed=True)
  1167. self.assertRaises(OverflowError, (-1).to_bytes, 2, 'big', signed=False)
  1168. self.assertRaises(OverflowError, (-1).to_bytes, 2, 'little', signed=False)
  1169. self.assertEqual((0).to_bytes(0, 'big'), b'')
  1170. self.assertEqual((1).to_bytes(5, 'big'), b'\x00\x00\x00\x00\x01')
  1171. self.assertEqual((0).to_bytes(5, 'big'), b'\x00\x00\x00\x00\x00')
  1172. self.assertEqual((-1).to_bytes(5, 'big', signed=True),
  1173. b'\xff\xff\xff\xff\xff')
  1174. self.assertRaises(OverflowError, (1).to_bytes, 0, 'big')
  1175. # gh-98783
  1176. class SubStr(str):
  1177. pass
  1178. self.assertEqual((0).to_bytes(1, SubStr('big')), b'\x00')
  1179. self.assertEqual((0).to_bytes(0, SubStr('little')), b'')
  1180. def test_from_bytes(self):
  1181. def check(tests, byteorder, signed=False):
  1182. def equivalent_python(byte_array, byteorder, signed=False):
  1183. if byteorder == 'little':
  1184. little_ordered = list(byte_array)
  1185. elif byteorder == 'big':
  1186. little_ordered = list(reversed(byte_array))
  1187. n = sum(b << i*8 for i, b in enumerate(little_ordered))
  1188. if signed and little_ordered and (little_ordered[-1] & 0x80):
  1189. n -= 1 << 8*len(little_ordered)
  1190. return n
  1191. for test, expected in tests.items():
  1192. try:
  1193. self.assertEqual(
  1194. int.from_bytes(test, byteorder, signed=signed),
  1195. expected)
  1196. except Exception as err:
  1197. raise AssertionError(
  1198. "failed to convert {} with byteorder={!r} and signed={}"
  1199. .format(test, byteorder, signed)) from err
  1200. # Test for all default arguments.
  1201. if byteorder == 'big' and not signed:
  1202. try:
  1203. self.assertEqual(
  1204. int.from_bytes(test),
  1205. expected)
  1206. except Exception as err:
  1207. raise AssertionError(
  1208. "failed to convert {} with default arguments"
  1209. .format(test)) from err
  1210. try:
  1211. self.assertEqual(
  1212. equivalent_python(test, byteorder, signed=signed),
  1213. expected
  1214. )
  1215. except Exception as err:
  1216. raise AssertionError(
  1217. "Code equivalent from docs is not equivalent for "
  1218. "conversion of {0} with byteorder={1!r} and signed={2}"
  1219. .format(test, byteorder, signed)) from err
  1220. # Convert signed big-endian byte arrays to integers.
  1221. tests1 = {
  1222. b'': 0,
  1223. b'\x00': 0,
  1224. b'\x00\x00': 0,
  1225. b'\x01': 1,
  1226. b'\x00\x01': 1,
  1227. b'\xff': -1,
  1228. b'\xff\xff': -1,
  1229. b'\x81': -127,
  1230. b'\x80': -128,
  1231. b'\xff\x7f': -129,
  1232. b'\x7f': 127,
  1233. b'\x00\x81': 129,
  1234. b'\xff\x01': -255,
  1235. b'\xff\x00': -256,
  1236. b'\x00\xff': 255,
  1237. b'\x01\x00': 256,
  1238. b'\x7f\xff': 32767,
  1239. b'\x80\x00': -32768,
  1240. b'\x00\xff\xff': 65535,
  1241. b'\xff\x00\x00': -65536,
  1242. b'\x80\x00\x00': -8388608
  1243. }
  1244. check(tests1, 'big', signed=True)
  1245. # Convert signed little-endian byte arrays to integers.
  1246. tests2 = {
  1247. b'': 0,
  1248. b'\x00': 0,
  1249. b'\x00\x00': 0,
  1250. b'\x01': 1,
  1251. b'\x00\x01': 256,
  1252. b'\xff': -1,
  1253. b'\xff\xff': -1,
  1254. b'\x81': -127,
  1255. b'\x80': -128,
  1256. b'\x7f\xff': -129,
  1257. b'\x7f': 127,
  1258. b'\x81\x00': 129,
  1259. b'\x01\xff': -255,
  1260. b'\x00\xff': -256,
  1261. b'\xff\x00': 255,
  1262. b'\x00\x01': 256,
  1263. b'\xff\x7f': 32767,
  1264. b'\x00\x80': -32768,
  1265. b'\xff\xff\x00': 65535,
  1266. b'\x00\x00\xff': -65536,
  1267. b'\x00\x00\x80': -8388608
  1268. }
  1269. check(tests2, 'little', signed=True)
  1270. # Convert unsigned big-endian byte arrays to integers.
  1271. tests3 = {
  1272. b'': 0,
  1273. b'\x00': 0,
  1274. b'\x01': 1,
  1275. b'\x7f': 127,
  1276. b'\x80': 128,
  1277. b'\xff': 255,
  1278. b'\x01\x00': 256,
  1279. b'\x7f\xff': 32767,
  1280. b'\x80\x00': 32768,
  1281. b'\xff\xff': 65535,
  1282. b'\x01\x00\x00': 65536,
  1283. }
  1284. check(tests3, 'big', signed=False)
  1285. # Convert integers to unsigned little-endian byte arrays.
  1286. tests4 = {
  1287. b'': 0,
  1288. b'\x00': 0,
  1289. b'\x01': 1,
  1290. b'\x7f': 127,
  1291. b'\x80': 128,
  1292. b'\xff': 255,
  1293. b'\x00\x01': 256,
  1294. b'\xff\x7f': 32767,
  1295. b'\x00\x80': 32768,
  1296. b'\xff\xff': 65535,
  1297. b'\x00\x00\x01': 65536,
  1298. }
  1299. check(tests4, 'little', signed=False)
  1300. class myint(int):
  1301. pass
  1302. self.assertIs(type(myint.from_bytes(b'\x00', 'big')), myint)
  1303. self.assertEqual(myint.from_bytes(b'\x01', 'big'), 1)
  1304. self.assertIs(
  1305. type(myint.from_bytes(b'\x00', 'big', signed=False)), myint)
  1306. self.assertEqual(myint.from_bytes(b'\x01', 'big', signed=False), 1)
  1307. self.assertIs(type(myint.from_bytes(b'\x00', 'little')), myint)
  1308. self.assertEqual(myint.from_bytes(b'\x01', 'little'), 1)
  1309. self.assertIs(type(myint.from_bytes(
  1310. b'\x00', 'little', signed=False)), myint)
  1311. self.assertEqual(myint.from_bytes(b'\x01', 'little', signed=False), 1)
  1312. self.assertEqual(
  1313. int.from_bytes([255, 0, 0], 'big', signed=True), -65536)
  1314. self.assertEqual(
  1315. int.from_bytes((255, 0, 0), 'big', signed=True), -65536)
  1316. self.assertEqual(int.from_bytes(
  1317. bytearray(b'\xff\x00\x00'), 'big', signed=True), -65536)
  1318. self.assertEqual(int.from_bytes(
  1319. bytearray(b'\xff\x00\x00'), 'big', signed=True), -65536)
  1320. self.assertEqual(int.from_bytes(
  1321. array.array('B', b'\xff\x00\x00'), 'big', signed=True), -65536)
  1322. self.assertEqual(int.from_bytes(
  1323. memoryview(b'\xff\x00\x00'), 'big', signed=True), -65536)
  1324. self.assertRaises(ValueError, int.from_bytes, [256], 'big')
  1325. self.assertRaises(ValueError, int.from_bytes, [0], 'big\x00')
  1326. self.assertRaises(ValueError, int.from_bytes, [0], 'little\x00')
  1327. self.assertRaises(TypeError, int.from_bytes, "", 'big')
  1328. self.assertRaises(TypeError, int.from_bytes, "\x00", 'big')
  1329. self.assertRaises(TypeError, int.from_bytes, 0, 'big')
  1330. self.assertRaises(TypeError, int.from_bytes, 0, 'big', True)
  1331. self.assertRaises(TypeError, myint.from_bytes, "", 'big')
  1332. self.assertRaises(TypeError, myint.from_bytes, "\x00", 'big')
  1333. self.assertRaises(TypeError, myint.from_bytes, 0, 'big')
  1334. self.assertRaises(TypeError, int.from_bytes, 0, 'big', True)
  1335. class myint2(int):
  1336. def __new__(cls, value):
  1337. return int.__new__(cls, value + 1)
  1338. i = myint2.from_bytes(b'\x01', 'big')
  1339. self.assertIs(type(i), myint2)
  1340. self.assertEqual(i, 2)
  1341. class myint3(int):
  1342. def __init__(self, value):
  1343. self.foo = 'bar'
  1344. i = myint3.from_bytes(b'\x01', 'big')
  1345. self.assertIs(type(i), myint3)
  1346. self.assertEqual(i, 1)
  1347. self.assertEqual(getattr(i, 'foo', 'none'), 'bar')
  1348. class ValidBytes:
  1349. def __bytes__(self):
  1350. return b'\x01'
  1351. class InvalidBytes:
  1352. def __bytes__(self):
  1353. return 'abc'
  1354. class MissingBytes: ...
  1355. class RaisingBytes:
  1356. def __bytes__(self):
  1357. 1 / 0
  1358. self.assertEqual(int.from_bytes(ValidBytes()), 1)
  1359. self.assertRaises(TypeError, int.from_bytes, InvalidBytes())
  1360. self.assertRaises(TypeError, int.from_bytes, MissingBytes())
  1361. self.assertRaises(ZeroDivisionError, int.from_bytes, RaisingBytes())
  1362. # gh-98783
  1363. class SubStr(str):
  1364. pass
  1365. self.assertEqual(int.from_bytes(b'', SubStr('big')), 0)
  1366. self.assertEqual(int.from_bytes(b'\x00', SubStr('little')), 0)
  1367. @support.cpython_only
  1368. def test_from_bytes_small(self):
  1369. # bpo-46361
  1370. for i in range(-5, 257):
  1371. b = i.to_bytes(2, signed=True)
  1372. self.assertIs(int.from_bytes(b, signed=True), i)
  1373. def test_access_to_nonexistent_digit_0(self):
  1374. # http://bugs.python.org/issue14630: A bug in _PyLong_Copy meant that
  1375. # ob_digit[0] was being incorrectly accessed for instances of a
  1376. # subclass of int, with value 0.
  1377. class Integer(int):
  1378. def __new__(cls, value=0):
  1379. self = int.__new__(cls, value)
  1380. self.foo = 'foo'
  1381. return self
  1382. integers = [Integer(0) for i in range(1000)]
  1383. for n in map(int, integers):
  1384. self.assertEqual(n, 0)
  1385. def test_shift_bool(self):
  1386. # Issue #21422: ensure that bool << int and bool >> int return int
  1387. for value in (True, False):
  1388. for shift in (0, 2):
  1389. self.assertEqual(type(value << shift), int)
  1390. self.assertEqual(type(value >> shift), int)
  1391. def test_as_integer_ratio(self):
  1392. class myint(int):
  1393. pass
  1394. tests = [10, 0, -10, 1, sys.maxsize + 1, True, False, myint(42)]
  1395. for value in tests:
  1396. numerator, denominator = value.as_integer_ratio()
  1397. self.assertEqual((numerator, denominator), (int(value), 1))
  1398. self.assertEqual(type(numerator), int)
  1399. self.assertEqual(type(denominator), int)
  1400. def test_square(self):
  1401. # Multiplication makes a special case of multiplying an int with
  1402. # itself, using a special, faster algorithm. This test is mostly
  1403. # to ensure that no asserts in the implementation trigger, in
  1404. # cases with a maximal amount of carries.
  1405. for bitlen in range(1, 400):
  1406. n = (1 << bitlen) - 1 # solid string of 1 bits
  1407. with self.subTest(bitlen=bitlen, n=n):
  1408. # (2**i - 1)**2 = 2**(2*i) - 2*2**i + 1
  1409. self.assertEqual(n**2,
  1410. (1 << (2 * bitlen)) - (1 << (bitlen + 1)) + 1)
  1411. if __name__ == "__main__":
  1412. unittest.main()