mini-gmp: mpz_congruent_p and mpz_probab_prime_p
bodrato at mail.dm.unipi.it
bodrato at mail.dm.unipi.it
Thu Mar 6 17:50:43 UTC 2014
Ciao,
Il Mer, 5 Marzo 2014 6:07 pm, Niels Möller ha scritto:
>> May I suggest to change the second comparison to mpz_cmp_ui (y, 1) <= 0?
>
> Hmm, why? y == 0 can happen if n divides a power of a (the base for the
> test), in which case n is composite (except in corner cases, see below),
> and in which case the test will eventually fail without explicitly
> checking for y == 0.
I agree with your analysis: with the suggested change the same answer is
given. Earlier, at no extra cost, I say.
The same comment applies to the current implementation in the main library.
Regards,
m
--
http://bodrato.it/
More information about the gmp-devel
mailing list