{{{id=3| K.=NumberField(x^2-x-1) import nosqlite db=nosqlite.Client('/home/psharaba/ECD/ECdb').db from psage.ellcurve.minmodel.sqrt5 import canonical_model /// }}} {{{id=4| def eta(alpha): return (alpha*alpha).trace() def find_small(alpha): while True: plus=a*alpha minus=alpha/a if eta(plus) 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'}) /// CPU time: 820.05 s, Wall time: 834.72 s }}} {{{id=12| g=open('/home/psharaba/ECD/list3.txt') /// }}} {{{id=11| %time for s in g.readlines(): b=s.split() h=f_height(eval(b[2])) n=n_height(eval(b[2])) db.N1k.update({'fh':float(h), 'nh': float(n), 'mtrx': b[3]}, ainv1=b[2]) /// Traceback (most recent call last): File "", line 1, in File "_sage_input_9.py", line 12, in exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8 -*-\\n" + _support_.preparse_worksheet_cell(base64.b64decode("Zm9yIHMgaW4gZy5yZWFkbGluZXMoKToKICAgIGI9cy5zcGxpdCgpCiAgICBoPWZfaGVpZ2h0KGV2YWwoYlsyXSkpCiAgICBuPW5faGVpZ2h0KGV2YWwoYlsyXSkpCiAgICBkYi5OMWsudXBkYXRlKHsnZmgnOmZsb2F0KGgpLCAnbmgnOiBmbG9hdChuKSwgJ210cngnOiBiWzNdfSwgYWludjE9YlsyXSk="),globals())+"\\n"); execfile(os.path.abspath("___code___.py")); print "CPU time: %.2f s, Wall time: %.2f s"%(cputime(__SAGE_t__), walltime(__SAGE_w__))' + '\n', '', 'single') File "", line 1, in File "/tmp/tmpxrZx_j/___code___.py", line 3, in exec compile(u"for s in g.readlines():\n b=s.split()\n h=f_height(eval(b[_sage_const_2 ]))\n n=n_height(eval(b[_sage_const_2 ]))\n db.N1k.update({'fh':float(h), 'nh': float(n), 'mtrx': b[_sage_const_3 ]}, ainv1=b[_sage_const_2 ])" + '\n', '', 'single') File "", line 3, in File "/tmp/tmpps2kWP/___code___.py", line 10, in f_height A2=E.period_lattice(phi2).basis() File "/usr/local/sage/local/lib/python2.6/site-packages/sage/schemes/elliptic_curves/ell_number_field.py", line 1835, in period_lattice return PeriodLattice_ell(self,embedding) File "/usr/local/sage/local/lib/python2.6/site-packages/sage/schemes/elliptic_curves/period_lattice.py", line 238, in __init__ e3 = AA(e3) File "parent.pyx", line 915, in sage.structure.parent.Parent.__call__ (sage/structure/parent.c:6820) File "coerce_maps.pyx", line 82, in sage.structure.coerce_maps.DefaultConvertMap_unique._call_ (sage/structure/coerce_maps.c:3254) File "coerce_maps.pyx", line 77, in sage.structure.coerce_maps.DefaultConvertMap_unique._call_ (sage/structure/coerce_maps.c:3157) File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 598, in _element_constructor_ if x.imag().is_zero(): File "element.pyx", line 776, in sage.structure.element.Element.is_zero (sage/structure/element.c:5916) File "element.pyx", line 764, in sage.structure.element.Element.__nonzero__ (sage/structure/element.c:5872) File "element.pyx", line 886, in sage.structure.element.Element.__richcmp__ (sage/structure/element.c:7306) File "element.pyx", line 858, in sage.structure.element.Element._richcmp (sage/structure/element.c:7183) File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 3413, in __cmp__ return self.sign() File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 3633, in sign return self.sign() File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 3636, in sign self.exactify() File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 2617, in exactify self._set_descr(self._descr.exactify()) File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 5500, in exactify arg.exactify() File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 2617, in exactify self._set_descr(self._descr.exactify()) File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 4996, in exactify root = ANRoot(QQx(red_pol), new_intv) File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 4417, in __init__ self._interval = self.refine_interval(interval, 64) File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 4528, in refine_interval return self._real_refine_interval(interval, prec) File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 4694, in _real_refine_interval uq = (center + interval.upper()) / 2 KeyboardInterrupt __SAGE__ }}} {{{id=14| temp=db("select N,eta,ainv1,ainv2 from N1k WHERE found='CL' ORDER BY N,eta,fh,nh") /// }}} {{{id=15| %time for v in range(len(temp)): h=f_height(eval(temp[v][2])) n=n_height(eval(temp[v][2])) db.N1k.update({'fh':float(h), 'nh':float(n)}, ainv1=temp[v][2]) /// Traceback (most recent call last): File "", line 1, in File "_sage_input_11.py", line 12, in exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8 -*-\\n" + _support_.preparse_worksheet_cell(base64.b64decode("Zm9yIHYgaW4gcmFuZ2UobGVuKHRlbXApKToKICAgIGg9Zl9oZWlnaHQoZXZhbCh0ZW1wW3ZdWzJdKSkKICAgIG49bl9oZWlnaHQoZXZhbCh0ZW1wW3ZdWzJdKSkKICAgIGRiLk4xay51cGRhdGUoeydmaCc6ZmxvYXQoaCksICduaCc6ZmxvYXQobil9LCBhaW52MT10ZW1wW3ZdWzJdKQ=="),globals())+"\\n"); execfile(os.path.abspath("___code___.py")); print "CPU time: %.2f s, Wall time: %.2f s"%(cputime(__SAGE_t__), walltime(__SAGE_w__))' + '\n', '', 'single') File "", line 1, in File "/tmp/tmpLFY7Av/___code___.py", line 3, in exec compile(u"for v in range(len(temp)):\n h=f_height(eval(temp[v][_sage_const_2 ]))\n n=n_height(eval(temp[v][_sage_const_2 ]))\n db.N1k.update({'fh':float(h), 'nh':float(n)}, ainv1=temp[v][_sage_const_2 ])" + '\n', '', 'single') File "", line 2, in File "/tmp/tmpps2kWP/___code___.py", line 6, in f_height A1=E.period_lattice(phi1).basis() File "/usr/local/sage/local/lib/python2.6/site-packages/sage/schemes/elliptic_curves/ell_number_field.py", line 1835, in period_lattice return PeriodLattice_ell(self,embedding) File "/usr/local/sage/local/lib/python2.6/site-packages/sage/schemes/elliptic_curves/period_lattice.py", line 238, in __init__ e3 = AA(e3) File "parent.pyx", line 915, in sage.structure.parent.Parent.__call__ (sage/structure/parent.c:6820) File "coerce_maps.pyx", line 82, in sage.structure.coerce_maps.DefaultConvertMap_unique._call_ (sage/structure/coerce_maps.c:3254) File "coerce_maps.pyx", line 77, in sage.structure.coerce_maps.DefaultConvertMap_unique._call_ (sage/structure/coerce_maps.c:3157) File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 598, in _element_constructor_ if x.imag().is_zero(): File "element.pyx", line 776, in sage.structure.element.Element.is_zero (sage/structure/element.c:5916) File "element.pyx", line 764, in sage.structure.element.Element.__nonzero__ (sage/structure/element.c:5872) File "element.pyx", line 886, in sage.structure.element.Element.__richcmp__ (sage/structure/element.c:7306) File "element.pyx", line 858, in sage.structure.element.Element._richcmp (sage/structure/element.c:7183) File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 3413, in __cmp__ return self.sign() File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 3633, in sign return self.sign() File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 3636, in sign self.exactify() File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 2617, in exactify self._set_descr(self._descr.exactify()) File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 5501, in exactify iv = QQbar_I * (arg.conjugate() - arg) / 2 File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 3228, in conjugate return AlgebraicNumber(self._descr.conjugate(self)) File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 2865, in __init__ AlgebraicNumber_base.__init__(self, QQbar, x) File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 2188, in __init__ self._value = self._descr._interval_fast(64) File "/usr/local/sage/local/lib/python2.6/site-packages/sage/rings/qqbar.py", line 5261, in _interval_fast v = self._value.polynomial()(gen_val) File "polynomial_rational_flint.pyx", line 456, in sage.rings.polynomial.polynomial_rational_flint.Polynomial_rational_flint.__call__ (sage/rings/polynomial/polynomial_rational_flint.cpp:6833) File "polynomial_element.pyx", line 626, in sage.rings.polynomial.polynomial_element.Polynomial.__call__ (sage/rings/polynomial/polynomial_element.c:7047) File "element.pyx", line 1459, in sage.structure.element.RingElement.__mul__ (sage/structure/element.c:12096) File "coerce.pyx", line 708, in sage.structure.coerce.CoercionModel_cache_maps.bin_op (sage/structure/coerce.c:6349) File "coerce.pyx", line 1179, in sage.structure.coerce.CoercionModel_cache_maps.get_action (sage/structure/coerce.c:11388) File "coerce.pyx", line 1302, in sage.structure.coerce.CoercionModel_cache_maps.discover_action (sage/structure/coerce.c:12568) File "parent.pyx", line 2109, in sage.structure.parent.Parent.get_action (sage/structure/parent.c:14536) File "parent.pyx", line 2183, in sage.structure.parent.Parent.discover_action (sage/structure/parent.c:15537) File "coerce_actions.pyx", line 186, in sage.structure.coerce_actions.detect_element_action (sage/structure/coerce_actions.c:4348) File "coerce_actions.pyx", line 28, in sage.structure.coerce_actions._record_exception (sage/structure/coerce_actions.c:2809) File "coerce.pyx", line 310, in sage.structure.coerce.CoercionModel_cache_maps._record_exception (sage/structure/coerce.c:3843) File "/usr/local/sage/local/lib/python/traceback.py", line 236, in format_exc return ''.join(format_exception(etype, value, tb, limit)) File "/usr/local/sage/local/lib/python/traceback.py", line 142, in format_exception list = list + format_tb(tb, limit) File "/usr/local/sage/local/lib/python/traceback.py", line 76, in format_tb return format_list(extract_tb(tb, limit)) File "/usr/local/sage/local/lib/python/traceback.py", line 101, in extract_tb line = linecache.getline(filename, lineno, f.f_globals) File "/usr/local/sage/local/lib/python/linecache.py", line 14, in getline lines = getlines(filename, module_globals) File "/usr/local/sage/local/lib/python/linecache.py", line 40, in getlines return updatecache(filename, module_globals) File "/usr/local/sage/local/lib/python/linecache.py", line 115, in updatecache fullname = os.path.join(dirname, basename) File "/usr/local/sage/local/lib/python/posixpath.py", line 59, in join def join(a, *p): KeyboardInterrupt __SAGE__ }}} {{{id=16| /// }}}