Bill Allombert on Thu, 12 Dec 2002 15:47:05 +0100


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

Preparing bugfix release 2.1.5


Hello developers,

As some of you may have gathered form the CVS, I am preparing behind the scene
a 2.1.5 bug fix release, so if you are aware of bugs in 2.1.4 that are not
yet fixed in the CVS, please report them to me.

I try to fix 2.2.4 F32:
 32- polredabs fails to reduce
   x^8-2*x^7-34*x^6+78*x^5+265*x^4-628*x^3-389*x^2+1237*x-449
   [typo in chk_gen_init: skipfirst not initialized properly]

I have some questions:
1) I do not thing the bracketed comment is accurate. I think the skipfirst
problem was fixed a long time ago in 2.2.3 F28 and 2.1.4 F7.

2) The patch responsible of the fix is as follow:

--- src/basemath/base1.c	2002/07/15 13:29:58	1.99
+++ src/basemath/base1.c	2002/08/03 15:37:29	1.100
@@ -1771,7 +1771,8 @@
     if (prev && degpol(prev) < N && !gegal(prev,P))
     {
       if (degpol(prev) * degpol(P) > 64) continue; /* too expensive */
-      P = (GEN)compositum(prev,P)[1];
+      P = compositum(prev,P);
+      P = (GEN)P[lg(P)-1];
       if (degpol(P) == N) continue;
       if (DEBUGLEVEL>2 && degpol(P) != degpol(prev))
         fprintferr("chk_gen_init: subfield %Z\n",P);

This assume the output of polcompositum is sorted by degree.
Is it true for 2.1.4 ? I remember something similar was changed
sometime ago.

Here the current list of bugs in 2.2 that may need to be investigated in 2.1

2.2.4
-----
8-  ? wrong rnfidealnorm[rel | abs] (wrong result if O_L not free O_K-module)
15- ? f(x=11,y=x)=local(t=ff(),u=t);1 --> parse error
28- ? ellsigma(ellinit([0,0,1,0,429]),matid(1))
  ***   bug in GP (Segmentation Fault), please report
32- polredabs fails to reduce
    x^8-2*x^7-34*x^6+78*x^5+265*x^4-628*x^3-389*x^2+1237*x-449
    [typo in chk_gen_init: skipfirst not initialized properly]
39- torsion group of [0,0,0,-6648,208633] reported as C2 instead of C6
? elltors(ellinit([0,0,0,-6648,208633]))
%5 = [2, [2], [[47, 0]]]
? elltors(ellinit([0,0,0,-6648,208633]),1)
%6 = [6, [6], [[26, 231]]]
56 ? polcoeff(1/x, -3)
    ***   bug in GP (Segmentation Fault), please report
57- (-2/x)/(-1/x) --> -2/-1

2.2.5
-----
18- ? ideallog(nfinit(y),2,idealstar(nfinit(y),4,1)) --> SEGV
23 ? Ser(a/b)
%1 = 1/b*x + O(x^17)
24 ? (x^2)^(1/2)
***   negative exponent in glog.
25 ? idealnorm(nfinit(y),matdiagonal([1+I])) --> SEGV
26- polylog(10,x^10) --> O(x^6)
30- ellisoncurve(ellinit([0,0,0,0,1]),[1/2-sqrt(3)/2*I,0]) --> 0
31- type 'gphelp' from command line: carriage return is missing
33- [HPUX:] should compile with cc -Ae, not -Aa
38- algdep(.1^5,1,10) --> 1
41- possible oo loop in ellap [ typo in appell1 ]

Comment appreciated.

Cheers,
Bill.