Yes, it is known. It no longer raises in Mateusz's new polys branch (see issue 1598):
In [3]: print factor( (1 - x**2 - y**2 + x**2 * y**2)**(S(1)/2)) (1 - x**2 - y**2 + x**2*y**2)**(1/2) Of course, it is still not factored there. Mateusz, is there some kind of deep option to factor? Also, you should use S(1)/2 (or simply sqrt). See http://docs.sympy.org/gotchas.html#python-numbers-vs-sympy-numbers Aaron Meurer On Feb 26, 2010, at 5:12 PM, Ben Goodrich wrote: > Is this a known issue? I see it on 0.7.0-git too. -- Ben > > Python 2.5.5 console for SymPy 0.6.6 > > These commands were executed: >>>> from __future__ import division >>>> from sympy import * >>>> x, y, z = symbols('xyz') >>>> k, m, n = symbols('kmn', integer=True) >>>> f, g, h = map(Function, 'fgh') > > Documentation can be found at http://sympy.org/ > > In [1]: factor( (1 - x**2 - y**2 + x**2 * y**2)**(1/2) ) > ERROR: An unexpected error occurred while tokenizing input > The following traceback may be corrupted or invalid > The error message is: ('EOF in multi-line statement', (33, 0)) > > ERROR: An unexpected error occurred while tokenizing input > The following traceback may be corrupted or invalid > The error message is: ('EOF in multi-line statement', (649, 0)) > > --------------------------------------------------------------------------- > PolynomialError Traceback (most recent call > last) > > /tmp/<ipython console> in <module>() > > /usr/lib/pymodules/python2.5/sympy/polys/factortools.pyc in factor(f, > *symbols, **flags) > 78 return > f > > 79 > ---> 80 coeff, factors = poly_factors(f, *symbols, > **flags) > > 81 > 82 result = 1 # XXX: don't include coeff in the leading > factor > > /usr/lib/pymodules/python2.5/sympy/polys/factortools.pyc in > poly_factors(f, *symbols, **flags) > 18 > """ > 19 if not isinstance(f, > Poly): > ---> 20 f = Poly(f, > *symbols) > 21 elif > symbols: > 22 raise SymbolsError("Redundant symbols were > given") > > /usr/lib/pymodules/python2.5/sympy/polys/polynomial.pyc in > __new__(cls, poly, *symbols, **flags) > 452 terms = Poly._permute(poly, > *symbols) > 453 > else: > --> 454 terms = Poly._decompose(poly, > *symbols) > > 455 > 456 if > terms: > > /usr/lib/pymodules/python2.5/sympy/polys/polynomial.pyc in > _decompose(poly, *symbols) > 596 > continue > > 597 > --> 598 raise PolynomialError("Can't decompose > %s" % factor) > 599 > else: > 600 coeff *= > factor > > PolynomialError: Can't decompose (1 - x**2 - y**2 + x**2*y**2)**0.5 > > -- > You received this message because you are subscribed to the Google Groups > "sympy" group. > To post to this group, send email to sy...@googlegroups.com. > To unsubscribe from this group, send email to > sympy+unsubscr...@googlegroups.com. > For more options, visit this group at > http://groups.google.com/group/sympy?hl=en. > -- You received this message because you are subscribed to the Google Groups "sympy" group. To post to this group, send email to sy...@googlegroups.com. To unsubscribe from this group, send email to sympy+unsubscr...@googlegroups.com. For more options, visit this group at http://groups.google.com/group/sympy?hl=en.
