{{{id=1| K.=NumberField(x^2-x-1) import nosqlite db=nosqlite.Client('/home/psharaba/ECdb').db from psage.ellcurve.minmodel.sqrt5 import canonical_model embs=K.embeddings(RR) /// }}} {{{id=2| #### Eta Function, galois function, falting height, naive height #### def eta(alpha): return (alpha*alpha).trace() def find_small(alpha): while True: plus=a*alpha minus=alpha/a if eta(plus)1: raise NotImplementedError if rank_bound == 0: return 0 Q = v[-1][0] for p in prime_range(sat_bound): if len(Q.division_points(p)) != 0: Q = Q.division_points(p)[0] return Q def regulator_new(E,sat_bound): Q = saturate_rank_01(E,sat_bound) if Q: return Q.height()/2, [Q] else: return int(1), [] def conjectural_sha(E, omega, reg=1, Lstar=1): M = E.tamagawa_product_bsd() sha = RR(sqrt(5))*Lstar*(E.torsion_order())^2/((omega)*reg*M) return sha def comp(weq,lstar): E=EllipticCurve(K,weq) p0,p1=period_stuff(E) omega=p0*p1 R=regulator_new(E,20) Sha=conjectural_sha(E,omega,R[0],lstar) return p0,p1,R[0],Sha /// }}} {{{id=6| def s_func(ainv1): E=EllipticCurve(K,ainv1) D=K(E.discriminant()) t='' for m in embs: e=m(D) if sgn(e)<0: t+='-' else: t+='+' t+=',' return t[:-1] /// }}} {{{id=3| ####Lstar function (with Parallel) #### @parallel(ncpus=6) def l_function(ainv): E=EllipticCurve(K,eval(ainv)) F=EllipticCurve(K,list(E.short_weierstrass_model().a_invariants())) v=F.lseries().dokchitser() f=v.taylor_series(1,6) r_an=0 while abs(f[r_an])<1e-10: r_an += 1 if r_an == 6: raise RuntimeError L=float(f[r_an]) db.N1k.update({'Lstar':L}, ainv1=ainv) print L /// }}} {{{id=20| import psage.modform.hilbert.sqrt5.sqrt5 as sq w = sq.primes_of_bounded_norm(K, 200) from psage.ellcurve.minmodel.sqrt5 import canonical_model def ap(E,p): return E.change_ring(p.residue_field()).trace_of_frobenius() def ap_comp(E): q=[] for i in range(len(w)): u=w[i] if E.has_good_reduction(w[i])==True: if not E.conductor().norm()%norm(w[i])==0: try: q.append(ap(E,w[i])) except: q.append('?') else: q.append('?') else: q.append('?') h=[] h.append(q[0]) h.append(q[2]) h.append(q[1]) for i in range(4,6): h.append(q[i]) for j in range(7,15): h.append(q[j]) h.append(q[3]) for k in range(15,35): h.append(q[k]) h.append(q[6]) for l in range(35,len(q)): h.append(q[l]) return h /// }}} {{{id=11| temp=db('select ainv1,Lstar from N1k ORDER BY N,eta') /// }}} {{{id=13| list(l_function(test)) /// WARNING: Output truncated! full_output.txt 0.359928959498 0.359928959498 0.359928959498 0.359928959498 0.359928959498 0.359928959498 0.396227861967 0.396227861967 0.396227861967 0.396227861967 0.422214159002 0.438646969129 0.422214159002 0.438646969129 0.438646969129 0.438646969129 0.438646969129 0.438646969129 0.438646969129 0.438646969129 0.467538645239 0.467538645239 0.438646969129 0.438646969129 0.493601465332 0.493601465332 0.493601465332 0.493601465332 0.493601465332 0.536078844311 0.493601465332 0.536078844311 0.493601465332 0.493601465332 0.556238984168 0.556238984168 0.536078844311 0.556238984168 0.536078844311 0.556238984168 0.585695878434 0.578909544516 0.585695878434 0.585695878434 0.585695878434 0.578909544516 0.594577551815 0.594577551815 0.60255817498 0.594577551815 0.60255817498 0.60255817498 0.60255817498 0.594577551815 0.594577551815 0.594577551815 0.586973216905 0.586973216905 0.586973216905 ... None), (((u'[a+1,-a-1,a,-28*a-21,-80*a-51]',), {}), None), (((u'[a+1,a+1,a,2*a-1,a-1]',), {}), None), (((u'[1,-a+1,0,-a-1,0]',), {}), None), (((u'[a+1,a+1,a,17*a-26,40*a-66]',), {}), None), (((u'[a,1,a,1797*a-2912,44574*a-72128]',), {}), None), (((u'[a+1,a+1,a,22*a-21,53*a-67]',), {}), None), (((u'[1,-a+1,0,-66*a-56,-333*a-241]',), {}), None), (((u'[1,-a+1,0,-71*a-46,-328*a-269]',), {}), None), (((u'[a,-1,a+1,-a-1,-a]',), {}), None), (((u'[1,a,0,-6*a+2,3*a]',), {}), None), (((u'[1,-a-1,0,-a-1,-a+1]',), {}), None), (((u'[a,1,1,-219*a-138,-2003*a-1237]',), {}), None), (((u'[a+1,-a+1,1,3*a+3,8*a+3]',), {}), None), (((u'[a+1,-a+1,1,-17*a-32,72*a-7]',), {}), None), (((u'[a,a+1,a+1,a,0]',), {}), None), (((u'[a,a+1,a+1,a-5,2*a-5]',), {}), None), (((u'[1,1,a+1,282*a-459,2862*a-4634]',), {}), None), (((u'[a+1,1,0,a+1,1]',), {}), None), (((u'[a+1,-1,1,-235*a-145,2010*a+1241]',), {}), None), (((u'[1,a+1,1,a,0]',), {}), None), (((u'[1,a+1,1,6*a-5,-2*a+6]',), {}), None), (((u'[a,1,1,0,0]',), {}), None), (((u'[0,-1,0,1,0]',), {}), None), (((u'[0,-a+1,0,-34*a-33,95*a+63]',), {}), None), (((u'[0,-a+1,0,-4*a,-4*a]',), {}), None), (((u'[0,1,0,a,-a]',), {}), None), (((u'[0,-1,0,-4,4]',), {}), None), (((u'[0,-a+1,0,-44*a-40,-196*a-112]',), {}), None), (((u'[0,-1,0,-24,-36]',), {}), None), (((u'[0,1,0,-119*a-100,-909*a-556]',), {}), None), (((u'[0,1,0,a-20,-29*a+12]',), {}), None), (((u'[0,-a-1,0,1024*a-1664,-19072*a+30876]',), {}), None), (((u'[0,-1,0,-384,-2772]',), {}), None), (((u'[0,-1,0,-64,220]',), {}), None), (((u'[0,-1,0,16,-180]',), {}), None), (((u'[0,-a+1,0,6*a+7,7*a+7]',), {}), None), (((u'[a,-a,a+1,0,-a-1]',), {}), None), (((u'[a+1,0,0,148*a-200,880*a-1384]',), {}), None), (((u'[a+1,0,0,-7*a+5,8*a-9]',), {}), None), (((u'[1,a-1,0,a,1]',), {}), None), (((u'[a+1,0,a+1,-3*a-3,-5*a-3]',), {}), None), (((u'[1,-a+1,a+1,130*a-216,963*a-1563]',), {}), None), (((u'[a+1,a,0,4*a-1,3*a]',), {}), None), (((u'[1,a+1,1,-a-3,-5*a-5]',), {}), None), (((u'[a+1,-1,a,-2*a-4,3]',), {}), None), (((u'[a+1,0,a+1,-193*a-168,-1862*a-1282]',), {}), None), (((u'[a,a+1,0,971*a-1565,17491*a-28313]',), {}), None), (((u'[a,a-1,a,0,0]',), {}), None), (((u'[0,a+1,a+1,a+1,0]',), {}), None), (((u'[1,1,1,a-1,0]',), {}), None), (((u'[a,a+1,1,-2637*a-1648,-81226*a-50237]',), {}), None), (((u'[a+1,a,a+1,-53*a-32,162*a+100]',), {}), None), (((u'[a,a-1,a+1,-a-1,-a]',), {}), None), (((u'[1,1,1,-9*a+4,-10*a-2]',), {}), None), (((u'[1,-1,0,1,0]',), {}), None), (((u'[1,-1,0,-4,3]',), {}), None), (((u'[a+1,a,0,1,0]',), {}), None), (((u'[1,-1,0,-59,190]',), {}), None), (((u'[1,-1,0,-29,-52]',), {}), None), (((u'[a+1,a,0,-35*a-54,-238*a-211]',), {}), None), (((u'[a+1,a,0,35*a-34,-84*a+119]',), {}), None), (((u'[1,a-1,a+1,-5574*a-3520,-243182*a-150513]',), {}), None), (((u'[a+1,a,0,-45*a-29,-207*a-327]',), {}), None), (((u'[a+1,-a-1,a+1,-a-2,0]',), {}), None), (((u'[a,a,0,8*a-15,11*a+17]',), {}), None), (((u'[a+1,-a-1,a+1,-6*a-17,13*a+28]',), {}), None), (((u'[a,-a-1,0,-1422*a-878,-29934*a-18511]',), {}), None), (((u'[1,1,0,-2,7*a+3]',), {}), None), (((u'[1,1,0,-35*a-52,210*a+51]',), {}), None), (((u'[a+1,a-1,a,1502*a-2457,33106*a-53599]',), {}), None), (((u'[1,1,0,-770*a-497,12257*a+7460]',), {}), None), (((u'[a+1,a-1,a,377*a-607,-4405*a+7107]',), {}), None), (((u'[1,1,0,-260*a-287,-3445*a-1696]',), {}), None), (((u'[1,1,0,-4785*a-3467,-204530*a-130554]',), {}), None), (((u'[1,1,0,-935*a-387,8620*a+7934]',), {}), None), (((u'[a+1,-a-1,1,-a-2,a+1]',), {}), None), (((u'[a,a,0,3*a-195,847*a-500]',), {}), None), (((u'[a+1,-a,1,-a,0]',), {}), None), (((u'[0,a+1,a,a,0]',), {}), None), (((u'[a+1,-a-1,1,-21*a-12,57*a+37]',), {}), None), (((u'[a,-a+1,1,1932*a-3156,50492*a-81759]',), {}), None), (((u'[0,1,a,-76*a-53,375*a+235]',), {}), None), (((u'[a,-a+1,1,-1,a-2]',), {}), None), (((u'[a+1,-a,1,-256*a-156,-2345*a-1447]',), {}), None), (((u'[a+1,-a-1,1,-a,0]',), {}), None), (((u'[0,-1,a,80172*a-129697,13036100*a-21092964]',), {}), None), (((u'[a,-a+1,1,45*a-81,-171*a+276]',), {}), None), (((u'[0,a+1,a+1,2*a-3,-4]',), {}), None), (((u'[a+1,0,a+1,-27*a-17,-92*a-57]',), {}), None), (((u'[a+1,-a,a+1,-2*a,-a]',), {}), None), (((u'[a,-a-1,a+1,-a-1,0]',), {}), None), (((u'[a,-a-1,a+1,-46*a-31,-159*a-88]',), {}), None), (((u'[a,-a-1,a+1,-a-6,-4*a+3]',), {}), None), (((u'[0,-1,a,a-1,-a+1]',), {}), None), (((u'[1,-1,a+1,271*a-438,-2634*a+4258]',), {}), None), (((u'[1,a+1,0,a+1,0]',), {}), None), (((u'[a,1,1,1,0]',), {}), None), (((u'[1,a+1,0,-89*a-64,-620*a-375]',), {}), None), (((u'[a+1,-a+1,0,74*a-123,422*a-683]',), {}), None), (((u'[a+1,-1,a+1,2*a-8,-7*a+10]',), {}), None), (((u'[a+1,0,a+1,a-2,a-1]',), {}), None), (((u'[0,a+1,0,a-2,-4*a-1]',), {}), None), (((u'[a+1,0,a+1,-79*a-32,295*a+191]',), {}), None), (((u'[a+1,0,a+1,-9*a-62,-199*a+39]',), {}), None), (((u'[a+1,0,a+1,-4*a-7,1]',), {}), None), (((u'[0,a+1,0,-159*a-82,-1220*a-753]',), {}), None), (((u'[1,0,1,0,0]',), {}), None), (((u'[a,a,1,0,0]',), {}), None), (((u'[1,-1,1,-3,3]',), {}), None), (((u'[a+1,1,0,-a+4,-2*a+2]',), {}), None), (((u'[1,0,1,-5,-8]',), {}), None), (((u'[1,0,1,-460,-3830]',), {}), None), (((u'[1,-1,1,-213,-1257]',), {}), None), (((u'[1,1,a,a-3,-a+1]',), {}), None), (((u'[a,a+1,a,-52*a-32,164*a+102]',), {}), None), (((u'[a,-a+1,a+1,-6*a-2,4*a+4]',), {}), None), (((u'[a,-a,a+1,-3*a-18,-18*a-3]',), {}), None), (((u'[a+1,a-1,a,-1024*a-651,-19873*a-12304]',), {}), None), (((u'[a,-a+1,a+1,-66*a-52,332*a+214]',), {}), None), (((u'[a,-a,a+1,-3*a+2,-2*a+1]',), {}), None), (((u'[a,-a,a+1,87*a-168,-552*a+879]',), {}), None), (((u'[a,-a,a,-2*a-1,a+1]',), {}), None), (((u'[a,a+1,a+1,6*a-7,-5*a+9]',), {}), None), (((u'[a+1,-a,a,a-4,2*a-4]',), {}), None), (((u'[a,a+1,a+1,-9*a-2,-44*a+20]',), {}), None), (((u'[1,0,a,-16*a-11,-39*a-24]',), {}), None), (((u'[a,0,a,-1771*a-1094,-43255*a-26733]',), {}), None), (((u'[a+1,-a,a,56*a-99,-255*a+409]',), {}), None), (((u'[1,0,a,-a,0]',), {}), None), (((u'[1,0,a,4*a,-6*a-7]',), {}), None), (((u'[0,a-1,0,a-4,-2]',), {}), None), (((u'[0,-a+1,0,-2*a,a+1]',), {}), None), (((u'[0,-a-1,0,0,a+1]',), {}), None), (((u'[0,0,0,-3*a-2,3*a+2]',), {}), None), (((u'[0,a-1,0,6*a-9,-11*a+13]',), {}), None), (((u'[0,a-1,0,36*a-64,120*a-192]',), {}), None), (((u'[0,-a-1,0,5*a-5,-3*a+8]',), {}), None), (((u'[0,a,0,-122*a-75,-843*a-522]',), {}), None), (((u'[0,a,0,44*a-68,148*a-244]',), {}), None), (((u'[0,-a+1,0,-107*a-85,-701*a-410]',), {}), None), (((u'[0,0,0,15*a-26,-36*a+58]',), {}), None), (((u'[0,0,0,-63*a-47,-312*a-194]',), {}), None), (((u'[0,a+1,0,-6*a-5,-18*a-12]',), {}), None), (((u'[0,0,0,255*a-406,-2304*a+3726]',), {}), None), (((u'[0,0,0,-63*a-7,-480*a-306]',), {}), None), (((u'[0,0,0,-1023*a-727,-19856*a-12626]',), {}), None), (((u'[0,1,0,68*a-112,352*a-572]',), {}), None), (((u'[0,a-1,a,-a+1,-1]',), {}), None), (((u'[a+1,-a-1,a,-3,-a+2]',), {}), None), (((u'[0,1,a,-5*a-4,5*a+4]',), {}), None), (((u'[0,a-1,a,3*a-4,-5*a+7]',), {}), None), (((u'[a+1,a+1,a,2*a+1,a+1]',), {}), None), (((u'[1,-a-1,a+1,-7*a-3,14*a+9]',), {}), None), (((u'[1,-a,a+1,4*a-8,-5*a+7]',), {}), None), (((u'[0,-a-1,a,2*a,-a]',), {}), None), (((u'[a+1,a,0,-4681*a-2895,179281*a+110801]',), {}), None), (((u'[1,-a,a+1,-26*a+27,-22*a+59]',), {}), None), (((u'[0,-1,a+1,-2*a-2,2*a+2]',), {}), None), (((u'[a,a-1,0,-42*a-27,142*a+89]',), {}), None), (((u'[0,a,a+1,2*a+2,-a]',), {}), None), (((u'[1,-1,a,2*a-2,a-3]',), {}), None), (((u'[a+1,1,1,3387*a-5493,114900*a-185944]',), {}), None), (((u'[1,-1,a,2*a-77,176*a-28]',), {}), None), (((u'[1,-1,a,27*a-52,91*a-148]',), {}), None), (((u'[1,-1,a,-123*a+173,776*a-703]',), {}), None), (((u'[a,-a-1,a,-3544*a-2272,125529*a+77866]',), {}), None), (((u'[1,-1,a,3452*a-6752,-138149*a+239247]',), {}), None), (((u'[0,-a,a+1,2*a-13,7*a-21]',), {}), None), (((u'[0,0,a+1,-55*a-25,-227*a-103]',), {}), None), (((u'[0,0,1,0,0]',), {}), None), (((u'[0,0,a+1,57415*a-92810,-7909437*a+12798927]',), {}), None), (((u'[1,-1,a+1,-2691333*a-1663335,2541889053*a+1570973833]',), {}), None), (((u'[0,0,a+1,0,-a]',), {}), None), (((u'[0,0,1,0,-7]',), {}), None), (((u'[0,0,1,-30,63]',), {}), None), (((u'[0,0,1,-270,-1708]',), {}), None), (((u'[0,0,a+1,0,6*a-14]',), {}), None), (((u'[1,-1,a,a-3,-2*a+3]',), {}), None), (((u'[a,-a-1,a+1,-4*a-2,5*a+3]',), {}), None), (((u'[a,-a-1,a,-32*a-17,-78*a-48]',), {}), None), (((u'[a,-a-1,a+1,4*a-8,-7*a+10]',), {}), None), (((u'[a+1,-a,0,-4*a-6,9*a+2]',), {}), None), (((u'[1,-1,a+1,13*a-29,33*a-60]',), {}), None), (((u'[a,-a+1,0,180*a-293,1508*a-2439]',), {}), None), (((u'[1,0,a+1,-a,-a]',), {}), None), (((u'[a+1,-1,a+1,5*a-12,-13*a+20]',), {}), None), (((u'[1,a,a+1,2*a-3,-2*a+2]',), {}), None), (((u'[a,-a+1,0,-54*a-29,-214*a-137]',), {}), None), (((u'[1,a,a+1,7*a-13,9*a-17]',), {}), None), (((u'[1,a,a+1,112*a-198,743*a-1188]',), {}), None), (((u'[1,-a+1,a+1,-3*a-2,2*a+1]',), {}), None), (((u'[a+1,0,0,a-1,0]',), {}), None), (((u'[1,-a+1,a+1,-43*a-22,150*a+97]',), {}), None), (((u'[a,-a,0,-2*a-2,3*a+2]',), {}), None), (((u'[1,-a+1,a,-12*a-5,-19*a-11]',), {}), None), (((u'[0,a+1,a,6*a-6,1]',), {}), None), (((u'[0,a,a,a,0]',), {}), None), (((u'[0,1,a,-3,a-1]',), {}), None), (((u'[0,a+1,a,206*a-356,1825*a-3024]',), {}), None), (((u'[a+1,a-1,1,-a-2,-3*a-2]',), {}), None), (((u'[0,a,a,11*a-30,-37*a+71]',), {}), None), (((u'[1,1,a+1,-3413*a-2108,-116137*a-71773]',), {}), None), (((u'[1,-a+1,0,61492*a-99493,-8737966*a+14138322]',), {}), None), (((u'[a,-a,0,556*a-910,-7512*a+12147]',), {}), None), (((u'[0,-a+1,a,-17*a-16,-45*a-30]',), {}), None), (((u'[a+1,a-1,1,-76*a-227,447*a+1348]',), {}), None), (((u'[a,a,1,107*a-173,620*a-1002]',), {}), None), (((u'[a+1,-a+1,1,-3*a+1,-2*a]',), {}), None), (((u'[a+1,-a+1,0,-62*a-76,411*a+171]',), {}), None), (((u'[a+1,-a+1,0,-12*a-251,-694*a+861]',), {}), None), (((u'[1,a+1,0,-8663*a-5338,455336*a+281393]',), {}), None), (((u'[1,a,1,-1,0]',), {}), None), (((u'[1,a+1,1,-143*a-90,-1140*a-705]',), {}), None), (((u'[0,a+1,a,6*a-7,10*a-15]',), {}), None), (((u'[0,a+1,a,436*a-707,5297*a-8596]',), {}), None), (((u'[1,0,a+1,-3*a-2,-3*a-2]',), {}), None), (((u'[1,0,a+1,-38*a-27,-143*a-91]',), {}), None), (((u'[a,-1,a,-2*a,a-1]',), {}), None), (((u'[0,0,0,-16*a-16,44*a+33]',), {}), None), (((u'[a,-1,a,88*a-130,401*a-661]',), {}), None), (((u'[a,-1,a,-2*a-50,-121*a+57]',), {}), None), (((u'[0,0,0,-11*a-11,72*a+54]',), {}), None), (((u'[0,0,0,296*a-487,-2984*a+4830]',), {}), None), (((u'[1,a-1,0,a-2,0]',), {}), None), (((u'[a,-a-1,a,-1,0]',), {}), None), (((u'[1,a+1,0,12*a-31,-34*a+63]',), {}), None), (((u'[a,1,a+1,-a,-a]',), {}), None), (((u'[a,0,a,-a-1,0]',), {}), None), (((u'[a,-a-1,a,5*a-6,4*a-7]',), {}), None), (((u'[1,a+1,0,2*a-1,1]',), {}), None), (((u'[a,0,a,4*a+4,-4*a-5]',), {}), None), (((u'[a+1,-1,0,-2*a,a]',), {}), None), (((u'[a,-a+1,a,23*a-39,-55*a+88]',), {}), None), (((u'[a+1,-a,a,-a,-a+1]',), {}), None), (((u'[a,0,1,-a,0]',), {}), None), (((u'[a+1,-a,a,-164*a-107,-1237*a-753]',), {}), None), (((u'[a,0,1,104*a-310,946*a-2218]',), {}), None), (((u'[1,a-1,0,-12*a-4,22*a+10]',), {}), None), (((u'[a+1,a,a+1,2015*a-3272,51754*a-84002]',), {}), None), (((u'[1,-a+1,0,3*a-7,7*a-10]',), {}), None), (((u'[1,-a+1,0,-7*a+23,-327*a+476]',), {}), None), (((u'[1,-a+1,0,373*a-757,-5395*a+8016]',), {}), None), (((u'[1,-a+1,0,23*a-47,-69*a+102]',), {}), None), (((u'[a,0,0,-23*a-29,77*a+69]',), {}), None), (((u'[a,0,0,-28*a-59,-a-78]',), {}), None), (((u'[1,0,0,721*a-1176,-11175*a+18093]',), {}), None), (((u'[1,0,0,9161*a-15731,534651*a-875694]',), {}), None), (((u'[1,0,0,-1,0]',), {}), None), (((u'[a,0,0,-1348*a-3959,-57589*a-102594]',), {}), None), (((u'[a,0,0,292*a-59,-769*a+882]',), {}), None), (((u'[1,0,0,4,1]',), {}), None), (((u'[a,0,a+1,-a-2,-a]',), {}), None), (((u'[a,0,a+1,-11*a-7,15*a+11]',), {}), None), (((u'[a,a,a+1,-a-3,-3*a-3]',), {}), None), (((u'[a,a-1,a+1,-2*a,-a-1]',), {}), None), (((u'[a+1,-1,a,a-3,-2*a+3]',), {}), None), (((u'[a+1,-a,a,4*a-8,6*a-11]',), {}), None), (((u'[a,a-1,a+1,-97*a-80,588*a+337]',), {}), None), (((u'[1,-a,0,39*a-63,-126*a+207]',), {}), None), (((u'[a+1,-a+1,a,-209*a-128,-1799*a-1112]',), {}), None), (((u'[a+1,-1,a,a-8,2*a-2]',), {}), None), (((u'[a+1,-1,a,46*a-98,209*a-380]',), {}), None), (((u'[1,a,a,-208*a-124,1576*a+971]',), {}), None), (((u'[1,a,a,-1033*a-654,-19938*a-12360]',), {}), None), (((u'[a,a+1,0,6224*a-10064,283403*a-458560]',), {}), None), (((u'[a+1,-a,a,-11*a+7,21*a-44]',), {}), None), (((u'[0,0,0,-a+1,a-2]',), {}), None), (((u'[0,1,0,2*a-2,-a+2]',), {}), None), (((u'[a+1,-a+1,1,a-3,1]',), {}), None), (((u'[0,0,0,-37*a-26,-136*a-86]',), {}), None), (((u'[0,1,0,-18*a+18,11*a-38]',), {}), None), (((u'[0,1,a,a,0]',), {}), None), (((u'[0,a+1,0,-566*a-373,-8475*a-5269]',), {}), None), (((u'[0,1,a,-9*a,-15*a-7]',), {}), None), (((u'[a,-a-1,1,4*a+3,-5*a+2]',), {}), None), (((u'[a,-a-1,1,4*a-42,-41*a+92]',), {}), None), (((u'[a,-a-1,1,-131*a-177,-1121*a-880]',), {}), None), (((u'[a+1,1,a,239756*a-387938,-67486715*a+109195812]',), {}), None), (((u'[a+1,1,a+1,-19683*a-12268,1576908*a+974960]',), {}), None), (((u'[1,-1,1,2160*a-3554,-59299*a+96135]',), {}), None), (((u'[0,a+1,a+1,a-4,-2*a+2]',), {}), None), (((u'[0,0,a+1,a,0]',), {}), None), (((u'[0,-a-1,a+1,15*a-26,-36*a+63]',), {}), None), (((u'[0,-a-1,a+1,35*a-26,-33*a+5]',), {}), None), (((u'[1,0,0,-3,-3]',), {}), None), (((u'[1,1,0,69*a-117,380*a-612]',), {}), None), (((u'[1,1,0,4*a-7,7*a-12]',), {}), None), (((u'[1,1,a,-3*a,-2*a-3]',), {}), None), (((u'[a+1,-a-1,1,-93*a-93,-525*a-394]',), {}), None), (((u'[a,1,1,6*a-9,-6*a+10]',), {}), None), (((u'[1,0,0,-53,-153]',), {}), None), (((u'[1,0,0,-28,272]',), {}), None), (((u'[1,0,0,-828,9072]',), {}), None), (((u'[a+1,-a-1,1,-1443*a-1443,-37245*a-27934]',), {}), None), (((u'[a+1,-a-1,1,7*a+7,-45*a-34]',), {}), None), (((u'[a+1,-a-1,1,-343*a-343,3875*a+2906]',), {}), None), (((u'[a+1,-a-1,1,-1668*a-1668,47355*a+35516]',), {}), None), (((u'[a+1,-a-1,1,-68*a-68,1275*a+956]',), {}), None), (((u'[a+1,-a-1,1,-2268*a-2268,10875*a+8156]',), {}), None), (((u'[a+1,-a-1,1,-26668*a-26668,3007355*a+2255516]',), {}), None), (((u'[0,-a+1,a+1,0,-a]',), {}), None), (((u'[0,-a,a+1,-7*a-3,12*a+10]',), {}), None), (((u'[a+1,-a+1,0,1,0]',), {}), None), (((u'[a+1,-a+1,0,-5*a-14,12*a+11]',), {}), None), (((u'[0,-a,a+1,213*a-463,1401*a-4842]',), {}), None), (((u'[a+1,-a+1,0,5*a-74,34*a-287]',), {}), None), (((u'[a,-a+1,1,0,0]',), {}), None), (((u'[a+1,-a,a,1,0]',), {}), None), (((u'[a,0,a+1,-a-1,-a]',), {}), None), (((u'[a,a,0,a,0]',), {}), None), (((u'[a,0,a+1,-a+4,7*a-13]',), {}), None), (((u'[a,a,0,-4*a,-6*a-5]',), {}), None), (((u'[0,-a-1,0,-a+1,0]',), {}), None), (((u'[0,-a-1,1,-4*a-2,10*a+6]',), {}), None), (((u'[0,-a-1,0,4*a-4,-4*a]',), {}), None), (((u'[a,1,1,-2*a-3,a+1]',), {}), None), (((u'[a+1,-a+1,a+1,3*a-8,-4*a+5]',), {}), None), (((u'[a,-a-1,1,-a-2,-1]',), {}), None), (((u'[0,a,1,2,a-2]',), {}), None), (((u'[0,-a-1,1,-1,2*a+1]',), {}), None), (((u'[1,-1,1,85*a-139,434*a-702]',), {}), None), (((u'[0,-a,1,-42*a-53,-192*a-140]',), {}), None), (((u'[a+1,-a,0,7*a-23,-24*a+20]',), {}), None), (((u'[1,0,0,-594*a-390,-8580*a-5345]',), {}), None), (((u'[a,0,0,-66265*a-40952,-9832128*a-6076598]',), {}), None), (((u'[a+1,-a,0,652*a-2048,27054*a-32629]',), {}), None), (((u'[1,a,a+1,834227*a-1349799,438442396*a-709414712]',), {}), None), (((u'[a+1,-a,0,1017*a-1853,21096*a-35088]',), {}), None), (((u'[a+1,-1,1,-2,-a+1]',), {}), None), (((u'[a+1,a-1,a,-4*a-1,a]',), {}), None), (((u'[a,a,a,-a-1,-2*a+1]',), {}), None), (((u'[a+1,0,a+1,9*a-17,-19*a+29]',), {}), None), (((u'[1,-1,1,2,-3]',), {}), None), (((u'[a,-a,a+1,2682*a-4688,86629*a-142543]',), {}), None), (((u'[a+1,a-1,a,-9*a-16,-32*a-29]',), {}), None), (((u'[1,-a+1,1,-84*a-52,-434*a-263]',), {}), None), (((u'[a+1,-a+1,a,-2*a-2,a-1]',), {}), None), (((u'[1,-1,1,-18,-19]',), {}), None), (((u'[1,-1,1,-268,-1619]',), {}), None), (((u'[1,-1,1,-88,317]',), {}), None), (((u'[a+1,1,1,0,0]',), {}), None), (((u'[a,1,a,-3*a-3,2*a]',), {}), None), (((u'[a,-1,1,5*a-16,16*a-18]',), {}), None), (((u'[1,-a-1,a+1,69*a-112,327*a-530]',), {}), None), (((u'[a,-1,1,-1,0]',), {}), None)] }}} {{{id=7| for r in range(len(temp)): A=temp[r] weq=eval(A[0]) f=f_height(weq) n=n_height(weq) s=s_func(weq) db.N1k.update({'fh':float(f), 'nh': float(n), 'sign': s},ainv1=A[0]) /// }}} {{{id=22| for r in range(len(temp)): print temp[r][1] /// WARNING: Output truncated! full_output.txt 0.359928959498 0.359928959498 0.359928959498 0.359928959498 0.359928959498 0.359928959498 0.396227861967 0.396227861967 0.396227861967 0.396227861967 0.422214159002 0.422214159002 0.438646969129 0.438646969129 0.438646969129 0.438646969129 0.438646969129 0.438646969129 0.438646969129 0.438646969129 0.438646969129 0.438646969129 0.467538645239 0.467538645239 0.493601465332 0.493601465332 0.493601465332 0.493601465332 0.493601465332 0.493601465332 0.493601465332 0.493601465332 0.536078844311 0.536078844311 0.536078844311 0.536078844311 0.556238984168 0.556238984168 0.556238984168 0.556238984168 0.585695878434 0.578909544516 0.585695878434 0.585695878434 0.585695878434 0.578909544516 0.60255817498 0.60255817498 0.60255817498 0.60255817498 0.594577551815 0.594577551815 0.594577551815 0.594577551815 0.594577551815 0.594577551815 0.586973216905 0.586973216905 0.586973216905 ... 1.26637346504 1.192564546 1.26637346504 1.192564546 1.192564546 1.33877798516 1.33877798516 1.51886547229 1.30306483839 1.37465043752 1.38530142966 1.44234314838 1.65222824075 1.44234314838 1.40058926314 1.39990300788 1.07349782684 1.39990300788 1.39990300788 1.39990300788 1.20197101283 1.31427048811 1.37348086034 1.20197101283 1.20197101283 1.37348086034 1.24990546388 1.40469137601 1.41253028804 1.70660373709 1.41253028804 1.70660373709 1.41446723275 1.39664575527 1.39664575527 1.56034162366 1.40546805065 1.46907610703 1.56034162366 1.44375302917 1.46590168838 1.36763058107 1.42805329888 1.19107564712 1.40853613745 1.46999753154 1.19107564712 1.19107564712 0.929961188229 0.929961188229 1.30232485775 1.24604293297 1.24604293297 1.24604293297 1.24604293297 1.416799465 1.64885551274 1.64885551274 1.48394747001 1.48394747001 }}} {{{id=29| for r in range(len(temp)): try: r0,r1,R,S=comp(eval(temp[r][0]),temp[r][1]) db.N1k.update({'real': (float(r0),float(r1)), 'reg': float(R), 'Sha': float(S)}, ainv1=temp[r][0]) print r print 'found sha' except NotImplementedError: E=EllipticCurve(K,eval(temp[r][0])) r0,r1=period_stuff(E) R='unknown' Sha='unknown' db.N1k.update({'real': (float(r0),float(r1)), 'reg': R, 'Sha': S}, ainv1=temp[r][0]) print r print 'cannot compute sha' except IndexError: E=EllipticCurve(K,eval(temp[r][0])) r0,r1=period_stuff(E) R='unknown' Sha='unknown' db.N1k.update({'real': (float(r0),float(r1)), 'reg': R, 'Sha': S}, ainv1=temp[r][0]) print r print r print 'has index error' /// WARNING: Output truncated! 0 found sha 1 found sha 2 found sha 3 found sha 4 found sha 5 found sha 6 found sha 7 found sha 8 found sha 9 found sha 10 found sha 11 found sha 12 found sha 13 found sha 14 found sha 15 found sha 16 found sha 17 found sha 18 cannot compute sha 19 found sha 20 cannot compute sha 21 cannot compute sha 22 found sha 23 found sha 24 found sha 25 found sha 26 found sha 27 found sha 28 found sha 29 ... 207 found sha 208 found sha 209 found sha 210 found sha 211 found sha 212 found sha 213 found sha 214 found sha 215 found sha 216 found sha 217 found sha 218 found sha 219 found sha 220 found sha 221 found sha 222 found sha 223 found sha 224 found sha 225 found sha 226 found sha 227 found sha 228 found sha 229 found sha 230 found sha 231 found sha 232 found sha 233 found sha 234 found sha 235 found sha 236 found sha }}} {{{id=38| temp[202] /// (u'[0,1,1,-131,-650]', 0.73638405757892988) }}} {{{id=37| comp(eval(temp[202][0]),temp[202][1]) /// Traceback (most recent call last): File "", line 1, in File "_sage_input_48.py", line 10, in exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8 -*-\\n" + _support_.preparse_worksheet_cell(base64.b64decode("Y29tcChldmFsKHRlbXBbMjAyXVswXSksdGVtcFsyMDJdWzFdKQ=="),globals())+"\\n"); execfile(os.path.abspath("___code___.py"))' + '\n', '', 'single') File "", line 1, in File "/tmp/tmpKIIuz_/___code___.py", line 3, in exec compile(u'comp(eval(temp[_sage_const_202 ][_sage_const_0 ]),temp[_sage_const_202 ][_sage_const_1 ])' + '\n', '', 'single') File "", line 1, in File "/tmp/tmpqx8RdZ/___code___.py", line 51, in comp R=regulator_new(E,_sage_const_20 ) File "/tmp/tmpqx8RdZ/___code___.py", line 36, in regulator_new Q = saturate_rank_01(E,sat_bound) File "/tmp/tmpqx8RdZ/___code___.py", line 28, in saturate_rank_01 Q = v[-_sage_const_1 ][_sage_const_0 ] IndexError: list index out of range }}} {{{id=36| comp((temp[200][0]), comp[200][1]) /// Traceback (most recent call last): File "", line 1, in File "_sage_input_45.py", line 10, in exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8 -*-\\n" + _support_.preparse_worksheet_cell(base64.b64decode("Y29tcCgodGVtcFsyMDBdWzBdKSwgY29tcFsyMDBdWzFdKQ=="),globals())+"\\n"); execfile(os.path.abspath("___code___.py"))' + '\n', '', 'single') File "", line 1, in File "/tmp/tmp0MRWlV/___code___.py", line 3, in exec compile(u'comp((temp[_sage_const_200 ][_sage_const_0 ]), comp[_sage_const_200 ][_sage_const_1 ])' + '\n', '', 'single') File "", line 1, in TypeError: 'function' object is unsubscriptable }}} {{{id=35| comp(temp[202][0],comp[202][1]) /// Traceback (most recent call last): File "", line 1, in File "_sage_input_34.py", line 10, in exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8 -*-\\n" + _support_.preparse_worksheet_cell(base64.b64decode("Y29tcCh0ZW1wWzIwMl1bMF0sY29tcFsyMDJdWzFdKQ=="),globals())+"\\n"); execfile(os.path.abspath("___code___.py"))' + '\n', '', 'single') File "", line 1, in File "/tmp/tmpQTYzW6/___code___.py", line 3, in exec compile(u'comp(temp[_sage_const_202 ][_sage_const_0 ],comp[_sage_const_202 ][_sage_const_1 ])' + '\n', '', 'single') File "", line 1, in TypeError: 'function' object is unsubscriptable }}} {{{id=34| EllipticCurve(K,eval(temp[141][0])).gens() /// [(-1 : 0 : 1)] }}} {{{id=32| sha_comp(temp[0]) /// Traceback (most recent call last): File "", line 1, in File "_sage_input_19.py", line 10, in exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8 -*-\\n" + _support_.preparse_worksheet_cell(base64.b64decode("c2hhX2NvbXAodGVtcFswXSk="),globals())+"\\n"); execfile(os.path.abspath("___code___.py"))' + '\n', '', 'single') File "", line 1, in File "/tmp/tmpI4UBvb/___code___.py", line 3, in exec compile(u'sha_comp(temp[_sage_const_0 ])' + '\n', '', 'single') File "", line 1, in File "/usr/local/sage/local/lib/python2.6/site-packages/sage/parallel/decorate.py", line 116, in g return f(*args, **kwds) TypeError: sha_comp() takes exactly 3 arguments (1 given) }}} {{{id=21| temp=db('select rlow,rhi,ainv1 from N1k ORDER BY N') for s in range(len(temp)): rl=temp[s][0] rh=temp[s][1] if rl==rh: db.N1k.update({'rank': rl},ainv1=temp[s][2]) if rl != rh: db.N1k.update({'rank': (rl,rh)},ainv1=temp[s][2]) /// (u'[1,1,1,-2160,-39540]', 0.4386469691285414) }}} {{{id=30| test=db('select N,ainv1 from N1k ORDER BY N') for t in range(len(test)): E=EllipticCurve(K,eval(test[t][1])) a=ap_comp(E) db.N1k.update({'a_p':str(a).replace(' ','')}, ainv1=test[t][1]) /// }}} {{{id=23| comp([1,1,1,-2160,-39540], 0.4386469691285414) /// Interrupting PARI/GP interpreter... Traceback (most recent call last): File "", line 1, in File "_sage_input_27.py", line 10, in exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8 -*-\\n" + _support_.preparse_worksheet_cell(base64.b64decode("Y29tcChbMSwxLDEsLTIxNjAsLTM5NTQwXSwgMC40Mzg2NDY5NjkxMjg1NDE0KQ=="),globals())+"\\n"); execfile(os.path.abspath("___code___.py"))' + '\n', '', 'single') File "", line 1, in File "/tmp/tmpJ2OBcW/___code___.py", line 3, in exec compile(u'comp([_sage_const_1 ,_sage_const_1 ,_sage_const_1 ,-_sage_const_2160 ,-_sage_const_39540 ], _sage_const_0p4386469691285414 )' + '\n', '', 'single') File "", line 1, in File "/tmp/tmpS3F0fM/___code___.py", line 51, in comp R=regulator_new(E,_sage_const_20 ) File "/tmp/tmpS3F0fM/___code___.py", line 36, in regulator_new Q = saturate_rank_01(E,sat_bound) File "/tmp/tmpS3F0fM/___code___.py", line 17, in saturate_rank_01 v = E.simon_two_descent() File "/usr/local/sage/local/lib/python2.6/site-packages/sage/schemes/elliptic_curves/ell_number_field.py", line 276, in simon_two_descent maxprob=maxprob, limbigprime=limbigprime) File "/usr/local/sage/local/lib/python2.6/site-packages/sage/schemes/elliptic_curves/gp_simon.py", line 100, in simon_two_descent s = gp.eval('ans=%s;'%cmd) File "/usr/local/sage/local/lib/python2.6/site-packages/sage/interfaces/expect.py", line 1054, in eval return '\n'.join([self._eval_line(L, **kwds) for L in code.split('\n') if L != '']) File "/usr/local/sage/local/lib/python2.6/site-packages/sage/interfaces/gp.py", line 393, in _eval_line wait_for_prompt=wait_for_prompt) File "/usr/local/sage/local/lib/python2.6/site-packages/sage/interfaces/expect.py", line 735, in _eval_line self._keyboard_interrupt() File "/usr/local/sage/local/lib/python2.6/site-packages/sage/interfaces/expect.py", line 754, in _keyboard_interrupt raise KeyboardInterrupt, "Ctrl-c pressed while running %s"%self KeyboardInterrupt: Ctrl-c pressed while running PARI/GP interpreter __SAGE__ }}} {{{id=28| /// }}}