bodr...@mail.dm.unipi.it wrote: > IIRC someone proposed to change the meaning of the second argument of > mpz_probab_prime_p. Now it is documented as "reps controls how many > [Miller-Rabin] tests are done". It should be something like "logprob > controls the likeliness of a wrong answer, the function will return > probably prime for a composite with a probability of 1 over 2^logprob".
What I think is unacceptable is that if a composite passes the test, it also passes the test when invoked a second time because the bases it's tested against are always the same. That's the why of a version that accepts a random state. So, a new function is recommendable in any case. That applies to mpz_millerrabin as well, of course, which is the one actually doing the PRNG calls. And while on that subject, there was also a request for a M-R test function accepting a specific base as parameter: http://gmplib.org/list-archives/gmp-devel/2002-December/000075.html And a suggestion to return the witness that proved the compositeness: http://gmplib.org/list-archives/gmp-devel/2008-January/000766.html In that message, Torbjörn also says that it'd be nice for a function called millerrabin to do a M-R test only, not also a Fermat test. _______________________________________________ gmp-devel mailing list gmp-devel@gmplib.org http://gmplib.org/mailman/listinfo/gmp-devel