Tested with mpz_probab_prime_p(p, 10) from GMP lib

To be completed...