Hello,
i work with sage for windows vmware, (my PC configuration is Intel
dual core 2, Ram=2GO),
i want to compute the trace of Hecke operator for a large prime
number.
My program is:
"
borneSup=103196184763
S=0
Bm=0
alpha=1
p=2
m1=2
while (p^alpha)<=borneSup:
while (p^alpha) <=borneSup:
m1=p^alpha
if gcd(m1,6)==1:
mangoldt=ln(p)/( m1^(3/2) )
x=PolynomialRing(QQ,'x').gen()
f=gen_laguerre(n-1,1,x)
j=log(m1)
Lg=float(f(j))
M = ModularSymbols(6,2)
if gcd(p,6)==1:
if alpha ==1:
T1 = M.hecke_operator(p)
Bm=T1.trace()
if alpha>=2:
T1 = M.hecke_operator(p^alpha)
trace1=T1.trace()
T2 = M.hecke_operator(p^(alpha-2))
trace2=T2.trace()
Bm=trace1-((p)*trace2)
S=S+(mangoldt*Bm*Lg)
alpha=alpha + 1
Bm=0
alpha=1
p=next_prime(p)
print "S=",float(S)
After a few minutes of program runnig i have the following error
message:
Traceback (most recent call last):
File "/home/sage/sage/local/lib/python2.5/site-packages/sage/
calculus/calculus.py", line 4955, in __float__
return self._convert(float)
File "/home/sage/sage/local/lib/python2.5/site-packages/sage/
calculus/calculus.py", line 4935, in _convert
fops = [typ(op) for op in self._operands]
can some one help me please.
--~--~---------~--~----~------------~-------~--~----~
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
-~----------~----~----~----~------~----~------~--~---
