numba psums

npsum.py

from numba import jit

@jit(nopython=True)
def psum(p):
    acc = 0
    for i in range(1,1000000):
        acc += pow(i, -p)
    return acc

import time
t0=time.time()
ps = [psum(p) for p in range(2,20)]
print('DT:', round(time.time()-t0,2))

DT=0.7 vs 4.5 on data.quant.cl

even faster 2nd time around!