generate.sage

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

973 bytes, as of the version from 2026-08-15 11:00 (current). Recorded here, not run.

import yaml
import os
from utils.utils import numbers_to_yaml
from utils.utils import real_interval_to_sage_string
from utils.utils import blur_real_interval

path = 'data/Series_expansions/Stieltjes_constants/'

prec10 = 100 #relative precision in base 10
n_max = 1000
n_range = [0..n_max]

print("n_max:",n_max)

prec2 = prec10 * 3.4 + 20
RIFprec = RealIntervalField(prec2)
RBFprec = RealBallField(prec2)


'''
numbers = {n:
	RIFprec(stieltjes(n).n(prec2))
	for n in n_range
}
'''
numbers = {}
for n in n_range:
	print("n:",n)
	numbers[n] = blur_real_interval(RIFprec(stieltjes(n).n(prec2)))
	#Blurring is done as stieltjes(n).n() might return an integer,
	#which is then interpreted by RIFprec as exact, 
	#which leads to display issues.

numbers = {str(n):
	real_interval_to_sage_string(number,prec10).replace('?','')
	for n, number in numbers.items()
}

filename = os.path.join(path, 'numbers.yaml')
yaml.dump(numbers, stream = open(filename, 'w'), sort_keys = False)