script/table_generation.py

def format(number):
    upper = number // (1 << 64)
    lower = number % (1 << 64)
    print("" + hex(upper) + "," + hex(lower) + ",")


for q in range(-342, 0):
    power5 = 5 ** -q
    z = 0
    while (1 << z) < power5:
        z += 1
    if q >= -27:
        b = z + 127
        c = 2 ** b // power5 + 1
        format(c)
    else:
        b = 2 * z + 2 * 64
        c = 2 ** b // power5 + 1
        # truncate
        while c >= (1 << 128):
            c //= 2
        format(c)

for q in range(0, 308 + 1):
    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)