check_program.py

back to table · edit · history · where entries came from · files · download

455 bytes, as of the version from 2026-11-02 01:06 (current). Recorded here, not run.

import numberdb.sage as numberdb
from pathlib import Path
import mpmath as mp,json
p=Path(__file__).parent
d=json.loads((p/'document.json').read_text());scope={}
exec(d['Programs']['program-exterior']['code'],scope)
rows=json.loads((p/'numerical-values.json').read_text())
with mp.workdps(210):
 for k in [13,14,15,16,-36]:
  err=abs(scope['mahler_q_exterior'](k)-mp.mpf(rows[str(k)]['value210']))
  assert err<mp.mpf('1e-115')
  print(k,mp.nstr(err,10))