# To run, load theta_data_22.sage, epsilon_data_22.sage, Multiple_Polylog_Good_Red_22.sage
# As the name suggests, we're currently limited to points of good reduction.
# Also, our implementation of 2-adic multiple zeta values is incomplete.
# Here are some sample commands:
word = (0,1,0,0,1)
point = 2
prime = 3
geom_prec = 11
padic_prec = 13
print(Li(word, point, prime, geom_prec, padic_prec))
word = (0,1,0,0,1)
prime = 3
padic_prec = 13
geom_prec = 11
print(zeta(word, prime, padic_prec, geom_prec))
# And here's one of many test we've run:
#Shuffle product tests
point = 3
prime = 5
geom_prec = 17
padic_prec = 12
def Li_a(word):
return Li(word, point, prime, geom_prec, padic_prec)
def zeta_p(word):
return zeta(word, prime, geom_prec, padic_prec)
print(3*Li_a((0,0,1,0,1))+6*Li_a((0,0,0,1,1))+Li_a((0,1,0,0,1))-Li_a((0,0,1))*Li_a((0,1)))