# How-to: sum

Suppose I have defined a function

f(x)=x^2


and a list

x=[0,1,2,3,4]


Now I can calculate the product

x[2]*f(x[2])


But: If I use the command

k=var('k')
sum(x[k]*f(x[k]),k,0,4)


I will get the message:

TypeError: unable to convert x (=k) to an integer


How can I realize such a summation?

edit retag close merge delete

Sort by » oldest newest most voted

Here are some ways to sum the squares of the numbers from 0 to 1000, with timings:

sage: var('x');f = x^2;
sage: timeit('sum(f,x,0,1000)')
5 loops, best of 3: 271 ms per loop
sage: timeit('sum([f(x=t) for t in [0..1000]])')
5 loops, best of 3: 39 ms per loop
sage: timeit('sum([t^2 for t in [0..1000]])')
625 loops, best of 3: 1.21 ms per loop
sage: timeit('reduce(lambda x,y:x+y^2,[0..1000],0)')
625 loops, best of 3: 1.13 ms per loop
sage: timeit('reduce(lambda x,y:x+y^2,xrange(0,1001),0)')
625 loops, best of 3: 945 µs per loop


So if speed matters to you, you should learn a little python!

For example, the syntax [EXPR for VAR in LIST] is very handy - it evaluates EXPR with VAR ranging over the list (or any iterable object) LIST. You can even use it with a filtering:

sum([x^2 for x in [0..1000] if is_prime(x)])


will sum the squares of primes in this range.

To understand the last two summation examples above you should read about lambda, reduce, and xrange in a python tutorial. It will pay off!

more

Given the syntax of the sum() function, which can be read by executing "sum?" (with the question mark) in the Sage notebook, you don't need to define a list to compute the sum. For example,

sage: var('x')
sage: f = x^2   # note that you don't write f(x) = x^2
sage: sum(f,x,0,4)
30


Also, be careful defining a variable, x, as well as a list with the same name. (You wrote f = x^2 as well as x=[0,1,2,3,4].) Hope this helps.

more

The first syntax is symbolic summation:

sage: var('a,b,i')
sage: sum(2^i, i, a, b)
-2^a + 2^(b + 1)


It will always try to evaluate the sum symbolically. If you use something that doesn't understand symbolic variables (like GF(8), which is implemented by polybori) then you get the "unable to convert to integer" error. Because polybori can only exponentiate by integers, and not by formal variables.

The second syntax is Python's list creation and summation of lists:

sage: [ 2^j for j in range(0,3) ]
[1, 2, 4]
sage: sum(_)
7


Note that the variable inside the list comprehension (which I called j here) is automatically created and will overwrite other variables. In your second example, calling the inner variable i overwrites the previously-defined symbolic variable var('i')

more

This is a great thread, if one can call this a thread. Many different points of view and followups!

( 2011-05-17 03:00:07 -0500 )edit

I met a similar problem. Take the following example.

i=var('i')
K.<a> = GF(8, 'a')
for k in K:
print sum(k^(2^i),i,0,3)


which prints:

Traceback (click to the left of this block for traceback)
...
TypeError: unable to convert x (=i) to an integer


Let us try another syntax.

i=var('i')
K.<a> = GF(8, 'a')
for k in K:
sum([k^(2^i) for i in [0..3]])


Now it is okay thanks to the previous answer. I wonder if it's possible to make it work with the first syntax.

more

f(x)=x^2
g(x)=x*f(x)
k=var('k')
sum(g(k),k,0,4)


This will return 100. The point is that you can use this syntax to sum the values of a proper function, but not arbitrary Python expressions.

more