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)[1].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 423 bits of accuracy, more than 127 decimal digits, and its entire ball lies inside the interval represented by the stored 100-digit decimal.
The attached independent checks use real balls at 512 bits, elementary theta series through $n=12$ with geometric tail bounds, and the theta-two and eta products through 80 factors with explicit omitted-product bounds. They check the quartic and eta identities, the elliptic relations against T295 and T296, the Ramanujan and Weber relations against T297, and the companion values against T504, T505 and T506. The square-point Gamma identity is a special-value check outside the theta family. Every comparison requires finite enclosures with radii less than $10^{-90}$ and overlap. Known-value controls and controls that must fail precede the comparisons. The checks are repeated on the stored values.