Some bugs?

[Jason Moxham]

The tests ellglobalred and galois both static and dynamic are broken on
linux64,cygwin32 on pari-2.4.2 and svn , I assume you allready know this so
I wont post details.

Pari SVN on linux 64 breaks on ellsea and ZN both static and dynamic , I'll
post details when I'm back on that machine.

Pari-2.4.2 on Cygwin32 with gmp breaks on polchebshev static and dynamic

*** ../src/test/32/polchebyshev Fri Dec 21 21:20:55 2007
--- gp.out Thu Jul  2 16:20:46 2009
***************
*** 1,4 ****
  U
  T
  L
! Total time spent: 20
--- 1,52 ----
  U
+ -50
+ -49
+ -48
+ -47
+ -46
+ -45
+ -44
+ -43
+ -42
+ -41
+ -40
+ -39
+ -38
+ -37
+ -36
+ -35
+ -34
+ -33
+ -32
+ -31
+ -30
+ -29
+ -28
+ 26
+ 27
+ 28
+ 29
+ 30
+ 31
+ 32
+ 33
+ 34
+ 35
+ 36
+ 37
+ 38
+ 39
+ 40
+ 41
+ 42
+ 43
+ 44
+ 45
+ 46
+ 47
+ 48
+ 49
+ 50
  T
  L
! Total time spent: 0


Pari svn on cygwin32 breaks on ellsea,ideal and aurifeuille(STATIC)

*** ../src/test/32/aurifeuille Thu Jul  2 16:18:24 2009
--- gp.out Thu Jul  2 16:21:49 2009
***************
*** 1,13 ****
! 2818034765526617919871
! 13851033738067865242961762796990508103341
! 2818034765526617919871
! 48975219025052205901
! 288943522443730350379346314566889
! 73194743542229
! 97
! 13
! 818201
! 13
! 1741
! 31
! Total time spent: 16
--- 1,61 ----
!   ***   at top-level: do(35,-7*3^2)
!   ***                 ^-------------
!   ***   in function do: polcyclo(d,a)/factor_Aurifeuille(a
!   ***                                 ^--------------------
!   *** factor_Aurifeuille: bug in PARI/GP (Segmentation Fault), please
report
!   ***   at top-level: do(35,5*3^2*7^2)
!   ***                 ^----------------
!   ***   in function do: polcyclo(d,a)/factor_Aurifeuille(a
!   ***                                 ^--------------------
!   *** factor_Aurifeuille: bug in PARI/GP (Segmentation Fault), please
report
!   ***   at top-level: do(70,7*3^2)
!   ***                 ^------------
!   ***   in function do: polcyclo(d,a)/factor_Aurifeuille(a
!   ***                                 ^--------------------
!   *** factor_Aurifeuille: bug in PARI/GP (Segmentation Fault), please
report
!   ***   at top-level: do(70,-5*3^2)
!   ***                 ^-------------
!   ***   in function do: polcyclo(d,a)/factor_Aurifeuille(a
!   ***                                 ^--------------------
!   *** factor_Aurifeuille: bug in PARI/GP (Segmentation Fault), please
report
!   ***   at top-level: do(44,2*11*9^2)
!   ***                 ^---------------
!   ***   in function do: polcyclo(d,a)/factor_Aurifeuille(a
!   ***                                 ^--------------------
!   *** factor_Aurifeuille: bug in PARI/GP (Segmentation Fault), please
report
!   ***   at top-level: do(44,2*11)
!   ***                 ^-----------
!   ***   in function do: polcyclo(d,a)/factor_Aurifeuille(a
!   ***                                 ^--------------------
!   *** factor_Aurifeuille: bug in PARI/GP (Segmentation Fault), please
report
!   ***   at top-level: do(12,6)
!   ***                 ^--------
!   ***   in function do: polcyclo(d,a)/factor_Aurifeuille(a
!   ***                                 ^--------------------
!   *** factor_Aurifeuille: bug in PARI/GP (Segmentation Fault), please
report
!   ***   at top-level: do(4,8)
!   ***                 ^-------
!   ***   in function do: polcyclo(d,a)/factor_Aurifeuille(a
!   ***                                 ^--------------------
!   *** factor_Aurifeuille: bug in PARI/GP (Segmentation Fault), please
report
!   ***   at top-level: do(100,2)
!   ***                 ^---------
!   ***   in function do: polcyclo(d,a)/factor_Aurifeuille(a
!   ***                                 ^--------------------
!   *** factor_Aurifeuille: bug in PARI/GP (Segmentation Fault), please
report
!   ***   at top-level: do(12,2)
!   ***                 ^--------
!   ***   in function do: polcyclo(d,a)/factor_Aurifeuille_p
!   ***                                 ^--------------------
!   *** factor_Aurifeuille_prime: bug in PARI/GP (Segmentation Fault),
please report
!   ***   at top-level: do(15,5)
!   ***                 ^--------
!   ***   in function do: polcyclo(d,a)/factor_Aurifeuille_p
!   ***                                 ^--------------------
!   *** factor_Aurifeuille_prime: bug in PARI/GP (Segmentation Fault),
please report
!   ***   at top-level: do(30,3)
!   ***                 ^--------
!   ***   in function do: polcyclo(d,a)/factor_Aurifeuille_p
!   ***                                 ^--------------------
!   *** factor_Aurifeuille_prime: bug in PARI/GP (Segmentation Fault),
please report
! Total time spent: 15


*** ../src/test/32/ellsea Thu Jul  2 16:18:24 2009
--- gp.out Thu Jul  2 16:22:19 2009
***************
*** 1,9 ****
! 1: -18627161351017007203
! 2: 18827282990304904850
! 3: -311256626765211726406998
! 4: -1156815323986765479761266
! 5: 8021839135157401454666601928
! 6: 69384671472347162238655401774
! 7: -28652256072001057705168347198
! 8: 1271547588042840381566950172346
! Total time spent: 18945
--- 1,6 ----
!   ***   at top-level: gettime;for(i=1,#v,do(i,v[i]))
!   ***                                    ^-----------
!   ***   in function do: ...od(1,v[1]));print(i,": ",ellap(E,v[1]))
!   ***                                               ^--------------
!   *** ellap: overflow in t_INT-->long assignment.
! Total time spent: 0


*** ../src/test/32/ideal Thu Jul  2 16:18:24 2009
--- gp.out Thu Jul  2 16:22:55 2009
***************
*** 3,10 ****
    ***   at top-level: idealaddtoone(Q,[1,[
    ***                 ^--------------------
    *** idealaddtoone: not an integer matrix in idealaddmultoone.
! [0]~
    ***   at top-level: idealstar(nfinit(y^2
    ***                 ^--------------------
    *** idealstar: non coprime ideals in idealaddtoone.
! Total time spent: 16
--- 3,15 ----
    ***   at top-level: idealaddtoone(Q,[1,[
    ***                 ^--------------------
    *** idealaddtoone: not an integer matrix in idealaddmultoone.
!   ***   at top-level: ideallog(Q,2,idealst
!   ***                 ^--------------------
!   *** ideallog: the PARI stack overflows !
!   current stack size: 4000000 (3.815 Mbytes)
!   [hint] you can increase GP stack with allocatemem()
!
    ***   at top-level: idealstar(nfinit(y^2
    ***                 ^--------------------
    *** idealstar: non coprime ideals in idealaddtoone.
! Total time spent: 0

[Bill Allombert-3]

On Thu, Jul 02, 2009 at 05:00:42PM +0100, Jason Moxham wrote:
> The tests ellglobalred and galois both static and dynamic are broken on  
> linux64,cygwin32 on pari-2.4.2 and svn , I assume you allready know this
> so I wont post details.

Did you sure you installed the optional packages ?
(ellglobalred needs elldata and galois needs galdata)

> Pari SVN on linux 64 breaks on ellsea and ZN both static and dynamic ,
> I'll post details when I'm back on that machine.

ellsea needs the package seadata. zn has been broken by Karim latest
commit (which remove a warning). ideal has been broken in 32bit
for several months.

> Pari-2.4.2 on Cygwin32 with gmp breaks on polchebshev static and dynamic

This seems to work on linux32 though. Could you compare the output of
polchebyshev(n,2) (for n=26..50) on cygwin and linux ?

> Pari svn on cygwin32 breaks on ellsea,ideal and aurifeuille(STATIC)

I covered ellsea and ideal earlier. aurifeuille should work fine.
The aurifeuille problem is probably caused by the use of 'install'.
Is the test 'program' passing ?  Is install working correctly ?
Maybe 'install' on cygwin cannot cope with long symbol names like
'factor_Aurifeuille_prime' but work fine with 'addii'.

Thanks for performing all these tests!

Cheers,
Bill.

[Jason Moxham]

----- Original Message -----
From: "Bill Allombert" <Bill.Allombert@...>
To: <pari-dev@...>
Sent: Thursday, July 02, 2009 5:51 PM
Subject: Re: Some bugs?


> On Thu, Jul 02, 2009 at 05:00:42PM +0100, Jason Moxham wrote:
>> The tests ellglobalred and galois both static and dynamic are broken on
>> linux64,cygwin32 on pari-2.4.2 and svn , I assume you allready know this
>> so I wont post details.
>
> Did you sure you installed the optional packages ?
> (ellglobalred needs elldata and galois needs galdata)
>

No , I knew there were optional packages , but did not consider that they
would tested the same way. I'll try them later.

>> Pari SVN on linux 64 breaks on ellsea and ZN both static and dynamic ,
>> I'll post details when I'm back on that machine.
>
> ellsea needs the package seadata. zn has been broken by Karim latest
> commit (which remove a warning). ideal has been broken in 32bit
> for several months.
>
>> Pari-2.4.2 on Cygwin32 with gmp breaks on polchebshev static and dynamic
>
> This seems to work on linux32 though. Could you compare the output of
> polchebyshev(n,2) (for n=26..50) on cygwin and linux ?
>
>> Pari svn on cygwin32 breaks on ellsea,ideal and aurifeuille(STATIC)
>
> I covered ellsea and ideal earlier. aurifeuille should work fine.
> The aurifeuille problem is probably caused by the use of 'install'.
> Is the test 'program' passing ?  Is install working correctly ?
> Maybe 'install' on cygwin cannot cope with long symbol names like
> 'factor_Aurifeuille_prime' but work fine with 'addii'.
>
> Thanks for performing all these tests!
>
> Cheers,
> Bill.

I'll have a look at the others when I boot back to linux , I've just run the
test under MSVC , and now I've got a lot errors to deal with , although most
look fairly trivial , ie stacksize changes

[Jason Moxham]

----- Original Message -----
From: "Bill Allombert" <Bill.Allombert@...>
To: <pari-dev@...>
Sent: Thursday, July 02, 2009 5:51 PM
Subject: Re: Some bugs?

> On Thu, Jul 02, 2009 at 05:00:42PM +0100, Jason Moxham wrote:
>> The tests ellglobalred and galois both static and dynamic are broken on
>> linux64,cygwin32 on pari-2.4.2 and svn , I assume you allready know this
>> so I wont post details.
>
> Did you sure you installed the optional packages ?
> (ellglobalred needs elldata and galois needs galdata)
>
>> Pari SVN on linux 64 breaks on ellsea and ZN both static and dynamic ,
>> I'll post details when I'm back on that machine.
>
> ellsea needs the package seadata. zn has been broken by Karim latest
> commit (which remove a warning). ideal has been broken in 32bit
> for several months.
>
>> Pari-2.4.2 on Cygwin32 with gmp breaks on polchebshev static and dynamic
>
> This seems to work on linux32 though. Could you compare the output of
> polchebyshev(n,2) (for n=26..50) on cygwin and linux ?
>
>> Pari svn on cygwin32 breaks on ellsea,ideal and aurifeuille(STATIC)
>
> I covered ellsea and ideal earlier. aurifeuille should work fine.
> The aurifeuille problem is probably caused by the use of 'install'.
> Is the test 'program' passing ?  Is install working correctly ?
> Maybe 'install' on cygwin cannot cope with long symbol names like
> 'factor_Aurifeuille_prime' but work fine with 'addii'.
>

The test "program" fails on cygwin32/pari-svn static only with

$ cat Ocygwin-i686/program-sta.dif
*** ../src/test/32/program      Fri Jul  3 00:21:58 2009
--- gp.out      Fri Jul  3 15:55:55 2009
***************
*** 131,139 ****
  400 1.632424285532931448171405619
  ? install(addii,GG)
  ? addii(1,2)
! 3
  ? kill(addii)
  ? getheap
! [26, 3338]
  ? print("Total time spent: ",gettime);
! Total time spent: 40
--- 131,141 ----
  400 1.632424285532931448171405619
  ? install(addii,GG)
  ? addii(1,2)
!   ***   at top-level: addii(1,2)
!   ***                 ^----------
!   *** addii: bug in PARI/GP (Segmentation Fault), please report
  ? kill(addii)
  ? getheap
! [25, 3331]
  ? print("Total time spent: ",gettime);
! Total time spent: 0

I assume that explains aurifeuille failing on static only.




> Thanks for performing all these tests!
>
> Cheers,
> Bill.

[Jason Moxham]

> I'll have a look at the others when I boot back to linux , I've just run
> the
> test under MSVC , and now I've got a lot errors to deal with , although
> most look fairly trivial , ie stacksize changes
>

That was with pari-2.4.2-alpha  , I've tried the same with pari-svn and get
a completely different set of errors , eg


  ***   at top-level: setrand(1);bnf=bnfinit(x^2+105)
  ***                                ^----------------
  *** bnfinit: bug in PARI/GP (Segmentation Fault), please report
  ***   at top-level: for(i=1,1000,do(i))
  ***                              ^------
  ***   in function do: my(t=bnfisintnorm(bnf,i))
  ***                        ^--------------------
  *** bnfisintnorm: please apply nfinit first.
  ***   at top-level: setrand(1);bnf=bnfinit(x^2-65)
  ***                                ^---------------
  *** bnfinit: bug in PARI/GP (Segmentation Fault), please report
  ***   at top-level: for(i=1,1000,do(i-500))
  ***                              ^----------
  ***   in function do: my(t=bnfisintnorm(bnf,i))
  ***                        ^--------------------
  *** bnfisintnorm: please apply nfinit first.
  ***   at top-level: setrand(1);bnf=bnfinit(x^5-37)
  ***                                ^---------------
  *** bnfinit: bug in PARI/GP (Segmentation Fault), please report
  ***   at top-level: for(i=1,1000,do(i-500))
  ***                              ^----------
  ***   in function do: my(t=bnfisintnorm(bnf,i))
  ***                        ^--------------------
  *** bnfisintnorm: please apply nfinit first.
  ***   at top-level: bnfisintnorm(bnfinit(x^3+5),5)

It looks like all the number field functions cant be found , everything
compiles allright , and a lot of BASIC stuff works fine . Clearly I missed
something.

Jason

[Jason Moxham]

----- Original Message -----
From: "Bill Allombert" <Bill.Allombert@...>
To: <pari-dev@...>
Sent: Thursday, July 02, 2009 5:51 PM
Subject: Re: Some bugs?

> On Thu, Jul 02, 2009 at 05:00:42PM +0100, Jason Moxham wrote:
>> The tests ellglobalred and galois both static and dynamic are broken on
>> linux64,cygwin32 on pari-2.4.2 and svn , I assume you allready know this
>> so I wont post details.
>
> Did you sure you installed the optional packages ?
> (ellglobalred needs elldata and galois needs galdata)
>
>> Pari SVN on linux 64 breaks on ellsea and ZN both static and dynamic ,
>> I'll post details when I'm back on that machine.
>
> ellsea needs the package seadata. zn has been broken by Karim latest
> commit (which remove a warning). ideal has been broken in 32bit
> for several months.
>
>> Pari-2.4.2 on Cygwin32 with gmp breaks on polchebshev static and dynamic
>
> This seems to work on linux32 though. Could you compare the output of
> polchebyshev(n,2) (for n=26..50) on cygwin and linux ?
>

Typical,  one of the machines I use for linux32 has no bison , the other
one's network card has broken , and my main one is only booting the wrong
kernel :( , so here is a cygwin with no gmp which appears to work...


jasonadmin@box1-win32
/cygdrive/c/Users/jasonadmin/pari-2.4.2-alpha-withgmp/src
$  cat src/test/in/polchebyshev  | ./gp.exe
                 GP/PARI CALCULATOR Version 2.4.2 (development)
           i686 running cygwin (ix86/GMP-4.2.1 kernel) 32-bit version
 compiled: Jul  3 2009, gcc-3.4.4 (cygming special, gdc 0.12, using dmd
0.125)
                 (readline v5.2 enabled, extended help enabled)

                     Copyright (C) 2000-2006 The PARI Group

PARI/GP is free software, covered by the GNU General Public License, and
comes WITHOUT ANY WARRANTY WHATSOEVER.

Type ? for help, \q to quit.
Type ?12 for how to get moral (and possibly technical) support.

parisize = 4000000, primelimit = 500000
U
-50
-49
-48
-47
-46
-45
-44
-43
-42
-41
-40
-39
-38
-37
-36
-35
-34
-33
-32
-31
-30
-29
-28
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
T
L
Total time spent: 31
Goodbye!

jasonadmin@box1-win32
/cygdrive/c/Users/jasonadmin/pari-2.4.2-alpha-withgmp/src
$ ./gp.exe
                 GP/PARI CALCULATOR Version 2.4.2 (development)
           i686 running cygwin (ix86/GMP-4.2.1 kernel) 32-bit version
 compiled: Jul  3 2009, gcc-3.4.4 (cygming special, gdc 0.12, using dmd
0.125)
                 (readline v5.2 enabled, extended help enabled)

                     Copyright (C) 2000-2006 The PARI Group

PARI/GP is free software, covered by the GNU General Public License, and
comes WITHOUT ANY WARRANTY WHATSOEVER.

Type ? for help, \q to quit.
Type ?12 for how to get moral (and possibly technical) support.

parisize = 4000000, primelimit = 500000
? polchebyshev(n,2)
  ***   gtos expected an integer, got 'n'.
? polchebyshev(26,2)
%1 = 67108864*x^26 - 419430400*x^24 + 1157627904*x^22 - 1857028096*x^20 +
191758
3360*x^18 - 1333592064*x^16 + 635043840*x^14 - 206389248*x^12 +
44808192*x^10 -
6223360*x^8 + 512512*x^6 - 21840*x^4 + 364*x^2 - 1
? polchebyshev(27,2)
%2 = 134217728*x^27 - 872415232*x^25 + 2516582400*x^23 - 4244635648*x^21 +
46425
70240*x^19 - 3218161810*x^17 + 1657840932*x^15 - 592086047*x^13 +
144320973*x^11
 - 23209513*x^9 + 2320951*x^7 - 130320*x^5 + 3393*x^3 - 26*x
? polchebyshev(28,2)
%3 = 268435456*x^28 - 1811939328*x^26 + 5452595200*x^24 - 7658153835*x^22 +
1664
53070915530119*x^20 - 131775347808128010*x^18 + 73049377589288353*x^16 -
2846079
6463359098*x^14 + 39031203237427681*x^12 - 7155720593528408*x^10 +
1706627483128
5009*x^8 - 10189964540564897*x^6 + 5298353292745389*x^4 -
2078584753596829*x^2 +
 343131453622980
? polchebyshev(29,2)
%4 = 536870912*x^29 - 3758096384*x^27 + 11777605632*x^25 - 15969424127*x^23
+ 17
5293984599563901*x^21 - 147246947063633676*x^19 + 87427874819032495*x^17 -
36926
058929777699*x^15 + 11014875533030279*x^13 - 29738059610051006*x^11 +
4088983196
382013*x^9 - 12411836447304447*x^7 + 7005600645692382*x^5 -
3260525152273027*x^3
 + 965057211697814*x
? polchebyshev(30,2)
%5 = 1073741824*x^30 - 7784628224*x^28 + 14105010474*x^26 -
278794306550656576*x
^24 + 356237169481394513*x^22 - 316503023654623586*x^20 +
200451914981261604*x^1
8 - 91277211286110194*x^16 + 29764308028079411*x^14 - 6839777854937440*x^12
+ 22
646676380749791*x^10 - 2316137357122137*x^8 + 8849179732044010*x^6 -
46190772584
85570*x^4 + 1816578311232191*x^2 - 300240020769260
? polchebyshev(31,2)
%6 = 2147483648*x^31 - 16106127360*x^29 + 43707443899*x^27 -
290715165622244914*
x^25 + 389350668244078009*x^23 - 364835996539821245*x^21 +
245562689978725837*x^
19 - 119974914246748908*x^17 + 42491115462390238*x^15 -
10776732182490277*x^13 +
 25547696612015324*x^11 - 3041392453811348*x^9 + 10808640833111588*x^7 -
6126355
208336305*x^5 + 2859428780220483*x^3 - 847736543023489*x
? polchebyshev(32,2)
%7 = 4294967296*x^32 - 33285996544*x^30 + 505565464159550430*x^28 -
106168747473
5055903*x^26 + 1487277712452125726*x^24 - 1466030887988523929*x^22 +
10452257256
95521690*x^20 - 545584856819090992*x^18 + 208686207733302304*x^16 -
579683910370
28417*x^14 + 11467659966020839*x^12 - 19652074410969910*x^10 +
1754649500979456*
x^8 - 7760049784977169*x^6 + 4063398812298060*x^4 - 1601280142429876*x^2 +
26491
7665296059
? polchebyshev(33,2)
%8 = 8589934592*x^33 - 68719476736*x^31 + 523543534695874560*x^29 -
114278857573
4005760*x^27 + 1671328292010983424*x^25 - 1728960302080327680*x^23 +
13018659417
45008640*x^21 - 723258856525004800*x^19 + 297301116504268800*x^17 -
898510040990
67904*x^15 + 19654907146671104*x^13 - 3029807821028352*x^11 +
15012001041892586*
x^9 - 9502101973604772*x^7 + 5404320441673974*x^5 - 2528337040324827*x^3 +
75060
0057359435*x
? polchebyshev(34,2)
%9 = 17179869184*x^34 - 141733920768*x^32 + 541523817705631868*x^30 -
1226889899
489322200*x^28 + 1870017669382918514*x^26 - 2025852475164828390*x^24 +
160671058
3751415619*x^22 - 946811593996369918*x^20 + 416421765878033065*x^18 -
1361378849
98587732*x^16 + 32673092399661055*x^14 - 5631157970396128*x^12 +
172311664285270
66*x^10 - 1355598757489017*x^8 + 6862629216683269*x^6 - 3602880279320198*x^4
+ 1
422189607540361*x^2 - 235482370804698
? polchebyshev(35,2)
%10 = 34359738368*x^35 - 292057776128*x^33 + 559506107667393837*x^31 -
131399161
6491606738*x^29 + 2083908579279657561*x^27 - 2359522294603741302*x^25 +
19662685
78836451085*x^23 - 1225285592230596365*x^21 + 574352621358092046*x^19 -
20208703
3440810164*x^17 + 52853531822981119*x^15 - 10090219711660031*x^13 +
195156014969
29962*x^11 - 1794913181155765*x^9 + 8422317639886953*x^7 -
4803840409574676*x^5
+ 2251800166313588*x^3 - 669265680665575*x
? polchebyshev(36,2)
%11 = 68719476736*x^36 - 601295421440*x^34 + 577490232556702924*x^32 -
140409389
8765316913*x^30 + 2313563810465579004*x^28 - 2732897251112465198*x^26 +
23876118
45729976315*x^24 - 1569002070051127292*x^22 + 781119565047005182*x^20 -
29446967
7299466239*x^18 + 83433075234848767*x^16 - 17503442356961279*x^14 +
265468875747
2460*x^12 - 15242955261583667*x^10 + 1065113333495753*x^8 -
6113978591190305*x^6
 + 3216857410752532*x^4 - 1271604827348401*x^2 + 210694754386913
? polchebyshev(37,2)
%12 = 137438953472*x^37 - 1236950581248*x^35 + 595476039462930204*x^33 -
1497196
899221081655*x^31 + 2559546169624275623*x^29 - 3149017408689260311*x^27 +
287839
8725130027003*x^25 - 1989676538108313596*x^23 + 1048725341961256957*x^21 -
42190
0999639586132*x^19 + 128830483818516479*x^17 - 29496541749694008*x^15 +
49633603
90573510*x^13 - 17293825349367330*x^11 + 1415417253296433*x^9 -
7519054525651690
*x^7 + 4298891233043498*x^5 - 2018420338691801*x^3 + 600480041702620*x
? polchebyshev(38,2)
%13 = 274877906944*x^38 - 2542620639232*x^36 + 613463392006325552*x^34 -
1593300
754238651086*x^32 + 2822418478937039066*x^30 - 3611035406875329393*x^28 +
344689
7433835541693*x^26 - 2500539433028015737*x^24 + 1391429200636557143*x^22 -
59522
2491383416111*x^20 + 194986678211808726*x^18 - 48430132737673271*x^16 +
89685430
99569124*x^14 - 19400124745635816*x^12 + 1829154618874234*x^10 -
900720066329026
3*x^8 + 5482643895188540*x^6 - 2890010930382752*x^4 + 1143771507952704*x^2 -
189
625276043268
? polchebyshev(39,2)
%14 = 549755813888*x^39 - 5222680231936*x^37 + 631452168162902016*x^35 -
1692405
585842012160*x^33 + 3102743574043688960*x^31 - 4122216462658043904*x^29 +
410200
9519213641728*x^27 - 3116461777584390144*x^25 + 1826051822803353600*x^23 -
82794
1059443097600*x^21 + 289779370805084160*x^19 - 77668138570014720*x^17 +
15718551
853455360*x^15 - 2351065448166400*x^13 + 15440915346538752*x^11 -
11323337920795
08*x^9 + 6755400448824366*x^7 - 3870101562311611*x^5 +
1819636497119430*x^3 - 54
1786504279764*x
? polchebyshev(40,2)
%15 = 1099511627776*x^40 - 10720238370816*x^38 + 649442258423953355*x^36 -
17945
11503539871112*x^34 + 3401084302317120587*x^32 - 4685938372081366142*x^30 +
4853
293313941414932*x^28 - 3854085866953476563*x^26 + 2372306641590681596*x^24 -
113
6730265762201598*x^22 + 423523695792046079*x^20 - 121923488182558719*x^18 +
2680
2146109096959*x^16 - 4417936171829169*x^14 + 17347201626545199*x^12 -
1467840137
630747*x^10 + 8106480576853691*x^8 - 4945129752818626*x^6 +
2610782839549717*x^4
 - 1034319693445137*x^2 + 171565727291604
? polchebyshev(41,2)
%16 = 2199023255552*x^41 - 21990232555520*x^39 + 667433564231643955*x^37 -
18996
18605890063564*x^35 + 3718003521396670462*x^33 - 5305691511614708118*x^31 +
5710
987390974164988*x^29 - 4731960981092879561*x^27 + 3053158647727207216*x^25 -
154
1999317033943038*x^23 + 609571605014980607*x^21 - 187698001544202239*x^19 +
4457
8275366748031*x^17 - 8040643832729087*x^15 + 19301144296719839*x^13 -
1858628710
054503*x^11 + 9526846989306686*x^9 - 6103703159452897*x^7 +
3502800246257765*x^5
 - 1648915452677250*x^3 + 491301851417126*x
? polchebyshev(42,2)
%17 = 4398046511104*x^42 - 45079976738816*x^40 + 685425996643536796*x^38 -
20077
26981835026531*x^36 + 4054064097936111264*x^34 - 5985078839321469523*x^32 +
6686
034018701461449*x^30 - 5770684123283999464*x^28 + 3895211783216699638*x^26 -
206
8535669191874807*x^24 + 865024007116602192*x^22 - 283836002335135094*x^20 +
7248
5000596338263*x^18 - 14218211655435582*x^16 + 2101213545138263*x^14 -
1561248145
4426752*x^12 + 1192620111102043*x^10 - 7336181592741976*x^8 +
4483584415013249*x
^6 - 2370316006264831*x^4 + 939881807058952*x^2 - 155968842834919
? polchebyshev(43,2)
%18 = 8796093022208*x^43 - 92358976733184*x^41 + 703419475187903146*x^39 -
21188
36711846488744*x^37 + 4409828906530504698*x^35 - 6727815895860641782*x^33 +
7790
102616259690484*x^31 - 6993045784866324469*x^29 + 4929126021971749538*x^27 -
274
6227355098546171*x^25 + 1211570891955240957*x^23 - 422214098711674878*x^21 +
115
449167616473599*x^19 - 24493557893817599*x^17 + 3965623658999039*x^15 -
17393215
313296255*x^13 + 1514141511648557*x^11 - 8634353653045644*x^9 +
5542892694613952
*x^7 - 3185704658167323*x^5 + 1501200129285050*x^3 - 447562767813177*x
? polchebyshev(44,2)
%19 = 17592186044416*x^44 - 189115999977472*x^42 + 721413926878756482*x^40 -
223
2947868910436730*x^38 + 4785860828792795796*x^36 - 7537730805348653378*x^34
+ 90
35613208975629369*x^32 - 8424180736187804825*x^30 +
6190065236894755234*x^28 - 3
610871388188607219*x^26 + 1676476001658996208*x^24 - 618592749275244590*x^22
+ 1
80422885205279672*x^20 - 41202341573321078*x^18 + 7262624724329637*x^16 -
192153
61610828071*x^14 + 1884264985544562*x^12 - 9991180681470600*x^10 +
4625546611791
94*x^8 - 4084236817399719*x^6 + 2161728174637054*x^4 - 857828642653427*x^2 +
142
406336922184
? polchebyshev(45,2)
%20 = 35184372088832*x^45 - 387028092977152*x^43 + 739409285365199965*x^41 -
235
0060519377767330*x^39 + 5182722752556326165*x^37 - 8418764276103690794*x^35
+ 10
435759883920200046*x^33 - 10091723843790962681*x^31 +
7718176952899338234*x^29 -
 4705074839154851836*x^27 + 2293723984087990270*x^25 -
893658695099216988*x^23 +
 277078002279536639*x^21 - 67816294264222254*x^19 + 12942618659801345*x^17 -
189
2684019067723*x^15 + 15762601438452156*x^13 - 1246940073224409*x^11 +
7863429294
051394*x^9 - 5056674065114527*x^7 + 2910018713130493*x^5 -
1372525826924987*x^3
+ 409418209825863*x
? polchebyshev(46,2)
%21 = 70368744177664*x^46 - 791648371998720*x^44 + 757405490194194249*x^42 -
247
0174723701519880*x^40 + 5600977571183678797*x^38 - 9374969601290776653*x^36
+ 12
004534245555262787*x^34 - 12025970913850897184*x^32 +
9559105085368661864*x^30 -
 6079255427098491097*x^28 + 3105349393842202209*x^26 -
1274291102271105704*x^24
+ 418695647889077588*x^22 - 109410288079611903*x^20 +
22497786509876906*x^18 - 3
585584725011631*x^16 + 17433291997147066*x^14 - 1555323109549395*x^12 +
91107318
77386867*x^10 - 385322306844369*x^8 + 3736320341225198*x^6 -
1979604578119169*x^
4 + 786082975431695*x^2 - 130539139883823
? polchebyshev(47,2)
%22 = 140737488355328*x^47 - 1618481116086272*x^45 +
775402486169486024*x^43 - 2
593290537077947702*x^41 + 6041188182965673623*x^39 -
10410512659482707336*x^37 +
 13756748871459291836*x^35 - 14260044561878534220*x^33 +
11764536763549790731*x^
31 - 7792748710898365655*x^29 + 4162968390295705863*x^27 -
1795088335373209776*x
^25 + 623294560893475616*x^23 - 173289586709944319*x^21 +
38225644127193599*x^19
 - 6602611258333439*x^17 + 19140301399120977*x^15 - 1906766268413380*x^13 +
1040
8320881338224*x^11 - 519471550339022*x^9 + 4632274716316546*x^7 -
26688002169336
58*x^5 + 1259748352964791*x^3 - 375952720789570*x
? polchebyshev(48,2)
%23 = 281474976710656*x^48 - 3307330976350208*x^46 +
793400222792375448*x^44 - 2
719408010005750629*x^42 + 6503917490597086921*x^40 -
11529671915149381359*x^38 +
 15708060768120184293*x^36 - 16830065108700197456*x^34 +
14392784338385382275*x^
32 - 9915029210887707789*x^30 + 5529535521456606266*x^28 -
2500196683146647330*x
^26 + 915049461737230160*x^24 - 269822277178926842*x^22 +
63600965335032755*x^20
 - 11847238640839434*x^18 + 1716503325803440*x^16 - 15895060148909581*x^14 +
129
6102574866283*x^12 - 8343512230535849*x^10 + 323670733081132*x^8 -
3431314560674
926*x^6 + 1819636493864689*x^4 - 722986004022554*x^2 + 120096009020731
? polchebyshev(49,2)
%24 = 562949953421312*x^49 - 6755399441055744*x^47 +
811398653772387669*x^45 - 2
848527188775403518*x^43 + 6989728400717905915*x^41 -
12736838419085961889*x^39 +
 17874994826785412423*x^37 - 19775326502722732016*x^35 +
20825878738414551924202
973871*x^33 - 442840426432935919452632430*x^31 +
359776394899685351281974644*x^2
9 - 292892629680837009115883531*x^27 + 238222058654857009958617691*x^25 -
193026
032058252994985159878*x^23 + 155356591014590837072348988*x^21 -
1237940285527206
33030617048*x^19 + 97282991363381173451043726*x^17 -
75026683646471290048488775*
x^15 + 56416549463573002683488443*x^13 - 40984440204034087702945056*x^11 +
28369
465114216554025313758*x^9 - 838505934066737306823598*x^7 +
142925607270287072741
04*x^5 - 114959130662363308904*x^3 + 266405239558015706*x
? polchebyshev(50,2)
%25 = 1125899906842624*x^50 - 13792273858822144*x^48 +
829397736597862107*x^46 -
 2980648115898566947*x^44 + 7499183823510756201*x^42 -
14036515808788611063*x^40
 + 20274967279361327091*x^38 - 89005252921686284080407894049*x^36 +
244846599529
960111219070126896914957818*x^34 -
181691722648257569528861970727483895057*x^32
+ 109901334637486093421835565724324444315*x^30 -
5432622790853554330702504222375
0212416*x^28 + 21939438187595227689119578028734304946*x^26 -
7216920454061229140
904518260735360226*x^24 + 1922654474846205830780576037518870734*x^22 -
411234426
233221352649087682301827326*x^20 +
69762980602827026331923086578336536*x^18 - 17
060267332849195295202813592753744097*x^16 +
172325931937787848819688405778161355
9*x^14 - 128961015191389280500991487493925293*x^12 +
156522972970309208863333123
86577923580*x^10 - 559010615983821584316742149697585133*x^8 +
122666890316042110
37213501181251082*x^6 - 2854038599704603513326911971349511228*x^4 +
132131414580
71075585336866474285658*x^2 - 120268396788654757119156692582581476
?



jasonadmin@box1-win32
/cygdrive/c/Users/jasonadmin/pari-2.4.2-alpha-nogmp/src
$ cat src/test/in/polchebyshev  | ./gp.exe
                 GP/PARI CALCULATOR Version 2.4.2 (development)
                i686 running cygwin (ix86 kernel) 32-bit version
 compiled: Jul  3 2009, gcc-3.4.4 (cygming special, gdc 0.12, using dmd
0.125)
                 (readline v5.2 enabled, extended help enabled)

                     Copyright (C) 2000-2006 The PARI Group

PARI/GP is free software, covered by the GNU General Public License, and
comes WITHOUT ANY WARRANTY WHATSOEVER.

Type ? for help, \q to quit.
Type ?12 for how to get moral (and possibly technical) support.

parisize = 4000000, primelimit = 500000
U
T
L
Total time spent: 0
Goodbye!


jasonadmin@box1-win32
/cygdrive/c/Users/jasonadmin/pari-2.4.2-alpha-nogmp/src
$ ./gp.exe
                 GP/PARI CALCULATOR Version 2.4.2 (development)
                i686 running cygwin (ix86 kernel) 32-bit version
 compiled: Jul  3 2009, gcc-3.4.4 (cygming special, gdc 0.12, using dmd
0.125)
                 (readline v5.2 enabled, extended help enabled)

                     Copyright (C) 2000-2006 The PARI Group

PARI/GP is free software, covered by the GNU General Public License, and
comes WITHOUT ANY WARRANTY WHATSOEVER.

Type ? for help, \q to quit.
Type ?12 for how to get moral (and possibly technical) support.

parisize = 4000000, primelimit = 500000
? polchebyshev(26,2)
%1 = 67108864*x^26 - 419430400*x^24 + 1157627904*x^22 - 1857028096*x^20 +
191758
3360*x^18 - 1333592064*x^16 + 635043840*x^14 - 206389248*x^12 +
44808192*x^10 -
6223360*x^8 + 512512*x^6 - 21840*x^4 + 364*x^2 - 1
? polchebyshev(27,2)
%2 = 134217728*x^27 - 872415232*x^25 + 2516582400*x^23 - 4244635648*x^21 +
46425
70240*x^19 - 3451650048*x^17 + 1778122752*x^15 - 635043840*x^13 +
154791936*x^11
 - 24893440*x^9 + 2489344*x^7 - 139776*x^5 + 3640*x^3 - 28*x
? polchebyshev(28,2)
%3 = 268435456*x^28 - 1811939328*x^26 + 5452595200*x^24 - 9646899200*x^22 +
1114
2168576*x^20 - 8820883456*x^18 + 4889837568*x^16 - 1905131520*x^14 +
515973120*x
^12 - 94595072*x^10 + 11202048*x^8 - 792064*x^6 + 29120*x^4 - 420*x^2 + 1
? polchebyshev(29,2)
%4 = 536870912*x^29 - 3758096384*x^27 + 11777605632*x^25 - 21810380800*x^23
+ 26
528972800*x^21 - 22284337152*x^19 + 13231325184*x^17 - 5588385792*x^15 +
1666990
080*x^13 - 343982080*x^11 + 47297536*x^9 - 4073472*x^7 + 198016*x^5 -
4480*x^3 +
 30*x
? polchebyshev(30,2)
%5 = 1073741824*x^30 - 7784628224*x^28 + 25367150592*x^26 - 49073356800*x^24
+ 6
2704844800*x^22 - 55710842880*x^20 + 35283533824*x^18 - 16066609152*x^16 +
52391
11680*x^14 - 1203937280*x^12 + 189190144*x^10 - 19348992*x^8 + 1188096*x^6 -
380
80*x^4 + 480*x^2 - 1
? polchebyshev(31,2)
%6 = 2147483648*x^31 - 16106127360*x^29 + 54492397568*x^27 -
109924319232*x^25 +
 147220070400*x^23 - 137950658560*x^21 + 92851404800*x^19 - 45364543488*x^17
+ 1
6066609152*x^15 - 4074864640*x^13 + 722362368*x^11 - 85995520*x^9 +
6449664*x^7
- 274176*x^5 + 5440*x^3 - 32*x
? polchebyshev(32,2)
%7 = 4294967296*x^32 - 33285996544*x^30 + 116769423360*x^28 -
245215789056*x^26
+ 343513497600*x^24 - 338606161920*x^22 + 241413652480*x^20 -
126012620800*x^18
+ 48199827456*x^16 - 13388840960*x^14 + 2648662016*x^12 - 361181184*x^10 +
32248
320*x^8 - 1736448*x^6 + 48960*x^4 - 544*x^2 + 1
? polchebyshev(33,2)
%8 = 8589934592*x^33 - 68719476736*x^31 + 249644974080*x^29 -
544923975680*x^27
+ 796951314432*x^25 - 824432394240*x^23 + 620777963520*x^21 -
344876646400*x^19
+ 141764198400*x^17 - 42844291072*x^15 + 9372188672*x^13 - 1444724736*x^11 +
150
492160*x^9 - 9922560*x^7 + 372096*x^5 - 6528*x^3 + 34*x
? polchebyshev(34,2)
%9 = 17179869184*x^34 - 141733920768*x^32 + 532575944704*x^30 -
1206617374720*x^
28 + 1839118417920*x^26 - 1992378286080*x^24 + 1580162088960*x^22 -
931166945280
*x^20 + 409541017600*x^18 - 133888409600*x^16 + 32133218304*x^14 -
5538111488*x^
12 + 662165504*x^10 - 52093440*x^8 + 2480640*x^6 - 62016*x^4 + 612*x^2 - 1
? polchebyshev(35,2)
%10 = 34359738368*x^35 - 292057776128*x^33 + 1133871366144*x^31 -
2662879723520*
x^29 + 4223160811520*x^27 - 4781707886592*x^25 + 3984756572160*x^23 -
2483111854
080*x^21 + 1163958681600*x^19 - 409541017600*x^17 + 107110727680*x^15 -
20448411
648*x^13 + 2769055744*x^11 - 254679040*x^9 + 14883840*x^7 - 496128*x^5 +
7752*x^
3 - 36*x
? polchebyshev(36,2)
%11 = 68719476736*x^36 - 601295421440*x^34 + 2409476653056*x^32 -
5858335391744*
x^30 + 9652938997760*x^28 - 11402534191104*x^26 + 9961891430400*x^24 -
654638579
7120*x^22 + 3259084308480*x^20 - 1228623052800*x^18 + 348109864960*x^16 -
730300
41600*x^14 + 11076222976*x^12 - 1171523584*x^10 + 81861120*x^8 - 3472896*x^6
+ 7
7520*x^4 - 684*x^2 + 1
? polchebyshev(37,2)
%12 = 137438953472*x^37 - 1236950581248*x^35 + 5111011082240*x^33 -
128505421496
32*x^31 + 21968757719040*x^29 - 27028229193728*x^27 + 24705490747392*x^25 -
1707
7528166400*x^23 + 9001280471040*x^21 - 3621204787200*x^19 +
1105760747520*x^17 -
 253170810880*x^15 + 42600857600*x^13 - 5112102912*x^11 + 418401280*x^9 -
218296
32*x^7 + 651168*x^5 - 9120*x^3 + 38*x
? polchebyshev(38,2)
%13 = 274877906944*x^38 - 2542620639232*x^36 + 10823317585920*x^34 -
28110560952
320*x^32 + 49795850829824*x^30 - 63709397385216*x^28 + 60813515685888*x^26 -
441
16947763200*x^24 + 24548946739200*x^22 - 10501493882880*x^20 +
3440144547840*x^1
8 - 854451486720*x^16 + 158231756800*x^14 - 21300428800*x^12 +
2008326144*x^10 -
 125520384*x^8 + 4775232*x^6 - 95760*x^4 + 760*x^2 - 1
? polchebyshev(39,2)
%14 = 549755813888*x^39 - 5222680231936*x^37 + 22883585753088*x^35 -
61332132986
880*x^33 + 112442243809280*x^31 - 149387552489472*x^29 +
148655260565504*x^27 -
112939386273792*x^25 + 66175421644800*x^23 - 30004268236800*x^21 +
1050149388288
0*x^19 - 2814663720960*x^17 + 569634324480*x^15 - 85201715200*x^13 +
9128755200*
x^11 - 669442048*x^9 + 31380096*x^7 - 842688*x^5 + 10640*x^3 - 40*x
? polchebyshev(40,2)
%15 = 1099511627776*x^40 - 10720238370816*x^38 + 48309792145408*x^36 -
133487583
559680*x^34 + 252995048570880*x^32 - 348570955808768*x^30 +
361019918516224*x^28
 - 286692288233472*x^26 + 176467791052800*x^24 - 84557483212800*x^22 +
315044816
48640*x^20 - 9069471989760*x^18 + 1993720135680*x^16 - 328635187200*x^14 +
39557
939200*x^12 - 3347210240*x^10 + 188280576*x^8 - 6460608*x^6 + 117040*x^4 -
840*x
^2 + 1
? polchebyshev(41,2)
%16 = 2199023255552*x^41 - 21990232555520*x^39 + 101842264522752*x^37 -
28985875
2872448*x^35 + 567322230128640*x^33 - 809584155426816*x^31 +
871427389521920*x^2
9 - 722039837032448*x^27 + 465874968379392*x^25 - 235290388070400*x^23 +
9301323
1534080*x^21 - 28640437862400*x^19 + 6802103992320*x^17 - 1226904698880*x^15
+ 1
64317593600*x^13 - 15823175680*x^11 + 1046003200*x^9 - 44301312*x^7 +
1076768*x^
5 - 12320*x^3 + 42*x
? polchebyshev(42,2)
%17 = 4398046511104*x^42 - 45079976738816*x^40 + 214404767416320*x^38 -
62802729
7890304*x^36 + 1268132043816960*x^34 - 1872163359424512*x^32 +
2091425734852608*
x^30 - 1805099592581120*x^28 + 1218442224992256*x^26 - 647048567193600*x^24
+ 27
0583946280960*x^22 - 88785357373440*x^20 + 22673679974400*x^18 -
4447529533440*x
^16 + 657270374400*x^14 - 71204290560*x^12 + 5439216640*x^10 - 276883200*x^8
+ 8
614144*x^6 - 141680*x^4 + 924*x^2 - 1
? polchebyshev(43,2)
%18 = 8796093022208*x^43 - 92358976733184*x^41 + 450799767388160*x^39 -
13578968
60303360*x^37 + 2826122840506368*x^35 - 4311648948977664*x^33 +
4992435625132032
*x^31 - 4481626574684160*x^29 + 3158924287016960*x^27 -
1759972102766592*x^25 +
776458280632320*x^23 - 270583946280960*x^21 + 73987797811200*x^19 -
156971630592
00*x^17 + 2541445447680*x^15 - 306726174720*x^13 + 26701608960*x^11 -
1599769600
*x^9 + 61529600*x^7 - 1360128*x^5 + 14168*x^3 - 44*x
? polchebyshev(44,2)
%19 = 17592186044416*x^44 - 189115999977472*x^42 + 946679511515136*x^40 -
293019
8488023040*x^38 + 6280272978903040*x^36 - 9891429941772288*x^34 +
11857034609688
576*x^32 - 11054678884220928*x^30 + 8122948166615040*x^28 -
4738386430525440*x^2
6 + 2199965128458240*x^24 - 811751838842880*x^22 + 236760952995840*x^20 -
540680
06092800*x^18 + 9530420428800*x^16 - 1270722723840*x^14 +
124607508480*x^12 - 86
38755840*x^10 + 399942400*x^8 - 11334400*x^6 + 170016*x^4 - 1012*x^2 + 1
? polchebyshev(45,2)
%20 = 35184372088832*x^45 - 387028092977152*x^43 + 1985717999763456*x^41 -
63111
96743434240*x^39 + 13918442818109440*x^37 - 22608982724050944*x^35 +
28025718168
354816*x^33 - 27101793393573888*x^31 + 20727522907914240*x^29 -
1263569714806784
0*x^27 + 6159902359683072*x^25 - 2399961958318080*x^23 +
744105852272640*x^21 -
182123809996800*x^19 + 34758003916800*x^17 - 5082890895360*x^15 +
555941191680*x
^13 - 43979120640*x^11 + 2399654400*x^9 - 84198400*x^7 + 1700160*x^5 -
16192*x^3
 + 46*x
? polchebyshev(46,2)
%21 = 70368744177664*x^46 - 791648371998720*x^44 + 4160551999504384*x^42 -
13569
072998383616*x^40 + 30767084124241920*x^38 - 51498238427004928*x^36 +
6594286627
8481920*x^34 - 66060621396836352*x^32 + 52509724700049408*x^30 -
333943424627507
20*x^28 + 17058191149891584*x^26 - 6999889045094400*x^24 +
2299963543388160*x^22
 - 601008572989440*x^20 + 123584013926400*x^18 - 19696202219520*x^16 +
238260510
7200*x^14 - 212565749760*x^12 + 13438064640*x^10 - 568339200*x^8 +
14734720*x^6
- 202400*x^4 + 1104*x^2 - 1
? polchebyshev(47,2)
%22 = 140737488355328*x^47 - 1618481116086272*x^45 + 8708132091985920*x^43 -
291
23863996530688*x^41 + 67845364991918080*x^39 - 116914919672119296*x^37 +
1544947
15281014784*x^35 - 160146960962027520*x^33 + 132121242793672704*x^31 -
875162078
33415680*x^29 + 46752079447851008*x^27 - 20159680449871872*x^25 +
69998890450944
00*x^23 - 1946122998251520*x^21 + 429291837849600*x^19 - 74150408355840*x^17
+ 9
848101109760*x^15 - 981072691200*x^13 + 70855249920*x^11 - 3536332800*x^9 +
1136
67840*x^7 - 2104960*x^5 + 18400*x^3 - 48*x
? polchebyshev(48,2)
%23 = 281474976710656*x^48 - 3307330976350208*x^46 +
18207912555970560*x^44 - 62
408279992565760*x^42 + 149259802982219776*x^40 - 264596923468480512*x^38 +
36048
7668989034496*x^36 - 386236788202536960*x^34 + 330303106984181760*x^32 -
2275421
40366880768*x^30 + 126898501358452736*x^28 - 57377552049635328*x^26 +
2099966713
5283200*x^24 - 6192209539891200*x^22 + 1459592248688640*x^20 -
271884830638080*x
^18 + 39392404439040*x^16 - 4344750489600*x^14 + 354276249600*x^12 -
20510730240
*x^10 + 795674880*x^8 - 18944640*x^6 + 239200*x^4 - 1200*x^2 + 1
? polchebyshev(49,2)
%24 = 562949953421312*x^49 - 6755399441055744*x^47 +
38034306228027392*x^45 - 13
3524692077117440*x^43 + 327643469960970240*x^41 - 597039211928879104*x^39 +
8378
90257650188288*x^37 - 926968291686088704*x^35 + 820753174930391040*x^33 -
587205
523527434240*x^31 + 341313210550321152*x^29 - 161507183547121664*x^27 +
62159014
720438272*x^25 - 19384308124876800*x^23 + 4865307495628800*x^21 -
97306149912576
0*x^19 + 152935217233920*x^17 - 18537602088960*x^15 + 1689625190400*x^13 -
11187
6710400*x^11 + 5127682560*x^9 - 151557120*x^7 + 2583360*x^5 - 20800*x^3 +
50*x
? polchebyshev(50,2)
%25 = 1125899906842624*x^50 - 13792273858822144*x^48 +
79375943432404992*x^46 -
285257296710205440*x^44 + 717695219914506240*x^42 - 1343338226839977984*x^40
+ 1
940377438768857088*x^38 - 2214424252361211904*x^36 +
2027743138063319040*x^34 -
1504714154039050240*x^32 + 910168561467523072*x^30 - 449912868452696064*x^28
+ 1
81695581490511872*x^26 - 59768283385036800*x^24 + 15922824531148800*x^22 -
34057
15246940160*x^20 + 577755265105920*x^18 - 76467608616960*x^16 +
7724000870400*x^
14 - 578029670400*x^12 + 30766095360*x^10 - 1098789120*x^8 + 24111360*x^6 -
2808
00*x^4 + 1300*x^2 - 1
?

>> Pari svn on cygwin32 breaks on ellsea,ideal and aurifeuille(STATIC)
>
> I covered ellsea and ideal earlier. aurifeuille should work fine.
> The aurifeuille problem is probably caused by the use of 'install'.
> Is the test 'program' passing ?  Is install working correctly ?
> Maybe 'install' on cygwin cannot cope with long symbol names like
> 'factor_Aurifeuille_prime' but work fine with 'addii'.
>
> Thanks for performing all these tests!
>
> Cheers,
> Bill.

[Bill Allombert-3]

On Fri, Jul 03, 2009 at 04:12:26PM +0100, Jason Moxham wrote:

>
>> I'll have a look at the others when I boot back to linux , I've just
>> run the
>> test under MSVC , and now I've got a lot errors to deal with , although
>> most look fairly trivial , ie stacksize changes
>>
>
> That was with pari-2.4.2-alpha  , I've tried the same with pari-svn and
> get a completely different set of errors , eg
>
>
>  ***   at top-level: setrand(1);bnf=bnfinit(x^2+105)
>  ***                                ^----------------
>  *** bnfinit: bug in PARI/GP (Segmentation Fault), please report
>  ***   at top-level: for(i=1,1000,do(i))
>  ***                              ^------
>  ***   in function do: my(t=bnfisintnorm(bnf,i))
>  ***                        ^--------------------
>  *** bnfisintnorm: please apply nfinit first.
>  ***   at top-level: setrand(1);bnf=bnfinit(x^2-65)
>  ***                                ^---------------
>  *** bnfinit: bug in PARI/GP (Segmentation Fault), please report
>  ***   at top-level: for(i=1,1000,do(i-500))
>  ***                              ^----------
>  ***   in function do: my(t=bnfisintnorm(bnf,i))
>  ***                        ^--------------------
>  *** bnfisintnorm: please apply nfinit first.
>  ***   at top-level: setrand(1);bnf=bnfinit(x^5-37)
>  ***                                ^---------------
>  *** bnfinit: bug in PARI/GP (Segmentation Fault), please report
>  ***   at top-level: for(i=1,1000,do(i-500))
>  ***                              ^----------
>  ***   in function do: my(t=bnfisintnorm(bnf,i))
>  ***                        ^--------------------
>  *** bnfisintnorm: please apply nfinit first.
>  ***   at top-level: bnfisintnorm(bnfinit(x^3+5),5)
>
> It looks like all the number field functions cant be found , everything  
> compiles allright , and a lot of BASIC stuff works fine . Clearly I
> missed something.

Looking at this excerpt, bnfinit fails with a SEGV, so all function that need
the output of bnfinit also fail.

Try to do
\g1
bnfinit(x^2+105);
and let us now what that outputs.

Cheers,
Bill.

[Bill Allombert-3]

On Fri, Jul 03, 2009 at 04:05:37PM +0100, Jason Moxham wrote:

> The test "program" fails on cygwin32/pari-svn static only with
>
> $ cat Ocygwin-i686/program-sta.dif
> *** ../src/test/32/program      Fri Jul  3 00:21:58 2009
> --- gp.out      Fri Jul  3 15:55:55 2009
> ***************
> *** 131,139 ****
>  400 1.632424285532931448171405619
>  ? install(addii,GG)
>  ? addii(1,2)
> ! 3
>  ? kill(addii)
>  ? getheap
> ! [26, 3338]
>  ? print("Total time spent: ",gettime);
> ! Total time spent: 40
> --- 131,141 ----
>  400 1.632424285532931448171405619
>  ? install(addii,GG)
>  ? addii(1,2)
> !   ***   at top-level: addii(1,2)
> !   ***                 ^----------
> !   *** addii: bug in PARI/GP (Segmentation Fault), please report
>  ? kill(addii)
>  ? getheap
> ! [25, 3331]
>  ? print("Total time spent: ",gettime);
> ! Total time spent: 0
>
> I assume that explains aurifeuille failing on static only.

Yes, this is bug #828. If you have any idea how to fix it, please tell us.

Cheers,
Bill.

[Bill Allombert-3]

On Fri, Jul 03, 2009 at 04:45:28PM +0100, Jason Moxham wrote:
>
>> This seems to work on linux32 though. Could you compare the output of
>> polchebyshev(n,2) (for n=26..50) on cygwin and linux ?
>>
>
> Typical,  one of the machines I use for linux32 has no bison , the other  

You can install bison if you need it.

> one's network card has broken , and my main one is only booting the wrong
> kernel :( , so here is a cygwin with no gmp which appears to work...

>                 GP/PARI CALCULATOR Version 2.4.2 (development)
>           i686 running cygwin (ix86/GMP-4.2.1 kernel) 32-bit version
>
> parisize = 4000000, primelimit = 500000
> ? polchebyshev(27,2)
> %2 = 134217728*x^27 - 872415232*x^25 + 2516582400*x^23 - 4244635648*x^21
> + 46425
> 70240*x^19 - 3218161810*x^17 + 1657840932*x^15 - 592086047*x^13 +  
> 144320973*x^11
> - 23209513*x^9 + 2320951*x^7 - 130320*x^5 + 3393*x^3 - 26*x

>                 GP/PARI CALCULATOR Version 2.4.2 (development)
>                i686 running cygwin (ix86 kernel) 32-bit version
> ? polchebyshev(27,2)
> %2 = 134217728*x^27 - 872415232*x^25 + 2516582400*x^23 - 4244635648*x^21
> + 46425
> 70240*x^19 - 3451650048*x^17 + 1778122752*x^15 - 635043840*x^13 +  
> 154791936*x^11
> - 24893440*x^9 + 2489344*x^7 - 139776*x^5 + 3640*x^3 - 28*x

Theses results are different (the second is correct)
However I see that revision 9957 fixed a memory corruption bug in
polchebyshev(,2). This might explain the discrepancy.
If you want to be sure apply the patch given by 'svn diff -c 9957'.
But it is more important to focus on the SVN HEAD branch, because
2.4.2 is not supported.

Cheers,
Bill.

[Jason Moxham]

----- Original Message -----
From: "Bill Allombert" <Bill.Allombert@...>
To: <pari-dev@...>
Sent: Friday, July 03, 2009 4:47 PM
Subject: Re: Some bugs?

> On Fri, Jul 03, 2009 at 04:12:26PM +0100, Jason Moxham wrote:
>>
>>> I'll have a look at the others when I boot back to linux , I've just
>>> run the
>>> test under MSVC , and now I've got a lot errors to deal with , although
>>> most look fairly trivial , ie stacksize changes
>>>
>>
>> That was with pari-2.4.2-alpha  , I've tried the same with pari-svn and
>> get a completely different set of errors , eg
>>
>>
>>  ***   at top-level: setrand(1);bnf=bnfinit(x^2+105)
>>  ***                                ^----------------
>>  *** bnfinit: bug in PARI/GP (Segmentation Fault), please report
>>  ***   at top-level: for(i=1,1000,do(i))
>>  ***                              ^------
>>  ***   in function do: my(t=bnfisintnorm(bnf,i))
>>  ***                        ^--------------------
>>  *** bnfisintnorm: please apply nfinit first.
>>  ***   at top-level: setrand(1);bnf=bnfinit(x^2-65)
>>  ***                                ^---------------
>>  *** bnfinit: bug in PARI/GP (Segmentation Fault), please report
>>  ***   at top-level: for(i=1,1000,do(i-500))
>>  ***                              ^----------
>>  ***   in function do: my(t=bnfisintnorm(bnf,i))
>>  ***                        ^--------------------
>>  *** bnfisintnorm: please apply nfinit first.
>>  ***   at top-level: setrand(1);bnf=bnfinit(x^5-37)
>>  ***                                ^---------------
>>  *** bnfinit: bug in PARI/GP (Segmentation Fault), please report
>>  ***   at top-level: for(i=1,1000,do(i-500))
>>  ***                              ^----------
>>  ***   in function do: my(t=bnfisintnorm(bnf,i))
>>  ***                        ^--------------------
>>  *** bnfisintnorm: please apply nfinit first.
>>  ***   at top-level: bnfisintnorm(bnfinit(x^3+5),5)
>>
>> It looks like all the number field functions cant be found , everything
>> compiles allright , and a lot of BASIC stuff works fine . Clearly I
>> missed something.
>
> Looking at this excerpt, bnfinit fails with a SEGV, so all function that
> need
> the output of bnfinit also fail.
>
> Try to do
> \g1
> bnfinit(x^2+105);
> and let us now what that outputs.
>
> Cheers,
> Bill.

This is MSVC build on pari-2.4.2 just to show it works....


C:\Users\jasonadmin\pari-2.4.2-alpha-msvc32nogmp\src\pari_patches>gp
         GP/PARI CALCULATOR Version 2.4.2 (development CHANGES-1.1969)
             ix86 running Windows 3.2 (ix86 kernel) 32-bit version
                        compiled: Jul  3 2009, MSVC-1500
             (readline not compiled in, extended help not enabled)

                     Copyright (C) 2000-2006 The PARI Group

PARI/GP is free software, covered by the GNU General Public License, and
comes WITHOUT ANY WARRANTY WHATSOEVER.

Type ? for help, \q to quit.
Type ?12 for how to get moral (and possibly technical) support.

parisize = 4000000, primelimit = 500000
? bnfinit(x^2+105)
%1 = [[2, 0, 0; 0, 2, 0; 0, 0, 2], [1, 1, 1, 0, 0, 0, 1; 1, 1, 0, 1, 0, 1,
0; 1,
 0, 1, 1, 1, 0, 0], [;], Mat([0.E-37 + 8.680440892848727768120953206*I,
0.E-37 +
 10.75052560806965122274788703*I, 0.E-37 + 11.42712797621331226724875481*I,
0.E-
37 + 12.28745458497605239059615414*I, 0.E-37 +
3.432297943279603018509133350*I,
0.E-37 + 3.770599127351433540759567240*I, 0.E-37 +
7.947234138551688506535677533
*I, 0, 0, 0.E-37]), [[11, [7, 1]~, 1, 1, [4, -105; 1, 4]], [13, [-5, 1]~, 1,
1,
[5, -105; 1, 5]], [19, [-3, 1]~, 1, 1, [3, -105; 1, 3]], [2, [1, 1]~, 2, 1,
[1,
-105; 1, 1]], [3, [0, 1]~, 2, 1, [0, -105; 1, 0]], [5, [0, 1]~, 2, 1,
[0, -105;
1, 0]], [7, [0, 1]~, 2, 1, [0, -105; 1, 0]], [19, [3, 1]~, 1, 1, [-3, -105;
1, -
3]], [13, [5, 1]~, 1, 1, [-5, -105; 1, -5]], [11, [15, 1]~, 1, 1, [-4, -105;
1,
-4]]], 0, [x^2 + 105, [0, 1], -420, 1, [Mat([1, 0.E-38 -
10.24695076595959838322
103868*I]), [1, -10.24695076595959838322103868; 1,
10.24695076595959838322103868
], 0, [2, 0; 0, -210], [210, 0; 0, 2], [105, 0; 0, -1], [105, [0, -105; 1,
0]]],
 [0.E-38 - 10.24695076595959838322103868*I], [1, x], [1, 0; 0, 1], [1, 0,
0, -10
5; 0, 1, 1, 0]], [[8, [2, 2, 2], [[11, 7; 0, 1], [10, 5; 0, 1], [6, 3; 0,
1]]],
1, 1, [2, -1], []], [[-1, 0, 0; 0, -1, 0; 0, 0, -1],
[[-4.7957905455967410881238
87156 + 2.397255585669141291195666440*I], [-4.867534450455582420071478896 +
2.23
3670150445032372911300130*I], [-4.736198448394495460821504790 +
2.57197133451686
2895161734020*I]], [[-9.591581091193482176247774311 -
3.885929721510445185729620
327*I], [-9.735068900911164840142957792 - 6.283185307179586476925286767*I],
[-9.
472396896788990921643009580 - 6.283185307179586476925286767*I]]], 0]
? \g1
   debug = 1
? bnfinit(x^2+105)
Time disc. factorisation: 0
Treating p^k = 2^2
Time round4: 0
get_red_G: starting LLL, prec = 4 (4 + 0)
Time LLL basis: 0
Time mult. table: 0
Time matrices: 0
smallvectors looking for norm < 2.000015259
Time initalg & rootsof1: 0
R1 = 0, R2 = 1
D = 420

*** Bach constant: 0.300000
LIMC = 20, LIMC2 = 20
Time factor base: 0
Time sub factorbase (3 elements): 0
KCZ = 7, KC = 10, n = 15
1 2 3 4 5 6 7
#### Looking for 15 relations (small norms)
8 9 10 11 12 13 14 15
Time small norm relations: 0
  small norms gave 15 relations.
  nb. fact./nb. small norm = 8/14 = 0.571
Time hnfspec [10 x 15] --> [3 x 8]: 0

#### Computing regulator multiple

#### Tentative class number: 16

#### Computing check
Time bestappr/regulator: 15

#### Tentative regulator : 1

 ***** check = 2.071754
Computing powers for subFB: Vecsmall([5, 7, 9])
Time powFBgen: 0

(more relations needed: 1)

++++ cglob = 19: new relation (need 16)
rel = 1^1 5^1 7^1 9^1
Time for this relation: 0
Time hnfadd (4 + 4): 0

#### Computing regulator multiple

#### Tentative class number: 8

#### Computing check
Time bestappr/regulator: 0

#### Tentative regulator : 1

 ***** check = 1.035877
Time cleanarch: 0

#### Computing fundamental units

#### Computing class group generators
Time classgroup generators: 0
%2 = [[2, 0, 0; 0, 2, 0; 0, 0, 2], [1, 1, 1, 0, 0, 0, 1; 1, 1, 0, 1, 0, 1,
0; 1,
 0, 1, 1, 1, 0, 0], [;], Mat([0.E-37 + 8.680440892848727768120953206*I,
0.E-37 +
 10.75052560806965122274788703*I, 0.E-37 + 11.42712797621331226724875481*I,
0.E-
46 + 12.28745458497605239059615414*I, 0.E-37 +
3.432297943279603018509133350*I,
0.E-38 + 3.770599127351433540759567240*I, 0.E-37 +
7.947234138551688506535677533
*I, 0, 0, 0.E-37]), [[11, [7, 1]~, 1, 1, [4, -105; 1, 4]], [13, [-5, 1]~, 1,
1,
[5, -105; 1, 5]], [19, [-3, 1]~, 1, 1, [3, -105; 1, 3]], [2, [1, 1]~, 2, 1,
[1,
-105; 1, 1]], [3, [0, 1]~, 2, 1, [0, -105; 1, 0]], [5, [0, 1]~, 2, 1,
[0, -105;
1, 0]], [7, [0, 1]~, 2, 1, [0, -105; 1, 0]], [19, [3, 1]~, 1, 1, [-3, -105;
1, -
3]], [13, [5, 1]~, 1, 1, [-5, -105; 1, -5]], [11, [15, 1]~, 1, 1, [-4, -105;
1,
-4]]], 0, [x^2 + 105, [0, 1], -420, 1, [Mat([1, 0.E-38 -
10.24695076595959838322
103868*I]), [1, -10.24695076595959838322103868; 1,
10.24695076595959838322103868
], 0, [2, 0; 0, -210], [210, 0; 0, 2], [105, 0; 0, -1], [105, [0, -105; 1,
0]]],
 [0.E-38 - 10.24695076595959838322103868*I], [1, x], [1, 0; 0, 1], [1, 0,
0, -10
5; 0, 1, 1, 0]], [[8, [2, 2, 2], [[11, 7; 0, 1], [10, 5; 0, 1], [6, 3; 0,
1]]],
1, 1, [2, -1], []], [[-1, 0, 0; 0, -1, 0; 0, 0, -1],
[[-4.7957905455967410881238
87156 + 2.397255585669141291195666440*I], [-4.867534450455582420071478896 +
2.23
3670150445032372911300130*I], [-4.736198448394495460821504790 +
2.57197133451686
2895161734020*I]], [[-9.591581091193482176247774311 -
3.885929721510445185729620
327*I], [-9.735068900911164840142957792 - 6.283185307179586476925286767*I],
[-9.
472396896788990921643009580 - 6.283185307179586476925286767*I]]], 0]
?


and here is parisvn with MSVC32 no GMP


C:\Users\jasonadmin\parisvn-msvc32nogmp>gp
            GP/PARI CALCULATOR Version 2.4.3 (development svn-11782)
             ix86 running Windows 3.2 (ix86 kernel) 32-bit version
                        compiled: Jul  3 2009, MSVC-1500
             (readline not compiled in, extended help not enabled)

                     Copyright (C) 2000-2008 The PARI Group

PARI/GP is free software, covered by the GNU General Public License, and
comes
WITHOUT ANY WARRANTY WHATSOEVER.

Type ? for help, \q to quit.
Type ?12 for how to get moral (and possibly technical) support.

parisize = 4000000, primelimit = 500509
? bnfinit(x^2+105)
  ***   at top-level: bnfinit(x^2+105)
  ***                 ^----------------
  *** bnfinit: bug in PARI/GP (Segmentation Fault), please report
  ***   Break loop: type <Return> three times, or Control-d, to go back to
GP)
break> \g1
   debug = 1
break> bnfinit(x^2+105)
Time disc. factorisation: 0
Treating p^k = 2^2
Time round4: 0
get_red_G: starting LLL, prec = 4 (4 + 0)
Time LLL basis: 0
Time mult. table: 0
Time matrices: 0
Time nfinit & rootsof1: 0
R1 = 0, R2 = 1
D = 420
LIMC = 20, LIMC2 = 20
Time factor base: 0
Time sub factorbase (3 elements): 0
KCZ = 7, KC = 10, n = 15
1 2 3 4 5 6 7
#### Looking for 15 relations (small norms)
  ***   at top-level: bnfinit(x^2+105)
  ***                 ^----------------
  *** bnfinit: bug in PARI/GP (Segmentation Fault), please report
  ***   Break loop: type <Return> three times, or Control-d, to go back to
GP)
break>


Thanks
Jason

[Bill Allombert-3]

On Fri, Jul 03, 2009 at 05:54:01PM +0100, Jason Moxham wrote:

> ----- Original Message ----- From: "Bill Allombert"
> <Bill.Allombert@...>
> To: <pari-dev@...>
> Sent: Friday, July 03, 2009 4:47 PM
> Subject: Re: Some bugs?
>
>
> and here is parisvn with MSVC32 no GMP
> break> bnfinit(x^2+105)
> Time disc. factorisation: 0
> Treating p^k = 2^2
> Time round4: 0
> get_red_G: starting LLL, prec = 4 (4 + 0)
> Time LLL basis: 0
> Time mult. table: 0
> Time matrices: 0
> Time nfinit & rootsof1: 0
> R1 = 0, R2 = 1
> D = 420
> LIMC = 20, LIMC2 = 20
> Time factor base: 0
> Time sub factorbase (3 elements): 0
> KCZ = 7, KC = 10, n = 15
> 1 2 3 4 5 6 7
> #### Looking for 15 relations (small norms)
>  ***   at top-level: bnfinit(x^2+105)
>  ***                 ^----------------
>  *** bnfinit: bug in PARI/GP (Segmentation Fault), please report
>  ***   Break loop: type <Return> three times, or Control-d, to go back to
> GP)

So this is a problem with the function small_norm.
Maybe you could use debugger to pin-point it ?
(there is a definite possibility that small_norm has a bug which does
manifest itself on Linux).

Cheers,
Bill.

[Jason Moxham]

----- Original Message -----
From: "Bill Allombert" <Bill.Allombert@...>
To: <pari-dev@...>
Sent: Friday, July 03, 2009 5:45 PM
Subject: Re: Some bugs?

> On Fri, Jul 03, 2009 at 04:45:28PM +0100, Jason Moxham wrote:
>>
>>> This seems to work on linux32 though. Could you compare the output of
>>> polchebyshev(n,2) (for n=26..50) on cygwin and linux ?
>>>
>>
>> Typical,  one of the machines I use for linux32 has no bison , the other
>
> You can install bison if you need it.
>
>> one's network card has broken , and my main one is only booting the wrong
>> kernel :( , so here is a cygwin with no gmp which appears to work...
>
>>                 GP/PARI CALCULATOR Version 2.4.2 (development)
>>           i686 running cygwin (ix86/GMP-4.2.1 kernel) 32-bit version
>>
>> parisize = 4000000, primelimit = 500000
>> ? polchebyshev(27,2)
>> %2 = 134217728*x^27 - 872415232*x^25 + 2516582400*x^23 - 4244635648*x^21
>> + 46425
>> 70240*x^19 - 3218161810*x^17 + 1657840932*x^15 - 592086047*x^13 +
>> 144320973*x^11
>> - 23209513*x^9 + 2320951*x^7 - 130320*x^5 + 3393*x^3 - 26*x
>
>>                 GP/PARI CALCULATOR Version 2.4.2 (development)
>>                i686 running cygwin (ix86 kernel) 32-bit version
>> ? polchebyshev(27,2)
>> %2 = 134217728*x^27 - 872415232*x^25 + 2516582400*x^23 - 4244635648*x^21
>> + 46425
>> 70240*x^19 - 3451650048*x^17 + 1778122752*x^15 - 635043840*x^13 +
>> 154791936*x^11
>> - 24893440*x^9 + 2489344*x^7 - 139776*x^5 + 3640*x^3 - 28*x
>
> Theses results are different (the second is correct)
> However I see that revision 9957 fixed a memory corruption bug in
> polchebyshev(,2). This might explain the discrepancy.
> If you want to be sure apply the patch given by 'svn diff -c 9957'.
> But it is more important to focus on the SVN HEAD branch, because
> 2.4.2 is not supported.
>
> Cheers,
> Bill.

Thanks
The Sage project will be using Pari SVN HEAD for the Windows port so I can
forget this as well :)

Jason

[Jason Moxham]

----- Original Message -----
From: "Bill Allombert" <Bill.Allombert@...>
To: <pari-dev@...>
Sent: Friday, July 03, 2009 11:41 PM
Subject: Re: Some bugs?

> On Fri, Jul 03, 2009 at 05:54:01PM +0100, Jason Moxham wrote:
>> ----- Original Message ----- From: "Bill Allombert"
>> <Bill.Allombert@...>
>> To: <pari-dev@...>
>> Sent: Friday, July 03, 2009 4:47 PM
>> Subject: Re: Some bugs?
>>
>>
>> and here is parisvn with MSVC32 no GMP
>> break> bnfinit(x^2+105)
>> Time disc. factorisation: 0
>> Treating p^k = 2^2
>> Time round4: 0
>> get_red_G: starting LLL, prec = 4 (4 + 0)
>> Time LLL basis: 0
>> Time mult. table: 0
>> Time matrices: 0
>> Time nfinit & rootsof1: 0
>> R1 = 0, R2 = 1
>> D = 420
>> LIMC = 20, LIMC2 = 20
>> Time factor base: 0
>> Time sub factorbase (3 elements): 0
>> KCZ = 7, KC = 10, n = 15
>> 1 2 3 4 5 6 7
>> #### Looking for 15 relations (small norms)
>>  ***   at top-level: bnfinit(x^2+105)
>>  ***                 ^----------------
>>  *** bnfinit: bug in PARI/GP (Segmentation Fault), please report
>>  ***   Break loop: type <Return> three times, or Control-d, to go back to
>> GP)
>
> So this is a problem with the function small_norm.
> Maybe you could use debugger to pin-point it ?
> (there is a definite possibility that small_norm has a bug which does
> manifest itself on Linux).
>
> Cheers,
> Bill.

I don't know how to use the windows CL debugger  , so doing it the old
fashioned way.
in buch2.c from line 2120 we have

    BOUND *= 1 + 1e-6;
    k = N; y[N] = z[N] = 0; x[N] = 0;
    for (av2 = avma;; avma = av2, step(x,y,inc,k))            // POSITION 1
    {                                                                        
         // POSITION 5
      do
      { /* look for primitive element of small norm, cf minim00 */
        int fl = 0;
        double p;
        if (k > 1)
        {
          long l = k-1;
          z[l] = 0;
          for (j=k; j<=N; j++) z[l] += q[l][j]*x[j];
          p = (double)x[k] + z[k];
          y[l] = y[k] + p*p*v[k];
          if (l <= skipfirst && !y[1]) fl = 1;
          x[l] = (long)floor(-z[l] + 0.5);
          k = l;
        }
        for(;; step(x,y,inc,k))
        {
          if (!fl)
          {
            p = (double)x[k] + z[k];
            if (y[k] + p*p*v[k] <= BOUND) break;

            step(x,y,inc,k);

            p = (double)x[k] + z[k];
            if (y[k] + p*p*v[k] <= BOUND) break;
          }
          fl = 0; inc[k] = 1;
          if (++k > N) goto ENDIDEAL;
        }
      } while (k > 1);

      /* element complete */                    //        POSITION 2
      if (zv_content(x) !=1) continue; /* not primitive */     // POSITION 3
      gx = ZM_zc_mul(IDEAL,x);                // POSITION 4
      if (ZV_isscalar(gx)) continue;


the code goes thru pos 1 thru to pos 2 then 3 does not go to 4 but goes
back to 1 , it never gets to pos 5
there is no exit condition on the for loop so it must crash in
step(x,y,inc,k)  , whether this the cause , I dont know ?

I'll keep looking
Jason

[Jason Moxham]

----- Original Message -----
From: "Jason Moxham" <jason@...>
To: <pari-dev@...>
Sent: Saturday, July 04, 2009 12:35 AM
Subject: Re: Some bugs?

> ----- Original Message -----
> From: "Bill Allombert" <Bill.Allombert@...>
> To: <pari-dev@...>
> Sent: Friday, July 03, 2009 11:41 PM
> Subject: Re: Some bugs?
>
>
>> On Fri, Jul 03, 2009 at 05:54:01PM +0100, Jason Moxham wrote:
>>> ----- Original Message ----- From: "Bill Allombert"
>>> <Bill.Allombert@...>
>>> To: <pari-dev@...>
>>> Sent: Friday, July 03, 2009 4:47 PM
>>> Subject: Re: Some bugs?
>>>
>>>
>>> and here is parisvn with MSVC32 no GMP
>>> break> bnfinit(x^2+105)
>>> Time disc. factorisation: 0
>>> Treating p^k = 2^2
>>> Time round4: 0
>>> get_red_G: starting LLL, prec = 4 (4 + 0)
>>> Time LLL basis: 0
>>> Time mult. table: 0
>>> Time matrices: 0
>>> Time nfinit & rootsof1: 0
>>> R1 = 0, R2 = 1
>>> D = 420
>>> LIMC = 20, LIMC2 = 20
>>> Time factor base: 0
>>> Time sub factorbase (3 elements): 0
>>> KCZ = 7, KC = 10, n = 15
>>> 1 2 3 4 5 6 7
>>> #### Looking for 15 relations (small norms)
>>>  ***   at top-level: bnfinit(x^2+105)
>>>  ***                 ^----------------
>>>  *** bnfinit: bug in PARI/GP (Segmentation Fault), please report
>>>  ***   Break loop: type <Return> three times, or Control-d, to go back
>>> to
>>> GP)
>>
>> So this is a problem with the function small_norm.
>> Maybe you could use debugger to pin-point it ?
>> (there is a definite possibility that small_norm has a bug which does
>> manifest itself on Linux).
>>
>> Cheers,
>> Bill.
>
>
> I don't know how to use the windows CL debugger  , so doing it the old
> fashioned way.
> in buch2.c from line 2120 we have
>
>    BOUND *= 1 + 1e-6;
>    k = N; y[N] = z[N] = 0; x[N] = 0;
>    for (av2 = avma;; avma = av2, step(x,y,inc,k))            // POSITION 1
>
> {
> // POSITION 5
>      do
>      { /* look for primitive element of small norm, cf minim00 */
>        int fl = 0;
>        double p;
>        if (k > 1)
>        {
>          long l = k-1;
>          z[l] = 0;
>          for (j=k; j<=N; j++) z[l] += q[l][j]*x[j];
>          p = (double)x[k] + z[k];
>          y[l] = y[k] + p*p*v[k];
>          if (l <= skipfirst && !y[1]) fl = 1;
>          x[l] = (long)floor(-z[l] + 0.5);
>          k = l;
>        }
>        for(;; step(x,y,inc,k))
>        {
>          if (!fl)
>          {
>            p = (double)x[k] + z[k];
>            if (y[k] + p*p*v[k] <= BOUND) break;
>
>            step(x,y,inc,k);
>
>            p = (double)x[k] + z[k];
>            if (y[k] + p*p*v[k] <= BOUND) break;
>          }
>          fl = 0; inc[k] = 1;
>          if (++k > N) goto ENDIDEAL;
>        }
>      } while (k > 1);
>
>      /* element complete */                    //        POSITION 2
>      if (zv_content(x) !=1) continue; /* not primitive */     // POSITION
> 3
>      gx = ZM_zc_mul(IDEAL,x);                // POSITION 4
>      if (ZV_isscalar(gx)) continue;
>
>
> the code goes thru pos 1 thru to pos 2 then 3 does not go to 4 but goes
> back to 1 , it never gets to pos 5
> there is no exit condition on the for loop so it must crash in
> step(x,y,inc,k)  , whether this the cause , I dont know ?
>
> I'll keep looking
> Jason
>
>

I've got a solution , although I'm not if this addresses the underlying
cause.

In buch2.c there is a inline function called step if we rename this to
step_buch2local and all the calls to it , then it works.

Why does this work?

There is another inline function called step in bibli1.c  , and it does't
like two functions with the same name , they are both inlined , but perhaps
the name scope is still global? , I'm getting a little out of my knowledge
range here :)

I try a web search , see if I can come up with something.

Jason

[Jason Moxham]

----- Original Message -----
From: "Jason Moxham" <jason@...>
To: <pari-dev@...>
Sent: Saturday, July 04, 2009 1:05 AM
Subject: Re: Some bugs?

> ----- Original Message -----
> From: "Jason Moxham" <jason@...>
> To: <pari-dev@...>
> Sent: Saturday, July 04, 2009 12:35 AM
> Subject: Re: Some bugs?
>
>
>> ----- Original Message -----
>> From: "Bill Allombert" <Bill.Allombert@...>
>> To: <pari-dev@...>
>> Sent: Friday, July 03, 2009 11:41 PM
>> Subject: Re: Some bugs?
>>
>>
>>> On Fri, Jul 03, 2009 at 05:54:01PM +0100, Jason Moxham wrote:
>>>> ----- Original Message ----- From: "Bill Allombert"
>>>> <Bill.Allombert@...>
>>>> To: <pari-dev@...>
>>>> Sent: Friday, July 03, 2009 4:47 PM
>>>> Subject: Re: Some bugs?
>>>>
>>>>
>>>> and here is parisvn with MSVC32 no GMP
>>>> break> bnfinit(x^2+105)
>>>> Time disc. factorisation: 0
>>>> Treating p^k = 2^2
>>>> Time round4: 0
>>>> get_red_G: starting LLL, prec = 4 (4 + 0)
>>>> Time LLL basis: 0
>>>> Time mult. table: 0
>>>> Time matrices: 0
>>>> Time nfinit & rootsof1: 0
>>>> R1 = 0, R2 = 1
>>>> D = 420
>>>> LIMC = 20, LIMC2 = 20
>>>> Time factor base: 0
>>>> Time sub factorbase (3 elements): 0
>>>> KCZ = 7, KC = 10, n = 15
>>>> 1 2 3 4 5 6 7
>>>> #### Looking for 15 relations (small norms)
>>>>  ***   at top-level: bnfinit(x^2+105)
>>>>  ***                 ^----------------
>>>>  *** bnfinit: bug in PARI/GP (Segmentation Fault), please report
>>>>  ***   Break loop: type <Return> three times, or Control-d, to go back
>>>> to
>>>> GP)
>>>
>>> So this is a problem with the function small_norm.
>>> Maybe you could use debugger to pin-point it ?
>>> (there is a definite possibility that small_norm has a bug which does
>>> manifest itself on Linux).
>>>
>>> Cheers,
>>> Bill.
>>
>>
>> I don't know how to use the windows CL debugger  , so doing it the old
>> fashioned way.
>> in buch2.c from line 2120 we have
>>
>>    BOUND *= 1 + 1e-6;
>>    k = N; y[N] = z[N] = 0; x[N] = 0;
>>    for (av2 = avma;; avma = av2, step(x,y,inc,k))            // POSITION
>> 1
>>
>> { // POSITION 5
>>      do
>>      { /* look for primitive element of small norm, cf minim00 */
>>        int fl = 0;
>>        double p;
>>        if (k > 1)
>>        {
>>          long l = k-1;
>>          z[l] = 0;
>>          for (j=k; j<=N; j++) z[l] += q[l][j]*x[j];
>>          p = (double)x[k] + z[k];
>>          y[l] = y[k] + p*p*v[k];
>>          if (l <= skipfirst && !y[1]) fl = 1;
>>          x[l] = (long)floor(-z[l] + 0.5);
>>          k = l;
>>        }
>>        for(;; step(x,y,inc,k))
>>        {
>>          if (!fl)
>>          {
>>            p = (double)x[k] + z[k];
>>            if (y[k] + p*p*v[k] <= BOUND) break;
>>
>>            step(x,y,inc,k);
>>
>>            p = (double)x[k] + z[k];
>>            if (y[k] + p*p*v[k] <= BOUND) break;
>>          }
>>          fl = 0; inc[k] = 1;
>>          if (++k > N) goto ENDIDEAL;
>>        }
>>      } while (k > 1);
>>
>>      /* element complete */                    //        POSITION 2
>>      if (zv_content(x) !=1) continue; /* not primitive */     // POSITION
>> 3
>>      gx = ZM_zc_mul(IDEAL,x);                // POSITION 4
>>      if (ZV_isscalar(gx)) continue;
>>
>>
>> the code goes thru pos 1 thru to pos 2 then 3 does not go to 4 but goes
>> back to 1 , it never gets to pos 5
>> there is no exit condition on the for loop so it must crash in
>> step(x,y,inc,k)  , whether this the cause , I dont know ?
>>
>> I'll keep looking
>> Jason
>>
>>
>
>
> I've got a solution , although I'm not if this addresses the underlying
> cause.
>
> In buch2.c there is a inline function called step if we rename this to
> step_buch2local and all the calls to it , then it works.
>
> Why does this work?
>
> There is another inline function called step in bibli1.c  , and it does't
> like two functions with the same name , they are both inlined , but
> perhaps the name scope is still global? , I'm getting a little out of my
> knowledge range here :)
>
> I try a web search , see if I can come up with something.
>
> Jason
>

They are both declaired inline static , so it should work , look like a
compiler error?

Jason
 

[Jason Moxham]

----- Original Message -----
From: "Jason Moxham" <jason@...>
To: <pari-dev@...>
Sent: Saturday, July 04, 2009 2:12 AM
Subject: Re: Some bugs?

>
> ----- Original Message -----
> From: "Jason Moxham" <jason@...>
> To: <pari-dev@...>
> Sent: Saturday, July 04, 2009 1:05 AM
> Subject: Re: Some bugs?
>
>
>> ----- Original Message -----
>> From: "Jason Moxham" <jason@...>
>> To: <pari-dev@...>
>> Sent: Saturday, July 04, 2009 12:35 AM
>> Subject: Re: Some bugs?
>>
>>
>>> ----- Original Message -----
>>> From: "Bill Allombert" <Bill.Allombert@...>
>>> To: <pari-dev@...>
>>> Sent: Friday, July 03, 2009 11:41 PM
>>> Subject: Re: Some bugs?
>>>
>>>
>>>> On Fri, Jul 03, 2009 at 05:54:01PM +0100, Jason Moxham wrote:
>>>>> ----- Original Message ----- From: "Bill Allombert"
>>>>> <Bill.Allombert@...>
>>>>> To: <pari-dev@...>
>>>>> Sent: Friday, July 03, 2009 4:47 PM
>>>>> Subject: Re: Some bugs?
>>>>>
>>>>>
>>>>> and here is parisvn with MSVC32 no GMP
>>>>> break> bnfinit(x^2+105)
>>>>> Time disc. factorisation: 0
>>>>> Treating p^k = 2^2
>>>>> Time round4: 0
>>>>> get_red_G: starting LLL, prec = 4 (4 + 0)
>>>>> Time LLL basis: 0
>>>>> Time mult. table: 0
>>>>> Time matrices: 0
>>>>> Time nfinit & rootsof1: 0
>>>>> R1 = 0, R2 = 1
>>>>> D = 420
>>>>> LIMC = 20, LIMC2 = 20
>>>>> Time factor base: 0
>>>>> Time sub factorbase (3 elements): 0
>>>>> KCZ = 7, KC = 10, n = 15
>>>>> 1 2 3 4 5 6 7
>>>>> #### Looking for 15 relations (small norms)
>>>>>  ***   at top-level: bnfinit(x^2+105)
>>>>>  ***                 ^----------------
>>>>>  *** bnfinit: bug in PARI/GP (Segmentation Fault), please report
>>>>>  ***   Break loop: type <Return> three times, or Control-d, to go back
>>>>> to
>>>>> GP)
>>>>
>>>> So this is a problem with the function small_norm.
>>>> Maybe you could use debugger to pin-point it ?
>>>> (there is a definite possibility that small_norm has a bug which does
>>>> manifest itself on Linux).
>>>>
>>>> Cheers,
>>>> Bill.
>>>
>>>
>>> I don't know how to use the windows CL debugger  , so doing it the old
>>> fashioned way.
>>> in buch2.c from line 2120 we have
>>>
>>>    BOUND *= 1 + 1e-6;
>>>    k = N; y[N] = z[N] = 0; x[N] = 0;
>>>    for (av2 = avma;; avma = av2, step(x,y,inc,k))            // POSITION
>>> 1
>>>
>>> { // POSITION 5
>>>      do
>>>      { /* look for primitive element of small norm, cf minim00 */
>>>        int fl = 0;
>>>        double p;
>>>        if (k > 1)
>>>        {
>>>          long l = k-1;
>>>          z[l] = 0;
>>>          for (j=k; j<=N; j++) z[l] += q[l][j]*x[j];
>>>          p = (double)x[k] + z[k];
>>>          y[l] = y[k] + p*p*v[k];
>>>          if (l <= skipfirst && !y[1]) fl = 1;
>>>          x[l] = (long)floor(-z[l] + 0.5);
>>>          k = l;
>>>        }
>>>        for(;; step(x,y,inc,k))
>>>        {
>>>          if (!fl)
>>>          {
>>>            p = (double)x[k] + z[k];
>>>            if (y[k] + p*p*v[k] <= BOUND) break;
>>>
>>>            step(x,y,inc,k);
>>>
>>>            p = (double)x[k] + z[k];
>>>            if (y[k] + p*p*v[k] <= BOUND) break;
>>>          }
>>>          fl = 0; inc[k] = 1;
>>>          if (++k > N) goto ENDIDEAL;
>>>        }
>>>      } while (k > 1);
>>>
>>>      /* element complete */                    //        POSITION 2
>>>      if (zv_content(x) !=1) continue; /* not primitive */     //
>>> POSITION 3
>>>      gx = ZM_zc_mul(IDEAL,x);                // POSITION 4
>>>      if (ZV_isscalar(gx)) continue;
>>>
>>>
>>> the code goes thru pos 1 thru to pos 2 then 3 does not go to 4 but goes
>>> back to 1 , it never gets to pos 5
>>> there is no exit condition on the for loop so it must crash in
>>> step(x,y,inc,k)  , whether this the cause , I dont know ?
>>>
>>> I'll keep looking
>>> Jason
>>>
>>>
>>
>>
>> I've got a solution , although I'm not if this addresses the underlying
>> cause.
>>
>> In buch2.c there is a inline function called step if we rename this to
>> step_buch2local and all the calls to it , then it works.
>>
>> Why does this work?
>>
>> There is another inline function called step in bibli1.c  , and it does't
>> like two functions with the same name , they are both inlined , but
>> perhaps the name scope is still global? , I'm getting a little out of my
>> knowledge range here :)
>>
>> I try a web search , see if I can come up with something.
>>
>> Jason
>>
>
> They are both declaired inline static , so it should work , look like a
> compiler error?
>
> Jason
>
>

On MSVC INLINE was only defined as __inline not __inline static , that
solves the bnfinit() ,
doesn't help with rest though.....

Jason

[Jason Moxham]

On Thursday 02 July 2009 23:11:31 Jason Moxham wrote:

> ----- Original Message -----
> From: "Bill Allombert" <Bill.Allombert@...>
> To: <pari-dev@...>
> Sent: Thursday, July 02, 2009 5:51 PM
> Subject: Re: Some bugs?
>
> > On Thu, Jul 02, 2009 at 05:00:42PM +0100, Jason Moxham wrote:
> >> The tests ellglobalred and galois both static and dynamic are broken on
> >> linux64,cygwin32 on pari-2.4.2 and svn , I assume you allready know this
> >> so I wont post details.
> >
> > Did you sure you installed the optional packages ?
> > (ellglobalred needs elldata and galois needs galdata)
>
> No , I knew there were optional packages , but did not consider that they
> would tested the same way. I'll try them later.
>
> >> Pari SVN on linux 64 breaks on ellsea and ZN both static and dynamic ,
> >> I'll post details when I'm back on that machine.
> >
> > ellsea needs the package seadata. zn has been broken by Karim latest
> > commit (which remove a warning). ideal has been broken in 32bit
> > for several months.
> >
> >> Pari-2.4.2 on Cygwin32 with gmp breaks on polchebshev static and dynamic
> >
> > This seems to work on linux32 though. Could you compare the output of
> > polchebyshev(n,2) (for n=26..50) on cygwin and linux ?
> >
> >> Pari svn on cygwin32 breaks on ellsea,ideal and aurifeuille(STATIC)
> >
> > I covered ellsea and ideal earlier. aurifeuille should work fine.
> > The aurifeuille problem is probably caused by the use of 'install'.
> > Is the test 'program' passing ?  Is install working correctly ?
> > Maybe 'install' on cygwin cannot cope with long symbol names like
> > 'factor_Aurifeuille_prime' but work fine with 'addii'.
> >
> > Thanks for performing all these tests!
> >
> > Cheers,
> > Bill.
>
> I'll have a look at the others when I boot back to linux , I've just run
> the test under MSVC , and now I've got a lot errors to deal with , although
> most look fairly trivial , ie stacksize changes

I manage to get rid of a few more bugs, and I'm left with these

polred,rnf,rnfkummer all have this error

x^7 + Mod(7*y - 7, y^2 - y - 1)*x^6 + Mod(-21*y + 28, y^2 - y - 1)*x^5 - 35*
x^4 + Mod(35*y - 49, y^2 - y - 1)*x^3 + Mod(-7*y + 84, y^2 - y - 1)*x^2 + Mo
d(-14*y + 21, y^2 - y - 1)*x + Mod(-y - 43, y^2 - y - 1)
  ***   at top-level: rnfpolredabs(nfinit(
  ***                 ^--------------------
  *** rnfpolredabs: could not open requested file ./MPQS.gpa/FREL.
x^3 + Mod(y^2 - 2, y^3 - y - 1)*x^2 + Mod(-y + 1, y^3 - y - 1)*x + Mod(y - 1
, y^3 - y - 1)
  *** nfinit: Warning: non-monic polynomial. Result of the form [nf,c].
x^2 - 3646554366
  ***   at top-level: ...*x^2+135345425000900625;#polredabs(p,4)
  ***                                             ^--------------
  *** polredabs: could not open requested file ./MPQS.gpa/FREL.


ffisom   has this error


-------------e=0--------------
[0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
  ***   at top-level: fptest(10007,Mod(1,1
  ***                 ^--------------------
  ***   in function fptest: ...,if(subst(P,x,C[i])==0,0,error("fptest("a","l
  ***                                                   ^--------------------
  ***   user error: fptest(a,10007,Mod(1, 10007)*x^30 + Mod(7812, 10007)*x^28
+ Mod(7090, 10007)*x^27 + Mod(7645, 10007)*x^26 + Mod(4110, 10007)*x^25 +
Mod(3307, 10007)*x^24 + Mod(5763, 10007)*x^23 + Mod(7900, 10007)*x^22 +
Mod(3872, 10007)*x^21 + Mod(8123, 10007)*x^20 + Mod(4076, 10007)*x^19 +
Mod(3265, 10007)*x^18 + Mod(3777, 10007)*x^17 + Mod(3398, 10007)*x^16 +
Mod(5674, 10007)*x^15 + Mod(4018, 10007)*x^14 + Mod(6820, 10007)*x^13 +
Mod(6479, 10007)*x^12 + Mod(984, 10007)*x^11 + Mod(5652, 10007)*x^10 +
Mod(1129, 10007)*x^9 + Mod(7573, 10007)*x^8 + Mod(1822, 10007)*x^7 + Mod(837,
10007)*x^6 + Mod(4169, 10007)*x^5 + Mod(4787, 10007)*x^4 + Mod(1616,
10007)*x^3 + Mod(5185, 10007)*x^2 + Mod(2649, 10007)*x + Mod(1483,
10007),Mod(1, 10007)*x^30 + Mod(1, 10007)*x + Mod(2, 10007))
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0]
-------------e=1--------------
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0]
-------------e=2--------------
[0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
-------------e=3--------------
[0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
-------------e>=4--------------
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
----------large p---------------
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0]

and nffactor


  ***   Warning: new stack size = 16000000 (15.259 Mbytes).
  ***   at top-level: ...17057741307681944498*x^48+1269570586472186440
  ***                                             ^--------------------
  ***   bug in PARI/GP (Segmentation Fault), please report
[x^48 + 12*x^46 + 948*x^44 + 7200*x^42 + 152361*x^40 + 815832*x^38 + 9475380
*x^36 + 44654004*x^34 + 299137536*x^32 + 1335241260*x^30 + 5029216452*x^28 +
 15282825984*x^26 + 37737671337*x^24 + 79579803672*x^22 + 143658877428*x^20
+ 222699104460*x^18 + 303698198961*x^16 + 348787956312*x^14 + 312863646960*x
^12 + 212893847424*x^10 + 111407984496*x^8 + 43762394880*x^6 + 11836253952*x
^4 + 1904684544*x^2 + 136048896]~
[x^48 - 104*x^46 + 4664*x^44 - 122476*x^42 + 2137838*x^40 - 26567700*x^38 +
245144964*x^36 - 1725955872*x^34 + 9441692003*x^32 - 40611588644*x^30 + 1383
56971048*x^28 - 374714866240*x^26 + 807289826646*x^24 - 1380693858220*x^22 +
 1866021172640*x^20 - 1978766780068*x^18 + 1630151673857*x^16 - 102950530102
4*x^14 + 489498952012*x^12 - 170832297056*x^10 + 42133382284*x^8 - 690450713
6*x^6 + 669868016*x^4 - 28899680*x^2 + 16]~
[x^64 - 6384*x^62 + 18261761*x^60 - 31231019568*x^58 + 35925400902280*x^56 -
 29635423138225800*x^54 + 18244443900381139917*x^52 - 8609789775431197305288
*x^50 + 3173715440318358526295493*x^48 - 926253189958924421506713024*x^46 +
216130107574547887816493973792*x^44 - 40600173780591579547211667354912*x^42
+ 6168621134132051706341715912515370*x^40 - 76013943427734813586799139295141
5744*x^38 + 76052955647065900426700689494084391434*x^36 - 617531184355441677
2593815034596603661344*x^34 + 406137362468726923132369140307868626461744*x^3
2 - 21559959594011538596615817038079394912524336*x^30 + 91912201022866770299
8458013093798368177514170*x^28 - 3125139901367629320825272825691018301064633
1888*x^26 + 840017792466500823059432352741453680599704357258*x^24 - 17651171
087877975905738341268934550365403253370816*x^22 + 28594069425611136664663397
3109729677687028435758880*x^20 - 3510138236669349338216717168753743720690200
440965920*x^18 + 31970659328447652254136627725804141367170891626180389*x^16
- 210577501276259853380917767547980363666145460503950096*x^14 + 972633665701
140285445722173171004563157651576027463941*x^12 - 30387893683293821798844628
51269146550248145196800679728*x^10 + 615892792889658839667059092889638564907
5890120481042824*x^8 - 76951350792700777906576936340585240393085356418600329
04*x^6 + 5510293734552561962521574495431679567021222445632508873*x^4 - 19987
78331104544904932086470347413669495129560426038280*x^2 + 2738927449953408337
77347939263771534786080723599733441]~
[x^8 + (-1531349/22619785297920*y^14 + 559416343/10977248747520*y^12 - 25247
684591/1995863408640*y^10 + 103982297321/80715064320*y^8 - 76714376249/12229
55520*y^6 + 111774992009/74186640*y^4 - 26292202739/1686060*y^2 + 1564021312
/54549)*x^6 + (19000645309/37322645741568*y^14 - 96552535399/249482926080*y^
12 + 536960987254321/5488624373760*y^10 - 4194944378809/407651840*y^8 + 5247
721708074653/10089383040*y^6 - 10860233882239/843030*y^4 + 630194805031381/4
636665*y^2 - 796961577488/2755)*x^4 + (-29172555558095/18661322870784*y^14 +
 5062445860141/4248161280*y^12 - 92318704732060083/304923576320*y^10 + 64819
3176741591037/20178766080*y^8 - 8268979844194150111/5044691520*y^6 + 1896413
90835405983/4636665*y^4 - 223283186783900366/515185*y^2 + 8946625664367872/9
405)*x^2 + (41018074761203773/23326653588480*y^14 - 409009053454624763/30492
3576320*y^12 + 117044119382956746521/343039023360*y^10 - 2927890424340695147
291/80715064320*y^8 + 4678762405288733527739/2522345760*y^6 - 38206315875334
762391/824296*y^4 + 2278645136881067296952/4636665*y^2 - 2967400235256172868
48/272745), x^8 + (-16146091/248817638277120*y^14 + 48843457/997931704320*y^
12 - 267727721489/21954497495040*y^10 + 9174250019/7337733120*y^8 - 82557283
8631/13452510720*y^6 + 10012035071/6744240*y^4 - 286258292621/18546660*y^2 +
 140187040/4959)*x^6 + (1723176289/3392967794688*y^14 - 1059662276341/274431
2186880*y^12 + 48715274000629/498965852160*y^10 - 69104247053959/6726255360*
y^8 + 476508837255737/917216640*y^6 - 238731150890387/18546660*y^4 + 5725141
8840649/421515*y^2 - 26269467866896/90915)*x^4 + (-29113232197835/1866132287
0784*y^14 + 1631934517436717/1372156093440*y^12 - 92145545361015407/30492357
6320*y^10 + 323550303468175019/10089383040*y^8 - 8257197132303274099/5044691
520*y^6 + 378862459377191059/9273330*y^4 - 223115572203411654/515185*y^2 + 2
59332507639972352/272745)*x^2 + (40511288842122127/23326653588480*y^14 - 454
626079701495571/343039023360*y^12 + 231451885839430725263/686078046720*y^10
- 60389017554715183063/1681563840*y^8 + 2319938341838169992923/1261172880*y^
6 - 42704610621363386060/927333*y^4 + 2268063476824046324248/4636665*y^2 - 9
8528590048656718144/90915), x^8 + (181633/186613228707840*y^14 - 6489/138601
62560*y^12 - 45386471/5488624373760*y^10 + 178250659/7337733120*y^8 - 259150
32257/10089383040*y^6 + 126060163/1498720*y^4 - 4457231129/4636665*y^2 + 945
4484/4959)*x^6 + (-11712673/1169988894720*y^14 + 2040142733/274431218688*y^1
2 - 93021416071/52522721280*y^10 + 1688263641205/10762008576*y^8 - 890698149
9371/1834433280*y^6 + 14859994827859/296746560*y^4 - 27539242721/421515*y^2
- 10412804654/90915)*x^4 + (-65739851389/2915831698560*y^14 + 8080067384671/
548862437376*y^12 - 615919067582923/228692682240*y^10 + 7139625968047/558581
76*y^8 - 12035658240536459/5044691520*y^6 + 2909281718459207/148373280*y^4 -
 114027481266257/1545555*y^2 + 26994839222348/272745)*x^2 + (351875661187594
13/746452914831360*y^14 - 175794408000830383/5488624373760*y^12 + 2827285081
9680893321/4390899499008*y^10 - 44274526240947916361/107620085760*y^8 + 4488
83263777097920391/40357532160*y^6 - 2698220779200958903/20465280*y^4 + 28100
59487221296958/4636665*y^2 - 81508587175798061/90915), x^8 + (3074279/248817
638277120*y^14 - 127157047/10977248747520*y^12 + 86499746437/21954497495040*
y^10 - 2843672497/4747944960*y^8 + 109487803351/2690502144*y^6 - 89990421617
/74186640*y^4 + 52396810109/3709332*y^2 - 13670576288/272745)*x^6 + (-361812
87281/186613228707840*y^14 + 446988241097/2744312186880*y^12 - 2791944232079
/57774993408*y^10 + 85916604196211/13452510720*y^8 - 3999057670537517/100893
83040*y^6 + 3070998321176/272745*y^4 - 593787426492949/4636665*y^2 + 3730615
4735728/90915)*x^4 + (64206977604683/93306614353920*y^14 - 781339956430123/1
372156093440*y^12 + 151234315443075961/914770728960*y^10 - 42955239518225908
7/20178766080*y^8 + 1307386034599317571/1008938304*y^6 - 168786176615470804/
4636665*y^4 + 127255299187877218/309111*y^2 - 351207138657683456/272745)*x^2
 + (-18530706211095341/23326653588480*y^14 + 1798847927551103773/27443121868
80*y^12 - 16252626552763832987/85759755840*y^10 + 653145277562478020663/2690
5021440*y^8 - 742591759365166575365/504469152*y^6 + 1530656313511014482581/3
7093320*y^4 - 432314981823305760184/927333*y^2 + 131680859590220167616/90915
), x^8 + (3692011/248817638277120*y^14 - 8644699/645720514560*y^12 + 9566281
9793/21954497495040*y^10 - 51309940321/80715064320*y^8 + 113348313275/269050
2144*y^6 - 91885783933/74186640*y^4 + 53072545513/3709332*y^2 - 13761686272/
272745)*x^6 + (-1260198851/6434938920960*y^14 + 450923192773/2744312186880*y
^12 - 53418556138459/1097724874752*y^10 + 43167277559347/6726255360*y^8 - 40
11960654104563/10089383040*y^6 + 209270128922867/18546660*y^4 - 594571618949
411/4636665*y^2 + 37329116221552/90915)*x^4 + (3385038072653/4910874439680*y
^14 - 782493705352717/1372156093440*y^12 + 151409707529271949/914770728960*y
^10 - 214941552982993939/10089383040*y^8 + 1308011461994334679/1008938304*y^
6 - 337682297591934347/9273330*y^4 + 127297072825004458/309111*y^2 - 3513766
27887928064/272745)*x^2 + (-18884591927033039/23326653588480*y^14 + 76223129
967189533/114346341120*y^12 - 131815203944675271599/686078046720*y^10 + 1236
67566697208858791/5044691520*y^8 - 373763922533488473487/252234576*y^6 + 116
278861122276976/2805*y^4 - 433850608352102566712/927333*y^2 + 39636262539108
2293312/272745), x^8 + (1192483/67859355893760*y^14 - 46386469/3659082915840
*y^12 + 5757671539/1995863408640*y^10 - 9816397597/40357532160*y^8 + 3150430
0903/3668866560*y^6 - 6499115329/49457760*y^4 + 1346529151/1686060*y^2 - 127
928104/54549)*x^6 + (-377107043/23326653588480*y^14 + 1193294195/99793170432
*y^12 - 3898646479517/1372156093440*y^10 + 83324546651/326121472*y^8 - 90270
807160297/10089383040*y^6 + 203979621263/1586880*y^4 - 6453521496733/9273330
*y^2 + 3740713928/2755)*x^4 + (-55838547409/2455437219840*y^14 + 81535888034
65/548862437376*y^12 - 1244323126435859/457385364480*y^10 + 2090484253417429
/16143012864*y^8 - 12134063030016491/5044691520*y^6 + 165453933807769/872784
0*y^4 - 93498888636323/1545555*y^2 + 14213849400992/272745)*x^2 + (118472423
79734429/248817638277120*y^14 - 177668994581955997/5488624373760*y^12 + 2861
1858755746067827/4390899499008*y^10 - 134874618655070911657/322860257280*y^8
 + 152404618518667558583/13452510720*y^6 - 4203829397104185047/31236480*y^4
+ 2871462472553380217/4636665*y^2 - 250004160368295232/272745), x^8 + (64582
73/186613228707840*y^14 - 3719191/152461788160*y^12 + 29024100473/5488624373
760*y^10 - 32135840999/80715064320*y^8 + 23297956603/2017876608*y^6 - 211350
0871/16485920*y^4 + 408266029/927333*y^2 - 316594868/272745)*x^6 + (-1110768
17459/373226457415680*y^14 + 18106275053/85759755840*y^12 - 509560146214411/
10977248747520*y^10 + 195109797344857/53810042880*y^8 - 2298235154532701/201
78766080*y^6 + 87655160117623/59349312*y^4 - 32564739888281/4636665*y^2 + 96
8008850258/90915)*x^4 + (83518124920963/93306614353920*y^14 - 91728776500189
/144437483520*y^12 + 42591591062192407/304923576320*y^10 - 88154696029372417
3/80715064320*y^8 + 86773202255301785/252234576*y^6 - 665429827885409833/148
373280*y^4 + 2218979482348247/103037*y^2 - 8878201815711796/272745)*x^2 + (-
246642305436057809/248817638277120*y^14 + 3860028638255806541/5488624373760*
y^12 - 3395672666468176433467/21954497495040*y^10 + 3904062052800146734661/3
22860257280*y^8 - 1024736920815948655471/2690502144*y^6 + 294713641271965414
9381/593493120*y^4 - 22118315020431326588/927333*y^2 + 9833627869746746873/2
72745), x^8 + (38996983/746452914831360*y^14 - 45610777/1219694305280*y^12 +
 184040442823/21954497495040*y^10 - 6929849477/10089383040*y^8 + 19138097395
1/8071506432*y^6 - 6002605199/16485920*y^4 + 8595505583/3709332*y^2 - 144096
9832/272745)*x^6 + (-56627432063/186613228707840*y^14 + 1182421824943/548862
4373760*y^12 - 260406121516747/5488624373760*y^10 + 200262331426543/53810042
880*y^8 - 297480952643123/2522345760*y^6 + 92330503999361/59349312*y^4 - 710
13869059163/9273330*y^2 + 1098811062152/90915)*x^4 + (83613652919587/9330661
4353920*y^14 - 1744863264395929/2744312186880*y^12 + 42642092214501583/30492
3576320*y^10 - 882655350615967327/80715064320*y^8 + 173789516969368021/50446
9152*y^6 - 666541614309814567/148373280*y^4 + 2223816735480437/103037*y^2 -
306949426473056/9405)*x^2 + (-742452471836109913/746452914831360*y^14 + 3873
208908206849579/5488624373760*y^12 - 3407275945042823422193/21954497495040*y
^10 + 1305807204358253484253/107620085760*y^8 - 3084740075192713038863/80715
06432*y^6 + 173955024027274456667/34911360*y^4 - 22194121937138006191/927333
*y^2 + 3289110492164463424/90915)]~
[x^2 + (-872560111/1750783970525184*y^18 - 103213549/291797328420864*y^16 -
4176139757/7204872306688*y^14 - 4416978655/10807308460032*y^12 - 13066402491
575/72949332105216*y^10 - 215588492387/1736888859648*y^8 - 892238763206647/1
09423998157824*y^6 - 89512294102405/18237333026304*y^4 - 18290521472141/2532
96292032*y^2 - 284080544581/14072016224), x^2 + (872560111/1750783970525184*
y^18 - 103213549/291797328420864*y^16 + 4176139757/7204872306688*y^14 - 4416
978655/10807308460032*y^12 + 13066402491575/72949332105216*y^10 - 2155884923
87/1736888859648*y^8 + 892238763206647/109423998157824*y^6 - 89512294102405/
18237333026304*y^4 + 18290521472141/253296292032*y^2 - 284080544581/14072016
224), x^2 + (-10955729699/63028222938906624*y^19 + 1635275/72949332105216*y^
17 - 471200483153/2334378627366912*y^15 + 68074817/2701827115008*y^13 - 1631
44665311371/2626175955787776*y^11 + 10788711995/1519777752192*y^9 - 10805585
507756651/3939263933681664*y^7 + 476476835441/4559333256576*y^5 - 6980744239
9435/3039555504384*y^3 - 4848316573/3518004056*y)*x + (-240314803/1313087977
893888*y^18 + 2177899/20842666315776*y^16 - 10352388961/48632888070144*y^14
+ 650809375/5403654230016*y^12 - 3599460371687/54711999078912*y^10 + 2204400
27581/6079111008768*y^8 - 35061053609941/11723999802624*y^6 + 11692976745493
/9118666513152*y^4 - 1447968833323/63324073008*y^2 - 24126350651/7036008112)
, x^2 + (-4768637809/63028222938906624*y^19 - 116295749/1750783970525184*y^1
7 - 205045552843/2334378627366912*y^15 - 4961577191/64843850760192*y^13 - 10
134809023919/375167993683968*y^11 - 1681738838629/72949332105216*y^9 - 46714
70958726025/3939263933681664*y^7 - 93324108785933/109423998157824*y^5 - 2896
7498100409/3039555504384*y^3 - 245294011997/84432097344*y)*x + (-496188853/1
750783970525184*y^18 + 8017631/48632888070144*y^16 - 338265737/1029267472384
*y^14 + 338938437/1801218076672*y^12 - 7349975872277/72949332105216*y^10 + 1
12959067435/2026370336256*y^8 - 471966038240125/109423998157824*y^6 + 558940
2408335/3039555504384*y^4 - 8241706316231/253296292032*y^2 + 1651761915/1005
144016), x^2 + (-550587305/9004031848415232*y^19 - 8520109/250111995789312*y
^17 - 165341932517/2334378627366912*y^15 - 2570761369/64843850760192*y^13 -
56851147336519/2626175955787776*y^11 - 895619999627/72949332105216*y^9 - 356
0831957116295/3939263933681664*y^7 - 62314525927891/109423998157824*y^5 - 18
560457963575/3039555504384*y^3 - 371119778467/84432097344*y)*x + (496188853/
1750783970525184*y^18 + 8017631/48632888070144*y^16 + 338265737/102926747238
4*y^14 + 338938437/1801218076672*y^12 + 7349975872277/72949332105216*y^10 +
112959067435/2026370336256*y^8 + 471966038240125/109423998157824*y^6 + 55894
02408335/3039555504384*y^4 + 8241706316231/253296292032*y^2 + 1651761915/100
5144016), x^2 + (-21512837/1313087977893888*y^19 - 102134281/175078397052518
4*y^17 - 932348591/48632888070144*y^15 - 4407863827/64843850760192*y^13 - 16
4654686001/27355999539456*y^11 - 1540975821641/72949332105216*y^9 - 25030796
719091/82067998618368*y^7 - 110013175606897/109423998157824*y^5 - 7018393859
/2261574036*y^3 - 296247471259/28144032448*y)*x + (240314803/131308797789388
8*y^18 + 2177899/20842666315776*y^16 + 10352388961/48632888070144*y^14 + 650
809375/5403654230016*y^12 + 3599460371687/54711999078912*y^10 + 220440027581
/6079111008768*y^8 + 35061053609941/11723999802624*y^6 + 11692976745493/9118
666513152*y^4 + 1447968833323/63324073008*y^2 - 24126350651/7036008112), x^2
 + (21512837/1313087977893888*y^19 + 102134281/1750783970525184*y^17 + 93234
8591/48632888070144*y^15 + 4407863827/64843850760192*y^13 + 164654686001/273
55999539456*y^11 + 1540975821641/72949332105216*y^9 + 25030796719091/8206799
8618368*y^7 + 110013175606897/109423998157824*y^5 + 7018393859/2261574036*y^
3 + 296247471259/28144032448*y)*x + (240314803/1313087977893888*y^18 + 21778
99/20842666315776*y^16 + 10352388961/48632888070144*y^14 + 650809375/5403654
230016*y^12 + 3599460371687/54711999078912*y^10 + 220440027581/6079111008768
*y^8 + 35061053609941/11723999802624*y^6 + 11692976745493/9118666513152*y^4
+ 1447968833323/63324073008*y^2 - 24126350651/7036008112), x^2 + (550587305/
9004031848415232*y^19 + 8520109/250111995789312*y^17 + 165341932517/23343786
27366912*y^15 + 2570761369/64843850760192*y^13 + 56851147336519/262617595578
7776*y^11 + 895619999627/72949332105216*y^9 + 3560831957116295/3939263933681
664*y^7 + 62314525927891/109423998157824*y^5 + 18560457963575/3039555504384*
y^3 + 371119778467/84432097344*y)*x + (496188853/1750783970525184*y^18 + 801
7631/48632888070144*y^16 + 338265737/1029267472384*y^14 + 338938437/18012180
76672*y^12 + 7349975872277/72949332105216*y^10 + 112959067435/2026370336256*
y^8 + 471966038240125/109423998157824*y^6 + 5589402408335/3039555504384*y^4
+ 8241706316231/253296292032*y^2 + 1651761915/1005144016), x^2 + (4768637809
/63028222938906624*y^19 + 116295749/1750783970525184*y^17 + 205045552843/233
4378627366912*y^15 + 4961577191/64843850760192*y^13 + 10134809023919/3751679
93683968*y^11 + 1681738838629/72949332105216*y^9 + 4671470958726025/39392639
33681664*y^7 + 93324108785933/109423998157824*y^5 + 28967498100409/303955550
4384*y^3 + 245294011997/84432097344*y)*x + (-496188853/1750783970525184*y^18
 + 8017631/48632888070144*y^16 - 338265737/1029267472384*y^14 + 338938437/18
01218076672*y^12 - 7349975872277/72949332105216*y^10 + 112959067435/20263703
36256*y^8 - 471966038240125/109423998157824*y^6 + 5589402408335/303955550438
4*y^4 - 8241706316231/253296292032*y^2 + 1651761915/1005144016), x^2 + (1095
5729699/63028222938906624*y^19 - 1635275/72949332105216*y^17 + 471200483153/
2334378627366912*y^15 - 68074817/2701827115008*y^13 + 163144665311371/262617
5955787776*y^11 - 10788711995/1519777752192*y^9 + 10805585507756651/39392639
33681664*y^7 - 476476835441/4559333256576*y^5 + 69807442399435/3039555504384
*y^3 + 4848316573/3518004056*y)*x + (-240314803/1313087977893888*y^18 + 2177
899/20842666315776*y^16 - 10352388961/48632888070144*y^14 + 650809375/540365
4230016*y^12 - 3599460371687/54711999078912*y^10 + 220440027581/607911100876
8*y^8 - 35061053609941/11723999802624*y^6 + 11692976745493/9118666513152*y^4
 - 1447968833323/63324073008*y^2 - 24126350651/7036008112)]~
[x^24 + 69]~
[x^7 - 2*y*x^6 + y^2*x^5 - 28*x^3 + 4*y^2*x + 16/7*y^3, x^7 + 2*y*x^6 + y^2*
x^5 - 28*x^3 + 4*y^2*x - 16/7*y^3, x^7 - 2/7*y^3*x^6 - y^2*x^5 - 28*x^3 - 4*
y^2*x + 16*y, x^7 + 2/7*y^3*x^6 - y^2*x^5 - 28*x^3 - 4*y^2*x - 16*y]~
[x^2 + (-12035/386*y^15 + 8337/386*y^14 + 566267/772*y^13 - 392327/772*y^12
- 4449119/772*y^11 + 3083327/772*y^10 + 7320431/386*y^9 - 2537889/193*y^8 -
9975849/386*y^7 + 3461954/193*y^6 + 10712085/772*y^5 - 7445851/772*y^4 - 190
5321/772*y^3 + 1322419/772*y^2 + 25961/386*y - 17735/386), x^2 + (-7909/772*
y^15 - 1879/772*y^14 + 93535/386*y^13 + 44921/772*y^12 - 371330/193*y^11 - 9
0974/193*y^10 + 4991389/772*y^9 + 317305/193*y^8 - 7127131/772*y^7 - 1947859
/772*y^6 + 2110371/386*y^5 + 1295397/772*y^4 - 224091/193*y^3 - 160649/386*y
^2 + 19739/772*y + 4029/386), x^2 + (7909/772*y^15 - 1879/772*y^14 - 93535/3
86*y^13 + 44921/772*y^12 + 371330/193*y^11 - 90974/193*y^10 - 4991389/772*y^
9 + 317305/193*y^8 + 7127131/772*y^7 - 1947859/772*y^6 - 2110371/386*y^5 + 1
295397/772*y^4 + 224091/193*y^3 - 160649/386*y^2 - 19739/772*y + 4029/386),
x^2 + (12035/386*y^15 + 8337/386*y^14 - 566267/772*y^13 - 392327/772*y^12 +
4449119/772*y^11 + 3083327/772*y^10 - 7320431/386*y^9 - 2537889/193*y^8 + 99
75849/386*y^7 + 3461954/193*y^6 - 10712085/772*y^5 - 7445851/772*y^4 + 19053
21/772*y^3 + 1322419/772*y^2 - 25961/386*y - 17735/386), x^4 + (541/386*y^14
 - 12373/386*y^12 + 46089/193*y^10 - 137477/193*y^8 + 311469/386*y^6 - 14604
9/386*y^4 + 12308/193*y^2 - 686/193)*x^2 + (-765/772*y^14 + 8659/386*y^12 -
126401/772*y^10 + 90322/193*y^8 - 369401/772*y^6 + 38430/193*y^4 - 23865/772
*y^2 + 553/386), x^4 + (2052/193*y^14 - 96741/386*y^12 + 763031/386*y^10 - 1
265475/193*y^8 + 1759047/193*y^6 - 2000657/386*y^4 + 430823/386*y^2 - 13920/
193)*x^2 + (-28293/772*y^14 + 166359/193*y^12 - 5225705/772*y^10 + 4295689/1
93*y^8 - 23400329/772*y^6 + 6320077/386*y^4 - 2369321/772*y^2 + 45727/386),
x^4 + (-5155/193*y^15 - 15405/772*y^14 + 241769/386*y^13 + 180325/386*y^12 -
 944044/193*y^11 - 2807115/772*y^10 + 3072303/193*y^9 + 2267740/193*y^8 - 41
02400/193*y^7 - 11896515/772*y^6 + 4334901/386*y^5 + 2988835/386*y^4 - 43236
3/193*y^3 - 1014845/772*y^2 + 23130/193*y + 3602/193)*x^2 + (14406/193*y^15
+ 38997/772*y^14 - 338423/193*y^13 - 455769/386*y^12 + 2652454/193*y^11 + 70
72893/772*y^10 - 8698762/193*y^9 - 11354279/386*y^8 + 23683093/386*y^7 + 293
19183/772*y^6 - 13169735/386*y^5 - 7082173/386*y^4 + 3014249/386*y^3 + 21677
71/772*y^2 - 244333/386*y + 2027/386), x^4 + (-4705/193*y^15 - 4035/772*y^14
 + 220691/386*y^13 + 47935/386*y^12 - 1722817/386*y^11 - 767525/772*y^10 + 5
588151/386*y^9 + 654320/193*y^8 - 3658160/193*y^7 - 3833005/772*y^6 + 341354
7/386*y^5 + 1159745/386*y^4 - 194425/386*y^3 - 441715/772*y^2 - 129025/386*y
 - 9518/193)*x^2 + (69045/772*y^15 + 17477/772*y^14 - 1634153/772*y^13 - 208
093/386*y^12 + 12990539/772*y^11 + 3347413/772*y^10 - 43749763/772*y^9 - 576
5109/386*y^8 + 62724993/772*y^7 + 17320883/772*y^6 - 37352513/772*y^5 - 5621
249/386*y^4 + 7969659/772*y^3 + 2742727/772*y^2 - 182407/772*y - 34007/386),
 x^4 + (4705/193*y^15 - 4035/772*y^14 - 220691/386*y^13 + 47935/386*y^12 + 1
722817/386*y^11 - 767525/772*y^10 - 5588151/386*y^9 + 654320/193*y^8 + 36581
60/193*y^7 - 3833005/772*y^6 - 3413547/386*y^5 + 1159745/386*y^4 + 194425/38
6*y^3 - 441715/772*y^2 + 129025/386*y - 9518/193)*x^2 + (-69045/772*y^15 + 1
7477/772*y^14 + 1634153/772*y^13 - 208093/386*y^12 - 12990539/772*y^11 + 334
7413/772*y^10 + 43749763/772*y^9 - 5765109/386*y^8 - 62724993/772*y^7 + 1732
0883/772*y^6 + 37352513/772*y^5 - 5621249/386*y^4 - 7969659/772*y^3 + 274272
7/772*y^2 + 182407/772*y - 34007/386), x^4 + (5155/193*y^15 - 15405/772*y^14
 - 241769/386*y^13 + 180325/386*y^12 + 944044/193*y^11 - 2807115/772*y^10 -
3072303/193*y^9 + 2267740/193*y^8 + 4102400/193*y^7 - 11896515/772*y^6 - 433
4901/386*y^5 + 2988835/386*y^4 + 432363/193*y^3 - 1014845/772*y^2 - 23130/19
3*y + 3602/193)*x^2 + (-14406/193*y^15 + 38997/772*y^14 + 338423/193*y^13 -
455769/386*y^12 - 2652454/193*y^11 + 7072893/772*y^10 + 8698762/193*y^9 - 11
354279/386*y^8 - 23683093/386*y^7 + 29319183/772*y^6 + 13169735/386*y^5 - 70
82173/386*y^4 - 3014249/386*y^3 + 2167771/772*y^2 + 244333/386*y + 2027/386)
]~
[x^4 - 4*y*x^3 + (-1/2*y^14 + 1/2*y^10 - 7/2*y^6 + 15/2*y^2)*x^2 + (y^15 - y
^11 + 7*y^7 - 7*y^3)*x + (1/2*y^12 - y^8 + 5/2*y^4 - 2), x^4 + 4*y*x^3 + (-1
/2*y^14 + 1/2*y^10 - 7/2*y^6 + 15/2*y^2)*x^2 + (-y^15 + y^11 - 7*y^7 + 7*y^3
)*x + (1/2*y^12 - y^8 + 5/2*y^4 - 2), x^4 + (-2*y^11 - 10*y^3)*x^3 + (9/2*y^
14 - 1/2*y^10 + 55/2*y^6 - 7/2*y^2)*x^2 + (3*y^13 + y^9 + 17*y^5 + 3*y)*x +
(1/2*y^12 + y^8 + 5/2*y^4 + 4), x^4 + (2*y^11 + 10*y^3)*x^3 + (9/2*y^14 - 1/
2*y^10 + 55/2*y^6 - 7/2*y^2)*x^2 + (-3*y^13 - y^9 - 17*y^5 - 3*y)*x + (1/2*y
^12 + y^8 + 5/2*y^4 + 4), x^4 + (-2*y^13 - 10*y^5)*x^3 + (1/2*y^14 - 1/2*y^1
0 + 7/2*y^6 - 15/2*y^2)*x^2 + (3*y^15 + y^11 + 17*y^7 + 3*y^3)*x + (1/2*y^12
 - y^8 + 5/2*y^4 - 2), x^4 + (-y^13 - y^9 - 7*y^5 - 3*y)*x^3 + (7/2*y^14 - 1
/2*y^10 + 37/2*y^6 - 3/2*y^2)*x^2 + (-y^15 + y^11 - 7*y^7 + 7*y^3)*x + (-1/2
*y^12 - y^8 - 5/2*y^4 - 2), x^4 + (-y^13 + y^9 - 7*y^5 + 3*y)*x^3 + (-7/2*y^
14 + 1/2*y^10 - 37/2*y^6 + 3/2*y^2)*x^2 + (-3*y^15 - y^11 - 17*y^7 - 3*y^3)*
x + (-1/2*y^12 - y^8 - 5/2*y^4 - 2), x^4 + (y^13 - y^9 + 7*y^5 - 3*y)*x^3 +
(-7/2*y^14 + 1/2*y^10 - 37/2*y^6 + 3/2*y^2)*x^2 + (3*y^15 + y^11 + 17*y^7 +
3*y^3)*x + (-1/2*y^12 - y^8 - 5/2*y^4 - 2), x^4 + (y^13 + y^9 + 7*y^5 + 3*y)
*x^3 + (7/2*y^14 - 1/2*y^10 + 37/2*y^6 - 3/2*y^2)*x^2 + (y^15 - y^11 + 7*y^7
 - 7*y^3)*x + (-1/2*y^12 - y^8 - 5/2*y^4 - 2), x^4 + (2*y^13 + 10*y^5)*x^3 +
 (1/2*y^14 - 1/2*y^10 + 7/2*y^6 - 15/2*y^2)*x^2 + (-3*y^15 - y^11 - 17*y^7 -
 3*y^3)*x + (1/2*y^12 - y^8 + 5/2*y^4 - 2), x^4 + (-4*y^15 - 24*y^7)*x^3 + (
-9/2*y^14 + 1/2*y^10 - 55/2*y^6 + 7/2*y^2)*x^2 + (-y^13 + y^9 - 7*y^5 + 7*y)
*x + (1/2*y^12 + y^8 + 5/2*y^4 + 4), x^4 + (-3*y^15 - y^11 - 17*y^7 - 7*y^3)
*x^3 + (3/2*y^14 - 5/2*y^10 + 17/2*y^6 - 23/2*y^2)*x^2 + (3*y^13 + y^9 + 17*
y^5 + 3*y)*x + (-1/2*y^12 + y^8 - 5/2*y^4 + 4), x^4 + (-3*y^15 + y^11 - 17*y
^7 + 7*y^3)*x^3 + (-3/2*y^14 + 5/2*y^10 - 17/2*y^6 + 23/2*y^2)*x^2 + (-y^13
+ y^9 - 7*y^5 + 7*y)*x + (-1/2*y^12 + y^8 - 5/2*y^4 + 4), x^4 + (3*y^15 - y^
11 + 17*y^7 - 7*y^3)*x^3 + (-3/2*y^14 + 5/2*y^10 - 17/2*y^6 + 23/2*y^2)*x^2
+ (y^13 - y^9 + 7*y^5 - 7*y)*x + (-1/2*y^12 + y^8 - 5/2*y^4 + 4), x^4 + (3*y
^15 + y^11 + 17*y^7 + 7*y^3)*x^3 + (3/2*y^14 - 5/2*y^10 + 17/2*y^6 - 23/2*y^
2)*x^2 + (-3*y^13 - y^9 - 17*y^5 - 3*y)*x + (-1/2*y^12 + y^8 - 5/2*y^4 + 4),
 x^4 + (4*y^15 + 24*y^7)*x^3 + (-9/2*y^14 + 1/2*y^10 - 55/2*y^6 + 7/2*y^2)*x
^2 + (y^13 - y^9 + 7*y^5 - 7*y)*x + (1/2*y^12 + y^8 + 5/2*y^4 + 4)]~
[x^64 + 192*x^62 + 17568*x^60 + 1019520*x^58 + 42131676*x^56 + 1319651424*x^
54 + 32559096528*x^52 + 649228312512*x^50 + 10651553826426*x^48 + 1456394385
52224*x^46 + 1674922821206832*x^44 + 16307859539653056*x^42 + 13502367773216
7696*x^40 + 953248899971965824*x^38 + 5745239175305568960*x^36 + 29556064271
185194240*x^34 + 129595725382952883843*x^32 + 483002100692576612640*x^30 + 1
523870714370199019760*x^28 + 4047489983524093705152*x^26 + 89858128286488620
19536*x^24 + 16525310345394167002752*x^22 + 24893927149975603242048*x^20 + 3
0294355815129821928192*x^18 + 29274561574319887883226*x^16 + 219878017711043
40121824*x^14 + 12494344840480632094992*x^12 + 5187763623118143696192*x^10 +
 1502211081063677383836*x^8 + 283567347515314680480*x^6 + 311461554388845258
72*x^4 + 1543354925530003776*x^2 + 8057044481403681]~
[x + (-7/16*y^29 - 5/32*y^25 - 97/8*y^21 - 139/32*y^17 - 435/16*y^13 - 323/3
2*y^9 - 21/4*y^5 - 77/32*y), x + (-11/32*y^29 - 1/8*y^25 - 309/32*y^21 - 55/
16*y^17 - 797/32*y^13 - 7*y^9 - 355/32*y^5 + 9/16*y), x + (11/32*y^29 + 1/8*
y^25 + 309/32*y^21 + 55/16*y^17 + 797/32*y^13 + 7*y^9 + 355/32*y^5 - 9/16*y)
, x + (7/16*y^29 + 5/32*y^25 + 97/8*y^21 + 139/32*y^17 + 435/16*y^13 + 323/3
2*y^9 + 21/4*y^5 + 77/32*y), x + (-21/32*y^31 - 587/32*y^23 + 1/16*y^19 - 14
43/32*y^15 + 13/8*y^11 - 541/32*y^7 + 21/16*y^3), x + (-5/16*y^31 - 5/32*y^2
7 - 35/4*y^23 - 139/32*y^19 - 349/16*y^15 - 323/32*y^11 - 57/8*y^7 - 109/32*
y^3), x + (5/16*y^31 + 5/32*y^27 + 35/4*y^23 + 139/32*y^19 + 349/16*y^15 + 3
23/32*y^11 + 57/8*y^7 + 109/32*y^3), x + (21/32*y^31 + 587/32*y^23 - 1/16*y^
19 + 1443/32*y^15 - 13/8*y^11 + 541/32*y^7 - 21/16*y^3), x^4 + (-13/4*y^24 -
 91*y^16 - 909/4*y^8 - 169/2), x^4 + (-1/4*y^24 - 7*y^16 - 57/4*y^8 - 1/2)]~
[8711099/70204123*y^14 - 3396450/70204123*y^13 - 230089978/70204123*y^12 + 7
1459644/70204123*y^11 + 2039293754/70204123*y^10 - 522502724/70204123*y^9 -
7578045032/70204123*y^8 + 136410216/6382193*y^7 + 11598831422/70204123*y^6 -
 1582740050/70204123*y^5 - 6466526698/70204123*y^4 + 712163508/70204123*y^3
+ 865017354/70204123*y^2 - 11706800/70204123*y + 7921687/70204123]
[y, -123209112482/559553426209*y^11 - 236161397417/559553426209*y^10 + 52225
05497467/559553426209*y^9 + 7627164004768/559553426209*y^8 - 68684785347690/
559553426209*y^7 - 98327585435469/559553426209*y^6 + 334508906676131/5595534
26209*y^5 + 508054669424553/559553426209*y^4 - 499853398148011/559553426209*
y^3 - 780815391953932/559553426209*y^2 + 222263541657120/559553426209*y + 69
71304961116/11905392047, -176690268281/1119106852418*y^11 - 641210922141/223
8213704836*y^10 + 3748626002639/559553426209*y^9 + 20339566453621/2238213704
836*y^8 - 49228932581896/559553426209*y^7 - 262748733927015/2238213704836*y^
6 + 958463738831775/2238213704836*y^5 + 1369627950576313/2238213704836*y^4 -
 1439017147702531/2238213704836*y^3 - 2121737405539667/2238213704836*y^2 + 6
35967216650047/2238213704836*y + 19194169373855/47621568188, -68393259315/11
19106852418*y^11 - 143769468519/1119106852418*y^10 + 2890095068571/111910685
2418*y^9 + 4773072512403/1119106852418*y^8 - 38008182942981/1119106852418*y^
7 - 30714902181882/559553426209*y^6 + 92097658832451/559553426209*y^5 + 1570
24332093039/559553426209*y^4 - 134021586074273/559553426209*y^3 - 2404025913
19962/559553426209*y^2 + 58771325675463/559553426209*y + 4293056883813/23810
784094, -26071422312/559553426209*y^11 - 45431912634/559553426209*y^10 + 110
4266993455/559553426209*y^9 + 1416464504385/559553426209*y^8 - 1440773326506
3/559553426209*y^7 - 18255285107310/559553426209*y^6 + 69227537719650/559553
426209*y^5 + 94877281869432/559553426209*y^4 - 100252280561328/559553426209*
y^3 - 143113043669604/559553426209*y^2 + 41073387952215/559553426209*y + 128
0223975324/11905392047, -16163557893/1119106852418*y^11 - 64908872441/223821
3704836*y^10 + 685522417391/1119106852418*y^9 + 2130089587017/2238213704836*
y^8 - 4531221497592/559553426209*y^7 - 27490869176393/2238213704836*y^6 + 89
319951302323/2238213704836*y^5 + 141812611520979/2238213704836*y^4 - 1390748
45426265/2238213704836*y^3 - 220351698213059/2238213704836*y^2 + 72725691407
453/2238213704836*y + 2005246545505/47621568188, 9609209575/2238213704836*y^
11 + 25387385845/559553426209*y^10 - 383645617715/2238213704836*y^9 - 208031
3216715/1119106852418*y^8 + 6033944605985/2238213704836*y^7 + 53081257123225
/2238213704836*y^6 - 33975296862719/2238213704836*y^5 - 255089221619815/2238
213704836*y^4 + 47525895430055/2238213704836*y^3 + 398412597852115/223821370
4836*y^2 - 30904572571305/2238213704836*y - 890343326655/11905392047, 522817
60559/2238213704836*y^11 + 14008471016/559553426209*y^10 - 2227004686929/223
8213704836*y^9 - 336992949231/559553426209*y^8 + 28726688292139/223821370483
6*y^7 + 17917872421753/2238213704836*y^6 - 137376604832131/2238213704836*y^5
 - 105226511459813/2238213704836*y^4 + 207364021022899/2238213704836*y^3 + 1
70780696189363/2238213704836*y^2 - 84647799904795/2238213704836*y - 44402858
6648/11905392047, 106577461171/2238213704836*y^11 + 236835106673/22382137048
36*y^10 - 4499294217001/2238213704836*y^9 - 8003166405117/2238213704836*y^8
+ 59306288158573/2238213704836*y^7 + 25836461923340/559553426209*y^6 - 72093
028502701/559553426209*y^5 - 264172433170095/1119106852418*y^4 + 21076054618
3675/1119106852418*y^3 + 203797855675948/559553426209*y^2 - 94758333876515/1
119106852418*y - 7259569285101/47621568188, 175225361355/2238213704836*y^11
+ 367043962553/2238213704836*y^10 - 7404910374729/2238213704836*y^9 - 121761
49738023/2238213704836*y^8 + 97370244941211/2238213704836*y^7 + 784051152149
35/1119106852418*y^6 - 117954333048705/559553426209*y^5 - 401317015113399/11
19106852418*y^4 + 171430955359538/559553426209*y^3 + 307321505664217/5595534
26209*y^2 - 148696095658675/1119106852418*y - 10923889360787/47621568188, 18
4746017659/1119106852418*y^11 + 182236626649/559553426209*y^10 - 78287971960
99/1119106852418*y^9 - 5941166792275/559553426209*y^8 + 103098612489213/1119
106852418*y^7 + 153130332934513/1119106852418*y^6 - 502641225832273/11191068
52418*y^5 - 789225051651289/1119106852418*y^4 + 749848149414089/111910685241
8*y^3 + 1216440499421323/1119106852418*y^2 - 334431034048089/1119106852418*y
 - 5475399545957/11905392047, 421608384/2321798449*y^11 + 1306271719/4643596
898*y^10 - 17925705039/2321798449*y^9 - 38874965891/4643596898*y^8 + 2342247
56928/2321798449*y^7 + 502105436649/4643596898*y^6 - 2270324281083/464359689
8*y^5 - 2658757184017/4643596898*y^4 + 3428983243105/4643596898*y^3 + 405157
9820087/4643596898*y^2 - 1510458181035/4643596898*y - 36134279083/98799934]

[x + (-y + 1) 1]

[x^2 + (y + 2)*x + (y^2 + y + 1) 1]

[x^2 + (y + 2)*x + (1/25*y^8 - 3/5*y^5 - 87/25*y^2 + y + 1) 1]

[x^2 + (-2/15*y^7 + 7/3*y^4 + 79/15*y + 2)*x + (1/25*y^8 - 2/15*y^7 - 3/5*y^
5 + 7/3*y^4 - 87/25*y^2 + 79/15*y + 1) 1]

[x^2 + (2/15*y^7 - 7/3*y^4 - 94/15*y + 2)*x + (1/25*y^8 + 2/15*y^7 - 3/5*y^5
 - 7/3*y^4 - 87/25*y^2 - 94/15*y + 1) 1]


[x - y 3]

[x^2 + y 4]

[x^3 - y*x + y 5]


[x + Mod(-y, y^2 + 1) 1]

[x + Mod(y, y^2 + 1) 1]


[x + Mod(-y, y^2 + 1) 1]

[x + Mod(y, y^2 + 1) 1]


[x + 1 3]

[2*x + 1 2]

9
[x - y, x + y, x + (-373/2372*y^15 - 69/2372*y^13 + 2205/593*y^11 - 56641/47
44*y^9 + 30823/1186*y^7 - 16782/593*y^5 + 11109/593*y^3 - 15285/2372*y), x +
 (373/2372*y^15 + 69/2372*y^13 - 2205/593*y^11 + 56641/4744*y^9 - 30823/1186
*y^7 + 16782/593*y^5 - 11109/593*y^3 + 15285/2372*y), x^4 + (-231/1186*y^14
- 313/1186*y^12 + 10535/2372*y^10 - 5568/593*y^8 + 10658/593*y^6 - 4133/593*
y^4 + 1731/1186*y^2 - 648/593)*x^2 + (81/1186*y^14 + 181/2372*y^12 - 937/593
*y^10 + 4313/1186*y^8 - 4207/593*y^6 + 5463/1186*y^4 - 1270/593*y^2 + 1883/5
93), x^4 + (-93/1186*y^14 - 49/1186*y^12 + 4457/2372*y^10 - 3058/593*y^8 + 6
170/593*y^6 - 6793/593*y^4 + 8429/1186*y^2 - 2140/593)*x^2 + (81/1186*y^14 +
 181/2372*y^12 - 937/593*y^10 + 4313/1186*y^8 - 4207/593*y^6 + 5463/1186*y^4
 - 1270/593*y^2 + 1883/593), x^4 + (159/2372*y^14 + 281/1186*y^12 - 3437/237
2*y^10 - 217/593*y^8 + 4085/1186*y^6 - 8339/593*y^4 + 12947/1186*y^2 + 146/5
93)*x^2 + (-81/1186*y^14 - 181/2372*y^12 + 937/593*y^10 - 4313/1186*y^8 + 42
07/593*y^6 - 5463/1186*y^4 + 1270/593*y^2 + 489/593)]~
36


I've attatched the files as well as my email screws up the spacing

Thanks
Jason

-------------e=0--------------
[0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
  ***   at top-level: fptest(10007,Mod(1,1
  ***                 ^--------------------
  ***   in function fptest: ...,if(subst(P,x,C[i])==0,0,error("fptest("a","l
  ***                                                   ^--------------------
  ***   user error: fptest(a,10007,Mod(1, 10007)*x^30 + Mod(7812, 10007)*x^28 + Mod(7090, 10007)*x^27 + Mod(7645, 10007)*x^26 + Mod(4110, 10007)*x^25 + Mod(3307, 10007)*x^24 + Mod(5763, 10007)*x^23 + Mod(7900, 10007)*x^22 + Mod(3872, 10007)*x^21 + Mod(8123, 10007)*x^20 + Mod(4076, 10007)*x^19 + Mod(3265, 10007)*x^18 + Mod(3777, 10007)*x^17 + Mod(3398, 10007)*x^16 + Mod(5674, 10007)*x^15 + Mod(4018, 10007)*x^14 + Mod(6820, 10007)*x^13 + Mod(6479, 10007)*x^12 + Mod(984, 10007)*x^11 + Mod(5652, 10007)*x^10 + Mod(1129, 10007)*x^9 + Mod(7573, 10007)*x^8 + Mod(1822, 10007)*x^7 + Mod(837, 10007)*x^6 + Mod(4169, 10007)*x^5 + Mod(4787, 10007)*x^4 + Mod(1616, 10007)*x^3 + Mod(5185, 10007)*x^2 + Mod(2649, 10007)*x + Mod(1483, 10007),Mod(1, 10007)*x^30 + Mod(1, 10007)*x + Mod(2, 10007))
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
-------------e=1--------------
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
-------------e=2--------------
[0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
-------------e=3--------------
[0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
-------------e>=4--------------
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
----------large p---------------
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]

  ***   Warning: new stack size = 16000000 (15.259 Mbytes).
  ***   at top-level: ...17057741307681944498*x^48+1269570586472186440
  ***                                             ^--------------------
  ***   bug in PARI/GP (Segmentation Fault), please report
[x^48 + 12*x^46 + 948*x^44 + 7200*x^42 + 152361*x^40 + 815832*x^38 + 9475380
*x^36 + 44654004*x^34 + 299137536*x^32 + 1335241260*x^30 + 5029216452*x^28 +
 15282825984*x^26 + 37737671337*x^24 + 79579803672*x^22 + 143658877428*x^20
+ 222699104460*x^18 + 303698198961*x^16 + 348787956312*x^14 + 312863646960*x
^12 + 212893847424*x^10 + 111407984496*x^8 + 43762394880*x^6 + 11836253952*x
^4 + 1904684544*x^2 + 136048896]~
[x^48 - 104*x^46 + 4664*x^44 - 122476*x^42 + 2137838*x^40 - 26567700*x^38 +
245144964*x^36 - 1725955872*x^34 + 9441692003*x^32 - 40611588644*x^30 + 1383
56971048*x^28 - 374714866240*x^26 + 807289826646*x^24 - 1380693858220*x^22 +
 1866021172640*x^20 - 1978766780068*x^18 + 1630151673857*x^16 - 102950530102
4*x^14 + 489498952012*x^12 - 170832297056*x^10 + 42133382284*x^8 - 690450713
6*x^6 + 669868016*x^4 - 28899680*x^2 + 16]~
[x^64 - 6384*x^62 + 18261761*x^60 - 31231019568*x^58 + 35925400902280*x^56 -
 29635423138225800*x^54 + 18244443900381139917*x^52 - 8609789775431197305288
*x^50 + 3173715440318358526295493*x^48 - 926253189958924421506713024*x^46 +
216130107574547887816493973792*x^44 - 40600173780591579547211667354912*x^42
+ 6168621134132051706341715912515370*x^40 - 76013943427734813586799139295141
5744*x^38 + 76052955647065900426700689494084391434*x^36 - 617531184355441677
2593815034596603661344*x^34 + 406137362468726923132369140307868626461744*x^3
2 - 21559959594011538596615817038079394912524336*x^30 + 91912201022866770299
8458013093798368177514170*x^28 - 3125139901367629320825272825691018301064633
1888*x^26 + 840017792466500823059432352741453680599704357258*x^24 - 17651171
087877975905738341268934550365403253370816*x^22 + 28594069425611136664663397
3109729677687028435758880*x^20 - 3510138236669349338216717168753743720690200
440965920*x^18 + 31970659328447652254136627725804141367170891626180389*x^16
- 210577501276259853380917767547980363666145460503950096*x^14 + 972633665701
140285445722173171004563157651576027463941*x^12 - 30387893683293821798844628
51269146550248145196800679728*x^10 + 615892792889658839667059092889638564907
5890120481042824*x^8 - 76951350792700777906576936340585240393085356418600329
04*x^6 + 5510293734552561962521574495431679567021222445632508873*x^4 - 19987
78331104544904932086470347413669495129560426038280*x^2 + 2738927449953408337
77347939263771534786080723599733441]~
[x^8 + (-1531349/22619785297920*y^14 + 559416343/10977248747520*y^12 - 25247
684591/1995863408640*y^10 + 103982297321/80715064320*y^8 - 76714376249/12229
55520*y^6 + 111774992009/74186640*y^4 - 26292202739/1686060*y^2 + 1564021312
/54549)*x^6 + (19000645309/37322645741568*y^14 - 96552535399/249482926080*y^
12 + 536960987254321/5488624373760*y^10 - 4194944378809/407651840*y^8 + 5247
721708074653/10089383040*y^6 - 10860233882239/843030*y^4 + 630194805031381/4
636665*y^2 - 796961577488/2755)*x^4 + (-29172555558095/18661322870784*y^14 +
 5062445860141/4248161280*y^12 - 92318704732060083/304923576320*y^10 + 64819
3176741591037/20178766080*y^8 - 8268979844194150111/5044691520*y^6 + 1896413
90835405983/4636665*y^4 - 223283186783900366/515185*y^2 + 8946625664367872/9
405)*x^2 + (41018074761203773/23326653588480*y^14 - 409009053454624763/30492
3576320*y^12 + 117044119382956746521/343039023360*y^10 - 2927890424340695147
291/80715064320*y^8 + 4678762405288733527739/2522345760*y^6 - 38206315875334
762391/824296*y^4 + 2278645136881067296952/4636665*y^2 - 2967400235256172868
48/272745), x^8 + (-16146091/248817638277120*y^14 + 48843457/997931704320*y^
12 - 267727721489/21954497495040*y^10 + 9174250019/7337733120*y^8 - 82557283
8631/13452510720*y^6 + 10012035071/6744240*y^4 - 286258292621/18546660*y^2 +
 140187040/4959)*x^6 + (1723176289/3392967794688*y^14 - 1059662276341/274431
2186880*y^12 + 48715274000629/498965852160*y^10 - 69104247053959/6726255360*
y^8 + 476508837255737/917216640*y^6 - 238731150890387/18546660*y^4 + 5725141
8840649/421515*y^2 - 26269467866896/90915)*x^4 + (-29113232197835/1866132287
0784*y^14 + 1631934517436717/1372156093440*y^12 - 92145545361015407/30492357
6320*y^10 + 323550303468175019/10089383040*y^8 - 8257197132303274099/5044691
520*y^6 + 378862459377191059/9273330*y^4 - 223115572203411654/515185*y^2 + 2
59332507639972352/272745)*x^2 + (40511288842122127/23326653588480*y^14 - 454
626079701495571/343039023360*y^12 + 231451885839430725263/686078046720*y^10
- 60389017554715183063/1681563840*y^8 + 2319938341838169992923/1261172880*y^
6 - 42704610621363386060/927333*y^4 + 2268063476824046324248/4636665*y^2 - 9
8528590048656718144/90915), x^8 + (181633/186613228707840*y^14 - 6489/138601
62560*y^12 - 45386471/5488624373760*y^10 + 178250659/7337733120*y^8 - 259150
32257/10089383040*y^6 + 126060163/1498720*y^4 - 4457231129/4636665*y^2 + 945
4484/4959)*x^6 + (-11712673/1169988894720*y^14 + 2040142733/274431218688*y^1
2 - 93021416071/52522721280*y^10 + 1688263641205/10762008576*y^8 - 890698149
9371/1834433280*y^6 + 14859994827859/296746560*y^4 - 27539242721/421515*y^2
- 10412804654/90915)*x^4 + (-65739851389/2915831698560*y^14 + 8080067384671/
548862437376*y^12 - 615919067582923/228692682240*y^10 + 7139625968047/558581
76*y^8 - 12035658240536459/5044691520*y^6 + 2909281718459207/148373280*y^4 -
 114027481266257/1545555*y^2 + 26994839222348/272745)*x^2 + (351875661187594
13/746452914831360*y^14 - 175794408000830383/5488624373760*y^12 + 2827285081
9680893321/4390899499008*y^10 - 44274526240947916361/107620085760*y^8 + 4488
83263777097920391/40357532160*y^6 - 2698220779200958903/20465280*y^4 + 28100
59487221296958/4636665*y^2 - 81508587175798061/90915), x^8 + (3074279/248817
638277120*y^14 - 127157047/10977248747520*y^12 + 86499746437/21954497495040*
y^10 - 2843672497/4747944960*y^8 + 109487803351/2690502144*y^6 - 89990421617
/74186640*y^4 + 52396810109/3709332*y^2 - 13670576288/272745)*x^6 + (-361812
87281/186613228707840*y^14 + 446988241097/2744312186880*y^12 - 2791944232079
/57774993408*y^10 + 85916604196211/13452510720*y^8 - 3999057670537517/100893
83040*y^6 + 3070998321176/272745*y^4 - 593787426492949/4636665*y^2 + 3730615
4735728/90915)*x^4 + (64206977604683/93306614353920*y^14 - 781339956430123/1
372156093440*y^12 + 151234315443075961/914770728960*y^10 - 42955239518225908
7/20178766080*y^8 + 1307386034599317571/1008938304*y^6 - 168786176615470804/
4636665*y^4 + 127255299187877218/309111*y^2 - 351207138657683456/272745)*x^2
 + (-18530706211095341/23326653588480*y^14 + 1798847927551103773/27443121868
80*y^12 - 16252626552763832987/85759755840*y^10 + 653145277562478020663/2690
5021440*y^8 - 742591759365166575365/504469152*y^6 + 1530656313511014482581/3
7093320*y^4 - 432314981823305760184/927333*y^2 + 131680859590220167616/90915
), x^8 + (3692011/248817638277120*y^14 - 8644699/645720514560*y^12 + 9566281
9793/21954497495040*y^10 - 51309940321/80715064320*y^8 + 113348313275/269050
2144*y^6 - 91885783933/74186640*y^4 + 53072545513/3709332*y^2 - 13761686272/
272745)*x^6 + (-1260198851/6434938920960*y^14 + 450923192773/2744312186880*y
^12 - 53418556138459/1097724874752*y^10 + 43167277559347/6726255360*y^8 - 40
11960654104563/10089383040*y^6 + 209270128922867/18546660*y^4 - 594571618949
411/4636665*y^2 + 37329116221552/90915)*x^4 + (3385038072653/4910874439680*y
^14 - 782493705352717/1372156093440*y^12 + 151409707529271949/914770728960*y
^10 - 214941552982993939/10089383040*y^8 + 1308011461994334679/1008938304*y^
6 - 337682297591934347/9273330*y^4 + 127297072825004458/309111*y^2 - 3513766
27887928064/272745)*x^2 + (-18884591927033039/23326653588480*y^14 + 76223129
967189533/114346341120*y^12 - 131815203944675271599/686078046720*y^10 + 1236
67566697208858791/5044691520*y^8 - 373763922533488473487/252234576*y^6 + 116
278861122276976/2805*y^4 - 433850608352102566712/927333*y^2 + 39636262539108
2293312/272745), x^8 + (1192483/67859355893760*y^14 - 46386469/3659082915840
*y^12 + 5757671539/1995863408640*y^10 - 9816397597/40357532160*y^8 + 3150430
0903/3668866560*y^6 - 6499115329/49457760*y^4 + 1346529151/1686060*y^2 - 127
928104/54549)*x^6 + (-377107043/23326653588480*y^14 + 1193294195/99793170432
*y^12 - 3898646479517/1372156093440*y^10 + 83324546651/326121472*y^8 - 90270
807160297/10089383040*y^6 + 203979621263/1586880*y^4 - 6453521496733/9273330
*y^2 + 3740713928/2755)*x^4 + (-55838547409/2455437219840*y^14 + 81535888034
65/548862437376*y^12 - 1244323126435859/457385364480*y^10 + 2090484253417429
/16143012864*y^8 - 12134063030016491/5044691520*y^6 + 165453933807769/872784
0*y^4 - 93498888636323/1545555*y^2 + 14213849400992/272745)*x^2 + (118472423
79734429/248817638277120*y^14 - 177668994581955997/5488624373760*y^12 + 2861
1858755746067827/4390899499008*y^10 - 134874618655070911657/322860257280*y^8
 + 152404618518667558583/13452510720*y^6 - 4203829397104185047/31236480*y^4
+ 2871462472553380217/4636665*y^2 - 250004160368295232/272745), x^8 + (64582
73/186613228707840*y^14 - 3719191/152461788160*y^12 + 29024100473/5488624373
760*y^10 - 32135840999/80715064320*y^8 + 23297956603/2017876608*y^6 - 211350
0871/16485920*y^4 + 408266029/927333*y^2 - 316594868/272745)*x^6 + (-1110768
17459/373226457415680*y^14 + 18106275053/85759755840*y^12 - 509560146214411/
10977248747520*y^10 + 195109797344857/53810042880*y^8 - 2298235154532701/201
78766080*y^6 + 87655160117623/59349312*y^4 - 32564739888281/4636665*y^2 + 96
8008850258/90915)*x^4 + (83518124920963/93306614353920*y^14 - 91728776500189
/144437483520*y^12 + 42591591062192407/304923576320*y^10 - 88154696029372417
3/80715064320*y^8 + 86773202255301785/252234576*y^6 - 665429827885409833/148
373280*y^4 + 2218979482348247/103037*y^2 - 8878201815711796/272745)*x^2 + (-
246642305436057809/248817638277120*y^14 + 3860028638255806541/5488624373760*
y^12 - 3395672666468176433467/21954497495040*y^10 + 3904062052800146734661/3
22860257280*y^8 - 1024736920815948655471/2690502144*y^6 + 294713641271965414
9381/593493120*y^4 - 22118315020431326588/927333*y^2 + 9833627869746746873/2
72745), x^8 + (38996983/746452914831360*y^14 - 45610777/1219694305280*y^12 +
 184040442823/21954497495040*y^10 - 6929849477/10089383040*y^8 + 19138097395
1/8071506432*y^6 - 6002605199/16485920*y^4 + 8595505583/3709332*y^2 - 144096
9832/272745)*x^6 + (-56627432063/186613228707840*y^14 + 1182421824943/548862
4373760*y^12 - 260406121516747/5488624373760*y^10 + 200262331426543/53810042
880*y^8 - 297480952643123/2522345760*y^6 + 92330503999361/59349312*y^4 - 710
13869059163/9273330*y^2 + 1098811062152/90915)*x^4 + (83613652919587/9330661
4353920*y^14 - 1744863264395929/2744312186880*y^12 + 42642092214501583/30492
3576320*y^10 - 882655350615967327/80715064320*y^8 + 173789516969368021/50446
9152*y^6 - 666541614309814567/148373280*y^4 + 2223816735480437/103037*y^2 -
306949426473056/9405)*x^2 + (-742452471836109913/746452914831360*y^14 + 3873
208908206849579/5488624373760*y^12 - 3407275945042823422193/21954497495040*y
^10 + 1305807204358253484253/107620085760*y^8 - 3084740075192713038863/80715
06432*y^6 + 173955024027274456667/34911360*y^4 - 22194121937138006191/927333
*y^2 + 3289110492164463424/90915)]~
[x^2 + (-872560111/1750783970525184*y^18 - 103213549/291797328420864*y^16 -
4176139757/7204872306688*y^14 - 4416978655/10807308460032*y^12 - 13066402491
575/72949332105216*y^10 - 215588492387/1736888859648*y^8 - 892238763206647/1
09423998157824*y^6 - 89512294102405/18237333026304*y^4 - 18290521472141/2532
96292032*y^2 - 284080544581/14072016224), x^2 + (872560111/1750783970525184*
y^18 - 103213549/291797328420864*y^16 + 4176139757/7204872306688*y^14 - 4416
978655/10807308460032*y^12 + 13066402491575/72949332105216*y^10 - 2155884923
87/1736888859648*y^8 + 892238763206647/109423998157824*y^6 - 89512294102405/
18237333026304*y^4 + 18290521472141/253296292032*y^2 - 284080544581/14072016
224), x^2 + (-10955729699/63028222938906624*y^19 + 1635275/72949332105216*y^
17 - 471200483153/2334378627366912*y^15 + 68074817/2701827115008*y^13 - 1631
44665311371/2626175955787776*y^11 + 10788711995/1519777752192*y^9 - 10805585
507756651/3939263933681664*y^7 + 476476835441/4559333256576*y^5 - 6980744239
9435/3039555504384*y^3 - 4848316573/3518004056*y)*x + (-240314803/1313087977
893888*y^18 + 2177899/20842666315776*y^16 - 10352388961/48632888070144*y^14
+ 650809375/5403654230016*y^12 - 3599460371687/54711999078912*y^10 + 2204400
27581/6079111008768*y^8 - 35061053609941/11723999802624*y^6 + 11692976745493
/9118666513152*y^4 - 1447968833323/63324073008*y^2 - 24126350651/7036008112)
, x^2 + (-4768637809/63028222938906624*y^19 - 116295749/1750783970525184*y^1
7 - 205045552843/2334378627366912*y^15 - 4961577191/64843850760192*y^13 - 10
134809023919/375167993683968*y^11 - 1681738838629/72949332105216*y^9 - 46714
70958726025/3939263933681664*y^7 - 93324108785933/109423998157824*y^5 - 2896
7498100409/3039555504384*y^3 - 245294011997/84432097344*y)*x + (-496188853/1
750783970525184*y^18 + 8017631/48632888070144*y^16 - 338265737/1029267472384
*y^14 + 338938437/1801218076672*y^12 - 7349975872277/72949332105216*y^10 + 1
12959067435/2026370336256*y^8 - 471966038240125/109423998157824*y^6 + 558940
2408335/3039555504384*y^4 - 8241706316231/253296292032*y^2 + 1651761915/1005
144016), x^2 + (-550587305/9004031848415232*y^19 - 8520109/250111995789312*y
^17 - 165341932517/2334378627366912*y^15 - 2570761369/64843850760192*y^13 -
56851147336519/2626175955787776*y^11 - 895619999627/72949332105216*y^9 - 356
0831957116295/3939263933681664*y^7 - 62314525927891/109423998157824*y^5 - 18
560457963575/3039555504384*y^3 - 371119778467/84432097344*y)*x + (496188853/
1750783970525184*y^18 + 8017631/48632888070144*y^16 + 338265737/102926747238
4*y^14 + 338938437/1801218076672*y^12 + 7349975872277/72949332105216*y^10 +
112959067435/2026370336256*y^8 + 471966038240125/109423998157824*y^6 + 55894
02408335/3039555504384*y^4 + 8241706316231/253296292032*y^2 + 1651761915/100
5144016), x^2 + (-21512837/1313087977893888*y^19 - 102134281/175078397052518
4*y^17 - 932348591/48632888070144*y^15 - 4407863827/64843850760192*y^13 - 16
4654686001/27355999539456*y^11 - 1540975821641/72949332105216*y^9 - 25030796
719091/82067998618368*y^7 - 110013175606897/109423998157824*y^5 - 7018393859
/2261574036*y^3 - 296247471259/28144032448*y)*x + (240314803/131308797789388
8*y^18 + 2177899/20842666315776*y^16 + 10352388961/48632888070144*y^14 + 650
809375/5403654230016*y^12 + 3599460371687/54711999078912*y^10 + 220440027581
/6079111008768*y^8 + 35061053609941/11723999802624*y^6 + 11692976745493/9118
666513152*y^4 + 1447968833323/63324073008*y^2 - 24126350651/7036008112), x^2
 + (21512837/1313087977893888*y^19 + 102134281/1750783970525184*y^17 + 93234
8591/48632888070144*y^15 + 4407863827/64843850760192*y^13 + 164654686001/273
55999539456*y^11 + 1540975821641/72949332105216*y^9 + 25030796719091/8206799
8618368*y^7 + 110013175606897/109423998157824*y^5 + 7018393859/2261574036*y^
3 + 296247471259/28144032448*y)*x + (240314803/1313087977893888*y^18 + 21778
99/20842666315776*y^16 + 10352388961/48632888070144*y^14 + 650809375/5403654
230016*y^12 + 3599460371687/54711999078912*y^10 + 220440027581/6079111008768
*y^8 + 35061053609941/11723999802624*y^6 + 11692976745493/9118666513152*y^4
+ 1447968833323/63324073008*y^2 - 24126350651/7036008112), x^2 + (550587305/
9004031848415232*y^19 + 8520109/250111995789312*y^17 + 165341932517/23343786
27366912*y^15 + 2570761369/64843850760192*y^13 + 56851147336519/262617595578
7776*y^11 + 895619999627/72949332105216*y^9 + 3560831957116295/3939263933681
664*y^7 + 62314525927891/109423998157824*y^5 + 18560457963575/3039555504384*
y^3 + 371119778467/84432097344*y)*x + (496188853/1750783970525184*y^18 + 801
7631/48632888070144*y^16 + 338265737/1029267472384*y^14 + 338938437/18012180
76672*y^12 + 7349975872277/72949332105216*y^10 + 112959067435/2026370336256*
y^8 + 471966038240125/109423998157824*y^6 + 5589402408335/3039555504384*y^4
+ 8241706316231/253296292032*y^2 + 1651761915/1005144016), x^2 + (4768637809
/63028222938906624*y^19 + 116295749/1750783970525184*y^17 + 205045552843/233
4378627366912*y^15 + 4961577191/64843850760192*y^13 + 10134809023919/3751679
93683968*y^11 + 1681738838629/72949332105216*y^9 + 4671470958726025/39392639
33681664*y^7 + 93324108785933/109423998157824*y^5 + 28967498100409/303955550
4384*y^3 + 245294011997/84432097344*y)*x + (-496188853/1750783970525184*y^18
 + 8017631/48632888070144*y^16 - 338265737/1029267472384*y^14 + 338938437/18
01218076672*y^12 - 7349975872277/72949332105216*y^10 + 112959067435/20263703
36256*y^8 - 471966038240125/109423998157824*y^6 + 5589402408335/303955550438
4*y^4 - 8241706316231/253296292032*y^2 + 1651761915/1005144016), x^2 + (1095
5729699/63028222938906624*y^19 - 1635275/72949332105216*y^17 + 471200483153/
2334378627366912*y^15 - 68074817/2701827115008*y^13 + 163144665311371/262617
5955787776*y^11 - 10788711995/1519777752192*y^9 + 10805585507756651/39392639
33681664*y^7 - 476476835441/4559333256576*y^5 + 69807442399435/3039555504384
*y^3 + 4848316573/3518004056*y)*x + (-240314803/1313087977893888*y^18 + 2177
899/20842666315776*y^16 - 10352388961/48632888070144*y^14 + 650809375/540365
4230016*y^12 - 3599460371687/54711999078912*y^10 + 220440027581/607911100876
8*y^8 - 35061053609941/11723999802624*y^6 + 11692976745493/9118666513152*y^4
 - 1447968833323/63324073008*y^2 - 24126350651/7036008112)]~
[x^24 + 69]~
[x^7 - 2*y*x^6 + y^2*x^5 - 28*x^3 + 4*y^2*x + 16/7*y^3, x^7 + 2*y*x^6 + y^2*
x^5 - 28*x^3 + 4*y^2*x - 16/7*y^3, x^7 - 2/7*y^3*x^6 - y^2*x^5 - 28*x^3 - 4*
y^2*x + 16*y, x^7 + 2/7*y^3*x^6 - y^2*x^5 - 28*x^3 - 4*y^2*x - 16*y]~
[x^2 + (-12035/386*y^15 + 8337/386*y^14 + 566267/772*y^13 - 392327/772*y^12
- 4449119/772*y^11 + 3083327/772*y^10 + 7320431/386*y^9 - 2537889/193*y^8 -
9975849/386*y^7 + 3461954/193*y^6 + 10712085/772*y^5 - 7445851/772*y^4 - 190
5321/772*y^3 + 1322419/772*y^2 + 25961/386*y - 17735/386), x^2 + (-7909/772*
y^15 - 1879/772*y^14 + 93535/386*y^13 + 44921/772*y^12 - 371330/193*y^11 - 9
0974/193*y^10 + 4991389/772*y^9 + 317305/193*y^8 - 7127131/772*y^7 - 1947859
/772*y^6 + 2110371/386*y^5 + 1295397/772*y^4 - 224091/193*y^3 - 160649/386*y
^2 + 19739/772*y + 4029/386), x^2 + (7909/772*y^15 - 1879/772*y^14 - 93535/3
86*y^13 + 44921/772*y^12 + 371330/193*y^11 - 90974/193*y^10 - 4991389/772*y^
9 + 317305/193*y^8 + 7127131/772*y^7 - 1947859/772*y^6 - 2110371/386*y^5 + 1
295397/772*y^4 + 224091/193*y^3 - 160649/386*y^2 - 19739/772*y + 4029/386),
x^2 + (12035/386*y^15 + 8337/386*y^14 - 566267/772*y^13 - 392327/772*y^12 +
4449119/772*y^11 + 3083327/772*y^10 - 7320431/386*y^9 - 2537889/193*y^8 + 99
75849/386*y^7 + 3461954/193*y^6 - 10712085/772*y^5 - 7445851/772*y^4 + 19053
21/772*y^3 + 1322419/772*y^2 - 25961/386*y - 17735/386), x^4 + (541/386*y^14
 - 12373/386*y^12 + 46089/193*y^10 - 137477/193*y^8 + 311469/386*y^6 - 14604
9/386*y^4 + 12308/193*y^2 - 686/193)*x^2 + (-765/772*y^14 + 8659/386*y^12 -
126401/772*y^10 + 90322/193*y^8 - 369401/772*y^6 + 38430/193*y^4 - 23865/772
*y^2 + 553/386), x^4 + (2052/193*y^14 - 96741/386*y^12 + 763031/386*y^10 - 1
265475/193*y^8 + 1759047/193*y^6 - 2000657/386*y^4 + 430823/386*y^2 - 13920/
193)*x^2 + (-28293/772*y^14 + 166359/193*y^12 - 5225705/772*y^10 + 4295689/1
93*y^8 - 23400329/772*y^6 + 6320077/386*y^4 - 2369321/772*y^2 + 45727/386),
x^4 + (-5155/193*y^15 - 15405/772*y^14 + 241769/386*y^13 + 180325/386*y^12 -
 944044/193*y^11 - 2807115/772*y^10 + 3072303/193*y^9 + 2267740/193*y^8 - 41
02400/193*y^7 - 11896515/772*y^6 + 4334901/386*y^5 + 2988835/386*y^4 - 43236
3/193*y^3 - 1014845/772*y^2 + 23130/193*y + 3602/193)*x^2 + (14406/193*y^15
+ 38997/772*y^14 - 338423/193*y^13 - 455769/386*y^12 + 2652454/193*y^11 + 70
72893/772*y^10 - 8698762/193*y^9 - 11354279/386*y^8 + 23683093/386*y^7 + 293
19183/772*y^6 - 13169735/386*y^5 - 7082173/386*y^4 + 3014249/386*y^3 + 21677
71/772*y^2 - 244333/386*y + 2027/386), x^4 + (-4705/193*y^15 - 4035/772*y^14
 + 220691/386*y^13 + 47935/386*y^12 - 1722817/386*y^11 - 767525/772*y^10 + 5
588151/386*y^9 + 654320/193*y^8 - 3658160/193*y^7 - 3833005/772*y^6 + 341354
7/386*y^5 + 1159745/386*y^4 - 194425/386*y^3 - 441715/772*y^2 - 129025/386*y
 - 9518/193)*x^2 + (69045/772*y^15 + 17477/772*y^14 - 1634153/772*y^13 - 208
093/386*y^12 + 12990539/772*y^11 + 3347413/772*y^10 - 43749763/772*y^9 - 576
5109/386*y^8 + 62724993/772*y^7 + 17320883/772*y^6 - 37352513/772*y^5 - 5621
249/386*y^4 + 7969659/772*y^3 + 2742727/772*y^2 - 182407/772*y - 34007/386),
 x^4 + (4705/193*y^15 - 4035/772*y^14 - 220691/386*y^13 + 47935/386*y^12 + 1
722817/386*y^11 - 767525/772*y^10 - 5588151/386*y^9 + 654320/193*y^8 + 36581
60/193*y^7 - 3833005/772*y^6 - 3413547/386*y^5 + 1159745/386*y^4 + 194425/38
6*y^3 - 441715/772*y^2 + 129025/386*y - 9518/193)*x^2 + (-69045/772*y^15 + 1
7477/772*y^14 + 1634153/772*y^13 - 208093/386*y^12 - 12990539/772*y^11 + 334
7413/772*y^10 + 43749763/772*y^9 - 5765109/386*y^8 - 62724993/772*y^7 + 1732
0883/772*y^6 + 37352513/772*y^5 - 5621249/386*y^4 - 7969659/772*y^3 + 274272
7/772*y^2 + 182407/772*y - 34007/386), x^4 + (5155/193*y^15 - 15405/772*y^14
 - 241769/386*y^13 + 180325/386*y^12 + 944044/193*y^11 - 2807115/772*y^10 -
3072303/193*y^9 + 2267740/193*y^8 + 4102400/193*y^7 - 11896515/772*y^6 - 433
4901/386*y^5 + 2988835/386*y^4 + 432363/193*y^3 - 1014845/772*y^2 - 23130/19
3*y + 3602/193)*x^2 + (-14406/193*y^15 + 38997/772*y^14 + 338423/193*y^13 -
455769/386*y^12 - 2652454/193*y^11 + 7072893/772*y^10 + 8698762/193*y^9 - 11
354279/386*y^8 - 23683093/386*y^7 + 29319183/772*y^6 + 13169735/386*y^5 - 70
82173/386*y^4 - 3014249/386*y^3 + 2167771/772*y^2 + 244333/386*y + 2027/386)
]~
[x^4 - 4*y*x^3 + (-1/2*y^14 + 1/2*y^10 - 7/2*y^6 + 15/2*y^2)*x^2 + (y^15 - y
^11 + 7*y^7 - 7*y^3)*x + (1/2*y^12 - y^8 + 5/2*y^4 - 2), x^4 + 4*y*x^3 + (-1
/2*y^14 + 1/2*y^10 - 7/2*y^6 + 15/2*y^2)*x^2 + (-y^15 + y^11 - 7*y^7 + 7*y^3
)*x + (1/2*y^12 - y^8 + 5/2*y^4 - 2), x^4 + (-2*y^11 - 10*y^3)*x^3 + (9/2*y^
14 - 1/2*y^10 + 55/2*y^6 - 7/2*y^2)*x^2 + (3*y^13 + y^9 + 17*y^5 + 3*y)*x +
(1/2*y^12 + y^8 + 5/2*y^4 + 4), x^4 + (2*y^11 + 10*y^3)*x^3 + (9/2*y^14 - 1/
2*y^10 + 55/2*y^6 - 7/2*y^2)*x^2 + (-3*y^13 - y^9 - 17*y^5 - 3*y)*x + (1/2*y
^12 + y^8 + 5/2*y^4 + 4), x^4 + (-2*y^13 - 10*y^5)*x^3 + (1/2*y^14 - 1/2*y^1
0 + 7/2*y^6 - 15/2*y^2)*x^2 + (3*y^15 + y^11 + 17*y^7 + 3*y^3)*x + (1/2*y^12
 - y^8 + 5/2*y^4 - 2), x^4 + (-y^13 - y^9 - 7*y^5 - 3*y)*x^3 + (7/2*y^14 - 1
/2*y^10 + 37/2*y^6 - 3/2*y^2)*x^2 + (-y^15 + y^11 - 7*y^7 + 7*y^3)*x + (-1/2
*y^12 - y^8 - 5/2*y^4 - 2), x^4 + (-y^13 + y^9 - 7*y^5 + 3*y)*x^3 + (-7/2*y^
14 + 1/2*y^10 - 37/2*y^6 + 3/2*y^2)*x^2 + (-3*y^15 - y^11 - 17*y^7 - 3*y^3)*
x + (-1/2*y^12 - y^8 - 5/2*y^4 - 2), x^4 + (y^13 - y^9 + 7*y^5 - 3*y)*x^3 +
(-7/2*y^14 + 1/2*y^10 - 37/2*y^6 + 3/2*y^2)*x^2 + (3*y^15 + y^11 + 17*y^7 +
3*y^3)*x + (-1/2*y^12 - y^8 - 5/2*y^4 - 2), x^4 + (y^13 + y^9 + 7*y^5 + 3*y)
*x^3 + (7/2*y^14 - 1/2*y^10 + 37/2*y^6 - 3/2*y^2)*x^2 + (y^15 - y^11 + 7*y^7
 - 7*y^3)*x + (-1/2*y^12 - y^8 - 5/2*y^4 - 2), x^4 + (2*y^13 + 10*y^5)*x^3 +
 (1/2*y^14 - 1/2*y^10 + 7/2*y^6 - 15/2*y^2)*x^2 + (-3*y^15 - y^11 - 17*y^7 -
 3*y^3)*x + (1/2*y^12 - y^8 + 5/2*y^4 - 2), x^4 + (-4*y^15 - 24*y^7)*x^3 + (
-9/2*y^14 + 1/2*y^10 - 55/2*y^6 + 7/2*y^2)*x^2 + (-y^13 + y^9 - 7*y^5 + 7*y)
*x + (1/2*y^12 + y^8 + 5/2*y^4 + 4), x^4 + (-3*y^15 - y^11 - 17*y^7 - 7*y^3)
*x^3 + (3/2*y^14 - 5/2*y^10 + 17/2*y^6 - 23/2*y^2)*x^2 + (3*y^13 + y^9 + 17*
y^5 + 3*y)*x + (-1/2*y^12 + y^8 - 5/2*y^4 + 4), x^4 + (-3*y^15 + y^11 - 17*y
^7 + 7*y^3)*x^3 + (-3/2*y^14 + 5/2*y^10 - 17/2*y^6 + 23/2*y^2)*x^2 + (-y^13
+ y^9 - 7*y^5 + 7*y)*x + (-1/2*y^12 + y^8 - 5/2*y^4 + 4), x^4 + (3*y^15 - y^
11 + 17*y^7 - 7*y^3)*x^3 + (-3/2*y^14 + 5/2*y^10 - 17/2*y^6 + 23/2*y^2)*x^2
+ (y^13 - y^9 + 7*y^5 - 7*y)*x + (-1/2*y^12 + y^8 - 5/2*y^4 + 4), x^4 + (3*y
^15 + y^11 + 17*y^7 + 7*y^3)*x^3 + (3/2*y^14 - 5/2*y^10 + 17/2*y^6 - 23/2*y^
2)*x^2 + (-3*y^13 - y^9 - 17*y^5 - 3*y)*x + (-1/2*y^12 + y^8 - 5/2*y^4 + 4),
 x^4 + (4*y^15 + 24*y^7)*x^3 + (-9/2*y^14 + 1/2*y^10 - 55/2*y^6 + 7/2*y^2)*x
^2 + (y^13 - y^9 + 7*y^5 - 7*y)*x + (1/2*y^12 + y^8 + 5/2*y^4 + 4)]~
[x^64 + 192*x^62 + 17568*x^60 + 1019520*x^58 + 42131676*x^56 + 1319651424*x^
54 + 32559096528*x^52 + 649228312512*x^50 + 10651553826426*x^48 + 1456394385
52224*x^46 + 1674922821206832*x^44 + 16307859539653056*x^42 + 13502367773216
7696*x^40 + 953248899971965824*x^38 + 5745239175305568960*x^36 + 29556064271
185194240*x^34 + 129595725382952883843*x^32 + 483002100692576612640*x^30 + 1
523870714370199019760*x^28 + 4047489983524093705152*x^26 + 89858128286488620
19536*x^24 + 16525310345394167002752*x^22 + 24893927149975603242048*x^20 + 3
0294355815129821928192*x^18 + 29274561574319887883226*x^16 + 219878017711043
40121824*x^14 + 12494344840480632094992*x^12 + 5187763623118143696192*x^10 +
 1502211081063677383836*x^8 + 283567347515314680480*x^6 + 311461554388845258
72*x^4 + 1543354925530003776*x^2 + 8057044481403681]~
[x + (-7/16*y^29 - 5/32*y^25 - 97/8*y^21 - 139/32*y^17 - 435/16*y^13 - 323/3
2*y^9 - 21/4*y^5 - 77/32*y), x + (-11/32*y^29 - 1/8*y^25 - 309/32*y^21 - 55/
16*y^17 - 797/32*y^13 - 7*y^9 - 355/32*y^5 + 9/16*y), x + (11/32*y^29 + 1/8*
y^25 + 309/32*y^21 + 55/16*y^17 + 797/32*y^13 + 7*y^9 + 355/32*y^5 - 9/16*y)
, x + (7/16*y^29 + 5/32*y^25 + 97/8*y^21 + 139/32*y^17 + 435/16*y^13 + 323/3
2*y^9 + 21/4*y^5 + 77/32*y), x + (-21/32*y^31 - 587/32*y^23 + 1/16*y^19 - 14
43/32*y^15 + 13/8*y^11 - 541/32*y^7 + 21/16*y^3), x + (-5/16*y^31 - 5/32*y^2
7 - 35/4*y^23 - 139/32*y^19 - 349/16*y^15 - 323/32*y^11 - 57/8*y^7 - 109/32*
y^3), x + (5/16*y^31 + 5/32*y^27 + 35/4*y^23 + 139/32*y^19 + 349/16*y^15 + 3
23/32*y^11 + 57/8*y^7 + 109/32*y^3), x + (21/32*y^31 + 587/32*y^23 - 1/16*y^
19 + 1443/32*y^15 - 13/8*y^11 + 541/32*y^7 - 21/16*y^3), x^4 + (-13/4*y^24 -
 91*y^16 - 909/4*y^8 - 169/2), x^4 + (-1/4*y^24 - 7*y^16 - 57/4*y^8 - 1/2)]~
[8711099/70204123*y^14 - 3396450/70204123*y^13 - 230089978/70204123*y^12 + 7
1459644/70204123*y^11 + 2039293754/70204123*y^10 - 522502724/70204123*y^9 -
7578045032/70204123*y^8 + 136410216/6382193*y^7 + 11598831422/70204123*y^6 -
 1582740050/70204123*y^5 - 6466526698/70204123*y^4 + 712163508/70204123*y^3
+ 865017354/70204123*y^2 - 11706800/70204123*y + 7921687/70204123]
[y, -123209112482/559553426209*y^11 - 236161397417/559553426209*y^10 + 52225
05497467/559553426209*y^9 + 7627164004768/559553426209*y^8 - 68684785347690/
559553426209*y^7 - 98327585435469/559553426209*y^6 + 334508906676131/5595534
26209*y^5 + 508054669424553/559553426209*y^4 - 499853398148011/559553426209*
y^3 - 780815391953932/559553426209*y^2 + 222263541657120/559553426209*y + 69
71304961116/11905392047, -176690268281/1119106852418*y^11 - 641210922141/223
8213704836*y^10 + 3748626002639/559553426209*y^9 + 20339566453621/2238213704
836*y^8 - 49228932581896/559553426209*y^7 - 262748733927015/2238213704836*y^
6 + 958463738831775/2238213704836*y^5 + 1369627950576313/2238213704836*y^4 -
 1439017147702531/2238213704836*y^3 - 2121737405539667/2238213704836*y^2 + 6
35967216650047/2238213704836*y + 19194169373855/47621568188, -68393259315/11
19106852418*y^11 - 143769468519/1119106852418*y^10 + 2890095068571/111910685
2418*y^9 + 4773072512403/1119106852418*y^8 - 38008182942981/1119106852418*y^
7 - 30714902181882/559553426209*y^6 + 92097658832451/559553426209*y^5 + 1570
24332093039/559553426209*y^4 - 134021586074273/559553426209*y^3 - 2404025913
19962/559553426209*y^2 + 58771325675463/559553426209*y + 4293056883813/23810
784094, -26071422312/559553426209*y^11 - 45431912634/559553426209*y^10 + 110
4266993455/559553426209*y^9 + 1416464504385/559553426209*y^8 - 1440773326506
3/559553426209*y^7 - 18255285107310/559553426209*y^6 + 69227537719650/559553
426209*y^5 + 94877281869432/559553426209*y^4 - 100252280561328/559553426209*
y^3 - 143113043669604/559553426209*y^2 + 41073387952215/559553426209*y + 128
0223975324/11905392047, -16163557893/1119106852418*y^11 - 64908872441/223821
3704836*y^10 + 685522417391/1119106852418*y^9 + 2130089587017/2238213704836*
y^8 - 4531221497592/559553426209*y^7 - 27490869176393/2238213704836*y^6 + 89
319951302323/2238213704836*y^5 + 141812611520979/2238213704836*y^4 - 1390748
45426265/2238213704836*y^3 - 220351698213059/2238213704836*y^2 + 72725691407
453/2238213704836*y + 2005246545505/47621568188, 9609209575/2238213704836*y^
11 + 25387385845/559553426209*y^10 - 383645617715/2238213704836*y^9 - 208031
3216715/1119106852418*y^8 + 6033944605985/2238213704836*y^7 + 53081257123225
/2238213704836*y^6 - 33975296862719/2238213704836*y^5 - 255089221619815/2238
213704836*y^4 + 47525895430055/2238213704836*y^3 + 398412597852115/223821370
4836*y^2 - 30904572571305/2238213704836*y - 890343326655/11905392047, 522817
60559/2238213704836*y^11 + 14008471016/559553426209*y^10 - 2227004686929/223
8213704836*y^9 - 336992949231/559553426209*y^8 + 28726688292139/223821370483
6*y^7 + 17917872421753/2238213704836*y^6 - 137376604832131/2238213704836*y^5
 - 105226511459813/2238213704836*y^4 + 207364021022899/2238213704836*y^3 + 1
70780696189363/2238213704836*y^2 - 84647799904795/2238213704836*y - 44402858
6648/11905392047, 106577461171/2238213704836*y^11 + 236835106673/22382137048
36*y^10 - 4499294217001/2238213704836*y^9 - 8003166405117/2238213704836*y^8
+ 59306288158573/2238213704836*y^7 + 25836461923340/559553426209*y^6 - 72093
028502701/559553426209*y^5 - 264172433170095/1119106852418*y^4 + 21076054618
3675/1119106852418*y^3 + 203797855675948/559553426209*y^2 - 94758333876515/1
119106852418*y - 7259569285101/47621568188, 175225361355/2238213704836*y^11
+ 367043962553/2238213704836*y^10 - 7404910374729/2238213704836*y^9 - 121761
49738023/2238213704836*y^8 + 97370244941211/2238213704836*y^7 + 784051152149
35/1119106852418*y^6 - 117954333048705/559553426209*y^5 - 401317015113399/11
19106852418*y^4 + 171430955359538/559553426209*y^3 + 307321505664217/5595534
26209*y^2 - 148696095658675/1119106852418*y - 10923889360787/47621568188, 18
4746017659/1119106852418*y^11 + 182236626649/559553426209*y^10 - 78287971960
99/1119106852418*y^9 - 5941166792275/559553426209*y^8 + 103098612489213/1119
106852418*y^7 + 153130332934513/1119106852418*y^6 - 502641225832273/11191068
52418*y^5 - 789225051651289/1119106852418*y^4 + 749848149414089/111910685241
8*y^3 + 1216440499421323/1119106852418*y^2 - 334431034048089/1119106852418*y
 - 5475399545957/11905392047, 421608384/2321798449*y^11 + 1306271719/4643596
898*y^10 - 17925705039/2321798449*y^9 - 38874965891/4643596898*y^8 + 2342247
56928/2321798449*y^7 + 502105436649/4643596898*y^6 - 2270324281083/464359689
8*y^5 - 2658757184017/4643596898*y^4 + 3428983243105/4643596898*y^3 + 405157
9820087/4643596898*y^2 - 1510458181035/4643596898*y - 36134279083/98799934]

[x + (-y + 1) 1]

[x^2 + (y + 2)*x + (y^2 + y + 1) 1]

[x^2 + (y + 2)*x + (1/25*y^8 - 3/5*y^5 - 87/25*y^2 + y + 1) 1]

[x^2 + (-2/15*y^7 + 7/3*y^4 + 79/15*y + 2)*x + (1/25*y^8 - 2/15*y^7 - 3/5*y^
5 + 7/3*y^4 - 87/25*y^2 + 79/15*y + 1) 1]

[x^2 + (2/15*y^7 - 7/3*y^4 - 94/15*y + 2)*x + (1/25*y^8 + 2/15*y^7 - 3/5*y^5
 - 7/3*y^4 - 87/25*y^2 - 94/15*y + 1) 1]


[x - y 3]

[x^2 + y 4]

[x^3 - y*x + y 5]


[x + Mod(-y, y^2 + 1) 1]

[x + Mod(y, y^2 + 1) 1]


[x + Mod(-y, y^2 + 1) 1]

[x + Mod(y, y^2 + 1) 1]


[x + 1 3]

[2*x + 1 2]

9
[x - y, x + y, x + (-373/2372*y^15 - 69/2372*y^13 + 2205/593*y^11 - 56641/47
44*y^9 + 30823/1186*y^7 - 16782/593*y^5 + 11109/593*y^3 - 15285/2372*y), x +
 (373/2372*y^15 + 69/2372*y^13 - 2205/593*y^11 + 56641/4744*y^9 - 30823/1186
*y^7 + 16782/593*y^5 - 11109/593*y^3 + 15285/2372*y), x^4 + (-231/1186*y^14
- 313/1186*y^12 + 10535/2372*y^10 - 5568/593*y^8 + 10658/593*y^6 - 4133/593*
y^4 + 1731/1186*y^2 - 648/593)*x^2 + (81/1186*y^14 + 181/2372*y^12 - 937/593
*y^10 + 4313/1186*y^8 - 4207/593*y^6 + 5463/1186*y^4 - 1270/593*y^2 + 1883/5
93), x^4 + (-93/1186*y^14 - 49/1186*y^12 + 4457/2372*y^10 - 3058/593*y^8 + 6
170/593*y^6 - 6793/593*y^4 + 8429/1186*y^2 - 2140/593)*x^2 + (81/1186*y^14 +
 181/2372*y^12 - 937/593*y^10 + 4313/1186*y^8 - 4207/593*y^6 + 5463/1186*y^4
 - 1270/593*y^2 + 1883/593), x^4 + (159/2372*y^14 + 281/1186*y^12 - 3437/237
2*y^10 - 217/593*y^8 + 4085/1186*y^6 - 8339/593*y^4 + 12947/1186*y^2 + 146/5
93)*x^2 + (-81/1186*y^14 - 181/2372*y^12 + 937/593*y^10 - 4313/1186*y^8 + 42
07/593*y^6 - 5463/1186*y^4 + 1270/593*y^2 + 489/593)]~
36

[Bill Allombert-3]

On Sat, Jul 04, 2009 at 12:24:00PM +0100, Jason Moxham wrote:

> On Thursday 02 July 2009 23:11:31 Jason Moxham wrote:
> I manage to get rid of a few more bugs, and I'm left with these
>
> polred,rnf,rnfkummer all have this error
>
> x^7 + Mod(7*y - 7, y^2 - y - 1)*x^6 + Mod(-21*y + 28, y^2 - y - 1)*x^5 - 35*
> x^4 + Mod(35*y - 49, y^2 - y - 1)*x^3 + Mod(-7*y + 84, y^2 - y - 1)*x^2 + Mo
> d(-14*y + 21, y^2 - y - 1)*x + Mod(-y - 43, y^2 - y - 1)
>   ***   at top-level: rnfpolredabs(nfinit(
>   ***                 ^--------------------
>   *** rnfpolredabs: could not open requested file ./MPQS.gpa/FREL.

There might a problem with the way MPQS creates temporary files.
Try factor(2^128+1)
The code to create temporary files is in src/language/es.c.
Maybe we use a wrong directory or a wrong filename.

>   ***   at top-level: fptest(10007,Mod(1,1
>   ***                 ^--------------------
>   ***   in function fptest: ...,if(subst(P,x,C[i])==0,0,error("fptest("a","l
>   ***                                                   ^--------------------
>   ***   user error: fptest(a,10007,Mod(1, 10007)*x^30 + Mod(7812, 10007)*x^28
> + Mod(7090, 10007)*x^27 + Mod(7645, 10007)*x^26 + Mod(4110, 10007)*x^25 +
> Mod(3307, 10007)*x^24 + Mod(5763, 10007)*x^23 + Mod(7900, 10007)*x^22 +
> Mod(3872, 10007)*x^21 + Mod(8123, 10007)*x^20 + Mod(4076, 10007)*x^19 +
> Mod(3265, 10007)*x^18 + Mod(3777, 10007)*x^17 + Mod(3398, 10007)*x^16 +
> Mod(5674, 10007)*x^15 + Mod(4018, 10007)*x^14 + Mod(6820, 10007)*x^13 +
> Mod(6479, 10007)*x^12 + Mod(984, 10007)*x^11 + Mod(5652, 10007)*x^10 +
> Mod(1129, 10007)*x^9 + Mod(7573, 10007)*x^8 + Mod(1822, 10007)*x^7 + Mod(837,
> 10007)*x^6 + Mod(4169, 10007)*x^5 + Mod(4787, 10007)*x^4 + Mod(1616,
> 10007)*x^3 + Mod(5185, 10007)*x^2 + Mod(2649, 10007)*x + Mod(1483,
> 10007),Mod(1, 10007)*x^30 + Mod(1, 10007)*x + Mod(2, 10007))

factorff returns a wrong result.

Try the following command:
{
factorff(x^30 + 7812*x^28 + 7090*x^27 + 7645*x^26 + 4110*x^25 + 3307*x^24 + 5763*x^23 + 7900*x^22 + 3872*x^21 + 8123*x^20 + 4076*x^19 + 3265*x^18 + 3777*x^17 + 3398*x^16 + 5674*x^15 + 4018*x^14 + 6820*x^13 + 6479*x^12 + 984*x^11 + 5652*x^10 + 1129*x^9 + 7573*x^8 + 1822*x^7 + 837*x^6 + 4169*x^5 + 4787*x^4 + 1616*x^3 + 5185*x^2 + 2649*x + 1483, 10007, a^30 + a + 2)
}

> and nffactor
>
>   ***   Warning: new stack size = 16000000 (15.259 Mbytes).
>   ***   at top-level: ...17057741307681944498*x^48+1269570586472186440
>   ***                                             ^--------------------
>   ***   bug in PARI/GP (Segmentation Fault), please report

This one might be a consequence of the previous problem. (nffactor
use factorff internally).

Cheers,
Bill.

[Bill Allombert-3]

On Sat, Jul 04, 2009 at 02:36:17AM +0100, Jason Moxham wrote:
>
> On MSVC INLINE was only defined as __inline not __inline static , that  
> solves the bnfinit() ,

Excellent!

> doesn't help with rest though.....

Be very careful with files in Odos because they were meant for DOS
and djgpp rather than Windows, so they might be completly off.

Cheers,
Bill.

[Jason Moxham]

----- Original Message -----
From: "Bill Allombert" <Bill.Allombert@...>
To: <pari-dev@...>
Sent: Saturday, July 04, 2009 1:19 PM
Subject: Re: Some bugs?

> On Sat, Jul 04, 2009 at 12:24:00PM +0100, Jason Moxham wrote:
>> On Thursday 02 July 2009 23:11:31 Jason Moxham wrote:
>> I manage to get rid of a few more bugs, and I'm left with these
>>
>> polred,rnf,rnfkummer all have this error
>>
>> x^7 + Mod(7*y - 7, y^2 - y - 1)*x^6 + Mod(-21*y + 28, y^2 - y - 1)*x^5 -
>> 35*
>> x^4 + Mod(35*y - 49, y^2 - y - 1)*x^3 + Mod(-7*y + 84, y^2 - y - 1)*x^2 +
>> Mo
>> d(-14*y + 21, y^2 - y - 1)*x + Mod(-y - 43, y^2 - y - 1)
>>   ***   at top-level: rnfpolredabs(nfinit(
>>   ***                 ^--------------------
>>   *** rnfpolredabs: could not open requested file ./MPQS.gpa/FREL.
>
> There might a problem with the way MPQS creates temporary files.
> Try factor(2^128+1)
> The code to create temporary files is in src/language/es.c.
> Maybe we use a wrong directory or a wrong filename.
>

Yes it is , mpqs assumes the function below creates the directory as well as
testing existence

line 4175 in language/es.c
static int
pari_dir_exists(const char *s) { return 0; }

change to

static int
pari_dir_exists(const char *s) { return mkdir(s); }

So this is the same as the UNIX one but without the mode
So I suppose we will need this for MSVC and MinGW

this fixes polred,rnf,rnfkummer  tests for Win32 MSVC

Jason