Status: New
Owner: ----
Labels: Type-Defect Priority-Medium
New issue 3575 by david.la...@wanadoo.fr: problem with method simplify()
http://code.google.com/p/sympy/issues/detail?id=3575
Hello,
I fall on the following error using the simplify() method on dFx :
Can you have a look to this ?
Thank you very much,
Cheers,
David Landriu
Python 2.6.5 (r265:79063, Apr 16 2010, 13:09:56)
[GCC 4.4.3] on linux2
Type "help", "copyright", "credits" or "license" for more information.
from sympy import exp,log,Symbol,pprint
x=Symbol('x')
dFx = - 0.00040033738827579*x*exp(x)**0.428344867337811 +
0.00250425613983611*x +
0.000357730249543611*exp(x)**0.428344867337811*log(0.011733044661259*exp(x))
- 0.00167981480713803*exp(x)**0.428344867337811 +
1.92510582333829e-5*exp(x)**0.856689734675622*log(0.011733044661259*exp(x))**2
dFx.simplify()
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "/home/landriu/sympy-0.7.2/sympy/core/expr.py", line 2733, in
simplify
return simplify(self)
File "/home/landriu/sympy-0.7.2/sympy/simplify/simplify.py", line 2911,
in simplify
expr1 = cancel(expr0)
File "/home/landriu/sympy-0.7.2/sympy/polys/polytools.py", line 5563, in
cancel
c, P, Q = F.cancel(G)
File "/home/landriu/sympy-0.7.2/sympy/polys/polytools.py", line 3178, in
cancel
result = F.cancel(G, include=include)
File "/home/landriu/sympy-0.7.2/sympy/polys/polyclasses.py", line 649, in
cancel
cF, cG, F, G = dmp_cancel(F, G, lev, dom, include=False)
File "/home/landriu/sympy-0.7.2/sympy/polys/euclidtools.py", line 1884,
in dmp_cancel
_, p, q = dmp_inner_gcd(f, g, u, K)
File "/home/landriu/sympy-0.7.2/sympy/polys/euclidtools.py", line 1588,
in dmp_inner_gcd
h, cff, cfg = _dmp_inner_gcd(f, g, u, K)
File "/home/landriu/sympy-0.7.2/sympy/polys/euclidtools.py", line 1553,
in _dmp_inner_gcd
return dmp_ff_prs_gcd(f, g, u, K)
File "/home/landriu/sympy-0.7.2/sympy/polys/euclidtools.py", line 1139,
in dmp_ff_prs_gcd
result = _dmp_ff_trivial_gcd(f, g, u, K)
File "/home/landriu/sympy-0.7.2/sympy/polys/euclidtools.py", line 961, in
_dmp_ff_trivial_gcd
return _dmp_simplify_gcd(f, g, u, K)
File "/home/landriu/sympy-0.7.2/sympy/polys/euclidtools.py", line 982, in
_dmp_simplify_gcd
F = dmp_content(f, u, K)
File "/home/landriu/sympy-0.7.2/sympy/polys/euclidtools.py", line 1804,
in dmp_content
cont = dmp_gcd(cont, c, v, K)
File "/home/landriu/sympy-0.7.2/sympy/polys/euclidtools.py", line 1631,
in dmp_gcd
return dmp_inner_gcd(f, g, u, K)[0]
File "/home/landriu/sympy-0.7.2/sympy/polys/euclidtools.py", line 1588,
in dmp_inner_gcd
h, cff, cfg = _dmp_inner_gcd(f, g, u, K)
File "/home/landriu/sympy-0.7.2/sympy/polys/euclidtools.py", line 1553,
in _dmp_inner_gcd
return dmp_ff_prs_gcd(f, g, u, K)
File "/home/landriu/sympy-0.7.2/sympy/polys/euclidtools.py", line 1147,
in dmp_ff_prs_gcd
h = dmp_subresultants(F, G, u, K)[-1]
File "/home/landriu/sympy-0.7.2/sympy/polys/euclidtools.py", line 561, in
dmp_subresultants
return dmp_inner_subresultants(f, g, u, K)[0]
File "/home/landriu/sympy-0.7.2/sympy/polys/euclidtools.py", line 535, in
dmp_inner_subresultants
h = dmp_prem(f, g, u, K)
File "/home/landriu/sympy-0.7.2/sympy/polys/densearith.py", line 1323, in
dmp_prem
c = dmp_pow(lc_g, N, u-1, K)
File "/home/landriu/sympy-0.7.2/sympy/polys/densearith.py", line 1041, 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 sympy-issues@googlegroups.com.
To unsubscribe from this group, send email to
sympy-issues+unsubscr...@googlegroups.com.
For more options, visit this group at
http://groups.google.com/group/sympy-issues?hl=en.
