#11975: Chow-Heegner points
-------------------------------+--------------------------------------------
Reporter: was | Owner: cremona
Type: enhancement | Status: needs_work
Priority: major | Milestone: sage-5.0
Component: elliptic curves | Keywords:
Work_issues: | Upstream: N/A
Reviewer: John Cremona | Author: William Stein
Merged: | Dependencies:
-------------------------------+--------------------------------------------
Comment(by jdemeyer):
There is numerical noise and a 32/64-bit issue on hawk (OpenSolaris
06.2009-32):
{{{
sage -t --long -force_lib
devel/sage/sage/schemes/elliptic_curves/chow_heegner_fast.pyx
**********************************************************************
File
"/export/home/buildbot/build/sage/hawk-1/hawk_full/build/sage-5.0.beta3/devel
/sage-main/sage/schemes/elliptic_curves/chow_heegner_fast.pyx", line 469:
sage: s.newton(-.6*I)
Expected:
[(-0.0720852069059 - 0.638326735148*I, 6, 1.3877787807814457e-17)]
Got:
[(-0.0720852069059 - 0.638326735148*I, 6, 0.0)]
**********************************************************************
File
"/export/home/buildbot/build/sage/hawk-1/hawk_full/build/sage-5.0.beta3/devel
/sage-main/sage/schemes/elliptic_curves/chow_heegner_fast.pyx", line 473:
sage: s.newton([0, -.6*I])
Expected:
[(-0.605829586188, 7, 0.0), (-0.0720852069059 - 0.638326735148*I, 6,
1.3877787807814457e-17)]
Got:
[(-0.605829586188, 7, 0.0), (-0.0720852069059 - 0.638326735148*I, 6,
0.0)]
**********************************************************************
sage -t --long -force_lib
devel/sage/sage/schemes/elliptic_curves/chow_heegner.py
**********************************************************************
File
"/export/home/buildbot/build/sage/hawk-1/hawk_full/build/sage-5.0.beta3/devel
/sage-main/sage/schemes/elliptic_curves/chow_heegner.py", line 86:
sage: zz
Expected:
0.341867712408 + 99999999.6769*I
Got:
0.34186771241 + 99999999.6769*I
**********************************************************************
File
"/export/home/buildbot/build/sage/hawk-1/hawk_full/build/sage-5.0.beta3/devel
/sage-main/sage/schemes/elliptic_curves/chow_heegner.py", line 207:
sage: sl2z_rep_in_fundom(CC(1+I), eps=1) # must be small
Expected:
Traceback (most recent call last):
...
ValueError: prec (=0) must be >= 2 and <= 9223372036854775807.
Got:
Traceback (most recent call last):
File
"/export/home/buildbot/build/sage/hawk-1/hawk_full/build/sage-5.0.beta3/local/bin/ncadoctest.py",
line 1231, in run_one_test
self.run_one_example(test, example, filename, compileflags)
File
"/export/home/buildbot/build/sage/hawk-1/hawk_full/build/sage-5.0.beta3/local/bin/sagedoctest.py",
line 38, in run_one_example
OrigDocTestRunner.run_one_example(self, test, example, filename,
compileflags)
File
"/export/home/buildbot/build/sage/hawk-1/hawk_full/build/sage-5.0.beta3/local/bin/ncadoctest.py",
line 1172, in run_one_example
compileflags, 1) in test.globs
File "<doctest __main__.example_2[22]>", line 1, in <module>
sl2z_rep_in_fundom(CC(Integer(1)+I), eps=Integer(1)) # must be
small###line 207:
sage: sl2z_rep_in_fundom(CC(1+I), eps=1) # must be small
File
"/export/home/buildbot/build/sage/hawk-1/hawk_full/build/sage-5.0.beta3/local/lib/python
/site-packages/sage/schemes/elliptic_curves/chow_heegner.py", line 265, in
sl2z_rep_in_fundom
w = ComplexIntervalField(prec)(z)
File
"/export/home/buildbot/build/sage/hawk-1/hawk_full/build/sage-5.0.beta3/local/lib/python
/site-packages/sage/rings/complex_interval_field.py", line 70, in
ComplexIntervalField
C = ComplexIntervalField_class(prec)
File
"/export/home/buildbot/build/sage/hawk-1/hawk_full/build/sage-5.0.beta3/local/lib/python
/site-packages/sage/rings/complex_interval_field.py", line 144, in
__init__
ParentWithGens.__init__(self, self._real_field(), ('I',), False,
category = Fields())
File
"/export/home/buildbot/build/sage/hawk-1/hawk_full/build/sage-5.0.beta3/local/lib/python
/site-packages/sage/rings/complex_interval_field.py", line 246, in
_real_field
self.__real_field = real_mpfi.RealIntervalField(self._prec)
File "real_mpfi.pyx", line 277, in
sage.rings.real_mpfi.RealIntervalField (sage/rings/real_mpfi.c:3097)
File "real_mpfi.pyx", line 443, in
sage.rings.real_mpfi.RealIntervalField_class.__init__
(sage/rings/real_mpfi.c:3302)
ValueError: prec (=0) must be >= 2 and <= 2147483647.
**********************************************************************
File
"/export/home/buildbot/build/sage/hawk-1/hawk_full/build/sage-5.0.beta3/devel
/sage-main/sage/schemes/elliptic_curves/chow_heegner.py", line 949:
sage: f(t[0][0])
Expected:
4.4408920985e-16 - 1.11022302463e-16*I
Got:
4.4408920985e-16 - 8.66548586359e-17*I
**********************************************************************
File
"/export/home/buildbot/build/sage/hawk-1/hawk_full/build/sage-5.0.beta3/devel
/sage-main/sage/schemes/elliptic_curves/chow_heegner.py", line 1524:
sage: phi(v[1])
Expected:
1.0
Got:
1.0 - 8.49472402142e-17*I
**********************************************************************
File
"/export/home/buildbot/build/sage/hawk-1/hawk_full/build/sage-5.0.beta3/devel
/sage-main/sage/schemes/elliptic_curves/chow_heegner.py", line 2011:
sage: P.numerical_approx()
Expected:
(1.44444444444440...: -1.03703703703... : 1.00000000000000)
Got:
(1.44444444444447 + 2.74248939530208e-15*I : -1.03703703703708 -
3.49918222885543e-15*I : 1.00000000000000)
**********************************************************************
}}}
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/11975#comment:21>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en.
