checked with gmplib mpz_probab_prime_p(k, 10)

To be completed...