# 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)))