cachedpower.go 3.9 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120
  1. package fast
  2. import "math"
  3. const (
  4. kCachedPowersOffset = 348 // -1 * the first decimal_exponent.
  5. kD_1_LOG2_10 = 0.30102999566398114 // 1 / lg(10)
  6. kDecimalExponentDistance = 8
  7. )
  8. type cachedPower struct {
  9. significand uint64
  10. binary_exponent int16
  11. decimal_exponent int16
  12. }
  13. var (
  14. cachedPowers = [...]cachedPower{
  15. {0xFA8FD5A0081C0288, -1220, -348},
  16. {0xBAAEE17FA23EBF76, -1193, -340},
  17. {0x8B16FB203055AC76, -1166, -332},
  18. {0xCF42894A5DCE35EA, -1140, -324},
  19. {0x9A6BB0AA55653B2D, -1113, -316},
  20. {0xE61ACF033D1A45DF, -1087, -308},
  21. {0xAB70FE17C79AC6CA, -1060, -300},
  22. {0xFF77B1FCBEBCDC4F, -1034, -292},
  23. {0xBE5691EF416BD60C, -1007, -284},
  24. {0x8DD01FAD907FFC3C, -980, -276},
  25. {0xD3515C2831559A83, -954, -268},
  26. {0x9D71AC8FADA6C9B5, -927, -260},
  27. {0xEA9C227723EE8BCB, -901, -252},
  28. {0xAECC49914078536D, -874, -244},
  29. {0x823C12795DB6CE57, -847, -236},
  30. {0xC21094364DFB5637, -821, -228},
  31. {0x9096EA6F3848984F, -794, -220},
  32. {0xD77485CB25823AC7, -768, -212},
  33. {0xA086CFCD97BF97F4, -741, -204},
  34. {0xEF340A98172AACE5, -715, -196},
  35. {0xB23867FB2A35B28E, -688, -188},
  36. {0x84C8D4DFD2C63F3B, -661, -180},
  37. {0xC5DD44271AD3CDBA, -635, -172},
  38. {0x936B9FCEBB25C996, -608, -164},
  39. {0xDBAC6C247D62A584, -582, -156},
  40. {0xA3AB66580D5FDAF6, -555, -148},
  41. {0xF3E2F893DEC3F126, -529, -140},
  42. {0xB5B5ADA8AAFF80B8, -502, -132},
  43. {0x87625F056C7C4A8B, -475, -124},
  44. {0xC9BCFF6034C13053, -449, -116},
  45. {0x964E858C91BA2655, -422, -108},
  46. {0xDFF9772470297EBD, -396, -100},
  47. {0xA6DFBD9FB8E5B88F, -369, -92},
  48. {0xF8A95FCF88747D94, -343, -84},
  49. {0xB94470938FA89BCF, -316, -76},
  50. {0x8A08F0F8BF0F156B, -289, -68},
  51. {0xCDB02555653131B6, -263, -60},
  52. {0x993FE2C6D07B7FAC, -236, -52},
  53. {0xE45C10C42A2B3B06, -210, -44},
  54. {0xAA242499697392D3, -183, -36},
  55. {0xFD87B5F28300CA0E, -157, -28},
  56. {0xBCE5086492111AEB, -130, -20},
  57. {0x8CBCCC096F5088CC, -103, -12},
  58. {0xD1B71758E219652C, -77, -4},
  59. {0x9C40000000000000, -50, 4},
  60. {0xE8D4A51000000000, -24, 12},
  61. {0xAD78EBC5AC620000, 3, 20},
  62. {0x813F3978F8940984, 30, 28},
  63. {0xC097CE7BC90715B3, 56, 36},
  64. {0x8F7E32CE7BEA5C70, 83, 44},
  65. {0xD5D238A4ABE98068, 109, 52},
  66. {0x9F4F2726179A2245, 136, 60},
  67. {0xED63A231D4C4FB27, 162, 68},
  68. {0xB0DE65388CC8ADA8, 189, 76},
  69. {0x83C7088E1AAB65DB, 216, 84},
  70. {0xC45D1DF942711D9A, 242, 92},
  71. {0x924D692CA61BE758, 269, 100},
  72. {0xDA01EE641A708DEA, 295, 108},
  73. {0xA26DA3999AEF774A, 322, 116},
  74. {0xF209787BB47D6B85, 348, 124},
  75. {0xB454E4A179DD1877, 375, 132},
  76. {0x865B86925B9BC5C2, 402, 140},
  77. {0xC83553C5C8965D3D, 428, 148},
  78. {0x952AB45CFA97A0B3, 455, 156},
  79. {0xDE469FBD99A05FE3, 481, 164},
  80. {0xA59BC234DB398C25, 508, 172},
  81. {0xF6C69A72A3989F5C, 534, 180},
  82. {0xB7DCBF5354E9BECE, 561, 188},
  83. {0x88FCF317F22241E2, 588, 196},
  84. {0xCC20CE9BD35C78A5, 614, 204},
  85. {0x98165AF37B2153DF, 641, 212},
  86. {0xE2A0B5DC971F303A, 667, 220},
  87. {0xA8D9D1535CE3B396, 694, 228},
  88. {0xFB9B7CD9A4A7443C, 720, 236},
  89. {0xBB764C4CA7A44410, 747, 244},
  90. {0x8BAB8EEFB6409C1A, 774, 252},
  91. {0xD01FEF10A657842C, 800, 260},
  92. {0x9B10A4E5E9913129, 827, 268},
  93. {0xE7109BFBA19C0C9D, 853, 276},
  94. {0xAC2820D9623BF429, 880, 284},
  95. {0x80444B5E7AA7CF85, 907, 292},
  96. {0xBF21E44003ACDD2D, 933, 300},
  97. {0x8E679C2F5E44FF8F, 960, 308},
  98. {0xD433179D9C8CB841, 986, 316},
  99. {0x9E19DB92B4E31BA9, 1013, 324},
  100. {0xEB96BF6EBADF77D9, 1039, 332},
  101. {0xAF87023B9BF0EE6B, 1066, 340},
  102. }
  103. )
  104. func getCachedPowerForBinaryExponentRange(min_exponent, max_exponent int) (power diyfp, decimal_exponent int) {
  105. kQ := diyFpKSignificandSize
  106. k := int(math.Ceil(float64(min_exponent+kQ-1) * kD_1_LOG2_10))
  107. index := (kCachedPowersOffset+k-1)/kDecimalExponentDistance + 1
  108. cached_power := cachedPowers[index]
  109. _DCHECK(min_exponent <= int(cached_power.binary_exponent))
  110. _DCHECK(int(cached_power.binary_exponent) <= max_exponent)
  111. decimal_exponent = int(cached_power.decimal_exponent)
  112. power = diyfp{f: cached_power.significand, e: int(cached_power.binary_exponent)}
  113. return
  114. }