z=np.zeros(10,np.int32) y=[gmpy.mpz[1] *10 ]
summation = gmp.xmpz(0) for i in range(0,10): product = gmpy2.xmpz(z[i]) product *= y[i] summation += product The bug prevents this from running. My workaround defeats the purpose of using gmpy2.xmpz typed variables : for i in range(0,10): product = gmpy2.xmpz(int(z[i]))