simonbyrne

import Base: TwicePrecision, significand_bits, significand_mask, exponent_mask, exponent_bias

# Bits of 1/2π
#   1/2π == sum(x / 0x1p64^i for i,x = enumerate(INV2PI))
# Can be obtained by:
#
#    setprecision(BigFloat, 4096)
#    I = 0.5/big(pi)
#    for i = 1:19
#        I *= 0x1p64