#5698: Sage 3.4.1.rc1: doctest failure in
schemes/elliptic_curves/monsky_washnitzer.py
----------------------+-----------------------------------------------------
Reporter: mabshoff | Owner: mabshoff
Type: defect | Status: new
Priority: blocker | Milestone: sage-3.4.1
Component: doctest | Keywords:
----------------------+-----------------------------------------------------
The following happens all the sudden when updating FLINT:
{{{
sage -t -long
"devel/sage/sage/schemes/elliptic_curves/monsky_washnitzer.py"
**********************************************************************
File
"/scratch/mabshoff/sage-3.4.1.rc2/devel/sage/sage/schemes/elliptic_curves/monsky_washnitzer.py",
line 1487:
sage: B = A.change_ring(Integers(p**prec)); B # long
time
Expected:
[74311982 57996908]
[95877067 25828133]
Got:
[54572546 87269244]
[10839343 79974813]
**********************************************************************
File
"/scratch/mabshoff/sage-3.4.1.rc2/devel/sage/sage/schemes/elliptic_curves/monsky_washnitzer.py",
line 1490:
sage: B.det() # long
time
Expected:
10007
Got:
29700776
**********************************************************************
File
"/scratch/mabshoff/sage-3.4.1.rc2/devel/sage/sage/schemes/elliptic_curves/monsky_washnitzer.py",
line 1492:
sage: B.trace() # long
time
Expected:
66
Got:
34407310
**********************************************************************
File
"/scratch/mabshoff/sage-3.4.1.rc2/devel/sage/sage/schemes/elliptic_curves/monsky_washnitzer.py",
line 1562:
sage: B.det() # long
time
Expected:
11 + 484*t^2 + 451*t^3 + O(t^4)
Got:
245 + 240*t + 724*t^2 + 808*t^3 + O(t^4)
**********************************************************************
}}}
But it is unclear if FLINT is the cause here. Downgrading FLINT makes it
disappear, but the cause could be and likely is in some other patch. Note
also the following issue valgrind found:
{{{
==30646== Conditional jump or move depends on uninitialised value(s)
==30646== at 0x8F302EC: gdiv (in
/scratch/mabshoff/sage-3.4.1.rc0/local/lib/libpari-gmp.so.2)
==30646== by 0x828FAFE: QQ_to_t_FRAC (in
/scratch/mabshoff/sage-3.4.1.rc0/devel/sage-main/c_lib/libcsage.so)
==30646== by 0x22EB14E7:
__pyx_f_4sage_6matrix_21matrix_rational_dense_pari_GEN
(matrix_rational_dense.c:25089)
==30646== by 0x22EE0D3C:
__pyx_pf_4sage_6matrix_21matrix_rational_dense_21Matrix_rational_dense__invert_pari
(matrix_rational_dense.c:24157)
==30646== by 0x417E92: PyObject_Call (abstract.c:1861)
==30646== by 0x22ED954D:
__pyx_pf_4sage_6matrix_21matrix_rational_dense_21Matrix_rational_dense___invert__main
(matrix_rational_dense.c:9542)
==30646== by 0x417E92: PyObject_Call (abstract.c:1861)
==30646== by 0x22EAF80A:
__pyx_pf_4sage_6matrix_21matrix_rational_dense_21Matrix_rational_dense___invert__
(matrix_rational_dense.c:8933)
==30646== by 0x216C1751:
__pyx_pf_4sage_6matrix_7matrix0_6Matrix___invert__ (matrix0.c:22763)
==30646== by 0xEDAA1DA:
__pyx_f_4sage_9structure_7element_generic_power_c (element.c:23281)
==30646== by 0xEDAC821:
__pyx_pf_4sage_9structure_7element_11RingElement___pow__ (element.c:10954)
==30646== by 0x457FBC: wrap_ternaryfunc (typeobject.c:3621)
==30646== by 0x417E92: PyObject_Call (abstract.c:1861)
==30646== by 0x4816E1: PyEval_CallObjectWithKeywords (ceval.c:3442)
==30646== by 0x4D0F9D: wrapperdescr_call (descrobject.c:304)
==30646== by 0x417E92: PyObject_Call (abstract.c:1861)
==30646== by 0x216C13C0:
__pyx_pf_4sage_6matrix_7matrix0_6Matrix___pow__ (matrix0.c:23157)
==30646== by 0x418AC4: ternary_op (abstract.c:518)
==30646== by 0x4850BD: PyEval_EvalFrameEx (ceval.c:1063)
==30646== by 0x489A95: PyEval_EvalCodeEx (ceval.c:2836)
==30646== by 0x4D40D7: function_call (funcobject.c:517)
==30646== by 0x417E92: PyObject_Call (abstract.c:1861)
==30646== by 0x41E6CE: instancemethod_call (classobject.c:2519)
==30646== by 0x417E92: PyObject_Call (abstract.c:1861)
==30646== by 0x45A257: slot_tp_init (typeobject.c:4943)
==30646==
==30646== Conditional jump or move depends on uninitialised value(s)
==30646== at 0x8F3004B: gdiv (in
/scratch/mabshoff/sage-3.4.1.rc0/local/lib/libpari-gmp.so.2)
==30646== by 0x8E67C57: gauss_intern (in
/scratch/mabshoff/sage-3.4.1.rc0/local/lib/libpari-gmp.so.2)
==30646== by 0x8E67F45: gauss (in
/scratch/mabshoff/sage-3.4.1.rc0/local/lib/libpari-gmp.so.2)
==30646== by 0x22EE0D9C:
__pyx_pf_4sage_6matrix_21matrix_rational_dense_21Matrix_rational_dense__invert_pari
(matrix_rational_dense.c:24218)
==30646== by 0x417E92: PyObject_Call (abstract.c:1861)
==30646== by 0x22ED954D:
__pyx_pf_4sage_6matrix_21matrix_rational_dense_21Matrix_rational_dense___invert__main
(matrix_rational_dense.c:9542)
==30646== by 0x417E92: PyObject_Call (abstract.c:1861)
==30646== by 0x22EAF80A:
__pyx_pf_4sage_6matrix_21matrix_rational_dense_21Matrix_rational_dense___invert__
(matrix_rational_dense.c:8933)
==30646== by 0x216C1751:
__pyx_pf_4sage_6matrix_7matrix0_6Matrix___invert__ (matrix0.c:22763)
==30646== by 0xEDAA1DA:
__pyx_f_4sage_9structure_7element_generic_power_c (element.c:23281)
==30646== by 0xEDAC821:
__pyx_pf_4sage_9structure_7element_11RingElement___pow__ (element.c:10954)
==30646== by 0x457FBC: wrap_ternaryfunc (typeobject.c:3621)
==30646== by 0x417E92: PyObject_Call (abstract.c:1861)
==30646== by 0x4816E1: PyEval_CallObjectWithKeywords (ceval.c:3442)
==30646== by 0x4D0F9D: wrapperdescr_call (descrobject.c:304)
==30646== by 0x417E92: PyObject_Call (abstract.c:1861)
==30646== by 0x216C13C0:
__pyx_pf_4sage_6matrix_7matrix0_6Matrix___pow__ (matrix0.c:23157)
==30646== by 0x418AC4: ternary_op (abstract.c:518)
==30646== by 0x4850BD: PyEval_EvalFrameEx (ceval.c:1063)
==30646== by 0x489A95: PyEval_EvalCodeEx (ceval.c:2836)
==30646== by 0x4D40D7: function_call (funcobject.c:517)
==30646== by 0x417E92: PyObject_Call (abstract.c:1861)
==30646== by 0x41E6CE: instancemethod_call (classobject.c:2519)
==30646== by 0x417E92: PyObject_Call (abstract.c:1861)
==30646== by 0x45A257: slot_tp_init (typeobject.c:4943)
}}}
This is likely caused by the Brandt module patch when William optimized
the "small LA" cases.
Cheers,
Michael
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/5698>
Sage <http://sagemath.org/>
Sage - Open Source Mathematical Software: Building the Car Instead of
Reinventing the Wheel
--~--~---------~--~----~------------~-------~--~----~
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
-~----------~----~----~----~------~----~------~--~---
