Hi,
I have a problem using Dokchitser's L-functions Calculator. The given
example for Delta function works well, but my own programming fails
all the time. I am getting a strange error message "Unable to create
L-series, due to precision or other limits in PARI" and cannot figure
out what to do with it. For example:
sage: pari_precode = '{tau(n)= return((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)));a(n)=return(sum(d=1,n,if((floor(n*1.0/d)*d)==n,d*tau(d)*tau(n/
d))));}'
sage: L2 = Dokchitser(conductor=1, gammaV=[-1,0,0,1], weight=13,
eps=1,prec=10)
sage: print "\nneed",L2.num_coeffs(),"coefficients"
sage: L2.init_coeffs('a(k)', pari_precode=pari_precode)
need 15 coefficients
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "/home/kvankov/.sage/sage_notebook/worksheets/admin/19/code/
3.py", line 9, in <module>
exec compile(ur'L2.init_coeffs(\u0027a(k)\u0027,
pari_precode=pari_precode)' + '\n', '', 'single')
File "/usr/local/sage-2.11-ubuntu32-intel-i686-Linux/local/lib/
python2.5/site-packages/sympy/plotting/", line 1, in <module>
File "/usr/local/sage-2.11-ubuntu32-intel-i686-Linux/local/lib/
python2.5/site-packages/sage/lfunctions/dokchitser.py", line 286, in
init_coeffs
self._gp_eval('initLdata("%s", %s)'%(v, cutoff))
File "/usr/local/sage-2.11-ubuntu32-intel-i686-Linux/local/lib/
python2.5/site-packages/sage/lfunctions/dokchitser.py", line 220, in
_gp_eval
raise RuntimeError, "Unable to create L-series, due to precision
or other limits in PARI."
RuntimeError: Unable to create L-series, due to precision or other
limits in PARI.
The same example works perfectly fine in PARI/GP :
GP/PARI CALCULATOR Version 2.4.1 (alpha)
i686 running linux (ix86 kernel) 32-bit version
compiled: Oct 23 2007, gcc-4.1.2 (Ubuntu 4.1.2-0ubuntu4)
(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
? read("computel");
? conductor = 1;
? gammaV = [-1,0,0,1];
? weight = 13;
? sgn = 1;
? tau(n)= return((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)));
? a(n)=return(sum(d=1,n,if((floor(n*1.0/d)*d)==n,d*tau(d)*tau(n/d))));
? initLdata("a(k)");
? print(" L(2) = ",lval = L(2));
L(2) = 0.005480450400989257021755917461
May be I need to add something in "pari_precode"? Any suggestions?
Kirill
PS. I use Ubuntu binary SAGE 2.11 on Dell Latitude D800 laptop with
512MB memory. The same issue was with SAGE 2.10.
--~--~---------~--~----~------------~-------~--~----~
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to [EMAIL PROTECTED]
For more options, visit this group at
http://groups.google.com/group/sage-support
URLs: http://www.sagemath.org
-~----------~----~----~----~------~----~------~--~---
