#12718: polynomial substitution overflow hell (libsingular bug?)
---------------------------------------+------------------------------------
Reporter: was | Owner: malb
Type: defect | Status: needs_work
Priority: critical | Milestone: sage-5.0
Component: commutative algebra | Resolution:
Keywords: | Work issues:
Report Upstream: N/A | Reviewers: Volker Braun, William
Stein
Authors: Martin Albrecht | Merged in:
Dependencies: | Stopgaps:
---------------------------------------+------------------------------------
Changes (by jdemeyer):
* status: positive_review => needs_work
Comment:
On moufang (OSX 10.4 ppc32):
{{{
**********************************************************************
File
"/Users/buildbot/build/sage/moufang-1/moufang_full/build/sage-5.0.beta15/devel
/sage-main/sage/rings/polynomial/multi_polynomial_libsingular.pyx", line
3154:
sage: n=1000; f = x^n; f.subs(x = x^n)
Exception raised:
Traceback (most recent call last):
File
"/Users/buildbot/build/sage/moufang-1/moufang_full/build/sage-5.0.beta15/local/bin/ncadoctest.py",
line 1231, in run_one_test
self.run_one_example(test, example, filename, compileflags)
File
"/Users/buildbot/build/sage/moufang-1/moufang_full/build/sage-5.0.beta15/local/bin/sagedoctest.py",
line 38, in run_one_example
OrigDocTestRunner.run_one_example(self, test, example, filename,
compileflags)
File
"/Users/buildbot/build/sage/moufang-1/moufang_full/build/sage-5.0.beta15/local/bin/ncadoctest.py",
line 1172, in run_one_example
compileflags, 1) in test.globs
File "<doctest __main__.example_61[28]>", line 1, in <module>
n=Integer(1000); f = x**n; f.subs(x = x**n)###line 3154:
sage: n=1000; f = x^n; f.subs(x = x^n)
File "multi_polynomial_libsingular.pyx", line 3240, in
sage.rings.polynomial.multi_polynomial_libsingular.MPolynomial_libsingular.subs
(sage/rings/polynomial/multi_polynomial_libsingular.cpp:21411)
raise OverflowError("Exponent overflow (%d)."%(degree))
OverflowError: Exponent overflow (1000000).
**********************************************************************
File
"/Users/buildbot/build/sage/moufang-1/moufang_full/build/sage-5.0.beta15/devel
/sage-main/sage/rings/polynomial/multi_polynomial_libsingular.pyx", line
3157:
sage: n=100000; f = x^n; f.subs(x = x^n)
Expected:
Traceback (most recent call last):
...
OverflowError: Exponent overflow (10000000000).
Got:
Traceback (most recent call last):
File
"/Users/buildbot/build/sage/moufang-1/moufang_full/build/sage-5.0.beta15/local/bin/ncadoctest.py",
line 1231, in run_one_test
self.run_one_example(test, example, filename, compileflags)
File
"/Users/buildbot/build/sage/moufang-1/moufang_full/build/sage-5.0.beta15/local/bin/sagedoctest.py",
line 38, in run_one_example
OrigDocTestRunner.run_one_example(self, test, example, filename,
compileflags)
File
"/Users/buildbot/build/sage/moufang-1/moufang_full/build/sage-5.0.beta15/local/bin/ncadoctest.py",
line 1172, in run_one_example
compileflags, 1) in test.globs
File "<doctest __main__.example_61[29]>", line 1, in <module>
n=Integer(100000); f = x**n; f.subs(x = x**n)###line 3157:
sage: n=100000; f = x^n; f.subs(x = x^n)
File "multi_polynomial_libsingular.pyx", line 2347, in
sage.rings.polynomial.multi_polynomial_libsingular.MPolynomial_libsingular.__pow__
(sage/rings/polynomial/multi_polynomial_libsingular.cpp:16619)
singular_polynomial_pow(&_p, self._poly, exp, _ring)
File "polynomial.pyx", line 322, in
sage.libs.singular.polynomial.singular_polynomial_pow
(sage/libs/singular/polynomial.cpp:3976)
File "singular.pyx", line 667, in
sage.libs.singular.singular.overflow_check
(sage/libs/singular/singular.cpp:6526)
OverflowError: Exponent overflow (100000).
**********************************************************************
}}}
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/12718#comment:14>
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.
