This is a bug. It looks like it's been fixed in the git version. Aaron Meurer
On Sun, Nov 10, 2013 at 1:49 PM, George Gerber <[email protected]> wrote: > Good day, > > Why does the following code give a valueerror (sympy 0.7.3 windows)? > > from sympy import symbols, integrate > x,y,C1,C2,C3 = symbols('x y C1 C2 C3') > aa = -24596.7477524977*C1*x - 24596.7477524977*C2*y - 24596.7477524977*C3 - > 11284722.2222222*(C1*x + C2*y + C3)**2 > print integrate(aa,x) > > Best regards, > George > > ~~~~~~~~~~~~~~~~~~~~~Output of last call~~~~~~~~~~~~~~~~ > File "C:\Python27\lib\site-packages\sympy\integrals\integrals.py", line > 1154, > in _eval_integral > return poly.integrate().as_expr() > File "C:\Python27\lib\site-packages\sympy\polys\polytools.py", line 2126, > in i > ntegrate > f = f.to_field() > File "C:\Python27\lib\site-packages\sympy\polys\polytools.py", line 677, > in to > _field > result = f.rep.to_field() > File "C:\Python27\lib\site-packages\sympy\polys\polyclasses.py", line 267, > in > to_field > return f.convert(f.dom.get_field()) > File "C:\Python27\lib\site-packages\sympy\polys\polyclasses.py", line 278, > in > convert > return DMP(dmp_convert(f.rep, f.lev, f.dom, dom), dom, f.lev) > File "C:\Python27\lib\site-packages\sympy\polys\densebasic.py", line 576, > in d > mp_convert > return dup_convert(f, K0, K1) > File "C:\Python27\lib\site-packages\sympy\polys\densebasic.py", line 550, > in d > up_convert > return dup_strip([ K1.convert(c, K0) for c in f ]) > File "C:\Python27\lib\site-packages\sympy\polys\domains\domain.py", line > 94, i > n convert > return self.convert_from(element, base) > File "C:\Python27\lib\site-packages\sympy\polys\domains\domain.py", line > 84, i > n convert_from > result = _convert(element, base) > File > "C:\Python27\lib\site-packages\sympy\polys\domains\old_fractionfield.py", > line 91, in from_GlobalPolynomialRing > return K1(a.rep) > File "C:\Python27\lib\site-packages\sympy\polys\domains\domain.py", line > 69, i > n __call__ > return self.new(*args) > File > "C:\Python27\lib\site-packages\sympy\polys\domains\old_fractionfield.py", > line 34, in new > return self.dtype(element, self.dom, len(self.gens) - 1, ring=self) > File "C:\Python27\lib\site-packages\sympy\polys\polyclasses.py", line > 1013, in > __init__ > num, den = dmp_cancel(num, den, lev, dom) > File "C:\Python27\lib\site-packages\sympy\polys\euclidtools.py", line > 1907, in > dmp_cancel > _, p, q = dmp_inner_gcd(f, g, u, K) > File "C:\Python27\lib\site-packages\sympy\polys\euclidtools.py", line > 1618, in > dmp_inner_gcd > h, cff, cfg = _dmp_inner_gcd(f, g, u, K) > File "C:\Python27\lib\site-packages\sympy\polys\euclidtools.py", line > 1582, in > _dmp_inner_gcd > return dmp_ff_prs_gcd(f, g, u, K) > File "C:\Python27\lib\site-packages\sympy\polys\euclidtools.py", line > 1166, in > dmp_ff_prs_gcd > result = _dmp_ff_trivial_gcd(f, g, u, K) > File "C:\Python27\lib\site-packages\sympy\polys\euclidtools.py", line 989, > in > _dmp_ff_trivial_gcd > return _dmp_simplify_gcd(f, g, u, K) > File "C:\Python27\lib\site-packages\sympy\polys\euclidtools.py", line > 1011, in > _dmp_simplify_gcd > F = dmp_content(f, u, K) > File "C:\Python27\lib\site-packages\sympy\polys\euclidtools.py", line > 1832, in > dmp_content > cont = dmp_gcd(cont, c, v, K) > File "C:\Python27\lib\site-packages\sympy\polys\euclidtools.py", line > 1660, in > dmp_gcd > return dmp_inner_gcd(f, g, u, K)[0] > File "C:\Python27\lib\site-packages\sympy\polys\euclidtools.py", line > 1618, in > dmp_inner_gcd > h, cff, cfg = _dmp_inner_gcd(f, g, u, K) > File "C:\Python27\lib\site-packages\sympy\polys\euclidtools.py", line > 1582, in > _dmp_inner_gcd > return dmp_ff_prs_gcd(f, g, u, K) > File "C:\Python27\lib\site-packages\sympy\polys\euclidtools.py", line > 1174, in > dmp_ff_prs_gcd > h = dmp_subresultants(F, G, u, K)[-1] > File "C:\Python27\lib\site-packages\sympy\polys\euclidtools.py", line 570, > in > dmp_subresultants > return dmp_inner_subresultants(f, g, u, K)[0] > File "C:\Python27\lib\site-packages\sympy\polys\euclidtools.py", line 517, > in > dmp_inner_subresultants > h = dmp_prem(f, g, u, K) > File "C:\Python27\lib\site-packages\sympy\polys\densearith.py", line 1252, > in > dmp_prem > c = dmp_pow(lc_g, N, u - 1, K) > File "C:\Python27\lib\site-packages\sympy\polys\densearith.py", line 988, > in d > mp_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" 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/sympy. > For more options, visit https://groups.google.com/groups/opt_out. -- You received this message because you are subscribed to the Google Groups "sympy" 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/sympy. For more options, visit https://groups.google.com/groups/opt_out.
