Status: New
Owner: ----
Labels: Type-Defect Priority-Medium
New issue 2831 by [email protected]: bug in simplify
http://code.google.com/p/sympy/issues/detail?id=2831
Dear sympy
simplify is not able to do its job on the simple eq that follows:
v
cat bug2.py
from sympy import symbols,simplify
x,y = symbols('x,y')
z =
0.0514453901238014*y*x**1.33333333333333*x**2.66666666666667/(x*(-0.0261211729852336*y
- 1.0*x**2.66666666666667)**2) +
0.0257226950619007*y*x**1.33333333333333/(x*(-0.0261211729852336*y -
1.0*x**2.66666666666667))
print simplify(z)
======
python2.7 bug2.py
Traceback (most recent call last):
File "bug2.py", line 6, in <module>
print simplify(z)
File "/usr/local/lib/python2.7/site-packages/sympy/simplify/simplify.py",
line 1627, in simplify
expr = together(cancel(powsimp(expr)).expand())
File "/usr/local/lib/python2.7/site-packages/sympy/polys/polytools.py",
line 5220, in cancel
c, P, Q = F.cancel(G)
File "/usr/local/lib/python2.7/site-packages/sympy/polys/polytools.py",
line 3022, in cancel
result = F.cancel(G, include=include)
File "/usr/local/lib/python2.7/site-packages/sympy/polys/polyclasses.py",
line 630, in cancel
cF, cG, F, G = dmp_cancel(F, G, lev, dom, include=False)
File "/usr/local/lib/python2.7/site-packages/sympy/polys/euclidtools.py",
line 1846, in dmp_cancel
_, p, q = dmp_inner_gcd(f, g, u, K)
File "/usr/local/lib/python2.7/site-packages/sympy/polys/euclidtools.py",
line 1556, in dmp_inner_gcd
h, cff, cfg = _dmp_inner_gcd(f, g, u, K)
File "/usr/local/lib/python2.7/site-packages/sympy/polys/euclidtools.py",
line 1522, in _dmp_inner_gcd
return dmp_ff_prs_gcd(f, g, u, K)
File "/usr/local/lib/python2.7/site-packages/sympy/polys/euclidtools.py",
line 1120, in dmp_ff_prs_gcd
h = dmp_subresultants(F, G, u, K)[-1]
File "/usr/local/lib/python2.7/site-packages/sympy/polys/euclidtools.py",
line 545, in dmp_subresultants
return dmp_inner_subresultants(f, g, u, K)[0]
File "/usr/local/lib/python2.7/site-packages/sympy/polys/euclidtools.py",
line 520, in dmp_inner_subresultants
h = dmp_prem(f, g, u, K)
File "/usr/local/lib/python2.7/site-packages/sympy/polys/densearith.py",
line 1281, in dmp_prem
c = dmp_pow(lc_g, N, u-1, K)
File "/usr/local/lib/python2.7/site-packages/sympy/polys/densearith.py",
line 1005, in dmp_pow
raise ValueError("can't raise polynomial to a negative power")
ValueError: can't raise polynomial to a negative power
<<<<<
--
You received this message because you are subscribed to the Google Groups
"sympy-issues" 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/sympy-issues?hl=en.
