Re: Frobenius test and GMP library
From: Cristiano (cristiano.pi_at_NSquipo.it)
Date: Fri, 24 Dec 2004 11:42:42 GMT
> [rob@localhost frobenius]$ ./frobenius.exe 10
> Performing 10 iterations of the frobenius probable prime test on
> Probable prime with parameters
I get the same!
I don't know whether you have the time, but it would be very helpful to try
you should change the lines
mpz_urandomm(a, rand_state, N);
mpz_urandomm(b, rand_state, N);
with: a=a+1 and b=b+1 which should be
mpz_add_ui(a, a, 1);
mpz_add_ui(b, b, 1);
Using that code, the primes 107, 32507 and 74471 are declared composites
(which seems very strange to me). Do you get the same?