Bill James

I tested my latest gmpy2 bernoulli code with the 10 millionth bernoulli number last night. It completed in 6.5 hours, not bad. When mathematica broke the bernoulli record with this...

I uploaded the entire code directory that I test from. I am trying to create a working environment around gmpy2 that maximizes multi-processor use. So far, I have implemented multi-processor...

Yep, that looks like the issue. Strange number to set as a cutoff! I will fall back to a python factorial algorithm with some bit slicing. Thanks for the quick...

Looks like the issue runs deeper than just gamma(). The below code computing a mpz factorial fails as well, with result _inf._ 44787927 still works. Is there some trick to...

Looks like I have a hard limit on size of mpz that will convert to mpfr. The below code works. Change 2 ** 29 to 2 ** 30 and it...