Re: [sympy] Re: Limit - Bug

2014-02-18 Thread Christophe Bal
Thanks for this reference. Christophe Le 19 févr. 2014 00:39, "someone" a écrit : > Hi, > > > > what is the origin of this name gruntz ? > > "On Computing Limits in a Symbolic Manipulation System": > > ftp://ftp.inf.ethz.ch/doc/dissertations/th11432.ps.gz > > A nice text to read. > > -- > You r

Re: [sympy] Re: Limit - Bug

2014-02-18 Thread Aaron Meurer
Sergey pointed this out on an issue. The problem is this test from test_limits.py: assert limit((x**1000/((x + 1)**1000 + exp(-x))), x, oo) == 1 If you replace "limit" with "gruntz" it hangs. If your change doesn't hang this expression, then it's probably fine. Aaron Meurer On Tue, Feb 18,

Re: [sympy] Re: Limit - Bug

2014-02-18 Thread someone
Hi, > what is the origin of this name gruntz ? "On Computing Limits in a Symbolic Manipulation System": ftp://ftp.inf.ethz.ch/doc/dissertations/th11432.ps.gz A nice text to read. -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from

[sympy] Re: Limit - Bug

2014-02-18 Thread Avichal Dayal
ok = lambda w: (z in w.free_symbols and any(a.is_polynomial(z) or any(z in m.free_symbols and m.is_polynomial(z) for m in Mul.make_args(a)) for a in Add.make_args(w))) This is the ok function which checks whether it should apply t

Re: [sympy] Re: Limit - Bug

2014-02-18 Thread Christophe Bal
Thanks. I've just tried 2014-02-18 13:20 GMT+01:00 Sergey Kirpichev : > On Tuesday, February 18, 2014 12:51:50 PM UTC+4, Christophe Bal wrote: >> >> Hello, >> what is the origin of this name gruntz ? >> > > From the name of algorithm author. > > >> 2014-02-18 1:57 GMT+01:00 Aaron Meurer : >> >>

Re: [sympy] Re: Limit - Bug

2014-02-18 Thread Sergey Kirpichev
On Tuesday, February 18, 2014 12:51:50 PM UTC+4, Christophe Bal wrote: > > Hello, > what is the origin of this name gruntz ? > >From the name of algorithm author. > 2014-02-18 1:57 GMT+01:00 Aaron Meurer >: > >> Ah, so it seems that these incorrect limit heuristics have still not >> been complet

Re: [sympy] Re: Limit - Bug

2014-02-18 Thread Christophe Bal
Hello, what is the origin of this name gruntz ? 2014-02-18 1:57 GMT+01:00 Aaron Meurer : > Ah, so it seems that these incorrect limit heuristics have still not > been completely abolished from the code base. I would recommend just > using gruntz() for now. limit() is not that much smarter, but i

Re: [sympy] Re: Limit - Bug

2014-02-17 Thread Aaron Meurer
Ah, so it seems that these incorrect limit heuristics have still not been completely abolished from the code base. I would recommend just using gruntz() for now. limit() is not that much smarter, but it does have a bad tendency to give wrong answers when gruntz() works just fine. Aaron Meurer On

Re: [sympy] Re: Limit - Bug

2014-02-17 Thread Christophe Bal
Thanks for this. 2014-02-17 19:21 GMT+01:00 Avichal Dayal : > However gruntz gives the correct result:- > gruntz((x*exp(x)) / (exp(x)-1), x, -oo) gives 0 > > The part of code that is going wrong is already labeled as XXX: todo > More specifically the following:- > if abs(z0) is S.Infinity: >

[sympy] Re: Limit - Bug

2014-02-17 Thread Avichal Dayal
However gruntz gives the correct result:- gruntz((x*exp(x)) / (exp(x)-1), x, -oo) gives 0 The part of code that is going wrong is already labeled as XXX: todo More specifically the following:- if abs(z0) is S.Infinity: # XXX todo: this should probably be stated in the # neg

[sympy] Re: Limit - Bug

2014-02-17 Thread Avichal Dayal
A bug. Answer should be 0 Check wolfram's result: http://www.wolframalpha.com/input/?i=limit%28%28x*exp%28x%29%29%2F%28exp%28x%29-1%29%2C+x%2C+-oo%29 On Monday, 17 February 2014 22:45:38 UTC+5:30, Christophe Bal wrote: > > Hello. > > limit((x*exp(x))/(exp(x)-1), x, -oo) > > gives > > -oo > > Is