A simple sum causes division by zero exception

7sDream
13 ●1 ●4

I'm very new to SageMath and did some small experiments yesterday. This is my Code:

x, n, r, i = var('x, n, r, i')
f(x, n, r) = x * sum((1 + r) ^ i, i, 1, n)
f(1, 30, 0.0)

It works well if r != 0, but will raise a exception of "division by zero" if equal.

My ipynb file: nbviewer.jupyter.org/gist/7sDream/5db3cfd153269fd1a1cacaf0b60f69bb

It seems some optimizations of sum did not consider the range of variable r. But there is also no change if I added assume(r >= 0)before define f.

So how can I disable this sum optimizations? Or what is wrong with the way I use it?

Comments

Curiously :

sage: def ef(x, n, r): return sum([(1+r)^u for u in (1..n)])
sage: ef(1, 30, 0)
30
sage: ef(1, 30, 0.0)
30.0000000000000

Emmanuel Charpentier ( 4 years ago )

Interestingly, if we define fc = fast_callable( x * sum((1 + r) ^ i, i, 1, n), vars=[x,n,r] ), then fc(1, 30, 0.0) evaluates to NaN.

Max Alekseyev ( 4 years ago )

Note that the "optimization" takes place at the creation of f :

sage: var("x, n, r, i")
(x, n, r, i)
sage: f(x,n,r,i)=x*sum((1+r)^i,i,1,n)
sage: f
(x, n, r, i) |--> ((r + 1)^(n + 1) - r - 1)*x/r

FWIW, sympy has an interesting translation :

>>> x, n, r, i = symbols("x, n, r, i")
>>> f = Lambda((x, n, r), x*summation((1+r)**i, (i, 1, n)))
>>> f
Lambda((x, n, r), x*Piecewise((n, Eq(r + 1, 1)), (-(r - (r + 1)**(n + 1) + 1)/r, True)))

Emmanuel Charpentier ( 4 years ago )

add a comment

sage: sum((1+r)^i, i, 1, n, algorithm="sympy") cases(((r + 1 == 1, n), (1, ((r + 1)^(n + 1) - r - 1)/r))) sage: f(x,n,r,i)=x*sum((1+r)^i,i,1,n, algorithm="sympy") sage: f (x, n, r, i) |--> x*cases(((r + 1 == 1, n), (1, ((r + 1)^(n + 1) - r - 1)/r))) sage: f(1, 30, 0) 30 sage: f(1, 30, 0.0) 30

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

A simple sum causes division by zero exception

Comments

1 Answer

Your Answer

Question Tools

Stats

Related questions

A simple sum causes division by zero exception savecancel

Comments

1 Answer

Your Answer

Question Tools

Stats

Related questions

A simple sum causes division by zero exception