import numberdb.sage as numberdb
from sage.rings.rational_field import QQ
from sage.rings.real_arb import RealBallField
R = RealBallField(numberdb.bits(100, losing=128))
pi = R.pi()
gamma = lambda a, b: R(QQ(a) / QQ(b)).gamma()
two_14_3 = (R(QQ(14) / QQ(3)) * R(2).log()).exp()
W_Z3 = (R(6).sqrt() * gamma(1, 24) * gamma(5, 24) *
gamma(7, 24) * gamma(11, 24) / (32 * pi**3))
W_A3 = 9 * gamma(1, 3)**6 / (two_14_3 * pi**4)
W_A3dual = gamma(1, 4)**4 / (4 * pi**3)Each value is a real ball at numberdb.bits(digits, losing=128) bits, computed from (2) by arb's gamma function. The generator writes 100 digits.
Before the draft was filled, each gamma-product value was compared with the corresponding complete-elliptic-integral form of Watson's triple integrals, using arb ball arithmetic.