organicmaps/fast_double_parser
13
Fork 0
mirror of https://github.com/lemire/fast_double_parser.git synced 2025-04-03 20:04:57 +00:00
Code Issues 1 Projects Releases Wiki Activity Actions

Update table_generation.py

Browse source
This commit is contained in:
Daniel Lemire 2021-01-13 09:17:53 -05:00 • committed by GitHub
parent 6a47c2a5a1
commit d52e3059bd
No known key found for this signature in database
GPG key ID: 4AEE18F83AFDEB23
1 changed files with 10 additions and 8 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

18
script/table_generation.py
View file

@ -1,10 +1,4 @@
def format(number): def format(number):
# move the most significant bit in position
while(number < (1<<127)):
number *= 2
# then *truncate*
while(number >= (1<<128)):
number //= 2
upper = number // (1<<64) upper = number // (1<<64)
lower = number % (1<<64) lower = number % (1<<64)
print(""+hex(upper)+","+hex(lower)+",") print(""+hex(upper)+","+hex(lower)+",")
@ -14,16 +8,24 @@ for q in range(-342,0):
z = 0 z = 0
while( (1<<z) < power5) : while( (1<<z) < power5) :
z += 1 z += 1
if( q >= -17 ): if(q >= -27):
b = z + 127 b = z + 127
c = 2 ** b // power5 + 1 c = 2 ** b // power5 + 1
assert c < (1<<128)
format(c) format(c)
else: else:
b = 2 * z + 2 * 64 b = 2 * z + 2 * 64
c = 2 ** b // power5 + 1 c = 2 ** b // power5 + 1
# truncate
while(c >= (1<<128)):
c //= 2
format(c) format(c)
for q in range(0,308+1): for q in range(0,308+1):
power5 = 5 ** q power5 = 5 ** q
# move the most significant bit in position
while(power5 < (1<<127)):
power5 *= 2
# *truncate*
while(power5 >= (1<<128)):
power5 //= 2
format(power5) format(power5)