No, you misunderstood my response. The bug has not been fixed yet (anywhere). I meant to say "we should fix the bug". Sorry for the ambiguity.
Aaron Meurer On Sun, Dec 1, 2013 at 6:07 AM, Leryra <[email protected]> wrote: > Well, it doesn't seem to work with the dev version from there : > https://github.com/sympy/sympy > Is there another version ? > > Here are my error's messages : > >>>> integrate(-3*sqrt(3)*(a**8 - 5*a**6 + 10*a**4 - 406*a**2 + >>>> 224)/(4*a**2*(a**6 + 6*a**5 + 15*a**4 + 20*a**3 + 15*a**2 + 6*a + 1)),a) > > Traceback (most recent call last): > File "<pyshell#0>", line 1, in <module> > integrate(-3*sqrt(3)*(a**8 - 5*a**6 + 10*a**4 - 406*a**2 + > 224)/(4*a**2*(a**6 + 6*a**5 + 15*a**4 + 20*a**3 + 15*a**2 + 6*a + 1)),a) > File "C:\Python27\lib\site-packages\sympy\utilities\decorator.py", line > 35, in threaded_func > return func(expr, *args, **kwargs) > File "C:\Python27\lib\site-packages\sympy\integrals\integrals.py", line > 1283, in integrate > risch=risch, manual=manual) > File "C:\Python27\lib\site-packages\sympy\integrals\integrals.py", line > 555, in doit > conds=conds) > File "C:\Python27\lib\site-packages\sympy\integrals\integrals.py", line > 811, in _eval_integral > result, i = risch_integrate(f, x, separate_integral=True, conds=conds) > File "C:\Python27\lib\site-packages\sympy\integrals\risch.py", line 1682, > in risch_integrate > ans = integrate(fa.as_expr()/fd.as_expr(), DE.x, risch=False) > File "C:\Python27\lib\site-packages\sympy\utilities\decorator.py", line > 35, in threaded_func > return func(expr, *args, **kwargs) > File "C:\Python27\lib\site-packages\sympy\integrals\integrals.py", line > 1283, in integrate > risch=risch, manual=manual) > File "C:\Python27\lib\site-packages\sympy\integrals\integrals.py", line > 555, in doit > conds=conds) > File "C:\Python27\lib\site-packages\sympy\integrals\integrals.py", line > 887, in _eval_integral > parts.append(coeff * ratint(g, x)) > File "C:\Python27\lib\site-packages\sympy\integrals\rationaltools.py", > line 71, in ratint > L = ratint_logpart(r, Q, x, t) > File "C:\Python27\lib\site-packages\sympy\integrals\rationaltools.py", > line 209, in ratint_logpart > res = Poly(res, t, composite=False) > File "C:\Python27\lib\site-packages\sympy\polys\polytools.py", line 87, in > __new__ > return cls._from_expr(rep, opt) > File "C:\Python27\lib\site-packages\sympy\polys\polytools.py", line 197, > in _from_expr > rep, opt = _dict_from_expr(rep, opt) > File "C:\Python27\lib\site-packages\sympy\polys\polyutils.py", line 331, > in _dict_from_expr > rep, gens = _dict_from_expr_if_gens(expr, opt) > File "C:\Python27\lib\site-packages\sympy\polys\polyutils.py", line 272, > in _dict_from_expr_if_gens > (poly,), gens = _parallel_dict_from_expr_if_gens((expr,), opt) > File "C:\Python27\lib\site-packages\sympy\polys\polyutils.py", line 177, > in _parallel_dict_from_expr_if_gens > raise PolynomialError("%s contains an element of the generators set" % > factor) > PolynomialError: 1/(16*_t**2 + 72*sqrt(3)*_t + 243) contains an element of > the generators set > > Or could it be that I wrongly installed the dev version ? > I just did "python setup.py install" on top of my current version... Is that > wrong ? > > Le dimanche 1 décembre 2013 12:40:21 UTC+1, Leryra a écrit : >> >> In which version should it be fixed ? In the dev one ? >> >> I'll try it. >> Thanks you. >> >> Le vendredi 29 novembre 2013 21:54:13 UTC+1, Aaron Meurer a écrit : >>> >>> It's a bug in ratint. The bug should be fixed, but regardless, didn't >>> we make it so that integrate automatically pull out common >>> coefficients before it calls the algorithms? >>> >>> Aaron Meurer >>> >>> On Fri, Nov 29, 2013 at 9:07 AM, Ondřej Čertík <[email protected]> >>> wrote: >>> > Hi Leryra, >>> > >>> > On Thu, Nov 28, 2013 at 3:16 PM, Leryra <[email protected]> wrote: >>> >> Hi there ! >>> >> >>> >> I'm using Sympy in an extensive way to do some work on slow blow up >>> >> solutions of the critical wave equation, >>> >> The thing is that I stumbled accross this : >>> >> integrate(-3*sqrt(3)*(a**8 - 5*a**6 + 10*a**4 - 406*a**2 + >>> >> 224)/(4*a**2*(a**6 + 6*a**5 + 15*a**4 + 20*a**3 + 15*a**2 + 6*a + >>> >> 1)),a) >>> >> >>> >> And I can't get why sympy throw me some "PolynomialError: 1/(16*_t**2 >>> >> + >>> >> 72*sqrt(3)*_t + 243) contains an element of the generators set" >>> >> >>> >> Could someone explain me this error ? Those are "just polynomials" and >>> >> the >>> >> integral doesn't look so ugly.. It has some log, but still... >>> > >>> > I think you found a bug. If you try the integral without the sqrt(3) >>> > at the beginning, it works: >>> > >>> >>>> integrate(-3*(a**8 - 5*a**6 + 10*a**4 - 406*a**2 + >>> >>>> 224)/(4*a**2*(a**6 + 6*a**5 + 15*a**4 + 20*a**3 + 15*a**2 + 6*a + >>> >>>> 1)),a) >>> > 5 4 3 2 >>> > 3⋅a 40620⋅a + 182640⋅a + 317590⋅a + 259775⋅a + 89869⋅a + 6720 >>> > - ─── + ───────────────────────────────────────────────────────────── + >>> > 1008⋅l >>> > 4 6 5 4 3 2 >>> > 40⋅a + 200⋅a + 400⋅a + 400⋅a + 200⋅a + 40⋅a >>> > >>> > >>> > 2007⋅log(a + 1) >>> > og(a) - ─────────────── >>> > 2 >>> > >>> > >>> > >>> > So we need to figure out what's going on and fix it. >>> > >>> > Ondrej >>> > >>> >> >>> >> -- >>> >> 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. > > -- > 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.
