passagemath-pari 10.6.32__cp314-cp314-musllinux_1_2_x86_64.whl

sage/ext_data/pari/dokchitser/ex-bsw

@@ -0,0 +1,43 @@
+/*** EXAMPLE: Birch-Swinnerton-Dyer for an elliptic curve over Q  ***/
+/*** v1.2, July 2013, questions to tim.dokchitser@bristol.ac.uk   ***/
+/***                                                              ***/
+/*** type \rex-bsw or read("ex-bsw") at pari prompt to run this   ***/
+
+elleq  = [0,0,1,-7,6];            \\ L-function of this elliptic curve over Q
+points = [[1,-1],[2,-1],[0,-3]];  \\ basis of Mordell-Weil modulo torsion
+                                  \\ May change this to any other E/Q
+
+read("computel");                 \\ read the ComputeL package
+                                  \\ and set the default values
+default(realprecision,20);        \\ set working precision; used throughout
+
+ellorig = ellinit(elleq);
+gred    = ellglobalred(ellorig);
+ell     = ellchangecurve(ellorig,gred[2]);    \\ global minimal model
+
+                            \\ initialize L-function parameters
+conductor = gred[1];        \\ conductor for the exponential factor
+gammaV    = [0,1];          \\ list of gamma-factors
+weight    = 2;              \\ L(s)=sgn*L(weight-s)
+sgn       = ellrootno(ell); \\ sign in the functional equation
+a(k)      = ellak(ell,k);   \\ L-series coefficients a(k)
+
+initLdata("a(k)",1);  \\ initialize the vector of coefficients
+                      \\ we won't call checkfeq(cutoff) or L(s,cutoff)
+                      \\ so set cutoff=1, saves some coefficients
+
+print("EXAMPLE: Birch-Swinnerton-Dyer on an elliptic curve over Q");
+print("         with ",default(realprecision)," digits precision");
+print("Elliptic curve : ", elleq);
+print("Points         : ", points);
+print("Conductor      = ", conductor);
+print("Root number    = ", sgn);
+print("Mord.-Weil rnk = ", rnk  = length(points));
+print("Lead.term L(1) = ", lval = L(1,,rnk)/rnk!); \\ leading (order rnk)
+                                                   \\ term of L(s) at s=1
+print("Real period    = ", om   = ell.omega[1]*((ell.disc>0)+1));
+print("Regulator      = ", reg  = matdet(ellheightmatrix(ellorig,points)));
+print("Product(c_p)   = ", cp   = ellglobalred(ell)[3]);
+print("Torsion size   = ", tors = elltors(ell)[1]);
+print("Analytic #Sha  = ", sha  = lval/om/reg/cp*tors^2);
+print("               = ", round(sha),", error=",errprint(sha-round(sha)));

sage/ext_data/pari/dokchitser/ex-chgen

@@ -0,0 +1,48 @@
+/*** EXAMPLE: Dirichlet L-function of a general character modulo odd P   ***/
+/*** (illustration of functional equations with 2 different L-functions) ***/
+/***                                                                     ***/
+/*** v1.2, July 2013, questions to tim.dokchitser@bristol.ac.uk          ***/
+/*** type \rex-chgen or read("ex-chgen") at pari prompt to run this      ***/
+
+read("computel");                 \\ read the ComputeL package
+                                  \\ and set the default values
+default(realprecision,28);        \\ set working precision; used throughout
+
+P = 37;                           \\ May change this to any odd prime
+
+\\ Take the character defined by the following modulo P (has order P-1)
+\\    (primitive root) -> e^(2 Pi I/(P-1))
+\\ a(n) = if(n%P,exp(2*Pi*I/(P-1)*znlog(n,prim)),0);
+\\ znlog needs pari 2.1, the following works with earlier versions
+
+prim = znprimroot(P);
+avec = vector(P,k,0);
+for (k=0,P-1,avec[lift(prim^k)+1]=exp(2*Pi*I*k/(P-1)));
+a(n)  = avec[n%P+1];
+
+                          \\ initialize L-function parameters
+conductor = P;            \\ exponential factor
+gammaV    = [1];          \\ [0] for even, [1] for odd
+weight    = 1;            \\ L(s)=sgn*L(weight-s)
+sgn       = X;            \\ unknown, to be solved from functional equation
+
+initLdata("a(k)",,"conj(a(k))");  \\ initialize L-series coefficients a(k)
+                                  \\ and the ones of the dual motive
+
+sgneq = Vec(checkfeq());       \\ checkfeq() returns c1+X*c2, should be 0
+sgn   = -sgneq[2]/sgneq[1];    \\ hence solve for the sign
+
+print("EXAMPLE: L(s)=L(Chi,s), Chi=character of order P-1 modulo P");
+print("         with ",default(realprecision)," digits precision");
+print("         take P=",P);
+print("Chi(",lift(prim),")      = exp(2*Pi*I/",P-1,\
+      ")    defines the character (mod ",P,")");
+print("Parity      = ", if(gammaV[1],"odd","even"));
+print("Sign          we will try to solve from the functional equation");
+print("            = ", sgn);
+print("|sign|      = ", abs(sgn));
+print("Verifying functional equation. Error: ",errprint(checkfeq(1.1)));
+print("L(Chi,2)    = ",L(2));
+print("  (check)   = ",L(2,1.1));
+print("L(Chi,I)    = ",L(I));
+print("  (check)   = ",L(I,1.1));

sage/ext_data/pari/dokchitser/ex-chqua

@@ -0,0 +1,37 @@
+/*** EXAMPLE: Dirichlet L-function of a quadratic character modulo odd P ***/
+/*** v1.2, July 2013, questions to tim.dokchitser@bristol.ac.uk          ***/
+/***                                                                     ***/
+/*** type \rex-chqua or read("ex-chqua") at pari prompt to run this      ***/
+
+read("computel");                 \\ read the ComputeL package
+                                  \\ and set the default values
+default(realprecision,28);        \\ set working precision; used throughout
+
+P = 37;                     \\ May change this to any other odd prime
+                            \\ Chi = Legendre character modulo this P
+
+                            \\ initialize L-function parameters
+conductor = P;              \\ exponential factor
+gammaV    = [P%4==3];       \\ [0] if P=1 mod 4 (even), [1] if 3 mod 4 (odd)
+weight    = 1;              \\ L(s)=sgn*L(weight-s)
+sgn       = 1;              \\ sign in the functional equation
+
+a(n) = if(n%P,(-1)^((P-1)/znorder(Mod(n,P))),0);
+  \\ coefficients: 1 for n square mod P, -1 for non-square, 0 for n=0 mod P
+
+initLdata("a(k)");              \\ initialize L-function parameters
+
+zinit     = zetakinit(nfinit(quadpoly(quaddisc(if(P%4==3,-P,P)))));
+LPari(x)  = zetak(zinit,x)/zeta(x);
+  \\ alternatively the same L-function is the Dedekind zeta function of
+  \\ the maximal order of Q(sqrt(P)) (p=1 mod 4) or Q(sqrt(-P)) (p=3 mod 4)
+  \\ divided by the Riemann zeta function
+
+print("EXAMPLE: L(s)=L(Chi(P),s), Chi(P)=(Legendre) quadratic character modulo P");
+print("         with ",default(realprecision)," digits precision");
+print("         take P=",P);
+print("Verifying functional equation. Error: ",errprint(checkfeq()));
+print("L(Chi(P),2) = ",L(2));
+print("  (pari)    = ",LPari(2));
+print("L(Chi(P),I) = ",L(I));
+print("  (pari)    = ",LPari(I));

sage/ext_data/pari/dokchitser/ex-delta

@@ -0,0 +1,35 @@
+/*** EXAMPLE: Modular form Delta                                    ***/
+/*** (illustration of unusual coefficient growth)                   ***/
+/***                                                                ***/
+/*** v1.2, July 2013, questions to tim.dokchitser@bristol.ac.uk     ***/
+/*** type \rex-delta or read("ex-delta") at pari prompt to run this ***/
+
+read("computel");                 \\ read the ComputeL package
+                                  \\ and set the default values
+default(realprecision,40);        \\ set working precision; used throughout
+
+\\ * Coefficients are given by Ramanujan's tau function
+\\ * Re-define the default bound on the coefficients:
+\\   Deligne's estimate on tau(n) is better than the default
+\\   coefgrow(n)=(4n)^(11/2)     (by a factor 1024),
+\\   so re-defining coefgrow() improves efficiency (slightly faster)
+
+coefgrow(n) = 2*n^(11/2);
+tau(n)      = (5*sigma(n,3)+7*sigma(n,5))*n/12\
+              -35*sum(k=1,n-1,(6*k-4*(n-k))*sigma(k,3)*sigma(n-k,5));
+
+                            \\ initialize L-function parameters
+conductor = 1;              \\ exponential factor
+gammaV    = [0,1];          \\ list of gamma-factors
+weight    = 12;             \\ L(s)=sgn*L(weight-s)
+sgn       = 1;              \\ sign in the functional equation
+
+
+initLdata("tau(k)");        \\ L-series coefficients a(k)
+
+print("EXAMPLE: L-function associated to the modular form Delta of weight 12");
+print("         coefficients = Ramanujan's tau function");
+print("         with ",default(realprecision)," digits precision");
+print("Verifying functional equation. Error: ",errprint(checkfeq()));
+print("L(1)       = ",lval = L(1));
+print(" (check)   = ",lval2 = L(1,1.1),"  (err=",errprint(lval-lval2),")");

sage/ext_data/pari/dokchitser/ex-eisen

@@ -0,0 +1,30 @@
+/*** EXAMPLE: Eisenstein series G_k, k>2 even                       ***/
+/*** v1.2, July 2013, questions to tim.dokchitser@bristol.ac.uk     ***/
+/***                                                                ***/
+/*** type \rex-eisen or read("ex-eisen") at gp prompt to run this   ***/
+
+read("computel");                 \\ read the ComputeL package
+                                  \\ and set the default values
+default(realprecision,40);        \\ set working precision; used throughout
+
+K = 16;    \\ Our modular form is G_K with this K
+           \\ may change this to any even K
+
+conductor = 1;              \\ exponential factor
+gammaV    = [0,1];          \\ list of gamma-factors
+weight    = K;              \\ L(s)=sgn*L(weight-s)
+sgn       = (-1)^(K/2);     \\ sign in the functional equation
+
+\\ It has a simple pole in s=K
+Lpoles    = [K];
+Lresidues = [(-1)^(K/2)*sqrt(Pi)*bernfrac(K)/K];
+
+initLdata("sigma(k,K-1)");  \\ initialize L-series
+                            \\ Coefficients given by the divisor function
+
+print("EXAMPLE: L-function associated to the modular form G_",K," of weight ",K);
+print("         coefficients = divisor function sigma(n,",K-1,")");
+print("         with ",default(realprecision)," digits precision");
+print("Verifying functional equation. Error: ",errprint(checkfeq()));
+print("L(1)       = ",lval = L(1));
+print(" (check)   = ",lval2 = L(1,1.1),"  (err=",errprint(lval-lval2),")");

sage/ext_data/pari/dokchitser/ex-gen2

@@ -0,0 +1,38 @@
+/*** EXAMPLE: L-function of a genus 2 curve over Q                  ***/
+/*** v1.3, December 2013, questions to tim.dokchitser@bristol.ac.uk ***/
+/***                                                                ***/
+/*** type \rex-gen2 or read("ex-gen2") at pari prompt to run this   ***/
+
+
+read("computel");                 \\ read the ComputeL package
+                                  \\ and set the default values
+default(realprecision,14);        \\ set working precision; used throughout
+                                  \\ larger precision needs more coefficients
+
+/*********** C: y^2+(x^3+x+1)y=x^5+x^4  (genus2, hyperelliptic) **********/
+cond = 169;  \\ conductor
+sgn  = 1;    \\ sign in the functional equation
+\\ Vector of coefficients
+coef = [1,-3,-2,4,0,6,0,-3,3,0,0,-8,-5,0,0,3,3,-9,-6,0,0,0,6,6,7,15,-10,0,\
+  -3,0,0,-6,0,-9,0,12,15,18,10,0,-9,0,-8,0,0,-18,0,-6,-7,-21,-6,-20,-6,30,\
+  0,0,12,9,12,0,-1,0,0,5,0,0,6,12,-12,0,6,-9,0,-45,-14,-24,0,-30,8,0,20,27,\
+  0,0,0,24,6,0,-12,0,0,24,0,0,0,12,12,21,0,28,3,18,-20,15,0,18,-6,-40,0,0,\
+  -30,0,15,-36,0,-12,-15,-36,0,0,-11,3,18,0,0,0,2,9,16,0,36,0,0,-18,0,-9,\
+  -27,36,4,0,0,-18,0,9,0,0];
+/*************************************************************************/
+
+                              \\ initialize L-function parameters
+conductor = cond;             \\ exponential factor
+gammaV    = [0,0,1,1];        \\ list of gamma-factors
+weight    = 2;                \\ L(s)=sgn*L(weight-s)
+initLdata("coef[k]",-146);    \\ initialize the L-series data
+
+print("EXAMPLE: L-function of a genus 2 curve over Q");
+print("         with 14 digits precision");
+print("C: y^2+(x^3+x+1)y=x^5+x^4 (hyperelliptic)");
+print("Conductor      = ", conductor);
+print("Functional eq. = ", errprint(checkfeq()));
+print("L(1)           = ", L(1));
+print(" (check)       = ", L(1,1.1));
+print("L(2)           = ", L(2));
+print(" (check)       = ", L(2,1.1));

sage/ext_data/pari/dokchitser/ex-gen3

@@ -0,0 +1,49 @@
+/*** EXAMPLE: L-function of a genus 3 curve over Q                   ***/
+/*** (illustration of precision issues when not enough coefficients) ***/
+/***                                                                 ***/
+/*** v1.2, July 2013, questions to tim.dokchitser@bristol.ac.uk      ***/
+/*** type \rex-gen3 or read("ex-gen3") at pari prompt to run this    ***/
+
+
+read("computel");                 \\ read the ComputeL package
+                                  \\ and set the default values
+default(realprecision,20);        \\ set working precision; used throughout
+                                  \\ actually need more coefficients for this
+
+/*** C: y^2+(x^3+x^2+x+1)y=x^7+2*x^6+2*x^5+x^4 (genus 3, hyperelliptic) ***/
+conductor = 13555; \\ conductor
+sgn       = 1;     \\ sign in the functional equation
+\\ Vector of 300 coefficients
+coef = [1,-1,-2,-1,-1,2,-2,1,4,1,-3,2,0,2,2,-1,8,-4,-3,1,4,3,4,-2,-3,0,-4,2,\
+  2,-2,-4,1,6,-8,2,-4,3,3,0,-1,-8,-4,-4,3,-4,-4,-5,2,1,3,-16,0,13,4,3,-2,6,\
+  -2,7,-2,0,4,-8,-1,0,-6,6,-8,-8,-2,-7,4,-9,-3,6,3,6,0,-10,1,0,8,-2,-4,-8,4,\
+  -4,-3,5,4,0,-4,8,5,3,-2,2,-1,-12,3,4,16,11,0,-4,-13,-9,4,2,-3,-6,2,7,-6,\
+  -4,-2,0,-7,-16,2,-14,0,16,4,3,8,13,9,8,0,8,-6,6,-6,4,8,-8,8,15,-2,10,7,0,\
+  -4,-2,9,-2,-3,4,-6,2,-3,32,-6,4,0,-20,10,-26,-1,-8,0,10,8,-6,2,14,4,-4,8,\
+  -12,4,-2,4,6,3,-14,-5,5,4,-9,0,0,4,-3,-8,-24,5,8,-3,6,2,7,-2,0,-1,-13,12,\
+  -9,-3,-12,-4,-4,16,8,-11,16,0,9,4,-11,-13,14,9,4,-4,8,-2,18,-3,0,6,-5,-2,\
+  -12,-7,13,-6,0,4,-12,2,17,0,5,-7,20,16,26,-2,-6,14,-2,0,-1,-16,0,-4,4,-3,\
+  -1,8,-12,-13,16,-1,16,-8,-6,0,8,-8,13,6,-13,-6,-10,-6,7,-4,-2,-8,0,8,9,8,\
+  8,-15,-16,2,-10,-10,16,7,-6,0,16,4,24,2,-4,9,-17,2,-7,3,12,-4,0,-6\
+];
+/*************************************************************************/
+
+                              \\ initialize other L-function parameters
+gammaV    = [0,0,0,1,1,1];    \\ list of gamma-factors
+weight    = 2;                \\ L(s)=sgn*L(weight-s)
+initLdata("coef[k]",-300);    \\ initialize the L-series data
+
+print("EXAMPLE: L-function of a genus 3 curve over Q");
+print("         with ",default(realprecision)," digits precision");
+print("C: y^2+(x^3+x^2+x+1)y=x^7+2*x^6+2*x^5+x^4 (hyperelliptic)");
+
+print("Actually need ",cflength()," coefficients to properly verify functional equation");
+print("          and ",cflength(1)," coefficients to compute L(s)");
+print("        given ",length(coef)," coefficients only, but let's try anyway:");
+print("Functional eq. = ", errprint(checkfeq(1.1))," <- indication of actual precision");
+print("Conductor      = ", conductor);
+print("Sign           = ", sgn);
+print("L(1)           = ", L(1));
+print(" (check)       = ", L(1,1.1));
+print("L(2)           = ", L(2));
+print(" (check)       = ", L(2,1.1));

sage/ext_data/pari/dokchitser/ex-gen4

@@ -0,0 +1,54 @@
+/*** EXAMPLE: L-function of a genus 4 curve over Q                   ***/
+/*** (illustration of precision issues when not enough coefficients) ***/
+/***                                                                 ***/
+/*** v1.2, July 2013, questions to tim.dokchitser@bristol.ac.uk      ***/
+/*** type \rex-gen4 or read("ex-gen4") at pari prompt to run this    ***/
+
+
+read("computel");                 \\ read the ComputeL package
+                                  \\ and set the default values
+default(realprecision,20);        \\ set working precision; used throughout
+                                  \\ actually need more coefficients for this
+
+/*** C: y^2+(x^5+x+1)y=x^7-x^6+x^4 (genus 4, hyperelliptic) ***/
+conductor = 3*23*7333; \\ conductor
+sgn       = 1;         \\ sign in the functional equation
+\\ Vector of 450 coefficients
+coef=[1,-1,-3,-1,0,3,-1,0,2,0,2,3,1,1,0,0,-1,-2,-6,0,3,-2,2,0,0,-1,3,1,5,0,6,\
+      2,-6,1,0,-2,-2,6,-3,0,-1,-3,-3,-2,0,-2,3,0,6,0,3,-1,2,-3,0,0,18,-5,9,0,\
+     -3,-6,-2,1,0,6,-4,1,-6,0,-4,0,-5,2,0,6,-2,3,-11,0,-9,1,17,-3,0,3,-15,0,\
+     0,0,-1,-2,-18,-3,0,-6,-5,-6,4,0,-4,-3,-4,0,0,-2,-3,-3,-20,0,6,0,8,-18,0,\
+     -5,2,-9,1,0,4,3,3,-6,0,2,-12,1,9,0,-3,6,6,4,0,0,13,6,-22,0,-9,4,2,0,0,5,\
+     -18,2,5,0,2,0,-2,2,0,3,7,11,-6,0,-2,9,4,1,0,-17,-15,0,7,0,-12,3,1,15,0,0,\
+     -27,0,-10,0,-14,1,9,0,0,18,-2,-3,-3,0,5,-3,2,5,0,-6,3,-4,19,0,12,4,-5,-3,\
+     0,4,4,0,-12,0,14,-2,12,3,0,0,-6,20,15,0,-1,-6,21,-2,0,-8,-11,-18,19,0,6,\
+     0,16,-2,0,-9,33,-1,-6,0,3,-4,15,3,0,-3,-6,0,-51,0,20,2,4,12,0,-9,24,-9,2,\
+     0,10,3,-25,0,0,-6,0,4,-29,0,-11,0,3,-13,0,6,15,22,12,0,-12,9,22,4,0,-2,1,\
+     4,9,0,15,5,1,18,0,0,6,-5,2,0,3,-2,12,0,0,2,-11,2,12,0,-12,0,3,-7,0,11,19,\
+     6,10,0,9,2,6,9,0,-4,60,0,-3,0,-4,-17,-4,15,0,0,27,-7,-24,0,12,12,-18,0,0,\
+     -1,-24,15,16,0,3,4,27,27,0,0,-3,10,5,0,-2,14,-12,1,0,-9,-20,0,-2,0,-2,18,\
+     43,2,0,0,5,3,-23,0,36,-5,3,-3,0,-2,-6,5,7,0,-2,0,9,-3,0,-4,-3,-19,-18,0,\
+     0,-12,6,4,0,5,-4,0,-2,0,-39,4,-9,-4,0,2,66,12,-17,0,10,-14,6,0,0,-12,3,3,\
+     -6,0,12,0,-14,6,0,20,-12,-15,6,0,12,1,-29,-6,0,-21,-15,-1,46,0\
+  ];
+/*************************************************************************/
+
+                               \\ initialize other L-function parameters
+gammaV    = [0,0,0,0,1,1,1,1]; \\ list of gamma-factors
+weight    = 2;                 \\ L(s)=sgn*L(weight-s)
+initLdata("coef[k]",-450);     \\ initialize the L-series data
+
+print("EXAMPLE: L-function of a genus 4 curve over Q");
+print("         with ",default(realprecision)," digits precision");
+print("C: y^2+(x^5+x+1)y=x^7-x^6+x^4 (genus 4, hyperelliptic)");
+
+print("Actually need ",cflength()," coefficients to properly verify functional equation");
+print("          and ",cflength(1)," coefficients to compute L(s)");
+print("        given ",length(coef)," coefficients only, but let's try anyway:");
+print("Functional eq. = ", errprint(checkfeq(1.1))," <- indication of actual precision");
+print("Conductor      = ", conductor);
+print("Sign           = ", sgn);
+print("L(1)           = ", L(1));
+print(" (check)       = ", L(1,1.1));
+print("L(2)           = ", L(2));
+print(" (check)       = ", L(2,1.1));

sage/ext_data/pari/dokchitser/ex-nf

@@ -0,0 +1,48 @@
+/*** Dedekind zeta-function of a number field F/Q                   ***/
+/*** v1.3, December 2013, questions to tim.dokchitser@bristol.ac.uk ***/
+/***                                                                ***/
+/*** type \rex-nf or read("ex-nf") at pari prompt to run this       ***/
+
+fpol = x^4-x^2-1;                 \\ polynomial which defines F/Q
+                                  \\ may change this to any other polynomial
+
+read("computel");                 \\ read the ComputeL package
+                                  \\ and set the default values
+default(realprecision,20);        \\ set working precision; used throughout
+
+nf        = nfinit(fpol);         \\ initialize the number field F/Q
+zinit     = zetakinit(nf);        \\ and the built-in Dedekind zeta-function
+LPari(x)  = zetak(zinit,x);
+
+print("EXAMPLE: Dedekind zeta-function L(s)=zeta_F(s) of number field F");
+print("         with ",default(realprecision)," digits precision");
+print("F                   = Q(root of ",fpol,")");
+print("[F:Q]               = ",Fdegree = poldegree(fpol));
+print("Discriminant        = ",disc    = nf.disc);
+print("r1 (real emb.)      = ",r1      = nf.sign[1]);
+print("r2 (complex emb.)   = ",r2      = nf.sign[2]);
+print("Gamma factor        = ",gammaV  = concat(vector(r1+r2,X,0),vector(r2,X,1)));
+
+conductor = abs(disc);                   \\ exponential factor
+weight    = 1;                           \\ L(s)=sgn*L(weight-s)
+sgn       = 1;                           \\ sign in the functional equation
+Lpoles    = [1];                         \\ pole at s=1, Lresidues=automatic
+dzk       = dirzetak(nf,cflength());     \\ coefficients a(k) in L(s)
+
+initLdata("dzk[k]");         \\ initialize L-series
+
+\\ Determine residue at s=1 and check the functional equation
+print("Error in func. eq.  = ",errprint(checkfeq()));
+
+\\ Compare L(2) and built-in LPari(2)
+print("L(2)                = ",L(2));
+print(" (or using pari)    = ",LPari(2));
+
+\\ Compare residue at 1 and the one given by class number formula
+print("Residue at s=1");
+print(" (automatically determined) = ",Lresidues[1]);
+
+\\ Determine the residue at s=1 using the class number formula
+bnf       = bnfinit(fpol);
+residue   = - 2^(r1+r2) * Pi^(r2/2) * bnf.reg * bnf.clgp.no / bnf.tu[1];
+print(" (class number formula)     = ",residue);

sage/ext_data/pari/dokchitser/ex-shin

@@ -0,0 +1,50 @@
+/*** EXAMPLE: Shintani's zeta-function                              ***/
+/*** v1.3, December 2013, questions to tim.dokchitser@bristol.ac.uk ***/
+/***                                                                ***/
+/*** type \rex-shin or read("ex-shin") at pari prompt to run this   ***/
+
+read("computel");                 \\ read the ComputeL package
+                                  \\ and set the default values
+default(realprecision,19);        \\ set working precision; used throughout
+                                  \\ larger precision needs more coefficients
+
+                            \\ initialize L-function parameters
+GlobalC   = 1/2/sqrt(Pi);   \\ Global coefficient in front of L(s)
+conductor = 432;            \\ exponential factor
+gammaV    = [0,1,-1/6,1/6]; \\ list of gamma-factors
+weight    = 1;              \\ L(s)=sgn*L(weight-s)
+sgn       = 1;              \\ sign in the functional equation
+Lpoles    = [1/6,1];        \\ ``half'' of the poles of L(s)
+{Lresidues =                \\ and residues in there
+  [zeta(2/3)*2^(3/2)*3^(-1/4)*Pi^(-1/3)*gamma(1/3)/GlobalC,   \\ in s=1/6
+   gamma(1/12)*gamma(-1/12)/sqrt(3)/12/GlobalC];}             \\ in s=1
+  \\ note: can also set Lresidues=automatic, checkfeq()
+  \\ will then determine them to decent precision
+
+\\ First 400 Coefficients
+
+{c1=[1/3,1,1,1,1,1,1,4/3,1,1,1,1,1,1,1,2,1,3,1,1,1,1,1,2,5/3,1,1,1,1,1,1,2,1,3,
+    1,1,1,1,1,2,2,1,1,1,1,1,1,2,1,3,1,1,1,2,1,2,1,1,3,1,1,1,2,2,1,3,1,1,1,1,1,
+    4,1,3,1,1,1,1,1,2,1,3,1,1,5/3,1,1,2,1,1,1,1,1,1,2,2,1,3,1,1,1,1,1,2,1,1,1,
+    2,1,1,1,2,3,3,3,1,1,1,1,2,1,1,1,1,1,3,1,10/3,3,3,1,1,1,1,1,4,1,1,1,1,1,1,
+    1,2,1,3,1,1,2,1,1,2,1,1,1,1,1,3,1,4,3,4,1,1,1,1,3,2,1,1,1,1,1,1,1,2,1,3,1,
+    3,5/3,1,1,2,1,1,1,1,1,1,1,2,1,3,1,3,1,1,1,4];
+c2=[1,1,1,1,1,1,1,2,1,1,3,1,2,3,3,2,1,1,1,1,1,3,1,2,1,1,1,1,3,3,1,2,1,1,1,3,1,
+    3,1,2,3,1,3,1,1,3,1,2,3,3,1,3,3,4,1,4,1,3,1,1,1,3,1,2,1,1,4,1,3,3,1,2,3,1,
+    1,3,1,3,1,2,1,1,1,1,1,3,3,4,1,3,1,1,1,3,1,4,1,1,3,3,1,3,3,2,3,3,1,4,1,3,1,
+    2,1,1,3,1,1,3,3,4,4,3,3,1,1,3,3,4,1,1,1,1,1,3,1,2,3,1,3,1,3,3,1,6,1,1,1,1,
+    1,3,1,4,3,1,1,1,3,3,1,2,1,4,3,1,3,3,3,2,3,1,2,1,1,3,4,2,1,3,1,1,1,3,3,6,1,
+    1,2,1,3,3,1,2,1,1,1,1,1,3,1,6];}
+a(n) = if(n%4>2,c2[n\2]/sqrt(3),if(n%4>1,0,if(n%4>0,c1[n\2+1],c1[n/2]+c2[n/2]/sqrt(3))));
+coefgrow(n) = n^(1/3);  \\ approx. growth of the coefficients in general
+
+initLdata("a(k)");      \\ Initialized L-series with coefficients a(k)
+                        \\ actually uses cflength()=352 coeffs
+
+print("EXAMPLE: L(s)=Shintani's zeta function");
+print("         with ",default(realprecision)," digits precision");
+print("Verifying functional equation. Error: ",errprint(checkfeq()));
+print("L(2)     = ",L(2));
+print(" (check) = ",L(2,1.1));
+print("L(3)     = ",L(3));
+print(" (check) = ",L(3,1.1));

sage/ext_data/pari/dokchitser/ex-tau2

@@ -0,0 +1,30 @@
+read("computel");              /* for Steve Finch */
+default(realprecision,40);
+
+coefgrow(n) = 2*n^20;
+tau(n)      = (5*sigma(n,3)+7*sigma(n,5))*n/12\
+              -35*sum(k=1,n-1,(6*k-4*(n-k))*sigma(k,3)*sigma(n-k,5));
+tau2(n)      = tau(n)^2;
+
+conductor = 1;
+gammaV    = [0,1,-10];
+weight    = 23;
+Lpoles    = [];
+Lresidues = automatic;
+
+len=cflength();
+v1=vector(len,n,tau2(n));
+v2=vector(len,n,n^11);
+v3=vector(len,n,if(issquare(n),n^11,0));
+cfs=concat(dirdiv(dirmul(v1,v3),v2),vector(1000,k,0));
+
+initLdata("cfs[k]");
+
+{
+forstep(s0=12.5,22.0,0.5,
+  s0=round(2*s0)/2;
+  print("  L       = ",lval = L(s0));
+  print("  (check) = ",lval2 = L(s0,1.1),"  (err=",errprint(lval-lval2),")");
+  print("f(",s0,") = ",lval/zeta(2*s0-22)*zeta(s0-11));
+);
+}

sage/ext_data/pari/dokchitser/ex-zeta

@@ -0,0 +1,27 @@
+/*** EXAMPLE: Riemann's zeta-function                             ***/
+/*** v1.2, July 2013, questions to tim.dokchitser@bristol.ac.uk   ***/
+/***                                                              ***/
+/*** type \rex-zeta or read("ex-zeta") at pari prompt to run this ***/
+
+read("computel");                 \\ read the ComputeL package
+                                  \\ and set the default values
+default(realprecision,28);        \\ set working precision; used throughout
+
+                            \\ initialize L-function parameters
+conductor = 1;              \\ exponential factor = 1/sqrt(Pi)
+gammaV    = [0];            \\ list of gamma-factors
+weight    = 1;              \\ L(s)=sgn*L(weight-s)
+sgn       = 1;              \\ sign in the functional equation
+Lpoles    = [1];            \\ poles of zeta(s) with Re(s)>0.5
+Lresidues = [-1];           \\ and residues in there
+
+initLdata("1");             \\ initialize the package; all coeffs equal to 1
+
+
+print("EXAMPLE: L(s)=zeta(s), Riemann zeta function");
+print("         with ",default(realprecision)," digits precision");
+print("Verifying functional equation. Error: ",errprint(checkfeq()));
+print("zeta(2)    = ",L(2),     "   using the package");
+print("           = ",L(2,1.1), "   using a different cutoff point");
+print("           = ",zeta(2),  "   built-in pari zeta(s) function");
+print("           = ",Pi^2/6,   "   Pi^2/6");

sage/ext_data/pari/dokchitser/ex-zeta2

@@ -0,0 +1,47 @@
+/*** EXAMPLE: Riemann's zeta-function for large imaginary s         ***/
+/*** (illustration of precision issues when Im(s) is large)         ***/
+/***                                                                ***/
+/*** v1.2, July 2013, questions to tim.dokchitser@bristol.ac.uk     ***/
+/*** type \rex-zeta2 or read("ex-zeta2") at pari prompt to run this ***/
+
+read("computel");                 \\ read the ComputeL package
+                                  \\ and set the default values
+default(realprecision,20);        \\ set working precision; used throughout
+
+                            \\ initialize L-function parameters
+conductor = 1;              \\ exponential factor
+gammaV    = [0];            \\ list of gamma-factors
+weight    = 1;              \\ L(s)=sgn*L(weight-s)
+sgn       = 1;              \\ sign in the functional equation
+Lpoles    = [1];            \\ poles of zeta(s) with Re(s)>0.5
+Lresidues = [-1];           \\ and residues in there
+
+initLdata("1");             \\ initialize the package; all coeffs equal to 1
+
+print("EXAMPLE: L(s)=zeta(s), Riemann zeta function");
+print("         for s with large imaginary part");
+print("         with ",default(realprecision)," digits precision");
+print("Verifying functional equation. Error: ",errprint(checkfeq()));
+print("For large imaginary s there is precision loss and a warning is printed");
+print("");
+print("zeta(1/2+40 I) = ",L(1/2+40*I));
+print("  (pari)       = ",zeta(1/2+40*I));
+print("");
+
+print("Changing MaxImaginaryPart to 40 and re-calculating");
+
+MaxImaginaryPart=40;
+initLdata("1");      \\ re-initialize the package
+                     \\ MaxImaginaryPart changed -> more coeffs necessary
+
+print("zeta(1/2+40 I) = ",L(1/2+40*I));
+print("  (pari)       = ",zeta(1/2+40*I));
+print("");
+
+print("Locating a nearby zero of zeta(s) using Newton-Raphson (5 iterations)");
+z0=1/2+40*I;
+for(k=1,5,z0=1/2+I*imag(z0-L(z0)/L(z0,,1)));
+
+header = concatstr("|zeta(",z0,"| = ");
+print(header,errprint(abs(L(z0))));
+print(StrTab("  (pari)",length(header)-3)," = ",errprint(abs(zeta(z0))));

sage/ext_data/pari/dokchitser/testall

@@ -0,0 +1,13 @@
+#!/usr/bin/env bash
+echo "\\r ex-bsw" | sage -gp
+echo "\\r ex-chgen" | sage -gp
+echo "\\r ex-chqua" | sage -gp
+echo "\\r ex-delta" | sage -gp
+echo "\\r ex-eisen" | sage -gp
+echo "\\r ex-gen2" | sage -gp
+echo "\\r ex-gen3" | sage -gp
+echo "\\r ex-gen4" | sage -gp
+echo "\\r ex-nf" | sage -gp
+echo "\\r ex-shin" | sage -gp
+echo "\\r ex-zeta" | sage -gp
+echo "\\r ex-zeta2" | sage -gp