Paul Underwood on Fri, 27 Sep 2019 15:53:43 +0200 |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Fw: qfbsolve (and ispseudoprime) |
The correct expression is: gettime();k=198;n=2^k-3;if(qfbsolve(Qfb(1,0,1),n),print([k,gettime()])) > Sent: Friday, September 27, 2019 at 1:58 PM > From: "Paul Underwood" <paulunderwood@mindless.com> > To: pari-users@pari.math.u-bordeaux.fr > Subject: qfbsolve (and ispseudoprime) > > It is a pity that the function qfbsolve is much slower in 2.12 alpha where it solves for composite n too. > > For example: > > gettime();k=198;n=2^p-3;if(qfbsolve(Qfb(1,0,1),n),print([k,gettime()])) > > overflows the default stack whereas it gives an almost instantaneous answer with version 2.11 and version 2.9. > > Also, qfbsolve(Qfb(1,0,1),n) *was* quicker than ispseudoprime for positive integers of the form 4*K+1, > > Best > > Paul