{{{id=1| from psage.ellcurve.minmodel.sqrt5 import canonical_model K.=NumberField(x^2-x-1) import nosqlite db=nosqlite.Client(DATA+'localdb').db /// }}} {{{id=2| def eta(alpha): return (alpha*alpha).trace() def find_small(alpha): while True: plus=a*alpha minus=alpha/a if eta(plus)", line 1, in File "_sage_input_16.py", line 10, in exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8 -*-\\n" + _support_.preparse_worksheet_cell(base64.b64decode("Zi5yZWFkbGluZSgpLnNwbGl0KCk="),globals())+"\\n"); execfile(os.path.abspath("___code___.py"))' + '\n', '', 'single') File "", line 1, in File "/tmp/tmp0KhSuQ/___code___.py", line 2, in exec compile(u'f.readline().split()' + '\n', '', 'single') File "", line 1, in File "element.pyx", line 328, in sage.structure.element.Element.__getattr__ (sage/structure/element.c:2790) File "parent.pyx", line 277, in sage.structure.parent.getattr_from_other_class (sage/structure/parent.c:2930) File "parent.pyx", line 175, in sage.structure.parent.raise_attribute_error (sage/structure/parent.c:2699) AttributeError: 'sage.rings.real_mpfr.RealNumber' object has no attribute 'readline' }}} {{{id=7| db.CLI.delete() /// }}} {{{id=5| for s in f.readlines(): a1=s.split() a2=eval(a1[2]) f=f_height(a2) n=n_height(a2) db.CLI.insert({'N':int(a1[0]), 'cond':a1[1], 'weq':a1[2], 'rlow':a1[3], 'rhi':a1[4], 'sign':a1[5], 'tor':a1[6], 'ordD':a1[7], 'ordj':a1[8], 'tama':a1[9], 'kod':a1[10], 'mtrx':a1[11], 'found': a1[12], 'fh':float(f), 'nh':float(n)}) /// }}} {{{id=6| test=db('select N,cond,weq from CLI ORDER BY N') /// }}} {{{id=10| len(test) /// 232 }}} {{{id=8| for s in range(len(test)): E=EllipticCurve(K,eval(test[s][2])) n=int(E.conductor().norm()) print (test[s][0],n), test[s][2] /// WARNING: Output truncated! full_output.txt (124, 124) [a+1,a-1,a+1,5*a-16,4*a-21] (124, 124) [a+1,a-1,a+1,-35*a-56,-260*a-117] (124, 124) [a+1,a-1,a+1,-1,-a+1] (124, 124) [a,-a,a,4880*a-7979,199505*a-323066] (124, 124) [a,a+1,1,-188901*a-116781,-47466544*a-29336020] (124, 124) [a,-a,a,85*a-139,-436*a+705] (144, 144) [0,-a+1,0,2*a+3,-13*a-9] (144, 144) [0,-a+1,0,-3*a-2,-2*a-2] (144, 144) [0,-a+1,0,-48*a-32,-188*a-116] (144, 144) [0,a,0,48*a-80,188*a-304] (145, 145) [0,a+1,a,-9*a-9,-31*a-23] (145, 145) [0,a+1,a,a+1,0] (179, 179) [a+1,a-1,0,20*a-33,42*a-71] (179, 179) [1,1,a,1,0] (180, 180) [1,0,1,-19,26] (180, 180) [1,0,1,-69,-194] (180, 180) [1,0,1,1,2] (180, 180) [1,0,1,-289,1862] (180, 180) [1,0,1,-334,-2368] (180, 180) [1,0,1,-5334,-150368] (180, 180) [1,0,1,-14,-64] (180, 180) [1,0,1,-454,-544] (196, 196) [1,0,1,4,-6] (196, 196) [1,0,1,-36,-70] (196, 196) [1,0,1,-171,-874] (196, 196) [1,0,1,-1,0] (196, 196) [1,0,1,-2731,-55146] (196, 196) [1,0,1,-11,12] (199, 199) [0,a+1,1,-9*a,-16*a-15] (199, 199) [0,-a-1,1,-28663*a-17718,-2768362*a-1710944] (199, 199) [0,a+1,1,a,0] (205, 205) [a+1,0,a,8*a+4,6*a+1] (205, 205) [a,0,a,104*a-180,656*a-1064] (205, 205) [a+1,0,a,-2*a-1,0] (205, 205) [a+1,0,a,-22*a-16,48*a+31] (209, 209) [a+1,-a+1,a+1,-2*a-5,1] (209, 209) [a+1,-a+1,a+1,3*a-15,14*a-26] (209, 209) [1,a+1,a+1,-250*a-165,2099*a+1308] (209, 209) [a+1,-a+1,a+1,-2*a,-a] (209, 209) [a+1,-a+1,a+1,-47*a,-67*a-159] (209, 209) [a+1,-a+1,a+1,133*a-190,851*a-1361] (225, 225) [a+1,a,a+1,-25*a-25,-119*a-83] (225, 225) [a+1,a,a+1,-400*a-400,-6044*a-4433] (225, 225) [a+1,a,a+1,0,-4*a-3] (225, 225) [a+1,a,a+1,-50*a-50,46*a+47] (225, 225) [a+1,a,a+1,175*a+175,1081*a+767] (225, 225) [a+1,a,a+1,-675*a-675,11171*a+8547] (225, 225) [a+1,a,a+1,-550*a-550,15946*a+12097] (225, 225) [a+1,a,a+1,-10800*a-10800,758396*a+571497] (225, 225) [a,1,1,65146314*a-105408963,-302615596781*a+489642321128] (225, 225) [a,a,0,-9504680*a-5874303,16858314825*a+10419011887] (239, 239) [1,-a,a,-6*a-5,-7*a-4] (239, 239) [1,-a,a,-91*a-85,-585*a-309] (239, 239) [1,-a,a,-a-5,-17*a-15] (239, 239) [1,-a,a,-a,0] (244, 244) [a,-a+1,1,4*a-9,8*a-20] (244, 244) [a,-a+1,1,-a+1,0] (244, 244) [1,a,a,3149*a-5143,102292*a-165671] (275, 275) [a,a,a,332*a-536,-3529*a+5713] ... (649, 649) [1,a+1,0,-89*a-64,-620*a-375] (649, 649) [1,a+1,0,a+1,0] (649, 649) [a+1,-a+1,0,74*a-123,422*a-683] (704, 704) [0,-a+1,0,-7*a-5,-9*a-6] (704, 704) [0,-a+1,0,-107*a-85,-701*a-410] (704, 704) [0,a,0,44*a-68,148*a-244] (704, 704) [0,-a+1,0,-2*a,a+1] (719, 719) [0,-a-1,a,2*a-10,-13] (719, 719) [0,-a-1,a,2*a,-a] (725, 725) [0,a-1,a+1,148*a-246,-1151*a+1868] (725, 725) [0,a,a+1,2*a+2,-a] (755, 755) [1,0,a+1,4*a-5,7*a-15] (755, 755) [1,0,a+1,-a,-a] (775, 775) [a+1,a-1,1,-26*a-27,-133*a-62] (775, 775) [1,1,a+1,-3413*a-2108,-116137*a-71773] (775, 775) [a,-a,0,556*a-910,-7512*a+12147] (775, 775) [a+1,a-1,1,-a-2,-3*a-2] (775, 775) [1,-a+1,0,61492*a-99493,-8737966*a+14138322] (775, 775) [a+1,a-1,1,-76*a-227,447*a+1348] (779, 779) [0,a+1,a,6*a-7,10*a-15] (779, 779) [0,a+1,a,436*a-707,5297*a-8596] (781, 781) [a,-1,a,3*a-10,4*a-8] (781, 781) [a,-1,a,-2*a-50,-121*a+57] (781, 781) [a,-1,a,-2*a,a-1] (781, 781) [a,-1,a,88*a-130,401*a-661] (801, 801) [1,a-1,0,-4*a+8,7*a-10] (801, 801) [1,a-1,0,a-2,0] (829, 829) [a+1,-1,0,3*a-20,4*a-31] (829, 829) [a+1,-1,0,-2*a,a] (836, 836) [a,0,1,4*a,2*a-6] (836, 836) [a,0,1,-a,0] (836, 836) [a,0,1,104*a-310,946*a-2218] (855, 855) [a,a-1,a+1,-7*a-5,6*a+4] (855, 855) [1,-a,0,39*a-63,-126*a+207] (855, 855) [a,a-1,a+1,-2*a,-a-1] (855, 855) [a,a-1,a+1,-97*a-80,588*a+337] (869, 869) [a+1,-a+1,1,-9*a-8,12*a+11] (869, 869) [a+1,-a+1,1,a-3,1] (880, 880) [0,1,0,27*a-47,-79*a+127] (880, 880) [0,1,0,2*a-2,-a+2] (880, 880) [0,a+1,0,-566*a-373,-8475*a-5269] (880, 880) [0,1,0,-18*a+18,11*a-38] (899, 899) [1,1,a,7*a-15,15*a-32] (899, 899) [1,1,a,-3*a,-2*a-3] (899, 899) [a+1,a+1,1,-5*a-4,-8*a-3] (899, 899) [a,1,1,6*a-9,-6*a+10] (900, 900) [a+1,-a-1,1,-93*a-93,-525*a-394] (900, 900) [a+1,-a-1,1,-1443*a-1443,-37245*a-27934] (900, 900) [a+1,-a-1,1,7*a+7,-45*a-34] (900, 900) [a+1,-a-1,1,-343*a-343,3875*a+2906] (900, 900) [a+1,-a-1,1,-1668*a-1668,47355*a+35516] (900, 900) [a+1,-a-1,1,-2268*a-2268,10875*a+8156] (900, 900) [a+1,-a-1,1,-68*a-68,1275*a+956] (900, 900) [a+1,-a-1,1,-26668*a-26668,3007355*a+2255516] (905, 905) [a+1,-a+1,0,-4,a-6] (905, 905) [a+1,-a+1,0,5*a-74,34*a-287] (905, 905) [a+1,-a+1,0,-5*a-14,12*a+11] (905, 905) [a+1,-a+1,0,1,0] (979, 979) [a+1,-a+1,a,3*a-7,6] (979, 979) [a+1,-a+1,a,-2*a-2,a-1] }}} {{{id=9| X=db('select weq from QTF ORDER BY N') @parallel(8) def comp(ainv): E=EllipticCurve(K,eval(ainv)) E=canonical_model(E.global_minimal_model()) 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 Lstar=float(f[r_an]) d = {'Lstar':Lstar} db.QTF.update(d, weq=ainv) print d for s in range(len(X)): A=X[s][0] comp(A) print A /// Traceback (most recent call last): v=F.lseries().dokchitser() File "", line 1, in File "/tmp/tmp6IWfM6/___code___.py", line 2, in X=db('select weq from QTF ORDER BY N') File "/home/wstein/nosqlite/nosqlite.py", line 592, in __call__ return self.client(cmds, t, file=self.name, many=many, coerce=coerce) File "/home/wstein/nosqlite/nosqlite.py", line 467, in __call__ return self.server.execute(cmd, t, file, many) File "/home/wstein/nosqlite/nosqlite.py", line 350, in execute o = cursor.execute(c) sqlite3.OperationalError: no such table: QTF }}} {{{id=12| E=EllipticCurve(K,[1,a+1,a,a,0]);E /// Elliptic Curve defined by y^2 + x*y + a*y = x^3 + (a+1)*x^2 + a*x over Number Field in a with defining polynomial x^2 - x - 1 }}} {{{id=13| E=canonical_model(E.global_minimal_model());E /// Elliptic Curve defined by y^2 + x*y + a*y = x^3 + (a+1)*x^2 + a*x over Number Field in a with defining polynomial x^2 - x - 1 }}} {{{id=14| F=EllipticCurve(K,list(E.short_weierstrass_model().a_invariants()));F /// Elliptic Curve defined by y^2 = x^3 + (432*a-1107)*x + (-4752*a+11502) over Number Field in a with defining polynomial x^2 - x - 1 }}} {{{id=11| @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.CLI.update({'Lstar':L}, weq=ainv) print L /// }}} {{{id=15| test=db('SELECT weq from CLI ORDER BY N') /// }}} {{{id=17| len(test) /// 232 }}} {{{id=16| list(l_function(test2)) /// 1.35749642859 1.35749642859 1.06466608709 1.06466608709 1.11478897665 1.06466608709 1.11478897665 1.17955284361 1.17955284361 1.30711887696 1.30711887696 1.30711887696 1.30711887696 1.30711887696 1.30711887696 1.215230822 1.35877336643 1.35877336643 1.215230822 1.215230822 1.37909433053 1.37909433053 1.05590429426 1.37909433053 1.37909433053 1.05590429426 1.19943679382 1.19943679382 1.19943679382 1.19943679382 1.3175102718 1.3175102718 1.3175102718 1.3175102718 1.5258567169 1.5258567169 1.24091429536 1.24091429536 1.28518138673 1.28518138673 1.50855362857 1.50855362857 1.50855362857 1.50855362857 1.50855362857 1.50855362857 1.581657216 1.581657216 1.03019431847 1.03019431847 1.03019431847 1.03019431847 1.28840002385 1.28840002385 1.27353814662 1.27353814662 1.51085620887 1.51085620887 1.51085620887 1.34816930722 1.34816930722 1.34816930722 1.34816930722 1.32299818988 1.32299818988 1.192564546 1.192564546 1.192564546 1.192564546 1.07349782684 1.65222824075 1.65222824075 1.07349782684 1.07349782684 1.40058926314 1.40058926314 1.07349782684 1.07349782684 1.07349782684 1.07349782684 1.07349782684 1.20197101283 1.20197101283 1.20197101283 1.30232485775 1.20197101283 1.30232485775 [(((u'[0,-a+1,0,-34*a-33,95*a+63]',), {}), None), (((u'[0,-a+1,0,6*a+7,7*a+7]',), {}), None), (((u'[1,-a+1,0,-a-1,0]',), {}), None), (((u'[1,-a+1,0,-66*a-56,-333*a-241]',), {}), None), (((u'[a+1,-a+1,1,-17*a-32,72*a-7]',), {}), None), (((u'[1,-a+1,0,-71*a-46,-328*a-269]',), {}), None), (((u'[a+1,-a+1,1,3*a+3,8*a+3]',), {}), None), (((u'[a+1,0,0,-7*a+5,8*a-9]',), {}), None), (((u'[a+1,0,0,148*a-200,880*a-1384]',), {}), None), (((u'[a+1,a,0,-4,-7*a-2]',), {}), 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'[a+1,a,0,1,0]',), {}), 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,0]',), {}), None), (((u'[a,a,0,8*a-15,11*a+17]',), {}), None), (((u'[a,a,0,3*a-195,847*a-500]',), {}), None), (((u'[a+1,-a,1,4*a-5,2*a-10]',), {}), None), (((u'[a,-a+1,1,1932*a-3156,50492*a-81759]',), {}), None), (((u'[a,-a+1,1,-6,-3*a]',), {}), None), (((u'[a+1,-a,1,-256*a-156,-2345*a-1447]',), {}), None), (((u'[a+1,-a,a+1,3*a-10,9*a-19]',), {}), None), (((u'[a,-a+1,1,-1,a-2]',), {}), None), (((u'[a,-a+1,1,45*a-81,-171*a+276]',), {}), None), (((u'[a+1,-a,a+1,-2*a,-a]',), {}), None), (((u'[1,a+1,0,-4*a-4,-13*a-9]',), {}), None), (((u'[1,a+1,0,-89*a-64,-620*a-375]',), {}), None), (((u'[1,a+1,0,a+1,0]',), {}), None), (((u'[a+1,-a+1,0,74*a-123,422*a-683]',), {}), None), (((u'[0,-a+1,0,-107*a-85,-701*a-410]',), {}), None), (((u'[0,-a+1,0,-7*a-5,-9*a-6]',), {}), None), (((u'[0,a,0,44*a-68,148*a-244]',), {}), None), (((u'[0,-a+1,0,-2*a,a+1]',), {}), None), (((u'[0,-a-1,a,2*a-10,-13]',), {}), None), (((u'[0,-a-1,a,2*a,-a]',), {}), None), (((u'[0,a,a+1,2*a+2,-a]',), {}), None), (((u'[0,a-1,a+1,148*a-246,-1151*a+1868]',), {}), None), (((u'[1,0,a+1,4*a-5,7*a-15]',), {}), None), (((u'[1,0,a+1,-a,-a]',), {}), None), (((u'[1,1,a+1,-3413*a-2108,-116137*a-71773]',), {}), None), (((u'[a+1,a-1,1,-26*a-27,-133*a-62]',), {}), None), (((u'[a+1,a-1,1,-a-2,-3*a-2]',), {}), None), (((u'[a,-a,0,556*a-910,-7512*a+12147]',), {}), None), (((u'[1,-a+1,0,61492*a-99493,-8737966*a+14138322]',), {}), None), (((u'[a+1,a-1,1,-76*a-227,447*a+1348]',), {}), 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'[a,-1,a,-2*a-50,-121*a+57]',), {}), None), (((u'[a,-1,a,3*a-10,4*a-8]',), {}), None), (((u'[a,-1,a,88*a-130,401*a-661]',), {}), None), (((u'[a,-1,a,-2*a,a-1]',), {}), None), (((u'[1,a-1,0,a-2,0]',), {}), None), (((u'[1,a-1,0,-4*a+8,7*a-10]',), {}), None), (((u'[a+1,-1,0,3*a-20,4*a-31]',), {}), None), (((u'[a+1,-1,0,-2*a,a]',), {}), None), (((u'[a,0,1,-a,0]',), {}), None), (((u'[a,0,1,4*a,2*a-6]',), {}), None), (((u'[a,0,1,104*a-310,946*a-2218]',), {}), None), (((u'[a,a-1,a+1,-7*a-5,6*a+4]',), {}), None), (((u'[1,-a,0,39*a-63,-126*a+207]',), {}), None), (((u'[a,a-1,a+1,-2*a,-a-1]',), {}), None), (((u'[a,a-1,a+1,-97*a-80,588*a+337]',), {}), None), (((u'[a+1,-a+1,1,a-3,1]',), {}), None), (((u'[a+1,-a+1,1,-9*a-8,12*a+11]',), {}), None), (((u'[0,1,0,2*a-2,-a+2]',), {}), None), (((u'[0,1,0,27*a-47,-79*a+127]',), {}), None), (((u'[0,a+1,0,-566*a-373,-8475*a-5269]',), {}), None), (((u'[0,1,0,-18*a+18,11*a-38]',), {}), None), (((u'[a+1,-a-1,1,-93*a-93,-525*a-394]',), {}), None), (((u'[1,1,a,7*a-15,15*a-32]',), {}), None), (((u'[1,1,a,-3*a,-2*a-3]',), {}), 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,-5*a-4,-8*a-3]',), {}), None), (((u'[a,1,1,6*a-9,-6*a+10]',), {}), 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,-2268*a-2268,10875*a+8156]',), {}), None), (((u'[a+1,-a-1,1,-68*a-68,1275*a+956]',), {}), None), (((u'[a+1,-a-1,1,-26668*a-26668,3007355*a+2255516]',), {}), None), (((u'[a+1,-a+1,0,-5*a-14,12*a+11]',), {}), None), (((u'[a+1,-a+1,0,5*a-74,34*a-287]',), {}), None), (((u'[a+1,-a+1,0,-4,a-6]',), {}), None), (((u'[a+1,-a+1,a,3*a-7,6]',), {}), None), (((u'[a+1,-a+1,0,1,0]',), {}), None), (((u'[a+1,-a+1,a,-2*a-2,a-1]',), {}), None)] }}} {{{id=18| test=db('select weq,Lstar from CLI ORDER BY N,Lstar') /// }}} {{{id=28| test2=db('SELECT found from CLI ORDER BY N,Lstar,fh,nh') /// }}} {{{id=35| for r in range(10): print test[r][15] /// '[a+1,a-1,a+1,-35*a-56,-260*a-117]', '[a+1,a-1,a+1,-35*a-56,-260*a-117]', '[a+1,a-1,a+1,-35*a-56,-260*a-117]', '[a+1,a-1,a+1,-35*a-56,-260*a-117]', '[a+1,a-1,a+1,-35*a-56,-260*a-117]', '[a+1,a-1,a+1,-35*a-56,-260*a-117]', '[0,-a+1,0,-3*a-2,-2*a-2]', '[0,-a+1,0,-3*a-2,-2*a-2]', '[0,-a+1,0,-3*a-2,-2*a-2]', '[0,-a+1,0,-3*a-2,-2*a-2]', }}} {{{id=29| test=db('SELECT N,cond,weq,rlow,rhi,tor,sign,ordD,ordj,tama,kod,Lstar,fh,nh,mtrx,found from CLI ORDER BY N,Lstar,fh,nh') /// }}} {{{id=30| len(test) /// 232 }}} {{{id=31| db2=nosqlite.Client('/home/psharaba/table/ECdb').db /// }}} {{{id=32| db2.N1k.count() /// (740,) }}} {{{id=34| test[0][14] /// u"['[a+1,a-1,a+1,5*a-16,4*a-21]'," }}} {{{id=36| k=0 g=0 /// }}} {{{id=37| print k,g /// 96 136 }}} {{{id=33| for t in range(6): c=test[t] G=list(db2.N1k.find(ainv1=c[2])) H=list(db2.N1k.find(ainv2=c[2])) L=len(G+H) if L != 0: print test[t][0],test[t][2],test[t][11], G,H g+=1 else: print test[t][0],test[t][2],test[t][11], 'unique CL' k+=1 /// 124 [a+1,a-1,a+1,-1,-a+1] 0.744306279976 [{u'nh': 44.627447581702725, u'rlow': 0, u'tama': u'2,1', u'cond1': u'10*a-4', u'ainv2': u'[a,a+1,1,1,0]', u'ainv1': u'[a+1,a-1,a+1,-1,-a+1]', u'cond2': u'-10*a+6', u'tor': 6, u'Lstar': 0.74430627997599996, u'ordD': u'2,1', u'kod': u'I2,I1', u'sign': u'+,-', u'rhi': 0, u'eta': 252, u'label': u'124a2a', u'fh': -1.0116113737408514, u'found': u'J', u'N': 124, u'mtrx': u'matrix(2,[0,2,2,0])', u'ordj': u'2,1'}] [] 124 [a,-a,a,85*a-139,-436*a+705] 0.744306279976 [{u'nh': 57.872648363088942, u'rlow': 0, u'tama': u'1,2', u'cond1': u'10*a-4', u'ainv2': u'[a+1,-1,a+1,-87*a-54,435*a+269]', u'ainv1': u'[a,-a,a,85*a-139,-436*a+705]', u'cond2': u'-10*a+6', u'tor': 6, u'Lstar': 0.74430627997599996, u'ordD': u'1,2', u'kod': u'I1,I2', u'sign': u'+,+', u'rhi': 0, u'eta': 252, u'label': u'124a2b', u'fh': -0.66503778346087872, u'found': u'JI', u'N': 124, u'mtrx': u'matrix(2,[0,2,2,0])', u'ordj': u'1,2'}] [] 124 [a+1,a-1,a+1,5*a-16,4*a-21] 0.744306279976 unique CL 124 [a+1,a-1,a+1,-35*a-56,-260*a-117] 0.744306279976 [{u'nh': 69.194685401364225, u'rlow': 0, u'tama': u'3,6', u'cond1': u'10*a-4', u'ainv2': u'[a,a+1,1,35*a-89,204*a-342]', u'ainv1': u'[a+1,a-1,a+1,-35*a-56,-260*a-117]', u'cond2': u'-10*a+6', u'tor': 6, u'Lstar': 0.74430627997599996, u'ordD': u'3,6', u'kod': u'I3,I6', u'sign': u'+,+', u'rhi': 0, u'eta': 252, u'label': u'124a2c', u'fh': -0.11573163912682383, u'found': u'JI', u'N': 124, u'mtrx': u'matrix(2,[0,2,2,0])', u'ordj': u'3,6'}] [] 124 [a,-a,a,4880*a-7979,199505*a-323066] 0.744306279976 [{u'nh': 85.524064509145447, u'rlow': 0, u'tama': u'2,1', u'cond1': u'10*a-4', u'ainv2': u'[a+1,-1,a+1,-4882*a-3099,-199506*a-123561]', u'ainv1': u'[a,-a,a,4880*a-7979,199505*a-323066]', u'cond2': u'-10*a+6', u'tor': 2, u'Lstar': 0.74430627997599996, u'ordD': u'2,1', u'kod': u'I2,I1', u'sign': u'+,-', u'rhi': 0, u'eta': 252, u'label': u'124a2d', u'fh': 0.087000914927258283, u'found': u'JI', u'N': 124, u'mtrx': u'matrix(2,[0,2,2,0])', u'ordj': u'2,1'}] [] 124 [a,a+1,1,-188901*a-116781,-47466544*a-29336020] 0.744306279976 [{u'nh': 93.841852902287656, u'rlow': 0, u'tama': u'1,2', u'cond1': u'10*a-4', u'ainv2': u'[a+1,a-1,a+1,188901*a-305684,47160860*a-76307979]', u'ainv1': u'[a,a+1,1,-188901*a-116781,-47466544*a-29336020]', u'cond2': u'-10*a+6', u'tor': 2, u'Lstar': 0.74430627997599996, u'ordD': u'1,2', u'kod': u'I1,I2', u'sign': u'+,+', u'rhi': 0, u'eta': 252, u'label': u'124a2e', u'fh': 0.43357450520723095, u'found': u'JI', u'N': 124, u'mtrx': u'matrix(2,[0,2,2,0])', u'ordj': u'1,2'}] [] }}} {{{id=24| %time for t in range(len(test)): c=test[t] E = EllipticCurve(K,eval(c[2])) F = E.global_minimal_model() F = canonical_model(F).a_invariants() G1 = str([F[0],F[1],F[2],F[3],F[4]]).replace(' ','') H = list(db2.N1k.find(ainv1=G1))+list(db2.N1k.find(ainv2=G1)) L = len(H) if L==0: E = EllipticCurve(K, [galois_conjugate(K(alpha)) for alpha in F]) F = E.global_minimal_model() F = canonical_model(F).a_invariants() G2 = str([F[0],F[1],F[2],F[3],F[4]]).replace(' ','') J1 = K.ideal(eval(c[1])) J2 = K.ideal(galois_conjugate(K(eval(c[1])))) if J1.integral_basis()[1][0] > J2.integral_basis()[1][0]: J1,J2 = J2,J1 G1,G2 = G2,G1 h=f_height(eval(G1)) n=n_height(eval(G1)) db.N1k.insert({'N': int(J1.norm()), 'eta': int(eta_ideal(J1)), 'cond1': str(J1.gen(0)).replace(' ',''), 'cond2': str(J2.gen(0)).replace(' ',''), 'ainv1': G1, 'ainv2': G2, 'rlow': int(c[3]), 'rhi': int(c[4]), 'tor': int(c[6]), 'sign': c[5], 'ordD': c[7], 'ordj': c[8], 'tama': c[9], 'kod': c[10],'fh': float(h), 'nh': float(n),'mtrx': c[11], 'found': 'JI'}) /// }}} {{{id=26| /// }}}