#14658: multivariate factorization errors
-----------------------------+----------------------------------------------
Reporter: burcin | Owner: tbd
Type: defect | Status: new
Priority: major | Milestone: sage-5.11
Component: factorization | Keywords: singular factory
Work issues: | Report Upstream: N/A
Reviewers: | Authors:
Merged in: | Dependencies:
Stopgaps: |
-----------------------------+----------------------------------------------
Reported by Manuel Kauers.
On Sage 5.10.beta4:
{{{
sage: R.<x1,x2,x3> = QQ['x1,x2,x3']
sage: f = (x1+x2+x3)*(2*x1-x2+x3+2)*(4*x1+x2+x3+2)*(8*x1-x2+x3+4)
sage: f.factor()
---------------------------------------------------------------------------
RuntimeError Traceback (most recent call
last)
<ipython-input-3-429743410e57> in <module>()
----> 1 f.factor()
/home/burcin/sage/sage-5.10.beta2/local/lib/python2.7/site-
packages/sage/rings/polynomial/multi_polynomial_libsingular.so in
sage.rings.polynomial.multi_polynomial_libsingular.MPolynomial_libsingular.factor
(build/cythonized/sage/rings/polynomial/multi_polynomial_libsingular.cpp:25997)()
RuntimeError: Floating point exception
sage: g = (x1+x2+x3)*(2*x1-x2+x3+2)*(4*x1+x2+x3+2)
sage: g.factor()
---------------------------------------------------------------------------
RuntimeError Traceback (most recent call
last)
<ipython-input-5-689da5069cee> in <module>()
----> 1 g.factor()
/home/burcin/sage/sage-5.10.beta2/local/lib/python2.7/site-
packages/sage/rings/polynomial/multi_polynomial_libsingular.so in
sage.rings.polynomial.multi_polynomial_libsingular.MPolynomial_libsingular.factor
(build/cythonized/sage/rings/polynomial/multi_polynomial_libsingular.cpp:25997)()
RuntimeError: Floating point exception
}}}
These examples work with Singular 3.1.6 (should be updated in #14333):
{{{
SINGULAR /
A Computer Algebra System for Polynomial Computations / version
3-1-6
0<
by: W. Decker, G.-M. Greuel, G. Pfister, H. Schoenemann \ Dec 2012
FB Mathematik der Universitaet, D-67653 Kaiserslautern \
> ring R = 0,(x1,x2,x3),dp;
> poly f = (x1+x2+x3)*(2*x1-x2+x3+2)*(4*x1+x2+x3+2)*(8*x1-x2+x3+4);
> factorize(f);
[1]:
_[1]=1
_[2]=2*x1-x2+x3+2
_[3]=x1+x2+x3
_[4]=8*x1-x2+x3+4
_[5]=4*x1+x2+x3+2
[2]:
1,1,1,1,1
> poly g = (x1+x2+x3)*(2*x1-x2+x3+2)*(4*x1+x2+x3+2);
> factorize(g);
[1]:
_[1]=1
_[2]=4*x1+x2+x3+2
_[3]=2*x1-x2+x3+2
_[4]=x1+x2+x3
[2]:
1,1,1,1
}}}
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/14658>
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 unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/sage-trac?hl=en.
For more options, visit https://groups.google.com/groups/opt_out.
