# To run, load theta_data.sage, epsilon_data.sage, Multiple_Polylog_Good_Red.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's 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))