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