import numberdb.sage as numberdb
from sage.rings.complex_arb import ComplexBallField
r = 1
C = ComplexBallField(numberdb.bits(100, losing=96))
tau = C.gen(0) * C(r).sqrt()
C(0).jacobi_theta(tau)[3].real()The generator evaluates the positive square root in $\tau_r$ and the full theta constant using Arb complex balls at 429 bits for 100 requested significant digits. Every value retains at least 426 bits of accuracy, more than 128 decimal digits, and its entire ball lies inside the interval represented by the stored 100-digit decimal.
The attached checks compare every row with elementary theta series and the theta-four product at 512 bits, with explicit bounds for the omitted terms and factors. They check the quartic and eta identities, the elliptic relations against T295 and T296, the Ramanujan and Weber relation against T297, and the companion theta-three values against T504. The square-point Gamma identity is an independent special-value check. Every comparison requires finite enclosures with radii less than $10^{-90}$ and overlap; known-value and failure controls precede the comparisons. The checks are repeated on the stored values.