back to table · edit · history · where entries came from · files · download
3598 bytes, as of the version from 2026-10-29 07:39 (current). Recorded here, not run.
"""Reproduce all numerical and arithmetic checks for numberdb.org/T603.
sage -python checks.py > checks-output.jsonl
Needs generate.py, independent.py and the campaign's unmodified check.py
alongside this file. This is build verification, not independent actor review.
"""
import json
import time
import os
from math import factorial
import numberdb.sage as numberdb
from sage.rings.rational_field import QQ
from sage.rings.real_arb import RealBallField
from sage.rings.complex_arb import ComplexBallField
from generate import HemisphereDirichlet,logdet
from independent import continuation
from check import exactness,measure,agrees_with,names_its_rings
def emit(x):
print(json.dumps(x),flush=True)
def require_agreement(x,y):
assert x.is_finite() and y.is_finite()
assert x.rad_as_ball()<x.parent()(10)**(-130)
assert y.rad_as_ball()<y.parent()(10)**(-130)
assert x.overlaps(y), (x,y)
if __name__=='__main__':
import sage.version,mpmath
emit({'versions':{'sage':sage.version.version,'numberdb':numberdb.__version__,'mpmath':mpmath.__version__},
'skill_manifest':'agents/runs/20261029T070407Z-build-skills/manifest.json',
'local_skill_sha256':'28ecbeeccb0c5e5d6eb541b5a5b0f88aec47b154374d909085116d34f5629199',
'live_skill_sha256':'28ecbeeccb0c5e5d6eb541b5a5b0f88aec47b154374d909085116d34f5629199'})
C=ComplexBallField(numberdb.bits(240));R=RealBallField(numberdb.bits(240))
# First a nonzero control for the independent method.
control,info=continuation(1)
require_agreement(control,(2*R.pi()).log())
emit({'control':'independent interval log determinant = log(2*pi)','passed':True})
# Different special-function identity tests all negative-even derivatives.
for n in range(1,10):
ref=QQ((-1)**n*factorial(2*n))/2*C(2*n+1).zeta()/(2*C.pi())**(2*n)
require_agreement(C(-2*n).zetaderiv(1).real(),ref.real())
emit({'control':'zeta functional equation at -2,-4,...,-18','passed':True})
values={};ind_values={};min_bits=10000
for d in range(1,21):
start=time.monotonic()
a=logdet(d,150);b=logdet(d,200)
require_agreement(a,b)
independent,info=continuation(d)
require_agreement(b,independent)
for name,x,y in [('det',a.exp(),independent.exp()),('logdet',a,independent)]:
values[(d,name)]=x;ind_values[(d,name)]=y
min_bits=min(min_bits,int(x.accuracy()))
# Explicit elementary low-dimensional relations, separate from the
# Stirling implementation, checked before inclusion in the document.
if d==2:
require_agreement(b,(-C(0).zetaderiv(1)*(-1)-2*C(-1).zetaderiv(1)+QQ(1)/4).real())
if d==3:
ref=-C(0).zetaderiv(1)+C(-1).zetaderiv(1)-C(-2).zetaderiv(1)-QQ(1)/2
require_agreement(b,ref.real())
emit({'d':d,'logdet_150':str(a),'det_150':str(a.exp()),
'logdet_200':str(b),'independent_logdet':str(independent),
'independent':info,'overlap':True,'seconds':round(time.monotonic()-start,3)})
assert not exactness(values)
assert not agrees_with(values,lambda k:ind_values[k])
assert not names_its_rings(os.path.join(os.path.dirname(__file__),'generate.py'))
for d in (2,20):
x,info=continuation(d,K=50)
require_agreement(x,ind_values[(d,'logdet')])
emit({'changed_split_control':d,'K':50,'passed':True})
emit({'status':'PASS','entries':40,'minimum_relative_accuracy_bits_150':min_bits,
'measure':measure(values),'exactness':[],'independent_agrees_with':[]})