Tested with mpz_probab_prime(p,10) from GMP lib

To be completed...