roed314 / OMS

Converting Rob Pollack's overconvergent modular symbols .sage code to .py and putting it in Sage
4 stars 4 forks source link

Tq_eigenvalue isn't working when the eigenvalue is 0 #41

Closed rpollack9974 closed 11 years ago

rpollack9974 commented 11 years ago

sage: E = EllipticCurve('11a') sage: phi = E.PS_modular_symbol() sage: Phi = phi.lift(11,5,eigensymbol=True) sage: E.ap(19) 0

sage: Phi.Tq_eigenvalue(19)

ValueError Traceback (most recent call last)

/Users/rpollack/git/OMS/sage/modular/pollack_stevens/ in ()

/Applications/sage-5.6/local/lib/python2.7/site-packages/sage/misc/cachefunc.so in sage.misc.cachefunc.CachedMethodCaller.call (sage/misc/cachefunc.c:6564)()

/Applications/sage-5.6/local/lib/python2.7/site-packages/sage/misc/cachefunc.so in sage.misc.cachefunc.CachedMethod._instance_call (sage/misc/cachefunc.c:9045)()

/Applications/sage-5.6/local/lib/python2.7/site-packages/sage/modular/pollack_stevens/modsym.pyc in Tq_eigenvalue(self, q, p, M, check) 546 except IndexError: 547 raise ValueError("self is zero") --> 548 aq = self._map[g].find_scalar(qhecke._map[g], p, M, check) 549 if check: 550 verbose("Checking that this is actually an eigensymbol")

/Applications/sage-5.6/local/lib/python2.7/site-packages/sage/modular/pollack_stevens/dist.so in sage.modular.pollack_stevens.dist.Dist.find_scalar (sage/modular/pollack_stevens/dist.c:5409)() 270 n = self.precision_relative() 271 if n != other.precision_relative(): --> 272 raise ValueError("other should have the same number of moments") 273 verbose("n = %s"%n) 274 verbose("moment 0")

ValueError: other should have the same number of moments