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)[2].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 427 bits of accuracy, more than 128 decimal digits, and its ball lies inside the interval represented by the stored 100-digit decimal.
The attached checks compare every row with elementary theta series at 512 bits, bounding the omitted terms geometrically. They check the quartic and eta identities, using an eta product with an explicit omitted-product bound, and compare the elliptic, Ramanujan and Weber relations with T295, T296 and T297. They also check the elliptic-alpha relation to T299 and the independent Gamma formula at $r=1$. 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.