float2int_intbig.py 2.6 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899
  1. # check cases converting float to int, relying only on single precision float
  2. try:
  3. import ustruct as struct
  4. except:
  5. import struct
  6. import sys
  7. maxsize_bits = 0
  8. maxsize = sys.maxsize
  9. while maxsize:
  10. maxsize >>= 1
  11. maxsize_bits += 1
  12. # work out configuration values
  13. is_64bit = maxsize_bits > 32
  14. # 0 = none, 1 = long long, 2 = mpz
  15. ll_type = None
  16. if is_64bit:
  17. if maxsize_bits < 63:
  18. ll_type = 0
  19. else:
  20. if maxsize_bits < 31:
  21. ll_type = 0
  22. if ll_type is None:
  23. one = 1
  24. if one << 65 < one << 62:
  25. ll_type = 1
  26. else:
  27. ll_type = 2
  28. # basic conversion
  29. print(int(14187745.))
  30. print("%d" % 14187745.)
  31. if ll_type == 2:
  32. print(int(2.**100))
  33. print("%d" % 2.**100)
  34. testpass = True
  35. p2_rng = ((30,63,127),(62,63,127))[is_64bit][ll_type]
  36. for i in range(0,p2_rng):
  37. bitcnt = len(bin(int(2.**i))) - 3;
  38. if i != bitcnt:
  39. print('fail: 2.**%u was %u bits long' % (i, bitcnt));
  40. testpass = False
  41. print("power of 2 test: %s" % (testpass and 'passed' or 'failed'))
  42. # TODO why does 10**12 fail this test for single precision float?
  43. testpass = True
  44. p10_rng = 9 if (ll_type == 0 and ~is_64bit) else 11
  45. for i in range(0,p10_rng):
  46. digcnt = len(str(int(10.**i))) - 1;
  47. if i != digcnt:
  48. print('fail: 10.**%u was %u digits long' % (i, digcnt));
  49. testpass = False
  50. print("power of 10 test: %s" % (testpass and 'passed' or 'failed'))
  51. def fp2int_test(num, name, should_fail):
  52. try:
  53. x = int(num)
  54. passed = ~should_fail
  55. except:
  56. passed = should_fail
  57. print('%s: %s' % (name, passed and 'passed' or 'failed'))
  58. if ll_type != 2:
  59. if ll_type == 0:
  60. if is_64bit:
  61. neg_bad_fp = -1.00000005*2.**62.
  62. pos_bad_fp = 2.**62.
  63. neg_good_fp = -2.**62.
  64. pos_good_fp = 0.99999993*2.**62.
  65. else:
  66. neg_bad_fp = -1.00000005*2.**30.
  67. pos_bad_fp = 2.**30.
  68. neg_good_fp = -2.**30.
  69. pos_good_fp = 0.9999999499*2.**30.
  70. else:
  71. neg_bad_fp = -0.51*2.**64.
  72. pos_bad_fp = 2.**63.
  73. neg_good_fp = -2.**63.
  74. pos_good_fp = 1.9999998*2.**62.
  75. fp2int_test(neg_bad_fp, 'neg bad', True)
  76. fp2int_test(pos_bad_fp, 'pos bad', True)
  77. fp2int_test(neg_good_fp, 'neg good', False)
  78. fp2int_test(pos_good_fp, 'pos good', False)
  79. else:
  80. fp2int_test(-1.999999879*2.**127., 'large neg', False)
  81. fp2int_test(1.999999879*2.**127., 'large pos', False)
  82. fp2int_test(float('inf'), 'inf test', True)
  83. fp2int_test(float('nan'), 'NaN test', True)
  84. # test numbers < 1 (this used to fail; see issue #1044)
  85. fp2int_test(0.0001, 'small num', False)
  86. struct.pack('I', int(1/2))