Comment #48 on issue 3809 by [email protected]: isprime can be faster
http://code.google.com/p/sympy/issues/detail?id=3809
Case, I tried r814 and it seems to work fine for me, both on the command
line and with sympy calling it. Using quotes now raises a TypeError, as do
most of the oddball error cases.
I was going through the prp calls, and am getting what I think are
incorrect results from is_extra_strong_lucas_prp. There are a lot more
pseudoprimes being output than there should be. It looks like the code and
comments have misread Grantham 2000 (the "Lucas pseudoprime" page on
Wikipedia is also good).
Grantham/MPU/Nicely: ( U_s=0 AND V_s = +/- 2 ) OR (V_{s*2^t} = 0 for t
in 0..r-1)
gmpy_mpz_prp.c: U_s=0 OR V_s = +/-2 OR V_{s*2^t} = 0 for t
in 0..r-1
It's in mpz_prp.c, so I'll drop David a note. This should have no impact
for GMPY2 outside of this one call.
--
You received this message because this project is configured to send all
issue notifications to this address.
You may adjust your notification preferences at:
https://code.google.com/hosting/settings
--
You received this message because you are subscribed to the Google Groups
"sympy-issues" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/sympy-issues.
For more options, visit https://groups.google.com/groups/opt_out.
