Bill Allombert on Tue, 06 Jan 2015 12:17:05 +0100


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: Problems in 'make bench' under Cygwin/64 with testing elliptic


On Tue, Jan 06, 2015 at 10:56:32AM +0100, Andreas Kellers wrote:
> I did some further investigation using a failed test from 'Testing compat'.
> As far as I understand, the following statements should give 66:
> 
>   default(compatible, 3)
>   acurve = initell([0,0,1,-1,0])
>   apell(acurve,10007)
> 
> Doing this with the GP debug version (via './Configure -g') yields -1591
> instead.
> 
> Doing this with the GP standard version results in an error:
> 
>   ***   at top-level: apell(acurve,10007)
>   ***                 ^-------------------
>   *** apell: impossible inverse in Fl_inv: mod(11266880367774138369, 10007).

Both are wrong, the correct result is 66.
I assume you get the same result with
ellap(ellinit([0,0,1,-1,0]),10007)

This is a progress from the state where isprime(131) was wrong.

The problem is probably in the function "Fp_ellcard_Shanks" in
src/basemath/FpE.c.

I will investigate.

Thanks for all you have done so far on this!

Cheers,
Bill.