tags 635507 + patch pending tags 635917 + patch pending thanks
Dear maintainer, I've prepared an NMU for genus2reduction (versioned as 0.3-2.1) and uploaded it to DELAYED/2. Please feel free to tell me if I should delay it longer. -- .''`. : :' : Luca Falavigna <dktrkr...@debian.org> `. `' `-
diff -u genus2reduction-0.3/debian/changelog genus2reduction-0.3/debian/changelog
--- genus2reduction-0.3/debian/changelog
+++ genus2reduction-0.3/debian/changelog
@@ -1,3 +1,11 @@
+genus2reduction (0.3-2.1) unstable; urgency=medium
+
+ * Non-maintainer upload.
+ * debian/patches/pari_2.5.patch:
+ - Port to pari 2.5 API (Closes: #635507, #635917).
+
+ -- Luca Falavigna <dktrkr...@debian.org> Sun, 14 Aug 2011 15:26:01 +0200
+
genus2reduction (0.3-2) unstable; urgency=low
* Improve description and copyright file.
diff -u genus2reduction-0.3/debian/patches/series genus2reduction-0.3/debian/patches/series
--- genus2reduction-0.3/debian/patches/series
+++ genus2reduction-0.3/debian/patches/series
@@ -1,0 +2 @@
+pari_2.5.patch
only in patch2:
unchanged:
--- genus2reduction-0.3.orig/debian/patches/pari_2.5.patch
+++ genus2reduction-0.3/debian/patches/pari_2.5.patch
@@ -0,0 +1,290 @@
+Index: genus2reduction-0.3/genus2reduction.c
+===================================================================
+--- genus2reduction-0.3.orig/genus2reduction.c 2011-08-13 23:10:02.622871337 +0000
++++ genus2reduction-0.3/genus2reduction.c 2011-08-13 23:13:56.302871205 +0000
+@@ -86,7 +86,7 @@
+ long tt,d,d1,d2,d3,f1,f2,g,h,v1,v2,j2m,r1,r2,r3,r4,dismin,alpha,lambda;
+ long indice,vc6,mm,nb,dism,maxc,condp,R,alpha1,comp,beta,vy,dd,pr,ip,j,temp;
+
+- init(10000000,500000);
++ pari_init(10000000,500000);
+ jpol2=lisseq("a0;a1;a2;a3;a4;a5;a6;2^-2*(-120*a0*a6+20*a1*a5-8*a2*a4+3*a3^2)");
+ jpol4=lisseq("a0;a1;a2;a3;a4;a5;a6;2^-7*(240*(a0*a3*a4*a5+a1*a2*a3*a6)-400*(a0*a2*a5^2+a1^2*a4*a6)-64*(a0*a4^3+a2^3*a6)+16*(a1*a3*a4^2+a2^2*a3*a5)-672*a0*a3^2*a6+240*a1^2*a5^2-112*a1*a2*a4*a5-8*a1*a3^2*a5+16*a2^2*a4^2-16*a2*a3^2*a4+3*a3^4+2640*a0^2*a6^2-880*a0*a1*a5*a6+1312*a0*a2*a4*a6)");
+ jpol6=lisseq("a0;a1;a2;a3;a4;a5;a6;2^-10*(1600*(a0^2*a4^2*a5^2+a1^2*a2^2*a6^2)+1600*(a0*a1*a2*a5^3+a1^3*a4*a5*a6)+640*(a0*a1*a3*a4*a5^2+a1^2*a2*a3*a5*a6)-4000*(a0^2*a3*a5^3+a1^3*a3*a6^2)-384*(a0*a1*a4^3*a5+a1*a2^3*a5*a6)-640*(a0*a2^2*a4*a5^2+a1^2*a2*a4^2*a6)+80*(a0*a2*a3^2*a5^2+a1^2*a3^2*a4*a6)+192*(a0*a2*a3*a4^2*a5+a1*a2^2*a3*a4*a6)-48*(a0*a3^3*a4*a5+a1*a2*a3^3*a6)-224*(a1^2*a3*a4^2*a5+a1*a2^2*a3*a5^2)+64*(a1^2*a4^4+a2^4*a5^2)-64*(a1*a2*a3*a4^3+a2^3*a3*a4*a5)+16*(a1*a3^3*a4^2+a2^2*a3^3*a5)-4096*(a0^2*a4^3*a6+a0*a2^3*a6^2)+6400*(a0^2*a2*a5^2*a6+a0*a1^2*a4*a6^2)+10560*(a0^2*a3*a4*a5*a6+a0*a1*a2*a3*a6^2)+2624*(a0*a1*a3*a4^2*a6+a0*a2^2*a3*a5*a6)-4432*a0*a1*a3^2*a5*a6-8*a2*a3^4*a4+a3^6-320*a1^3*a5^3+64*a1^2*a2*a4*a5^2+176*a1^2*a3^2*a5^2+128*a1*a2^2*a4^2*a5+112*a1*a2*a3^2*a4*a5-28*a1*a3^4*a5+16*a2^2*a3^2*a4^2+5120*a0^3*a6^3-2544*a0^2*a3^2*a6^2+312*a0*a3^4*a6-14336*a0^2*a2*a4*a6^2+1024*a0*a2^2*a4^2*a6-2560*a0^2*a1*a5*a6^2-2240*a0*a1^2*a5^2*a6-6528*a0*a1*a2*a4*a5*a6-1568*a0*a2*a3^2*a4*a6)");
+@@ -115,7 +115,7 @@
+ pr=1; if(!cmpis((GEN)factp[1],2)) pr=2;
+ matmin=cgetg(6,19);
+ for(j=1;j<=5;j++)
+- matmin[j]=lgetg(lg(factp)-pr+1,18);
++ matmin[j]=(long)cgetg(lg(factp)-pr+1,18);
+ for(i=pr;i<lg(factp);i++)
+ {
+ p=(GEN)factp[i];
+@@ -400,7 +400,7 @@
+ {
+ if(myval(c3,p)>1)
+ {printf("bug c3");fflush(stdout);exit(1);}
+- mm=min(min(3*myval(c4,p)-4,3*myval(c5,p)-5),3*vc6-6);
++ mm=minss(minss(3*myval(c4,p)-4,3*myval(c5,p)-5),3*vc6-6);
+ if(alpha)
+ {
+ condp=dismin-mm-16;
+@@ -473,7 +473,7 @@
+ else
+ {
+ list=factorpadicnonun(polh,p,dismin-5);
+- nb=lg(list);prod=polun[varn(polh)];
++ nb=lg(list);prod=pol_1(varn(polh));
+ for(i=1;i<nb;i++)
+ {
+ facti=(GEN)list[i];
+@@ -524,7 +524,7 @@
+ }
+ if(cmpis(p,3)&&!cmpis(denom(theta),3)) goto modere;
+ list=factorpadicnonun(polh,p,dismin-10*alpha);
+- nb=lg(list);prod=polun[varn(polh)];
++ nb=lg(list);prod=pol_1(varn(polh));
+ for(i=1;i<nb;i++)
+ {
+ facti=(GEN)list[i];
+@@ -597,7 +597,7 @@
+ polh1=gcopy(polh);
+ theta1=gcopy(theta);
+ alpha1=alpha;
+- polh=gmul(gsubst(polh,0,ginv(polx[0])),gpuigs(polx[0],6));
++ polh=gmul(gsubst(polh,0,ginv(pol_x(0))),gpuigs(pol_x(0),6));
+ polf=polymini(gmul(polh,gpuigs(p,alpha)),p);
+ polh=(GEN)polf[1];
+ lambda=itos((GEN)polf[2]);
+@@ -655,7 +655,7 @@
+ d=val[6]-val[7]+(val[eps2]/eps);
+ if(r1&&alpha1==0)
+ {
+- polh1=gdiv(gsubst(polh1,0,gmul(p,polx[0])),gpuigs(p,3));
++ polh1=gdiv(gsubst(polh1,0,gmul(p,pol_x(0))),gpuigs(p,3));
+ }
+ polhp=gmul(unmodp,polh1);
+ if(!gcmp0(discsr(polhp)))
+@@ -671,10 +671,10 @@
+ {
+ if(gcmp1(theta1))
+ {
+- polh1=gdiv(gsubst(polh1,0,gmul(p,polx[0])),gpuigs(p,3));
++ polh1=gdiv(gsubst(polh1,0,gmul(p,pol_x(0))),gpuigs(p,3));
+ }
+ d=val[6]-3*val[3]+(val[eps2]/eps);
+- d1=min(val[7]-3*val[3],d/2);
++ d1=minss(val[7]-3*val[3],d/2);
+ if(d==2*d1) indice=d1;
+ else
+ {
+@@ -720,7 +720,7 @@
+ {
+ if(r1&&alpha1==0)
+ {
+- polh1=gdiv(gsubst(polh1,0,gmul(p,polx[0])),gpuigs(p, 3));
++ polh1=gdiv(gsubst(polh1,0,gmul(p,pol_x(0))),gpuigs(p, 3));
+ }
+ polhp=gmul(unmodp,polh1);
+ if(!gcmp0(discsr(polhp))) d=0;
+@@ -749,7 +749,7 @@
+ {
+ if(r2&&alpha==0)
+ {
+- polh=gdiv(gsubst(polh,0,gmul(p,polx[0])),gpuigs(p,3));
++ polh=gdiv(gsubst(polh,0,gmul(p,pol_x(0))),gpuigs(p,3));
+ }
+ polhp=gmul(unmodp,polh);
+ if(!gcmp0(discsr(polhp))) indice=0;
+@@ -1647,7 +1647,7 @@
+ else
+ {
+ v=myval((GEN)pol[2],p);
+- for(i=3;i<lx;i++) v=min(myval((GEN)pol[i],p),v);
++ for(i=3;i<lx;i++) v=minss(myval((GEN)pol[i],p),v);
+ }
+ return v;
+ }
+@@ -1673,8 +1673,8 @@
+ if(contpol) {printf("contenu non trivial, erreur\n");fflush(stdout);exit(1);}
+ pol=gdiv(pol,gpuigs(p,contpol));
+ r=(lx-4)*ggval(td,p)+r;
+- upol=gsubst(gmul(pol,gpuigs(td,lx-4)),vx,gdiv(polx[vx],td));
+- fact=gsubst(factorpadic4(upol,p,r),vx,gmul(polx[vx],td));
++ upol=gsubst(gmul(pol,gpuigs(td,lx-4)),vx,gdiv(pol_x(vx),td));
++ fact=gsubst(factorpadic4(upol,p,r),vx,gmul(pol_x(vx),td));
+ list=(GEN)fact[1];ep=(GEN)fact[2];nb=lg(list)-1;
+ for(i=1;i<=nb;i++)
+ if(!gcmp1((GEN)ep[i]))
+@@ -1698,7 +1698,7 @@
+ rec1:
+ alpha=polval(pol,p);polh=gdiv(pol,gpuigs(p,alpha));
+ if(gdivise(truecoeff(polh,6),p)&&gdivise(truecoeff(polh,5),p)&&gdivise(truecoeff(polh,4),p)&&gdivise(truecoeff(polh,3),p))
+- {polh=gmul(gpuigs(polx[0],6),gsubst(polh,0,ginv(polx[0])));}
++ {polh=gmul(gpuigs(pol_x(0),6),gsubst(polh,0,ginv(pol_x(0))));}
+ alpha=alpha&1;
+ beta=0;
+ lambda=0;
+@@ -1709,11 +1709,11 @@
+
+ labelm6:
+ theta=caltheta(polh,p,lambda);
+- polf[5]=zero; /* 1 si racine d'ordre 3 non dans F_p, 0 sinon */
++ polf[5]=(long)gen_0; /* 1 si racine d'ordre 3 non dans F_p, 0 sinon */
+ if(gcmp(theta,gun)>=0)
+ {
+ ent=gfloor(theta);pent=gpui(p,ent,0);
+- polh=gdiv(gsubst(polh,0,gmul(pent,polx[0])),gpuigs(pent,6-lambda));
++ polh=gdiv(gsubst(polh,0,gmul(pent,pol_x(0))),gpuigs(pent,6-lambda));
+ alpha=(alpha+lambda*itos(ent))&1;
+ beta=beta+itos(ent);
+ theta=gsub(theta,ent);
+@@ -1728,11 +1728,11 @@
+ rac=(GEN)fac[2];
+ if(lgef(rac)==5)
+ {
+- polf[5]=un;
++ polf[5]=(long)gen_1;
+ goto sortie;
+ }
+ rac=lift(gdiv(gneg(truecoeff(rac,0)),truecoeff(rac,1)));
+- polh=gsubst(polh,0,gadd(polx[0],rac));
++ polh=gsubst(polh,0,gadd(pol_x(0),rac));
+ lambda=6-maxord; goto labelm6;
+ }
+ if(maxord<=2&&lgef(polhp)>=7) goto sortie;
+@@ -1741,7 +1741,7 @@
+ {
+ if((myval(truecoeff(polh,2),p)>1-alpha)&&(myval(truecoeff(polh,1),p)>2-alpha)&&(myval(truecoeff(polh,0),p)>3-alpha))
+ {
+- pol=gsubst(polh,0,gmul(p,polx[0]));
++ pol=gsubst(polh,0,gmul(p,pol_x(0)));
+ if(alpha) pol=gmul(p,pol);
+ goto rec1;
+ }
+@@ -1753,7 +1753,7 @@
+ {
+ if(myval(truecoeff(polh,6),p)>=3&&myval(truecoeff(polh,5),p)>=2)
+ {
+- polh=gmul(gsubst(polh,0,gdiv(polx[0],p)),gpuigs(p,3));
++ polh=gmul(gsubst(polh,0,gdiv(pol_x(0),p)),gpuigs(p,3));
+ theta=gadd(theta,gun);
+ alpha=0;
+ beta=beta-1;
+@@ -1762,13 +1762,13 @@
+ }
+ if(lgef(polhp)==9&&!gcmp0(theta))
+ {
+- pro=gdiv(polhp,gpuigs(polx[0],3));
++ pro=gdiv(polhp,gpuigs(pol_x(0),3));
+ fac=factmz(pro,p);
+ if(!cmpis((GEN)fac[1],3))
+ {
+ rac=(GEN)fac[2];
+ rac=lift(gdiv(gneg(truecoeff(rac,0)),truecoeff(rac,1)));
+- pro=gsubst(polh,0,gadd(rac,gmul(p,polx[0])));
++ pro=gsubst(polh,0,gadd(rac,gmul(p,pol_x(0))));
+ if(polval(pro,p)>=3)
+ {
+ polh=gdiv(pro,gpuigs(p,3));
+@@ -1815,7 +1815,7 @@
+ if(lgef(polhp)!=6)
+ {printf("bug, on ne doit pas arriver dans discpart");fflush(stdout);exit(1);}
+ list=factorpadicnonun(polh,p,ord);
+- prod=polun[varn(polh)];
++ prod=pol_1(varn(polh));
+ for(i=1;i<lg(list);i++)
+ {
+ facti=(GEN)list[i];
+@@ -1837,7 +1837,7 @@
+ long
+ myvalzi(GEN b) /* valuation par rapport a 3 dans Z[i] */
+ {
+- return min(myval(greal(b),stoi(3)),myval(gimag(b),stoi(3)));
++ return minss(myval(greal(b),stoi(3)),myval(gimag(b),stoi(3)));
+ }
+
+ GEN
+@@ -1853,15 +1853,15 @@
+ polh=pol;
+ if(alpha) polh=gdiv(pol,p);
+ beta=0;
+- rac=gadd(lift(gdiv(gmul(unmodp,truecoeff(polh,3)),gmul(unmodp,truecoeff(polh,6)))),gi);
++ rac=gadd(lift(gdiv(gmul(unmodp,truecoeff(polh,3)),gmul(unmodp,truecoeff(polh,6)))),gen_I());
+
+ labelm6:
+- polh=gsubst(polh,0,gadd(polx[0],rac));
++ polh=gsubst(polh,0,gadd(pol_x(0),rac));
+ theta=calthetazi(polh);
+ if(gcmp(theta,gun)>=0)
+ {
+ ent=gfloor(theta);pent=gpui(p,ent,0);
+- polh=gdiv(gsubst(polh,0,gmul(pent,polx[0])),gpuigs(pent,3));
++ polh=gdiv(gsubst(polh,0,gmul(pent,pol_x(0))),gpuigs(pent,3));
+ alpha=(alpha+itos(ent))&1;
+ beta=beta+itos(ent);
+ theta=gsub(theta,ent);
+@@ -1879,7 +1879,7 @@
+ {
+ if(myvalzi(truecoeff(polh,6))>=3&&myvalzi(truecoeff(polh,5))>=2&&myvalzi(truecoeff(polh,4))>=1)
+ {
+- polh=gmul(gsubst(polh,0,gdiv(polx[0],p)),gpuigs(p,3));
++ polh=gmul(gsubst(polh,0,gdiv(pol_x(0),p)),gpuigs(p,3));
+ theta=gadd(theta,gun);
+ beta=beta-1;
+ alpha=0;
+@@ -1907,7 +1907,7 @@
+ /* on calcule la valuation de b dans Z[i, Y]/(Y^2-3) */
+ {
+ b=lift(b);
+- return min(2*myvalzi(truecoeff(b,0)),2*myvalzi(truecoeff(b,1))+1);
++ return minss(2*myvalzi(truecoeff(b,0)),2*myvalzi(truecoeff(b,1))+1);
+ }
+
+ GEN
+@@ -1936,10 +1936,10 @@
+ {printf("bug, le polynome n'est pas minimal\n");fflush(stdout);exit(1);}
+ polh=pol;alpha=0;
+ beta=0;
+- polh=gdivgs(gsubst(polh,0,gmul(polx[0],y)),27);
++ polh=gdivgs(gsubst(polh,0,gmul(pol_x(0),y)),27);
+ if(myvalzi2(truecoeff(polh,4))>0&&myvalzi2(truecoeff(polh,2))>0)
+ {
+- if(myvalzi2(gsub(truecoeff(polh,6),truecoeff(polh,0)))>0) rac=gi;
++ if(myvalzi2(gsub(truecoeff(polh,6),truecoeff(polh,0)))>0) rac=gen_I();
+ else rac=gun;
+ }
+ else
+@@ -1949,12 +1949,12 @@
+ }
+
+ labelm6:
+- polh=gsubst(polh,0,gadd(polx[0],rac));
++ polh=gsubst(polh,0,gadd(pol_x(0),rac));
+ theta=calthetazi2(polh);
+ if(gcmp(theta,gun)>=0)
+ {
+ ent=gfloor(theta);pent=gpui(y,ent,0);
+- polh=gdiv(gsubst(polh,0,gmul(pent,polx[0])),gpuigs(pent,3));
++ polh=gdiv(gsubst(polh,0,gmul(pent,pol_x(0))),gpuigs(pent,3));
+ alpha=(alpha+itos(ent))&1;
+ beta=beta+itos(ent);
+ theta=gsub(theta,ent);
+@@ -1998,7 +1998,7 @@
+ u=truecoeff(lift(u),0);
+ a=gmul(unmodp,greal(u));
+ b=gmul(unmodp,gimag(u));
+- u=gadd(a,gmul(b,gi));
++ u=gadd(a,gmul(b,gen_I()));
+ return u;
+ }
+
+@@ -2020,7 +2020,7 @@
+ {
+ fac=(GEN)factmod(polhp,p);
+ for(i=1;i<lg((GEN)fac[2]);i++)
+- m=max(m,itos(gcoeff(fac,i,2)));
++ m=maxss(m,itos(gcoeff(fac,i,2)));
+ for(i=1;i<lg((GEN)fac[1]);i++)
+ if(!cmpis(gcoeff(fac,i,2),m)) l=i;
+ z=gcoeff(fac,l,1);
signature.asc
Description: OpenPGP digital signature
