cheng/mpir
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Fix mpz_nextprime so it runs 25 rounds of miller-rabin in total (same as GMP).

Browse Source
This commit is contained in:
William Hart 2014-03-26 12:39:08 +00:00
parent 01404cb523
commit 0d8b59ca28
1 changed files with 14 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
mpz/nextprime.c
View File

@ -25,11 +25,25 @@ Boston, MA 02110-1301, USA.
/* This function is Obsolete 17/8/2009 */ /* This function is Obsolete 17/8/2009 */
/*
But people use it anyway!
FIXME: This function should prove the primality of x using
ECPP or APR-CL.
*/
void mpz_nextprime(mpz_ptr x, mpz_srcptr y) void mpz_nextprime(mpz_ptr x, mpz_srcptr y)
{ {
gmp_randstate_t rnd; gmp_randstate_t rnd;
gmp_randinit_default(rnd); gmp_randinit_default(rnd);
mpz_next_prime_candidate(x, y, rnd); mpz_next_prime_candidate(x, y, rnd);
while (!mpz_miller_rabin (x, 23, rnd)) /* we've done 2 rounds already, do another 23 */
{
mpz_add_ui(x, x, 2);
mpz_next_prime_candidate(x, x, rnd);
}
gmp_randclear(rnd); gmp_randclear(rnd);
} }