{{{id=1| K. = NumberField(x^2-x-1) import nosqlite db=nosqlite.Client('/home/psharaba/ECdb').db /// }}} {{{id=3| def saturate_rank_01(E, sat_bound): v = E.simon_two_descent() t2 = E.two_torsion_rank() sel2 = v[1] rank_bound = sel2 - t2 if rank_bound >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, reg=1, Lstar=1): M = E.tamagawa_product_bsd() print M sha = RR(sqrt(5))*Lstar*(E.torsion_order())^2/((real_periods(E))*reg*M) return sha /// }}} {{{id=4| def printer(): temp = db('select N,weq1,weq2,L,omega from N1k ORDER BY N,eta,R,V,U') J=open('/home/psharaba/RegList.txt','w') for y in temp: S=[] for entry in y: s=str(entry) s.replace(' ','') S.append(s) J.write(' '.join(S)+'\n') /// }}} {{{id=5| printer() /// }}} {{{id=6| /// }}}